Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “tracking errors”

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 1,783 records · Page 99
All sourcesNASA NTRSDOE OSTI
NASA NTRS2000

Nozomi Cis-Lunar Phase Orbit Determination

Japan's Institute of Space and Astronautical Science (ISAS) launched Nozomi, its first mission to the planet Mars using the newly developed M-V launch vehicle on July 3, 1998. Scientific objectives of the mission are to study the structure and dynamics of the Martian upper atmosphere and its interaction with the solar wind. Nozomi is a cooperative mission between ISAS and the National Aeronautics and Space Administration (NASA). The NASA contribution includes navigation and tracking services provided by the Jet Propulsion Laboratory (JPL). The spacecraft also serves as an engineering demonstration of basic technology for planetary exploration. One of the new technologies was a unique trajectory, developed by ISAS, which used solar gravitational perturbations at the weak stability boundary as an aid to achieve an Earth-Mars transfer orbit. This trajectory saves approximately 120 m/s of Delta V compared to direct hyperbolic insertion and is considered an enabling technology for the mission. Nozomi was the first spacecraft to employ this trajectory and provided on-orbit validation of the technique. The trajectory was achieved by initially placing the spacecraft in a highly elliptical cis-lunar phasing orbit. Six maneuvers were performed during this period to correct injection errors and target an outbound lunar swingby in September 1998. The gravity assist from the lunar swingby raised apogee to the vicinity of the weak stability boundary. After three more targeting maneuvers, Nozomi performed an inbound lunar swingby followed immediately by a powered Earth swingby in late December 1998. A 420 m/s Trans Mars Insertion (TMI) burn at the final Earth periapsis was intended to place the spacecraft on a heliocentric trajectory leading to Mars orbit insertion in October 1999. Orbit determination for Nozomi is performed in parallel by both ISAS and the Multi-Mission Navigation (MMNAV) group at JPL. This was an advantage for the mission because each group would generate solutions based on data collected from their respective tracking networks. Spacecraft events, such as sequence uplinks and maneuvers, were generally scheduled during passes at the Usuda tracking station in Japan. As a result, maneuver design and reconstruction was derived from MMNAV solutions based on JPL tracking data obtained immediately prior to or following maneuvers. Data was also exchanged between ISAS and MMNAV so orbit determination could be performed on joint data sets in support of critical targeting late in the cis-lunar phase. In this paper, information regarding the MMNAV orbit determination effort for the first six months of the mission is presented. The spacecraft trajectory is characterized first, followed by a discussion of the orbit determination estimation procedure and models. Results from selected orbit solutions are presented and compared against reconstructed trajectories. One area of emphasis in this paper is orbit determination in the vicinity of the weak stability boundary. Precise navigation was necessary to target the second lunar swingby and the powered Earth swingby. Delivery accuracy of 150 m was required for these critical encounters, but a number of factors contributed to the general degradation of orbit determination accuracy. This included the fact that the spacecraft was at apogee, at a range of 1.7 million km and moving at less than I km/sec perpendicular to the line of sight. Nozomi was also close to zero degrees declination where there are known limitations on orbit determination performance. Finally, S-band tracking data was acquired through the Nozomi backup low gain antenna. This antenna is offset from the axis of this spin stabilized spacecraft and superimposed large signatures in the Doppler and range data. These difficulties were overcome by combining long data arcs, spanning several maneuvers, with a high fidelity solar pressure model. The model included a physically accurate representation of the spacecraft structure and a high time resolution orientation model. Observation modeling included the removal of the spin induced Doppler bias, spin signature and per pass correction of range calibration errors applied for data leading up to critical events. As a result, all orbit determination goals were met. A second area of emphasis in this paper is the JPL tracking and orbit determination effort in support of the TMI maneuver. TMI occurred out of contact with ground stations and the JPL Goldstone tracking complex had the first pass following the bum. As a result, MMNAV had the responsibility to make a rapid assessment of the maneuver performance. MMNAV made the determination that a 100 m/s under bum had occurred and promptly informed ISAS via voice lines. ISAS immediately began preparations for a correction maneuver (TMIc), which had to be performed during the next Usuda pass. The near real time assessment by MMNAV provided accurate antenna frequency and pointing updates for the spacecraft acquisition at Usuda and the close coordination between the two agencies enabled the design and successful execution of the TMc maneuver. Propellant consumption during the correction burn dictated that the mission be redesigned. ISAS developed a new plan which adds 3 full solar orbits, two Earth swingbys and one lunar swingby with arrival at Mars in January 2004. The final Mars orbit will still enable the mission to achieve all of its science objectives.

