SEARCH · Engineering Papers

Results for “Related References”

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.

SourceResearch connectionFrom yearThrough year

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

At least 451 records · Page 25

All sources NASA NTRS DOE OSTI

NASA NTRS2019

Quantifying the Evolution of Molecular Production Rates of Comet 21P/Giacobini–Zinner with iSHELL/NASA-IRTF

This paper presents results from our molecular characterization of comet 21P/Giacobini–Zinner during its 2018 apparition. We followed the comet during four observing runs, for a total of eight nights. The high-resolution spectra of 21P were acquired using iSHELL—the near-infrared high-resolution immersion echelle spectrograph on NASA/IRTF (Maunakea, Hawaii). We detected many cometary emission lines across four customized iSHELL settings in the (2.9–5.2) μm range. CO abundances relative to water displayed a relatively constant behavior during the observing campaign, suggesting a possible correlation between the outgassing of water and CO. While CO mixing ratios in 21P are consistent with the reference median value for Jupiter family comets (JFCs), ethane was depleted relative to the reference median value for JFCs, except for one observing run (August 8 and 9) in which ethane appeared overabundant. We consider several possible causes: an ethane outburst, decreasing outgassing of water, or a seasonal effect owing to nucleus rotation. In the pre-perihelion runs, methanol was overabundant compared with the reference median value for JFCs, however it decreased steadily to a depleted value during 21P’s post-perihelion phase, suggesting a possible seasonal effect. We report 3σ upper limits for acetylene, formaldehyde, ammonia, and methane. The ratios of our measured 3σ upper limits for acetylene and detections of hydrogen cyanide are consistent with the depletion of C2 relative to CN in 21P observed at optical wavelengths. This result conﬁrms that 21P is depleted in the carbon-chain primary volatile responsible for C2.

Faggi, S.↗

NASA NTRS2019

Quantifying the Evolution of Molecular Production Rates of Comet 21P/Giacobini–Zinner with iSHELL/NASA-IRTF

This paper presents results from our molecular characterization of comet 21P/Giacobini-Zinner during its 2018 apparition. We followed the comet during four observing runs, for a total of 8 nights. The high-resolution spectra of 21P were acquired using iSHELL - the near-IR high resolution immersion echelle spectrograph on NASA/IRTF (Mauna Kea, Hawaii). We detected many cometary emission lines across four customized iSHELL settings in the (2.9 – 5.2) μm range. CO abundances relative to water displayed a relatively constant behavior during the observing campaign, suggesting a possible correlation between the outgassing of water and CO. While CO mixing ratios in 21P are consistent with the reference median value for Jupiter Family comets (JFCs), ethane was depleted relative to the reference median value for JFCs, excepting one observing run (8th and 9th Aug) in which ethane appeared overabundant. We consider several possible causes: an ethane outburst, decreasing outgassing of water, or a seasonal effect owing to nucleus rotation. In the pre-perihelion runs, methanol was overabundant compared with the reference median value for JFCs, however it decreased steadily to a depleted value during 21P’s post-perihelion phase, suggesting a possible seasonal effect. We report 3-s upper limits for acetylene, formaldehyde, ammonia and methane. The ratios of our measured 3-s upper limits for acetylene and detections of hydrogen cyanide are consistent with the depletion of C2 relative to CN in 21P observed at optical wavelengths. This result confirms that 21P is depleted in the carbon chain primary volatile responsible for C2 production.

Faggi S.↗

NASA NTRS1974

Relative motion of orbiting particles under the influence of perturbing forces. Volume 2: Analytical results

The mathematical developments carried out for this investigation are reported. In addition to describing and discussing the solutions which were acquired, there are compendia of data presented herein which summarize the equations and describe them as representative trace geometries. In this analysis the relative motion problems have been referred to two particular frames of reference; one which is inertially aligned, and one which is (local) horizon oriented. In addition to obtaining the classical initial values solutions, there are results which describe cases having applied specific forces serving as forcing functions. Also, in order to provide a complete state representation the speed components, as well as the displacements, have been described. These coordinates are traced on representative planes analogous to the displacement geometries. By this procedure a complete description of a relative motion is developed; and, as a consequence range rate as well as range information is obtained.

Eades, J. B., Jr.↗

DOE OSTI2021

SMC 2021 : Analyzing Resource Utilization and User Behavior on Titan Supercomputer

