SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Prior to this update, ORION was capable of processing cases and data that spanned a single UTM zone and assumed that all spatial inputs were given in units of meters. During the 2024 SMART bridge funding period, we completed upgrades to the code that permit general spatial projections, basin-scale datasets, and units.
Abstract. Robust projections and predictions of climate variability and change, particularly at regional scales, rely on the driving processes being represented with fidelity in model simulations. Consequently, the role of enhanced horizontal resolution in improved process representation in all components of the climate system continues to be of great interest. Recent simulations suggest the possibility of significant changes in both large-scale aspects of the ocean and atmospheric circulations and in the regional responses to climate change, as well as improvements in representations of small-scale processes and extremes, when resolution is enhanced. The first phase of the High-Resolution Model Intercomparison Project (HighResMIP1) was successful at producing a baseline multi-model assessment of global simulations with model grid spacings of 25–50 km in the atmosphere and 10–25 km in the ocean, a significant increase when compared to models with standard resolutions on the order of 1° that are typically used as part of the Coupled Model Intercomparison Project (CMIP) experiments. In addition to over 250 peer-reviewed manuscripts using the published HighResMIP1 datasets, the results were widely cited in the Intergovernmental Panel on Climate Change report and were the basis of a variety of derived datasets, including tracked cyclones (both tropical and extratropical), river discharge, storm surge, and impact studies. There were also suggestions from the few ocean eddy-rich coupled simulations that aspects of climate variability and change might be significantly influenced by improved process representation in such models. The compromises that HighResMIP1 made should now be revisited, given the recent major advances in modelling and computing resources. Aspects that will be reconsidered include experimental design and simulation length, complexity, and resolution. In addition, larger ensemble sizes and a wider range of future scenarios would enhance the applicability of HighResMIP. Therefore, we propose the High-Resolution Model Intercomparison Project phase 2 (HighResMIP2) to improve and extend the previous work, to address new science questions, and to further advance our understanding of the role of horizontal resolution (and hence process representation) in state-of-the-art climate simulations. With further increases in high-performance computing resources and modelling advances, along with the ability to take full advantage of these computational resources, an enhanced investigation of the drivers and consequences of variability and change in both large- and synoptic-scale weather and climate is now possible. With the arrival of global cloud-resolving models (currently run for relatively short timescales), there is also an opportunity to improve links between such models and more traditional CMIP models, with HighResMIP providing a bridge to link understanding between these domains. HighResMIP also aims to link to other CMIP projects and international efforts such as the World Climate Research Program lighthouse activities and various digital twin initiatives. It also has the potential to be used as training and validation data for the fast-evolving machine learning climate models.
This collaborative project leverages advanced characterization tools, catalyst treatment and electrochemical testing capabilities at LBNL to help commercialize a bimetallic nanosegrated catalysts for oxygen evolution reaction in proton-exchange membrane water electrolyzers. This project also has three industrial partners who are involved in catalyst scale-up synthesis, large electrode fabrication and large electrode electrochemical testing. This project investigates a broad range of topics including electrocatalysis, ink rheology and electrode fabrication and testing. This knowledge could help bridge the performance and durability gap from laboratory-developed materials (such as electrocatalysts) to applications in commercial devices. The research results also bring fundamental understanding of possible degradation mechanisms that occur under extremely oxidative potentials including material and morphological changes. The project will help accelerate the deployment of proton-exchange membrane water electrolyzers. The public benefits could include creating jobs and increasing public awareness of potential clean technologies for green hydrogen production. Ultimately, the project will help national wide effort of flighting climate change and achieving carbon neutrality by 2050.
The need for dedicated, accurate, expertly curated weather data is increasingly important as the share of variable renewable energy increases on the power system. Projections for futures with very high (50+% annual energy) shares of variable generation require ongoing assessment of data requirements from industry stakeholders in their power system operation and planning contexts. In March 2024, NREL organized a workshop entitled "Bridging the Gap Between Atmospheric Science and Grid Integration Workshop", which brought atmospheric scientists and power system experts together to refine the requirements of atmospheric datasets for grid integration, and to describe a holistic approach to creating new and regularly updated national scale wind datasets for power system planning and operations. The results of this workshop are being used to inform the near-term development and a longer-term strategy for DOE to produce relevant wind resource datasets and inform wider use of wind/solar/load data sets in power system planning. This presentation provides an overview of a preworkshop survey, an assessment of current state of the art of national-scale datasets for wind resource assessment and grid integration, insights on appropriate uses of the WTK-LED, power system perspectives on data needs, as well as recommended next steps as discussed in the workshop and how these steps support longer-term strategies.
Electrically exploded wires find uses throughout high-energy physics. For example, they are commonly used as high-temperature sources, X-ray generators, and in precision timing detonators. However, the detailed and complete physics that occurs is complex and still poorly understood. A full mechanistic description of these complex phenomena is beyond the scope of a single paper. Instead, we focus on the formation of metal vapor and its transition to plasma. This single transition is commonly assumed to comprise “bridge-burst”. We use a suite of diagnostics including a novel, fiber-based, high-speed, optical pyrometer to better characterize this transition. The primary finding from this project is that peak light output from an exploding wire does not temporally match the peak temperature. Additionally, it is found that peak light does not align with peak bridge-burst voltage and that the peak temperature is not voltage-dependent. These findings are non-intuitive and will allow for the correction of false assumptions previously made about this topic.
Generation 3 (Gen) concentrating solar power (CSP) systems requires a high operation temperature ≥750°C to increase the power-cycle efficiency towards ≥50%. However, such a high operation temperature poses a notable challenge to high temperature materials. On the receiver side, a critical challenge is the lack of solar selective absorbers that demonstrate both high optical-to-thermal energy conversion efficiency η therm approaching 95% and long-term thermal stability at >750°C in air. Existing solar absorbers either have limited η therm ≤89% or deteriorate significantly within 500 h at 750°C; some of these also require costly vacuum deposition for stringent thickness control. Therefore, it is highly desirable to simultaneously achieve η therm ~95% AND high thermal stability at 750°C for future generations of CSP receivers. This project has investigated and developed low-cost, highly scalable spray-coated transition metal oxide nanoparticle (NP) pigmented solar selective coatings on various types of Inconel alloy tube sections that are thermodynamically stable at 750-800°C in air, maintaining η therm >94.3%(93.2%) under a solar concentration ratio of C=1000 after 60 simulated day-night cycles at 750ºC (800ºC) (1 cycle=12 h at 750 or 800°C and 12h ramping down to 25°C). We have also coated up to 48 inches long receiver tubes for preliminary solar testing under Norwich Technology’s parabolic trough systems, demonstrating notably improved performance in solar heating compared to benchmark Pyromark 2500 coatings. Two key innovations have been developed in this project: (1) We achieved an unprecedented high thermal efficiency >94% by optimizing the d-band optical absorption spectra of transition metal ions, engineering their valences and stoichiometry; (2) We were able to maintain or even slightly increase the efficiency when operating at 750°C in air by engineering the interdiffusion of transition metal ions between the coating and the Inconel substrate to our advantage. Featuring high-temperature solar spectral selectivity and thermodynamic stability in air, as well as low-cost, highly scalable solution-chemical synthesis and spray coating techniques, this innovation in solar selective absorber technology alone accounts for nearly 40% of the targeted reduction in the levelized cost of electricity (LCOE) from the receiver, thermal energy storage, and operation & maintenance combined by 2030, as proposed by the U.S. Department of Energy. For a 110 MWe CSP power plant, this LCOE reduction from the solar selective coating alone transfers to ~$\$ $1.9 M increase in annual profit based on a sales price of 10¢/kWh (or an annual sale of $50 M/year). The high solar absorptance (~98%) NP pigment materials developed in this project are equally applicable to volumetric receivers either as a coating or as bigger ceramic microspheres after sintering, potentially benefiting Gen3 falling particle CSP technologies. All these contributions facilitate the larger scale deployment of CSP systems with energy storage capability to address the intermittency issue of solar energy towards dispatchable solar electricity, bridging the temporal gap between peak solar electricity production and peak electricity consumption. The collaboration with Norwich Technologies and Brayton Energy in this project will also facilitate the future commercialization of this new solar selective coating technology developed in this SIPS project.
Thermal transport in complex solids is governed by local structure, defects, and anisotropy, yet most continuum models still rely on oversimplified and homogenized conductivities. Here, we bridge atomistic and continuum descriptions by building finite element (FE) models directly from the site-projected thermal conductivity (SPTC), an atomic-level decomposition of the Green–Kubo thermal conductivity. We introduce a toolkit, the “Simulator Collection for Atomic-to-Continuum Scales (SCACS)”, which uses a graph neural network to predict SPTC on large atomic structures, coarse-grains these fields into anisotropic conductivity tensors, and embeds them into the heat-flow FE equation with a customized, anisotropy-aware adaptive mesh refinement scheme. Applied to silicon nanostructures, the resulting FE models act as representative volume elements, reproduce bulk conductivities, and capture interfacial and defect-driven anisotropy while maintaining thermodynamic consistency. Additionally, SCACS predicts experimental conductance trends and fields. This work demonstrates a general route for transferring atomistic transport information into device-scale thermal simulations with physics-based approximations.
This dataset contains aerodynamic quantities - including flow field values (momentum, energy, and vorticity) and summary values (coefficients of lift, drag, and momentum) - for 1,830 airfoil shapes computed using the HAM2D CFD (computational fluid dynamics) model. The airfoil shapes were designed using the separable shape tensor parameterization that encodes two-dimensional shapes as elements of the Grassmann manifold. This data-driven approach learns two independent spaces of parameter from a collection of sample airfoils. The first captures large-scale, linear perturbations, and the second defines small-scale, higher-order perturbations. For this dataset, we used the G2Aero database of over 19,000 airfoil shapes to learn a parameter space that captured a wide array of shape characteristics. We sampled airfoil designs over both parameter spaces to explore the full range of possible shape variations. The aerodynamic quantities for the generated airfoil were obtained using the HAM2D code, which is a finite-volume Reynolds-averaged Navier-Stokes (RANS) flow solver. We employ a fifth-order WENO scheme for spatial reconstruction with Roe's flux difference scheme for inviscid flux and second-order central differencing for viscous flux. A preconditioned GMRES method is applied for implicit integration. The Spalart-Allmaras 1-eq turbulence model is used for the turbulence closure, and the Medida-Baeder 2-eq transition model is applied to account for the effects of laminar turbulent transition. The airfoil grid is generated with a total of 400 points on the airfoil surface, the initial wall-normal spacing of y+ = 1, and an outer boundary located at 300 chord lengths away from the wall. The CFD simulations are performed at a freestream Mach number of 0.1, for or three different Reynolds' numbers (3M, 6M, and 9M), and for 25 angles of attack from -4 deg. to 20 deg. with 1 degree increments. Across all these various parameters, this dataset includes the results from over 250,000 CFD simulations. The simulations were performed using the Bridges-2 system at the Pittsburgh Supercomputing Center in February 2023 as part of the INTEGRATE project funded by the Advanced Research Projects Agency - Energy, in the U.S. Department of Energy. The data was collected, reformatted, and preprocessed for this OEDI submission in July 2023 under the Foundational AI for Wind Energy project funded by the U.S. Department of Energy Wind Energy Technologies Office. This dataset is intended to serve as a benchmark against which new artificial intelligence (AI) or machine learning (ML) tools may be tested. Baseline AI/ML methods for analyzing this dataset have been implemented, and a link to their repository containing those models has been provided. The .h5 data file structure can be found in the GitHub Repository resource under explore_airfoil_2k_data.ipynb.
This dataset contains aerodynamic quantities - including flow field values (momentum, energy, and vorticity) and summary values (coefficients of lift, drag, and momentum) - for 8,996 airfoil shapes, computed using the HAM2D CFD (computational fluid dynamics) model. The airfoil shapes were designed using the separable shape tensor parameterization that encodes two-dimensional shapes as elements of the Grassmann manifold. This data-driven approach learns two independent spaces of parameter from a collection of sample airfoils. The first captures large-scale, linear perturbations, and the second defines small-scale, higher-order perturbations. For this data, we used the G2Aero database of over 19,000 airfoil shapes to learn a parameter space that captured a wide array of shape characteristics. We fixed the linear deformations to be the mean over the database and sampled new shapes over a four-dimensional parameter space of higher-order perturbation. This sampling approaches allows for isolated analysis of non-linear airfoil shape deformations while holding other aspects (e.g., airfoil thickness) approximately constant. The aerodynamic quantities for the generated airfoil were obtained using the HAM2D code, which is a finite-volume Reynolds-averaged Navier-Stokes (RANS) flow solver. We employ a fifth-order WENO scheme for spatial reconstruction with Roe's flux difference scheme for inviscid flux and second-order central differencing for viscous flux. A preconditioned GMRES method is applied for implicit integration. The Spalart-Allmaras 1-eq turbulence model is used for the turbulence closure, and the Medida-Baeder 2-eq transition model is applied to account for the effects of laminar turbulent transition. The airfoil grid is generated with a total of 400 points on the airfoil surface, the initial wall-normal spacing of y+ = 1, and an outer boundary located at 300 chord lengths away from the wall. The CFD simulations are performed at a freestream Mach number of 0.1, Reynolds number of 9M, and at two angles of attack, 4 deg. and 12 deg. The simulations were performed using the Bridges-2 system at the Pittsburgh Supercomputing Center in February 2023 as part of the INTEGRATE project funded by the Advanced Research Projects Agency - Energy in the U.S. Department of Energy. The data was collected, reformatted, and preprocessed for this OEDI submission in July 2023 under the Foundational AI for Wind Energy project funded by the U.S. Department of Energy Wind Energy Technologies Office. This dataset is intended to serve as a benchmark against which new artificial intelligence (AI) or machine learning (ML) tools may be tested. Baseline AI/ML methods for analyzing this dataset have been implemented, and a link to their repository containing those models has been provided. The .h5 data file structure can be found in the GitHub Repository resource under explore_airfoil_9k_data.ipynb.
The Rucio Data Management System [1] is an important tool used by High Energy Physics experiments, including those at Fermi National Accelerator Laboratory, to store and manage exabyte-scale scientific datasets. Despite its central role in coordinating data across globally distributed storage sites, Rucio's command line interface (CLI) presents a steep learning curve, and makes it difficult for scientists to navigate through. To solve this issue, a containerized Model Context Protocol (MCP) [2] server was built that connects Large Language Models directly to Rucio, allowing AI agents to handle data tasks by using simple, natural language rather than memorized terminal commands. The core engineering focus of this project was moving the server away from slow terminal commands that require text parsing and replacing them with a native Python Client API toolset and a planned REST API framework. Moving to the Python API handles data operations directly in memory, which helps clear up formatting errors, provides the AI with clean, structured JSON data and speeds up tool execution. To prove that the system actually works, a benchmarking pipeline was also built with various questions to test the AI across four different model configurations. The questions included finding data scopes, tracking down specific datasets, and checking replication rules. Through benchmarking, early runs showed that with raw terminal text, the model would get confused and stuck, whereas switching to the Python API to feed the AI clean, structured data yielded massive improvement. By creating an intelligent and autonomous bridge to a storage network, this project shows how AI can be implemented in scientific data management, which ultimately helps scientists at Fermilab spend less time sorting through data and more time focusing on their experiments and analysis.
PPL Electric Utilities’ Keystone Solar Future Project (KSFP) is part of the U.S. Department of Energy’s Office of Energy Efficiency and Renewables, Energy Enabling Extreme Real-time Grid Integration of Solar Energy (ENERGISE) grant. The three-year project, beginning in October, 2017, funded an innovative pilot program to develop and evaluate a technology platform and framework for managing high levels of distributed energy resources (DER), and improving future DER interconnection. The KSFP aimed to bridge gaps between existing and future technology, including Distributed Energy Resource Management Systems (DERMS), by monitoring, managing, and optimizing a system with high penetration of solar and other DER. Specifically, the project worked to make Demand Management Systems (DMS) applications, devices, and real-time communications technology DER-aware to allow for advanced applications for utility operations. Here, by installing and implementing a DERMS, project designers transformed how DERs interact with PPL’s distribution system, offering learnings and results worth sharing industry-wide.
The goal of this project is to model and control a novel unidirectional cascaded multilevel bridgeless rectifier as an active front end in medium and high voltage applications. This topology has many advantages over a conventional cascaded H-bridge rectifier, such as lower implementation cost, higher reliability, and greater flexibility with similar power quality.The steady-state mathematical model is used to develop a method for the voltage balancing of dc cells. Power factor analysis is discussed to achieve unity power factor using fully controlled hbridge cells. Power loss, efficiency, and cost comparison studies between the traditional cascaded H-Bridge converter and the proposed bridgeless converter demonstrate the advantages. A novel control strategy is proposed to achieve dc voltage balancing, fast and robust grid synchronization and power factor correction under unbalanced load conditions. Simulation and experimental results validate the models and control method.
DHS’s Regional Resiliency Assessment Program (RRAP) undertook the Oregon State Transportation Systems RRAP project between 2018-2021 to assess the seismic vulnerabilities of the state’s transportation system to a Cascadia Subduction Zone (CSZ) earthquake (CISA 2021). The project primarily focused on the effects of a CSZ earthquake on the mobility of emergency response efforts into the region. In the event that bridges located at crucial river crossings fail during a CSZ earthquake, rivers will become natural barriers affecting the ability to support response and recovery efforts. The characteristics of rivers (e.g., width, depth, and flow) are critical data to inform planning for the recovery of transportation systems across river barriers and determine reopening times following a CSZ earthquake. As part of the transportation system assessment, an Argonne National Laboratory (Argonne) project team developed a River Characterization Tool (RCT) to quantify the physical characteristics of the river system in Oregon, which can be used to inform a system-level analysis of transportation recovery and reopening times. This report mainly discusses approaches, methods, development, and applications of the RCT. The use of the river data generated from this tool for a system-level analysis of transportation recovery is discussed in the Oregon Transportation Systems RRAP Resiliency Assessment report (CISA 2021). The primary goal of the RCT is to generate statewide dataset of widths, depths, and flows for all rivers across Oregon to assess the effects of river barriers on emergency response and recovery efforts following a CSZ earthquake. The RCT leverages remote sensing data that multispectral sensors collect from satellites, employs a deep learning model to characterize river predictions, and computes river channel geometric dimensions with hydraulic relationships and river predictions. The entire process of river characterization is automated and implemented through the RCT.
During fissionable material processing, all normal and credible abnormal conditions must remain safely subcritical. Nuclear Criticality Safety (NCS) uses a number of methods to determine subcriticality, one of which is the use of neutron transport codes such as MCNP6. In order to create models for use with MCNP6, both the geometry and materials in fissionable material processes must be known, or assumptions must be made and quantified for the impact to bias. One of the systems with a significant amount of bias due to material modeling assumptions is in the area of aqueous plutonium processing. These solutions are typically plutonium nitrate solutions or plutonium chloride solutions, which are modeled as fictitious plutonium metal-water mixtures because little is known about the actual density of the solution and there is no current predictive capability approved for use at Los Alamos National Laboratory (LANL) for modeling them. This research is currently underway to fill the gap and develop an algorithm for use with MCNP6 to model the density of plutonium chloride solutions. The method is to be validated with experimental data for density, and also validated with critical experiments using MCNP6. Note that the Chlorine Worth Study (CWS) was performed in December 2021 to help bridge the gap in chlorine data for critical experiments, and is currently awaiting International Criticality Safety Benchmark Evaluation Project (ICSBEP) review. This study was performed by LANL at the National Criticality Experiments Research Center (NCERC) at the Nevada National Security Site (NNSS). Additional information regarding this experiment may currently be found in LA-UR- 22-29180. Additionally, the Chemistry-Actinide Analytical Chemistry (C-AAC) at LANL has performed a number of solution density measurements for PuCl 3 -HC 1 -H 2 O, allowing for such data be used to create a semi-empirical density via the Pitzer method. The published dataset for the measurements is documented in LA-UR-22-25454. This method has already been tested successfully for aqueous plutonium nitrate solutions in SCALE. Current solution density measurements exist of plutonium concentrations of 0-~142g/L, all at 2M HC1, at temperatures 20-40°C. Additional data was taken for HC1-corrected density values, which essentially mimics the data for a pure PuCl x -water solution. The calculations in this report aim to support the current research by demonstrating the difference in system reactivity for the current modeling method when compared to the new proposed modeling with a density law implementation, which is being written as a Python tool to be used with MCNP6.
During the period August 24-28, 2016, a delayed component the GO-Amazon Project, we launched nine altitude-controlled free balloons (CMET). Smaller than typical rawinsondes, CMET are equipped with altitude control, global communication via Iridium satellite, and aspirated sensors. The aims of our effort were to examine the interactions among convective boundary layer and dual river breeze circulations near the confluence of the Tapajos and Amazon Rivers in the eastern Basin. The week-long field campaign was timed to examine the reestablishment of the breeze circulations shortly after the passage of a strong instability line on August 22. Nine CMET were launched at the Curua-Una hydroelectric dam (2.8S; 54.3W), timed to encounter the Tapajos river breeze front by late afternoon. Soundings were made to establish the thickness of interface between the easterly trade and westerly Tapajos breeze circulation. Careful use of sounding strategies allowed these free balloons to track along the northerly channeled flow in the lowest 300 m above the River. Following the river encounter, balloons tracked to the west, sounding to describe the diurnal course of boundary layer in the forest west of the Tapajos River. The longest flight traveled more than 770 km over three days and twice rested overnight in the rain forest canopy. Ancillary data from surface climate and flux stations as well as the Santarem radiosonde, satellite images will be used to illustrate how the breeze circulations are seen near the surface and how they were disrupted by larger-scale events. Comparisons with HYSPLIT trajectories illustrate how sensitive real trajectories are to the refraction that the encounter with the breeze effects. A second effort was intensive analyses of data from the eddy flux tower at the Tapajós National Forest near Santarém Pará Brazil, with special emphasis on the role of clouds and canopy presence on light absorbance and photosynthetic activity.
Abstract 2D transition metal dichalcogenides (TMDCs) with intense and tunable photoluminescence (PL) have opened up new opportunities for optoelectronic and photonic applications such as light‐emitting diodes, photodetectors, and single‐photon emitters. Among the standard characterization tools for 2D materials, Raman spectroscopy stands out as a fast and non‐destructive technique capable of probing material's crystallinity and perturbations such as doping and strain. However, a comprehensive understanding of the correlation between photoluminescence and Raman spectra in monolayer MoS 2 remains elusive due to its highly nonlinear nature. Here, the connections between PL signatures and Raman modes are systematically explored, providing comprehensive insights into the physical mechanisms correlating PL and Raman features. This study's analysis further disentangles the strain and doping contributions from the Raman spectra through machine‐learning models. First, a dense convolutional network (DenseNet) to predict PL maps by spatial Raman maps is deployed. Moreover, a gradient boosted trees model (XGBoost) with Shapley additive explanation (SHAP) to bridge the impact of individual Raman features in PL features is applied. Last, a support vector machine (SVM) to project PL features on Raman frequencies is adopted. This work may serve as a methodology for applying machine learning to characterizations of 2D materials.
Abstract Process modeling of Aerosol‐cloud interaction (ACI) is essential to bridging gaps between observational analysis and climate modeling of aerosol effects in the Earth system and eventually reducing climate projection uncertainties. In this study, we examine ACI in summertime precipitating shallow cumuli observed during the Aerosol Cloud meTeorology Interactions oVer the western ATlantic Experiment (ACTIVATE). Aerosols and precipitating shallow cumuli were extensively observed with in‐situ and remote‐sensing instruments during two research flight cases on 02 June and 07 June, respectively, during the ACTIVATE summer 2021 deployment phase. We perform observational analysis and large‐eddy simulation (LES) of aerosol effect on precipitating cumulus in these two cases. Given the measured aerosol size distributions and meteorological conditions, LES is able to reproduce the observed cloud properties by aircraft such as liquid water content (LWC), cloud droplet number concentration ( N c ) and effective radius r eff . However, it produces smaller liquid water path (LWP) and larger N c compared to the satellite retrievals. Both 02 and 07 June cases are over warm waters of the Gulf Stream and have a cloud top height over 3 km, but the 07 June case is more polluted and has larger LWC. We find that the N a ‐induced LWP adjustment is dominated by precipitation feedback for the 2 June precipitating case and there is no clear entrainment feedback in both cases. An increase of cloud fraction due to a decrease of aerosol number concentration is also shown in the simulations for the 02 June case.
Radiation imaging capabilities, essential in the nuclear nonproliferation regime, facilitate source localization and, in certain cases, spectroscopy. Scatter-based neutron cameras, which can measure the neutron signatures from special nuclear material, hold particular interest. Systems incorporating organic scintillators can extract neutron energy spectra, potentially distinguishing fission neutron sources from others, such as alpha-neutron sources. The development and testing of a scatter-based neutron imager, however, can be challenging without having an accurate simulation model or first constructing a prototype. This work describes a simulation pipeline that takes output from MCNPX-PoliMi simulations and creates the expected back-projection neutron images and neutron energy spectra. This pipeline was developed to improve the modeling of fast neutron imagers and bridge the current gap in literature, which predominantly focuses on gamma-ray Compton imager models. This work also reports on the significance of various real-world system considerations and their effects on the simulated detector responses. The pipeline was verified and validated with experimental data collected using a 252 Cf spontaneous fission source using a fast neutron scattering imager developed at the University of Michigan.