Ryne, Mark↗
DOE OSTI2013

Calibration of time-of-flight cameras for accurate intraoperative surface reconstruction

Purpose: In image-guided surgery (IGS) intraoperative image acquisition of tissue shape, motion, and morphology is one of the main challenges. Recently, time-of-flight (ToF) cameras have emerged as a new means for fast range image acquisition that can be used for multimodal registration of the patient anatomy during surgery. The major drawbacks of ToF cameras are systematic errors in the image acquisition technique that compromise the quality of the measured range images. In this paper, we propose a calibration concept that, for the first time, accounts for all known systematic errors affecting the quality of ToF range images. Laboratory and in vitro experiments assess its performance in the context of IGS.Methods: For calibration the camera-related error sources depending on the sensor, the sensor temperature and the set integration time are corrected first, followed by the scene-specific errors, which are modeled as function of the measured distance, the amplitude and the radial distance to the principal point of the camera. Accounting for the high accuracy demands in IGS, we use a custom-made calibration device to provide reference distance data, the cameras are calibrated too. To evaluate the mitigation of the error, the remaining residual error after ToF depth calibration was compared with that arising from using the manufacturer routines for several state-of-the-art ToF cameras. The accuracy of reconstructed ToF surfaces was investigated after multimodal registration with computed tomography (CT) data of liver models by assessment of the target registration error (TRE) of markers introduced in the livers.Results: For the inspected distance range of up to 2 m, our calibration approach yielded a mean residual error to reference data ranging from 1.5 ± 4.3 mm for the best camera to 7.2 ± 11.0 mm. When compared to the data obtained from the manufacturer routines, the residual error was reduced by at least 78% from worst calibration result to most accurate manufacturer data. After registration of the CT data with the ToF data, the mean TRE ranged from 3.7 ± 2.1 mm for point-based and 5.7 ± 1.9 mm for surface-based registration for the best camera to 6.2 ± 3.4 and 11.1 ± 2.8 mm, respectively. Compared to data provided by the manufacturer, the mean TRE decreased by 8%–60% for point-based and by 18%–74% for surface-based registration.Conclusions: Using the proposed calibration approach improved the measurement accuracy of all investigated ToF cameras. Although evaluated in the context of intraoperative image acquisition, the proposed calibration procedure can easily be applied to other medical applications using ToF cameras, such as patient positioning or respiratory motion tracking in radiotherapy.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2017

Solar Field Optical Characterization at Stillwater Geothermal/Solar Hybrid Plant

Concentrating solar power (CSP) can provide additional thermal energy to boost geothermal plant power generation. For a newly constructed solar field at a geothermal power plant site, it is critical to properly characterize its performance so that the prediction of thermal power generation can be derived to develop an optimum operating strategy for a hybrid system. In the past, laboratory characterization of a solar collector has often extended into the solar field performance model and has been used to predict the actual solar field performance, disregarding realistic impacting factors. In this work, an extensive measurement on mirror slope error and receiver position error has been performed in the field by using the optical characterization tool called Distant Observer (DO). Combining a solar reflectance sampling procedure, a newly developed solar characterization program called FirstOPTIC and public software for annual performance modeling called System Advisor Model (SAM), a comprehensive solar field optical characterization has been conducted, thus allowing for an informed prediction of solar field annual performance. The paper illustrates this detailed solar field optical characterization procedure and demonstrates how the results help to quantify an appropriate tracking-correction strategy to improve solar field performance. In particular, it is found that an appropriate tracking-offset algorithm can improve the solar field performance by about 15%. The work here provides a valuable reference for the growing CSP industry.