Resource utilization statistics of submitted jobs on a supercomputer can help us understand how users from various scientific domains use HPC platforms and better design a job scheduler. We explore to generate insight regarding workload distribution and usage pattern domains from job scheduler trace, GPU failure information, and project-specific information collected from Titan supercomputer. Furthermore, we want to know how the scheduler performance varies over time and how the users' scheduling behavior changes following a system failure. These observations have the potential to provide valuable insight, which is helpful to prepare for system failures. These practices will help us develop and apply novel machine learning algorithms in understanding system behavior, requirement, and better scheduling of HPC systems. There are two datasets, RUR and GPU. RUR: This dataset is the job scheduler traces collected from the Titan supercomputerfrom 01/01/2015 to 07/31/2019 (2015.csv - 2019.csv). These were collected usingResource Utilization Report (RUR), a Cray-developed resource-usage data collectionand reporting system. It contains the usage information of its critical resources (CPU,Memory, GPU, and I/O) of each running job on Titan during that period [2]. ProjectAreas: Every job is associated with a project ID. TheProjectAreas.csvdatasetprovides a mapping of the project ID to its domain science. GPU: There have been some hardware-related issues in the GPUs in Titan that caused some GPUs to fail, sometimes irrecoverably during some job runs. This dataset provides information regarding these failures during the execution of the submitted jobs. GPUs on Titan are uniquely identified by a serial number (SN), and they are installed in a location. A GPU can be installed in a location, then removed from that location following a failure, and then re-installed in a different location after fixing the problem. If the failure can't be recovered, the GPU might be removed entirely from Titan. There are two prominent types of failures that resulted in the removal of GPUs from Titan: Double Bit Error (DBE) and Out of the Bus (OTB). The dataset (gc_full.csv) has the following fields: 1. SN : Serial number of a GPU 2. location : The location where it is installed 3. insert : The time when it was inserted into that location 4. remove : The time when it was removed from that location 5. duration : Amount of time the GPU spent in this location 6. out : If the device was taken out entirely w/o a re-installment into a new location. 7. event : If the GPU was taken out entirely, the reason for its removal. To learn more about this dataset, please refer to the git repositoryhttps://github.com/olcf/TitanGPULifeand the related publication [1]. References [1] George Ostrouchov, Don Maxwell, Rizwan A Ashraf, Christian Engelmann, MallikarjunShankar, and James H Rogers. Gpu lifetimes on titan supercomputer: Survival analysisand reliability. InSC20: International Conference for High Performance Computing,Networking, Storage and Analysis, pages 1-14. IEEE, 2020. [2] Feiyi Wang, Sarp Oral, Satyabrata Sen, and Neena Imam. Learning from five-yearresource-utilization data of titan system. In2019 IEEE International Conference onCluster Computing (CLUSTER), pages 1-6. IEEE, 2019.

42 ENGINEERING↗

NASA NTRS1967

Device measures static friction of magnetic tape

Device measures the coefficient of static friction of magnetic tape over a range of temperatures and relative humidities. It uses a strain gage to measure the force of friction between a reference surface and the tape drawn at a constant velocity of approximately 0.0001 inch per second relative to the reference surface.

Cole, P. T.↗

DOE OSTI2021

Revisiting AMOC Transport Estimates From Observations and Models

Abstract Reference level assumptions used to calculate the Atlantic meridional overturning circulation transports at the RAPID and Meridional Overturning Variability Experiment (MOVE) observing arrays are revisited in an eddying ocean model. Observational transport calculation methods are complemented by several alternative approaches. At RAPID, the model transports from the observational method and the model truth (based on the actual model velocities) agree well in their mean and variability. There are substantial differences among the transport estimates obtained with various methods at the MOVE site. These differences result from relatively large and time‐varying reference velocities at depth in the model, not supporting a level‐of‐no‐motion. The methods that account for these reference velocities properly at MOVE produce transports that are in good agreement with the model truth. In contrast with the observational estimates, the model transport trends at MOVE and RAPID largely agree with each other on pentadal to multi‐decadal time scales.

Danabasoglu, Gokhan↗

NASA NTRS1987

Role of orientation reference selection in motion sickness, supplement 2S

Previous experiments with moving platform posturography have shown that different people have varying abilities to resolve conflicts among vestibular, visual, and proprioceptive sensory signals. The conceptual basis of the present proposal hinges on the similarities between the space motion sickness problem and the sensory orientation reference selection problems associated with benign paroxysmal positional vertigo (BPPV) syndrome. These similarities include both etiology related to abnormal vertical canal-otolith function, and motion sickness initiating events provoked by pitch and roll head movements. The objectives are to explore and quantify the orientation reference selection abilities of subjects and the relation of this selection to motion sickness in humans. The overall objectives are to determine: if motion sickness susceptibility is related to sensory orientation reference selection abilities of subjects; if abnormal vertical canal-otolith function is the source of abnormal posture control strategies and if it can be quantified by vestibular and oculomotor reflex measurements, and if it can be quantified by vestibular and oculomotor reflex measurements; and quantifiable measures of perception of vestibular and visual motion cues can be related to motion sickness susceptibility and to orientation reference selection ability.

Peterka, Robert J.↗

NASA NTRS1989

Corotation lag of the Jovian atmosphere, ionosphere, and magnetosphere

The Jovian ionosphere-magnetosphere coupling model of Hill (1979) was modified to include the rotational slippage of the neutral atmosphere at ionospheric heights, relative to a frame of reference corotating rigidly with Jupiter. Equations were derived for the altitude distributions of ionospheric neutral and ion velocities, and a generalized expression was obtained for the corotation lag of the magnetosphere. The results of calculations provide independent support for the expectation that vertical mixing in Jupiter's atmosphere is much more vigorous at high latitudes than near the equator. They also indicate that the observed corotation lag in the magnetosphere and the Io torus is largely attributable to the slippage of the neutral atmosphere itself, rather than to the slippage of ionospheric ions relative to ionospheric neutrons, as previously suggested.

Huang, T. S.↗

NASA NTRS1994

Interpretation of waverider performance data using computational fluid dynamics

A computational study was conducted to better understand experimental results obtained from wind-tunnel tests of a Mach 4 conical-flow-derived waverider and a comparative reference configuration, which showed that the aerodynamic performance of the reference configuration was slightly better than that of the waverider. The computational results showed that the predicted surface pressure values and the integrated lift and drag coefficients were much lower for the reference model because the reference model bottom is an expansion surface. However, the lift-drag ratios for the reference model were higher due to a relatively low drag for a comparable amount of lift. The results also showed that the reference model exhibited the same shock attachment characteristics as the conical-flow-derived waverider, and is therefore also a waverider. The shock attachment characteristic gives the waverider a performance advantage over conventional hypersonic vehicles, and the results suggest that altering the bottom surface does not cause significant performance degradation. Flowfield solutions also show that the conical-flow waverider model has better propulsion/airframe integration characteristics than the reference configuration. The results also suggest that generating flowfields other than conical ones may be used to design waveriders with improved aerodynamic performance.

COMPUTATIONAL FLUID DYNAMICS↗

DOE OSTI2025

Laser Confocal Microscopy Uncertainty Quantification Study

At Los Alamos National Laboratory (LANL), the Storage Safety and Engineering (SSE) team completes annual surveillance on a subset of in-use interim nuclear material storage containers in fulfilment of requirements outlined in DOE Manual M 441.1-1. The containers are selected through several methods, such as subject matter expert judgement, random selection, and trending items. Following these selections, the SSE team has the capacity to complete surveillance on 15-20 containers each fiscal year, composed of a combination of SAVY-4000 and Hagan storage containers. Through previous work, the stainless-steel components of the containers have been identified as life limiting components, with an emphasis on the thin-walled bodies. The team is focused on understanding the extent of general and pitting corrosion, due to observations of extensive corrosion from stored contents and bag-out-bag degradation. Quantifying corrosion effects on the thin-walled stainless steel container bodies, and understanding potential impacts to the respective design release rates and design qualification release rates is paramount to the team. To date, destructive examination (DE) has proven to be the most insightful method for developing an understanding on the extent of corrosion on used containers. To standardize this process, the SSE team developed a destructive examination guide for analyzing stainless steel components of the containers. Corroded containers of interest are identified during surveillance activities and set aside for sectioning and characterization. Following sectioning, a major step in the DE workflow is the utilization of laser confocal microscopy for scanning corroded samples of interest and extracting data on pits, such as count, depth, and equivalent diameter. Adhering to the techniques outlined in the DE guide, analysis has been completed on two Hagans and one SAVY-4000 container, with the maximum pit depth recorded as 139.1 ± 22.82 μm on a 17.5 year old Hagan. The findings from the completed destructive examinations will be utilized to support lifetime extension efforts of the SAVY-4000 as the team can better estimate corrosion rates and effects over time based on stored contents and age. Due to the implications of observing extreme pit depths that approach the nominal container body thickness of .0299 inches (0.759 mm) or minimum container thickness of 0.236” (0.6 mm), high confidence in the LCM measurements is desired. Through testing outlined in, it was concluded that the total error ascribed to the 20x objective when conducting large image mapping on the Keyence VK-X3050 laser confocal microscope (LCM) relative to a 50x objective (reference) is 16.4% (± 8.73%). For shallow features on the order of pristine SAVY surface defects (i.e. 5 μm), this uncertainty is appropriate. However, this conservative estimate of total error poses a fundamental concern for pit depths that approach the thickness of the measured samples. That is, with the measurement uncertainty currently employed on all measurements, the LCM would be unable to resolve if a pit with a depth of 515 μm is through wall. Standard step height samples were procured and used in the present study to assess the resolution and repeatability of height measurements. Understanding the resolution and repeatability of height measurements was the first focus of the team as it relates directly to pit depth, which is of primary concern. Calibration gratings were procured to evaluate the resolution and repeatability of measurements in the X and Y axes of the LCM stage. The results of the depth uncertainty study were conducted first and presented in the subsequent sections. The planar uncertainty study is appended to the depth study with conclusions from both summarized at the end of the report.