14 SOLAR ENERGY↗
DOE OSTI2021

Quantifying the Solar Energy Resource for Puerto Rico

After Hurricane Maria, multiple U.S. Department of Energy laboratories studied the state of the electric grid in Puerto Rico and analyzed grid resilience and grid integration of renewable energy. As part of the work done at the National Renewable Energy Laboratory, researchers created new solar resource data, conducted a technical potential and supply curve analysis, and studied the interannual variability of the solar resource. A new methodology was developed to downscale solar resource data from the National Solar Radiation Data Base (NSRDB) from a 4-km x 4-km spatial and 30-minute temporal resolution to a 2 km x 2 km and 5-minute resolution. This methodology primarily used simple physical principles to develop high-resolution cloud properties which were then used to compute solar radiation. The high-resolution datasets were validated against ground measurements and the error metrics were found to be similar to the original lower resolution dataset. Using 20 years of downscaled data from the NSRDB multi-year capacity factors for photovoltaics (PV) were developed for both single-axis tracking and fixed latitude-tilt configurations. Use of the multi-year data provides the ability to understand variability in capacity factors due to variability in weather over a long period of time. For Puerto Rico the coastal regions were found to have significant higher capacity factors than inland. Using land-use and terrain information a technical potential analysis was conducted for Puerto Rico. This analysis restricted single PV plant development to a maximum of 100 MW nameplate capacity. The nameplate capacity for each municipality were then determined. Based on our assumptions, 56 of the 78 total municipalities of Puerto Rico contain some level of solar capacity. Most of the interior municipalities did not have any capacity because of the geographic exclusions used in this study. The lowest capacity for a PV plant observed in a municipality was 10 MW. The maximum capacity within a county was 2,000 MW. Further a supply curve analysis was conducted by taking the results of the technical potential and quantifying system and transmission costs. The levelized cost of energy (LCOE) was calculated for each theoretical PV plant site, and the levelized cost of transmission was added to the LCOE to produce a total cost estimate for each site. The results of the supply curve analysis allow for a relative comparison of the cost for integrating new PV capacity into the grid. This analysis indicates that cheaper total LCOE sites tend to be larger in capacity. The total capacity in this study was found to be far beyond the maximum peak load for the island. However, this study does not consider the economic and market potential for development. The cumulative capacity presented in this study assumes that the best locations are developed first and ignores the complex decision paths for new power plant development. Therefore, this analysis can only be treated as illustrative. Finally this study investigates the impact of inter-annual variability of resource using a variety of metrices including probability of exceedance and variation in capacity factor and LCOE. This study demonstrates that the capacity factor or LCOE could vary by over 10% year to year. This clearly indicates the risks involved in using any particular year of data and clearly points to the use of multi-year data to reduce some of the risks related to variability in weather.

14 SOLAR ENERGY↗
DOE OSTI2020

An Optimal Power Control Strategy for Grid-Following Inverters in a Synchronous Frame

This work proposes a power control strategy based on the linear quadratic regulator with optimal reference tracking (LQR-ORT) for a three-phase inverter-based generator (IBG) using an LCL filter. The use of an LQR-ORT controller increases robustness margins and reduces the quadratic value of the power error and control inputs during transient response. A model in a synchronous reference frame that integrates power sharing and voltage–current (V–I) dynamics is also proposed. This model allows for analyzing closed-loop eigenvalue location and robustness margins. The proposed controller was compared against a classical droop approach using proportional-resonant controllers for the inner loops. Mathematical analysis and hardware-in-the-loop (HIL) experiments under variations in the LCL filter components demonstrate fulfillment of robustness and performance bounds of the LQR-ORT controller. Experimental results demonstrate accuracy of the proposed model and the effectiveness of the LQR-ORT controller in improving transient response, robustness, and power decoupling.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS1980