36 MATERIALS SCIENCE↗

NASA NTRS2002

Geoscience Laser Altimeter System (GLAS) for the ICESat Mission

The Geoscience Laser Altimeter System (GLAS) is a new generation lidar and is the primary science payload for NASA's ICESat Mission. The GLAS design combines a 10 cm precision surface lidar with a sensitive dual wavelength cloud and aerosol lidar. GLAS will precisely measure the heights of the Earth's polar ice sheets, establish a grid of accurate height profiles of the Earth's land topography, and profile the vertical distribution of clouds and aerosols on a global scale. GLAS will be integrated onto a small spacecraft built by Ball Aerospace, and will be launched into a polar orbit with a 590-630 km altitude at an inclination of 94 degrees. ICESat is is currently planned to launch in winter 2002/03 and GLAS is designed to operate continuously in space for a minimum of 3 years. GLAS will measure the vertical distance from orbit to the Earth's surface with pulses from a ND:YAG laser at a 40 Hz rate. Each 6 nsec wide 1064 nm laser pulse is used to produce a single range measurement. On the surface, the laser footprints have 66 m diameter and approx. 170 m center-center spacings. The GLAS receiver uses a I m diameter telescope to detect laser backscatter and a Si APD to detect the 1064 nm signals. The detector's output is sampled by a digital ranging receiver, which records each transmitted pulse and surface echo waveform with 1 nsec (15 cm) resolution. Each echo pulse is digitized and is reported to ground with a record length of from 200 to 544 samples, depending on the spacecraft's location . The GLAS location and epoch times are measured by a precision GPS receiver carried on the ICESat spacecraft. Initial processing of the echo waveforms within GLAS permits discrimination between cloud and surface echoes for selecting appropriate waveform samples. This selection is guided by an on-board DEM which is used to set the boundaries for the echo pulse search algorithm. Subsequent ground-based echo pulse analysis, along with GPS-based clock frequency estimates, permit final determination of the range to the surface, degree of pulse spreading, and vertical distribution of any vegetation illuminated by the laser. Accurate knowledge of the laser beam's pointing angle is needed to prevent height biases when measuring over tilted surfaces, such as near the boundaries of ice sheets. For surfaces with 2 deg. slopes, knowledge of pointing angle of the beam's centroid angle to better than 10 urad is needed. GLAS uses a stellar reference system (SRS) to measure the pointing angle of each laser firing relative to inertial space. The SRS uses a high precision star camera oriented toward local zenith and a gyroscope to determine the inertial orientation of the SRS optical bench. The far field pattern of each laser is measured pulse relative to the star camera with a laser reference system (LRS). GLAS will also measure the vertical distributions of clouds and aerosols by recording the vertical profiles of laser pulse backscatter at both 1064 and 532 nm. The 1064 rim measurements use the Si APD detector and will be used to measure the height and echo pulse shape from thicker clouds. The lidar receiver at 532 nm uses a narrow bandwidth etalon filter and highly sensitive photon counting detectors. The 532 nm backscatter profiles will be used to measure the vertical extent of thinner clouds and the atmospheric boundary layer. The GLAS instrument component development is complete and the instrument is undergoing final testing and qualification at NASA-Goddard. The GLAS "as-built" characteristics and its expected measurement performance will be discussed.

Abshire, James B.↗

NASA NTRS1974

SPAR reference manual

SPAR is a system of related programs which may be operated either in batch or demand (teletype) mode. Information exchange between programs is automatically accomplished through one or more direct access libraries, known collectively as the data complex. Card input is command-oriented, in free-field form. Capabilities available in the first production release of the system are fully documented, and include linear stress analysis, linear bifurcation buckling analysis, and linear vibrational analysis.

Whetstone, W. D.↗

NASA NTRS1950

Contribution to the problem of buckling of orthotropic plates, with special reference to plywood

Planar stress-strain relations and bending stress-strain relations are presented for elastic orthotropic plates and specialized to plywood. These relations are used to derive the differential equation and energy expression for the buckling of orthotropic rectangular plates whose principal stiffness directions are not parallel to the plate edges. Buckling analyses are made for the case of pure compression and pure shear of a long plate-strip.

WOODS↗

NASA NTRSDate not supplied

Effects of Communication Modality on Pilot-Controller Coordination during a Simulated m:N Operation

The last decade or so has seen growing interest in new control paradigms and concepts of operation for uncrewed aircraft systems (UAS) in which multiple aircraft are piloted remotely by a single or relatively small number of people. Referred to as “one-to-many” and “many-to- many” (alternatively, “multi-operator, multi-vehicle”)—and frequently expressed as the corresponding ratios, 1:N and m:N—such novel configurations of aircraft and the people who manage them are seen as critical to the path to future operations involving UAS. Examples of industry domains interested in these control paradigms are small package delivery services utilizing small UAS and passenger-carrying, short-range “Urban Air Mobility” (UAM) operations. Stakeholders in such operations have identified communication and coordination of flight activity with air traffic controllers (ATC) as a barrier to operations. In contrast to present-day flight operations, in which a pilot communicates with one ATC on one radio frequency for one aircraft, multi-vehicle operations potentially entail a significant increase in pilot task load for management of comms. New concepts, such as UAS Service Suppliers (USSs) and Providers of Services to UAM (PSUs), have been proposed to address the known bottleneck for Air Traffic Management (ATM) presented by multi-vehicle operations. While progress has been steadily made over years developing USSs and PSUs, it is generally expected that initial UAM operations will rely on traditional voice-over-radio communication with ATC for purposes of ATM. The current study was a human-in-the-loop simulation that had participants, each possessing a Private Pilot License, act as the ground-based pilot-in- command for multiple vehicles in a hypothetical UAM service in the San Francisco Bay Area. The experiment utilized a 2-by-3, within-subjects design in which the pilot’s Vehicle Load (4 vs. 12) and Comm System (Voice, Datalink, and a Hybrid) were manipulated. The task given to pilots was to use the Comm System to coordinate flight activity for all aircraft with appropriate controllers, having to obtain departure and arrival clearances at “vertiport” facilities and transition clearances for any intermediate airspaces along the route. Pilots were additionally responsible for compliance with vectoring instructions issued by ATC. Subjective workload questionnaires (NASA-TLX) were administered following each experimental trial. Screen recordings of the pilot’s Ground Control Station (GCS) and audio recordings of trials were subsequently coded to obtain performance metrics: response times and error rates. Presented in this paper are results related to pilot responses to vectoring instructions issued by ATC. Workload was found to be significantly higher in the 12-Vehicle condition compared to the 4-Vehicle condition, nearly maxing out the NASA-TLX overall workload scale. There was no significant difference made by the Comm System on workload ratings. Pilots’ response times to communications were fastest in the Voice condition, although overall “service time” for compliance was shorter in Datalink and Hybrid conditions in most cases. Errors by pilots were frequent in both Vehicle Load conditions, most perniciously when using the Voice system. The results of this study suggest tradeoffs in advantages and disadvantages of the three comm systems. Recommendations for communication system design are provided taking the tradeoffs into account.

Garrett G Sadler↗

NASA NTRSDate not supplied

Effects of Communication Modality on Pilot-Controller Coordination during a Simulated m:N Operation

urban air mobility↗

DOE OSTI2021