Calibration validation for the GEOS-3 altimeter

The absolute bias calibration for the GEOS-3 intensive mode altimeter was measured using two satellite passes whose groundtracks were within 1 km of the Bermuda laser station. The Bermuda laser tracked on the two passes, and was supported by two other NASA lasers on one pass and by the NASA Spacecraft Tracking and Data Network on the other pass. For each pass, the altimeter data around Bermuda was smoothed and extrapolated to the point closest to overhead at the laser site. After correcting for tide heights and sea state effects, the two passes give calibration biases which are in agreement to within 26 cm and have a weighted mean of -5.69 + or - 0.16m for correcting altimeter measurements to the center-of-mass of the spacecraft (i.e., including the antenna tracking point correction). It was found impossible to reconcile the two calibration passes, as well as a set of altimeter crossovers in the middle of the GEOS-3 calibration area, without allowing for a data time tag error. On the bias of a selected set of four crossovers, and an assessment of probable sources of timing error, it was concluded that one interpulse period (10.24 msec) should be added to the data time tags.

Martin, C. F.↗
NASA NTRS1984

Satellite tracking and Earth dynamics research programs

Following an upgrading program, ranging performance capabilities of a satellite-tracking pulsed laser system were assessed in terms of range accuracy, range noise, data yield, and reliability. With a shorter laser pulse duration (2.5 to 3.0 NSEC) and a new analog pulse processing system, the systematic range errors were reduced to 3 to 5 cm and range noise was reduced to 5 to 16 cm and range noise was reduced to 5 to 15 cm on Starlette and BE-C, and 10 to 18 cm on LAGEOS. Maximum pulse repetition rate was increased to 30 pulses per minute and significant improvement was made in signal to noise ratio by installing a 3 A interference filter and by reducing the range gate window to 200 to 400 nsec. The solution to a problem involving leakage of a fraction of the laser oscillator pulse through the pulse chopper was outlined.

Pearlman, M. R.↗
NASA NTRS1993

Sidereal-Rate Generator

Innovative circuit derives sidereal seconds from precise 60-Hz source. Sidereal-rate generator puts out sidereal clock signals with error of only about 1.8 seconds per year. Output pulse periods uniform, solar and sidereal outputs available from single source, circuit compatible with normal digital clocks, and circuit driven by readily available reference frequency. Applications in stellar tracking, celestial navigation, and celestial photography.

Franke, John M.↗
NASA NTRS1994

Mean sea surface and geoid gradient comparisons with TOPEX altimeter data

Cycles 4 to 54 of TOPEX data have been analyzed through comparisons with the mean sea surface given on the disturbed geophysical data record (GDR). Two inverted barometer correction procedures were considered for the data reduction. One used a constant atmospheric pressure for all data while the one adopted for use, for most computations, introduced a cycle average pressure. The maximum difference between the two estimates was 3.0 cm with a clear annual signal. With the modified correction the TOPEX sea surface was compared to The Ohio State University (OSU) mean sea surface, given on the GDR, to estimate three translations ( delta x = -2.3 cm; delta y = 25.0 cm; delta z = -0.3 cm) and a bias (43.3 cm) between the two surfaces. The only significant translation is delta y which indicates the reference frame of the TOPEX system differs from that used in the OSU mean sea surface system. The bias between the TOPEX mean sea surface and the OSU mean sea surface was used to estimate an equatorial radius of 6,378,136.55 m based on an 18-cm biased estimate of the TOPEX altimeter. Examination of the average difference, by cycle, between the TOPEX sea surface and the OSU mean sea surface suggested a bias change of 3.1 +/- 2.2 mm/yr with a positive sign indicating the average ocean surface is rising or the altimeter measured distance is decreasing. Models were implemented that solved directly for a bias, bias rate annual/semiannual, and tide correction terms. The computations indicated that a simultaneous solution for this bias, bias rate, and annual/semiannual terms gave the most accurate results. Nonsimultaneous solutions led to slightly different bias rate values. The root mean square difference between the TOPEX sea surface and OSU sea surface, after translation and bias correction, was +/- 17 cm for a typical cycle. Some locations were indentified where the difference could reach 2.3 cm and were repeated over several cycles indicating errors in the mean sea surface. Most of the large differences occur in regions lacking altimeter data prior to the TOPEX/POSEIDON mission and/or areas of significant bathymetric signature. Geoid gradients are needed for the reduction of individual track data to a reference track. The accuracy of the determination of such gradients was determined through the comparison of predicted along-track gradients to the observed gradients. Among four mean sea surfaces tested the best agreement was found with the OSU mean sea surface placed on the TOPEX geophysical data record.