Effects of knock intensity measurement technique and fuel chemical composition on the research octane number (RON) of FACE gasolines: Part 1 – Lambda and knock characterization

The Research and Motor Octane Number (RON and MON) rate the knock propensity of gasoline in the Cooperative Fuel Research (CFR) engine by comparing the knock intensity of sample fuels relative to that of primary reference fuels (PRF), a binary blend of iso-octane and n-heptane. Important differences exist between standard octane testing and automotive spark ignition (SI) engine knock testing including speed, load, air-to-fuel equivalence ratio (lambda), and knock characterization, which lead to a discrepancy between a fuel’s RON rating and its knock resistance characterized on an automotive SI engine based on knock-limited spark advance. This publication is the first of a set of three publications which modify operating parameters of the RON test method (ASTM D2699) to investigate the effects of these differences with automotive SI engine knock-limited spark advance testing. A fuel’s standard RON is evaluated at the lambda of the highest knock intensity, whereas automotive SI engines typically operate at stoichiometry. Here, we analyze the effects of a stoichiometric RON rating methodology. Additionally, the knock intensity response from the standard CFR knockmeter system is compared to a cylinder pressure oscillation-based knock intensity at several lambda settings. All experiments were performed with a set of seven Coordinating Research Council (CRC) Fuels for Advanced Combustion Engines (FACE) gasolines with approximately 95 RON. The fuel chemical composition impacted the lambda of the highest knock intensity, which resulted in fuel-specific offsets between the standard and stoichiometric RON ratings. The knock system comparison showed significant offsets between cylinder pressure-based and knockmeter-based knock intensity levels.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Effects of knock intensity measurement technique and fuel chemical composition on the research octane number (RON) of FACE gasolines: Part 2 – Effects of spark timing

The Research and Motor Octane Number (RON, MON) characterize a fuel’s knock resistance by rating the knock intensity of a sample fuel relative to that of Primary Reference Fuels (PRF) in a Cooperative Fuel Research (CFR) Engine. A fuel’s octane number is regulated to prevent damage from autoignition leading to knocking combustion in spark-ignition engines. The operational differences between the standard RON rating and modern engine operation are explored in a three-part publication series. The previous study focused on the effects of lambda and knock characterization. This second study primarily focuses on the effects of spark timing on RON determination. Following the findings from the first publication, the knock intensity was captured by the knockmeter and by the maximum amplitude of pressure oscillations (MAPO) at the lambda of peak knock intensity and stoichiometry. Knock-limited spark advance tests were conducted for a set of seven Fuels for Advanced Combustion Engines (FACE) from the Coordinating Research Council (CRC) with varying chemical composition, PRFs, and Toluene Standardization Fuels (TSFs). For retarded spark timings, pre-spark low-temperature heat release was found for low RON PRFs. Low RON PRFs also showed knocking characteristics before reaching the center of combustion suggesting that the use of knock-limited spark advance (KLSA) was preferred over the knock-limited combustion phasing. Primarily paraffinic fuels tended towards increased pressure oscillations while dominantly aromatic fuels experienced higher pressure rise rates. A MAPO-based KLSA correlated best to Octane Index at a negative K-factor suggesting beyond RON operation despite being at otherwise RON conditions. At stoichiometry, the MAPO-based KLSA did neither correlate to RON nor Octane Index. Finally, good agreement was found between KLSA-based effective RON from this study to the MAPO-based effective RON from the first study.

10 SYNTHETIC FUELS↗

DOE OSTI2025

Computing material volume fractions on a superimposed mesh as applied to Monte Carlo particle transport simulations

Here, we present a newly implemented ray tracing algorithm in OpenMC for efficiently computing material volume fractions on superimposed meshes in complex geometries. By firing rays along each coordinate direction through the geometry, the approach accumulates track-length data in each mesh element, thereby determining the fractional composition of each material. Scaling studies on three different models—a random tetrahedra configuration, the Frascati Neutron Generator ITER dose rate benchmark, and a stellarator design—show excellent parallel performance, with nearly linear speedup on modern multi-threaded and distributed-memory systems. An analysis of the residual error relative to high-resolution reference solutions demonstrated that under optimal conditions it decreases as 1/R, where R is the number of rays fired, making it straightforward to achieve user-prescribed accuracy. This new functionality enables practical, mesh-based approaches for detailed nuclear analyses in production Monte Carlo workflows without resorting to expensive, fully conformal or unstructured meshing.

Monte Carlo↗