Rapp, Richard H.↗
NASA NTRS2005

PlanWorks: A Debugging Environment for Constraint Based Planning Systems

Numerous planning and scheduling systems employ underlying constraint reasoning systems. Debugging such systems involves the search for errors in model rules, constraint reasoning algorithms, search heuristics, and the problem instance (initial state and goals). In order to effectively find such problems, users must see why each state or action is in a plan by tracking causal chains back to part of the initial problem instance. They must be able to visualize complex relationships among many different entities and distinguish between those entities easily. For example, a variable can be in the scope of several constraints, as well as part of a state or activity in a plan; the activity can arise as a consequence of another activity and a model rule. Finally, they must be able to track each logical inference made during planning. We have developed PlanWorks, a comprehensive system for debugging constraint-based planning and scheduling systems. PlanWorks assumes a strong transaction model of the entire planning process, including adding and removing parts of the constraint network, variable assignment, and constraint propagation. A planner logs all transactions to a relational database that is tailored to support queries for of specialized views to display different forms of data (e.g. constraints, activities, resources, and causal links). PlanWorks was specifically developed for the Extensible Universal Remote Operations Planning Architecture (EUROPA(sub 2)) developed at NASA, but the underlying principles behind PlanWorks make it useful for many constraint-based planning systems. The paper is organized as follows. We first describe some fundamentals of EUROPA(sub 2). We then describe PlanWorks' principal components. We then discuss each component in detail, and then describe inter-component navigation features. We close with a discussion of how PlanWorks is used to find model flaws.

Daley, Patrick↗
NASA NTRS2011

Achieving Exact and Constant Turnaround Ratio in a DDS-Based Coherent Transponder

A report describes a non-standard direct digital synthesizer (DDS) implementation that can be used as part of a coherent transponder so as to allow any rational turnaround ratio to be exactly achieved and maintained while the received frequency varies. (A coherent transponder is a receiver-transmitter in which the transmitted carrier is locked to a pre-determined multiple of the received carrier's frequency and phase. That multiple is called the turnaround ratio.) The report also describes a general model for coherent transponders that are partly digital. A partially digital transponder is one in which analog signal processing is used to convert the signals between high frequencies at which they are radiated and relatively low frequencies at which they are converted to or from digital form, with most of the complex processing performed digitally. There is a variety of possible architectures for such a transponder, and different ones can be selected by choosing different parameter values in the general model. Such a transponder uses a DDS to create a low-frequency quasi-sinusoidal signal that tracks the received carrier s phase, and another DDS to generate an IF or near-baseband version of the transmitted carrier. With conventional DDS implementations, a given turnaround ratio can be achieved only approximately, and the error varies slightly as the received frequency changes. The non-conventional implementation employed here allows any rational turnaround ratio to be exactly maintained.

D'Addario, Larry R.↗
NASA NTRS1993

Ionospheric calibration for single frequency altimeter measurements

This report investigates the potential of using Global Positioning System (GPS) data and a model of the ionosphere to supply a measure of the sub-satellite Total Electron Current (TEC) of the required accuracy (10 TECU rms) for the purpose of calibrating single frequency radar altimeter measurements. Since climatological (monthly mean) models are known to be in error by as much as 50 percent, this work focused on the Parameterized Real-Time Ionospheric Specification Model (PRISM) which has the capability to improve model accuracy by ingesting (adjusting to) in situ ionospheric measurements. A set of globally distributed TEC measurements were generated using GPS data and were used as input to improve the accuracy of the PRISM model. The adjusted PRISM TEC values were compared to TOPEX dual frequency TEC measurements (which are considered truth) for a number of TOPEX sub-satellite tracks. The adjusted PRISM values generally compared to the TOPEX measurements within the 10 TECU accuracy requirements when the sub-satellite track passed within 300 to 400 km of the GPS TEC data or when the track passed through a night time ionosphere. However, when the sub-satellite points were greater than 300 to 400 km away from the GPS TEC data or when a local noon ionosphere was sampled, the adjusted PRISM values generally differed by greater than 10 TECU rms with data excursions from the TOPEX TEC measurements of as much as 40 TECU (an 8 cm path delay error at K band). Therefore, it can be concluded from this analysis that an unrealistically large number of GPS stations would be needed to predict sub-satellite TEC at the 10 TECU level in the day time ionosphere using a model such as PRISM. However, a technique currently being studied at the Jet Propulsion Laboratory (JPL) may provide a means of supplying adequate TEC data to meet the 10 TECU ionospheric correction accuracy when using a realistic number of ionospheric stations. This method involves using global GPS TEC data to estimate a global grid of vertical ionospheric TEC as a function of time (i.e. every one half hour) in a sun-fixed longitude frame. Working in a sun-fixed longitude frame, one is not limited by the spatial decorrelation distance of the ionosphere, but instead is limited more by the temporal correlations of the ionosphere in the sun-fixed frame which are a smaller effect. It is the opinion of the authors that using the global sun-fixed TEC grid data, in particular, ingesting it into PRISM, offers the best possibility of meeting the 10 TECU ionospheric correction accuracy requirement, and should be the subject of further study.

Schreiner, William S.↗
NASA NTRS1986

Cloud-top meridional momentum transports on Saturn and Jupiter

Cloud-tracked wind measurements reported by Sromovsky et al. were analyzed to determine meridional momentum transports in Saturn's northern middle latitudes. Results are expressed in terms of eastward and northward velocity components (u and v), and eddy components u and v. At most latitudes between 13 and 44 deg N (planetocentric), the transport by the mean flow ( ) is measurably southward, tending to support Saturn's large equatorial jet, and completely dominating the eddy transport. Meridional velocities are near zero at the peak of the relatively weak westward jet; along the flanks of that jet, measurements indicate divergent flow out of the jet. In this region the dominant eddy transport ( ) is northward on the north side of the jet, but not resolvable on the south side. Eddy transports at most other latitudes are not significantly different from measurement error. The conversion of eddy kinetic energy to mean kinetic energy, indicated by the correlation between and d /dy (where y is meridional distance) is clearly smaller than various values reported for Jupiter, and not significantly different from zero. Both Jovian and Saturnian results may be biased by the tendency for cloud tracking to favor high contrast features, and thus may not be entirely representative of the cloud level motions as a whole.

Stromovsky, L. A.↗
NASA NTRS1992

Orbiter-orbiter and orbiter-lander tracking using same-beam interferometry

Two spacecraft orbiting Mars will subtend a small angle as viewed from Earth. This angle will usually be smaller than the beam width of a single radio antenna. Thus the two spacecraft may be tracked simultaneously by a single Earth-based antenna. The same-beam interferometry (SBI) technique involves using two widely separated antennas, each observing the two spacecraft, to produce a measurement of the angular separation of the two spacecraft in the plane of the sky. The information content of SBI data is thus complementary to the line-of-sight information provided by conventional Doppler data. The inclusion of SBI data with the Doppler data in a joint orbit estimation procedure can desensitize the solution to gravity mismodeling and result in improved orbit determination accuracy. This article presents an overview of the SBI technique, a measurement error analysis, and an error covariance analysis of some examples of the application of SBI to orbit determination. For hypothetical scenarios involving the Mars Observer and the Russian Mars '94 spacecraft, orbit determination accuracy improvements of up to an order of magnitude are predicted, relative to the accuracy that can be obtained by using only Doppler data acquired separately from each spacecraft. Relative tracking between a Mars orbiter and a lander fixed on the surface of Mars is also studied. Results indicate that the lander location may be determined to a few meters, while the orbiter ephemeris may be determined with accuracy similar to the orbiter-orbiter case.

Folkner, W. M.↗
DOE OSTI2022

Thermodynamically Driven Synthetic Optimization for Cation-Disordered Rock Salt Cathodes

Relating the synthesis conditions of materials to their functional performance has long been an experience-based trial-and-error process. However, this methodology is not always efficient in identifying an appropriate protocol and can lead to overlooked opportunities for the performance optimization of materials through simple modifications of the synthesis process. In this work, we systematically track the structural evolution in the synthesis of a representative disordered rocksalt (a promising next-generation Li-ion cathode material) at the scale of both the long-range crystal structure and the short-range atomic structure using various in situ and ex situ techniques, including transmission electron microscopy, X-ray diffraction, and pair distribution function analysis. Further, an optimization strategy is proposed for the synthesis protocol, leading to a remarkably enhanced capacity (specific energy) of 313 mAh g –1 (987 Wh kg –1 ) at a low rate (20 mA g –1 ), with a capacity of more than 140 mAh g –1 retained even at a very high cycling rate of 2,000 mA g –1 . This strategy is further rationalized using ab initio calculations, and important opportunities for synthetic optimization demonstrated in this study are highlighted.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Using CUBIT to Create Unstructured Mesh Models for MCNP Simulations

The Monte Carlo N-Particle (MCNP) transport code version 6 (also known as MCNP6) has the capability for tracking particles on unstructured mesh (UM) geometry models embedded into constructive solid geometry (CSG) cells. This feature has been developed for performing calculations of complex geometry models because creating CSG models is a time-consuming and error-prone process as the complexities of geometries increase. An MCNP UM calculation requires UM geometry input files. The UM capability was originally designed to work with UM models created with the Abaqus/CAE software and ASCII input files that it generates. The Abaqus-formatted input files needed for MCNP UM calculations must have the correct Abaqus syntax and meet the additional MCNP requirements. Several software packages can generate a UM model formatted as an Abaqus input file. Cubit, the Sandia National Laboratory automated mesh generation toolkit, can generate a UM model formatted as an Abaqus input file. However, the Abaqus input files created by Cubit cannot be used for MCNP simulations. A Python script has been developed to convert an Abaqus file created by Cubit to an Abaqus file that the MCNP code can process. This report describes the process of using Cubit to create UM models for MCNP calculations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Generating MCNP Input Files for Unstructured Mesh Geometries

The Los Alamos National Laboratory’s (LANL) Monte Carlo N-Particle (MCNP)1 transport code version 6.3 (also known as MCNP6.3) has the capability for tracking particles on unstructured mesh (UM) geometry models embedded into constructive solid geometry (CSG) cells. This feature has been developed for performing calculations of complex geometry models because manually creating CSG models is time-consuming and error-prone as the complexities of geometries increase. A UM geometry model is a collection of finite elements representing a solid geometry. The first step of the MCNP UM calculation is using other software packages to create a finite element mesh representation of a solid 3D geometry because the MCNP code cannot be used to generate a UM model. Computer-aided design (CAD) software is typically used to create a solid geometry model, which is later imported into mesh generation software to create a UM model. Some mesh generation software packages may also be used to create solid geometries and thus CAD files are not needed. The MCNP UM feature was originally designed for models generated by the Abaqus/CAE software suite. The MCNP code version 6.0 and later can process UM models formatted as Abaqus input files. Starting with a 6.3 version, the MCNP code can process HDF5 mesh input files. We only focus on the UM models formatted as Abaqus input files in this report since currently no external software can be used to generate HDF5 mesh input files for MCNP UM calculations. The MCNP code version 6.3 can be used to convert the Abaqus mesh input files into the HDF5 mesh input files, but this option is typically used by the MCNP code development team to test the HDF5 mesh input file feature. Several software packages (such as Abaqus, Attila4MC, or Cubit) can be used to create the Abaqus input files for MCNP UM calculations. An MCNP UM calculation using an Abaqus model requires two input file types: MCNP and Abaqus input files. The Abaqus input files needed for MCNP UM calcu lations must have the correct Abaqus syntax and meet the additional requirements by the MCNP code. The MCNP code can process only Abaqus input files that make use of part and assembly definitions, where elements in each part must be grouped into one or more element sets (i.e., elset) using *Elset keyword lines with specified naming formats. The MCNP and Abaqus input files required for MCNP UM simulations must be related; pseudo-cells in an MCNP input file must be constructed from mesh model data from an Abaqus input file. For large complex UM models, it is tedious to manually create MCNP UM input files. The um pre op (unstructured mesh pre operations) program with the -m option can be used to create a skeleton MCNP input file from an Abaqus input file [6]. Since the um pre op program was written in Fortran and was not written for optimized performance, this program is a deprecated feature in the MCNP code version 6.3 and may be removed in the next release of the code. To improve calculation flow of multiphysics calculations, a Python3 code called write mcnp um input has been developed to generate an MCNP input file instead of using the um_pre_op -m option. This Python code was initially released to the public in 2020. We have updated this Python code for MCNP6.3 and it was used to generate the MCNP input files used to verify the MCNP6.3 code. The write_mcnp_um_input code is included with the MCNP6.3 code package which will be released to the public through the Radiation Safety Information Computational Center (RSICC) at Oak Ridge National Laboratory. This report is a revision of LA-UR-20-27139 report.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

FY25 MOOSE Usability Improvements: 3D Meshing Capabilities, Initiation of Geometry Support for Monte Carlo Tools, and Enhancement of MOOSE/Workbench User Input Interactions

Usability improvements have been made to MOOSE and Workbench in FY25 to enhance usability and user workflows. Assorted enhancement have been made to MOOSE’s intrinsic meshing capabilities in order to enable more flexible and complex meshing of nuclear reactor systems, in particular for 3D applications. Mesh generators have been added to perform operations such as batch mesh generation, surface mesh generation, and creation of 3D transition layers. These mesh generation capabilities make it much easier to generate high quality non-extruded 3D meshes. Additionally, work to integrate Monte Carlo reactor physics simulations into MOOSE-based multi-physics workflows has reached another milestone with the implementation of the Constructive Solid Geometry (CSG) base framework. This framework lays the foundation for mesh generators to offer the user a generic CSG output option (as opposed to a finite element mesh). To support users, workshop on the MOOSE Reactor Module was delivered which featured hands-on examples using the NEAMS Workbench on INL’s High Performance Computing system. Recent updates to the NEAMS Workbench, WASP, and the MOOSE language server have introduced several improvements aimed at making MOOSE-based simulation setup and input management faster, more accurate, and easier to use. Key capabilities that have been added include multi-tab-stop autocompletion, visual input diagnostics, developer-directed data visualizations, upgraded ParaView integration, and Workspace-level file tracking. Together, these changes make it easier for users to build, validate, and manage complex MOOSE-based simulation models — especially those involving reusable components, included files, and datasets. The improvements are designed to save time, reduce input errors, and help users get to a successful simulation run faster, with more confidence in the results.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