Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “flux analysis”

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,135 records · Page 63
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Implementation and evaluation of the unified stomatal optimization approach in the Functionally Assembled Terrestrial Ecosystem Simulator (FATES)

Abstract. Stomata play a central role in regulating the exchange of carbon dioxide and water vapor between ecosystems and the atmosphere. Their function is represented in land surface models (LSMs) by conductance models. The Functionally Assembled Terrestrial Ecosystem Simulator (FATES) is a dynamic vegetation demography model that can simulate both detailed plant demographic and physiological dynamics. To evaluate the effect of stomatal conductance model formulation on forest water and carbon fluxes in FATES, we implemented an optimality-based stomatal conductance model – the Medlyn (MED) model – that simulates the relationship between photosynthesis (A) and stomatal conductance to water vapor (gsw) as an alternative to the FATES default Ball–Woodrow–Berry (BWB) model. To evaluate how the behavior of FATES is affected by stomatal model choice, we conducted a model sensitivity analysis to explore the response of gsw to climate forcing, including atmospheric CO2 concentration, air temperature, radiation, and vapor pressure deficit in the air (VPDa). We found that modeled gsw values varied greatly between the BWB and MED formulations due to the different default stomatal slope parameters (g1). After harmonizing g1 and holding the stomatal intercept parameter (g0) constant for both model formulations, we found that the divergence in modeled gsw was limited to conditions when the VPDa exceeded 1.5 kPa. We then evaluated model simulation results against measurements from a wet evergreen forest in Panama. Results showed that both the MED and BWB model formulations were able to capture the magnitude and diurnal changes of measured gsw and A but underestimated both by about 30 % when the soil was predicted to be very dry. Comparison of modeled soil water content from FATES to a reanalysis product showed that FATES captured soil drying well, but translation of drying soil to modeled physiology reduced the models' ability to match observations. Our study suggests that the parameterization of stomatal conductance models and current model response to drought are the critical areas for improving model simulation of CO2 and water fluxes in tropical forests.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Modeling ExB drift transport in conceptual slot divertor configurations

At DIII-D, a slot divertor concept, called small-angle-slot (SAS), is under development, aiming to enable detachment at relatively low plasma edge density. We report on simulations using the SOLPS-ITER two-dimensional edge code to examine the performance of conceptual “SAS 2” slot configurations. The focus of the analysis is on E × B drift effects on upstream density at detachment (UDD), with detachment marked by electron temperature T e ≤ 3 eV at the outer strike point (OSP). With toroidal field such that radial E × B drift carries particles from the OSP towards the private flux region (PFR), placing the OSP near the inner slot wall gives ≈20% lower UDD than having the OSP near the outer wall. The inner wall effectively traps the radial E × B drift flux, resulting in low T e and associated radial electric field in the PFR, and thus small losses from the slot to the inner target via poloidal E × B drift flux. With toroidal field reversed such that radial E × B drift is reversed, OSP placement near the inner wall gives ≈10% lower UDD than OSP placement near the outer wall. Although radial E × B flux is from the OSP towards the outer wall, this flux largely escapes the slot, raising the UDD. Finally, a change in the slot shaping is suggested with the goal of eliminating such E × B-driven particle losses from the slot.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2023

The role of plasma–atom and molecule interactions on power & particle balance during detachment on the MAST Upgrade Super-X divertor

This paper shows first quantitative analysis of the detachment processes in the MAST Upgrade Super-X divertor (SXD). We identify an unprecedented impact of plasma-molecular interactions involving molecular ions (likely $D^+_2$), resulting in strong ion sinks (Molecular Activated Recombination—MAR), leading to a reduction of ion target flux. The MAR ion sinks exceed the divertor ion sources before electron-ion recombination (EIR) starts to occur, suggesting that significant ionisation occurs outside of the divertor chamber. In the EIR region, $T_e \ll 0.2$ eV is observed and MAR remains significant in these deep detached phases. The total ion sink strength demonstrates the capability for particle (ion) exhaust in the Super-X Configuration. Molecular Activated Dissociation is the dominant volumetric neutral atom creation process can lead to an electron cooling of 20% of $P_\mathrm{SOL}$. The measured total radiative power losses in the divertor chamber are consistent with inferred hydrogenic radiative power losses. This suggests that intrinsic divertor impurity radiation, despite the carbon walls, is minor in the divertor chamber. This contrasts previous TCV results, which may be associated with enhanced plasma-neutral interactions and reduced chemical erosion in the detached, tightly baffled SXD. The above observations have also been observed in higher heat flux (narrower SOL width) type I ELMy H-mode discharges. This provides evidence that the characterisation in this paper may be general.

plasma-molecular interactions↗
DOE OSTI2025

Development of an atomic layer deposition system for deposition of alumina as a hydrogen permeation barrier

Tritium permeation into and through materials poses a critical challenge for the development of nuclear fusion reactors. Minimizing tritium permeation is essential for the safe and efficient use of available fuel supplies. In this work, we present the design, construction, and validation of custom atomic layer deposition (ALD) and deuterium permeation measurement systems aimed at developing thin-film hydrogen permeation barriers. Using the ALD system, we deposited conformal $\mathrm{Al}_{2}\mathrm{O}_{3}$ films on copper foil substrates and characterized their growth behavior, morphology, and composition. ALD growth rates of ∼1.1 Å/cycle were achieved for temperatures between 100 ∘ C and 210 ∘ C. Permeation measurements on bare and alumina coated copper foils revealed a significant reduction in deuterium flux with the addition of a ∼10 nm $\mathrm{Al}_{2}\mathrm{O}_{3}$ layer. While bare copper followed diffusion-limited transport consistent with Sievert’s law, the alumina-coated samples exhibited surface-limited, pore-mediated transport with linear pressure dependence. Arrhenius analysis showed distinct differences in activation energy for the two transport regimes, and permeation reduction factors exceeding an order of magnitude were observed. These results demonstrate the potential of ALD-grown $\mathrm{Al}_{2}\mathrm{O}_{3}$ films as effective hydrogen isotope barriers and provide a foundation for future studies on film optimization and integration into fusion-relevant components.

atomic layer deposition↗
DOE OSTI2021

How does vertical wind shear influence entrainment in squall lines?

The influence of vertical wind shear on updraft entrainment in squall lines is not well understood. To address this knowledge gap, a suite of high-resolution idealized numerical model simulations of squall lines were run in various vertical wind shear (hereafter “shear”) environments to study the effects of shear on entrainment in deep convective updrafts. Low-level horizontal mass flux into the leading edge of the cold pool was strongest in the simulations with the strongest low-level shear. These simulations consequently displayed wider updrafts, less entrainment-driven dilution, and larger buoyancy than the simulations with comparatively weak low-level shear. An analysis of vertical accelerations along trajectories that passed through updrafts showed larger net accelerations from buoyancy in the simulations with stronger low-level shear, which demonstrates how less entrainment-driven dilution equated to stronger updrafts. Additionally, the effects of upper-level shear on entrainment and updraft vertical velocities were generally less pronounced than the effects of low-level shear. We argue that in addition to the outflow boundary-shear interactions and their effect on updraft tilt established by previous authors, decreased entrainment-driven dilution is yet another beneficial effect of strong low-level shear on squall line updraft intensity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Linking resource availability to pantropical forest canopy resistance and resilience to cyclone disturbance

Statement of purpose: Tropical cyclones are intensifying and occurring at higher latitudes in recent decades, but the mechanisms underpinning the resistance (ability to withstand disturbance-induced change) and resilience (pace of return to pre-disturbance reference values) of tropical forests to cyclones remains largely unexplored at the pantropical scale. We conducted a meta-analysis to investigate the role of soil resource availability (i.e., total soil phosphorus concentration) in mediating site-level forest canopy resistance and resilience to cyclones pan-tropically. We evaluated cyclone-induced and post-cyclone litterfall mass (g/m2/day), phosphorus (P) and nitrogen (N) fluxes (mg/m2/day), as well as concentrations (mg/g) across 73 case studies in Australia, Guadeloupe, Hawaii, Mexico, Puerto Rico, and Taiwan. The dataset zip file includes three data and two metadata files: - The compiled Litterfall Mass Flux data from tropical forests across the globe prior to and after varying tropical cyclone disturbances are provided in Litterfall_Mass.csv. This data file also includes site location, geographical characteristics, elevation, soil phosphorus concentration, geology, and several variables related to each tropical cyclone disturbance. - The compiled Litterfall Nitrogen and Phosphorus Flux data from tropical forests across the globe prior to and after varying tropical cyclone disturbances are provided in Litterfall_Nutrients.csv. This data file also includes site location, geographical characteristics, elevation, soil phosphorus concentration, geology, and several variables related to each tropical cyclone disturbance. - Tropical cyclone track data compiled from HURDAT2 and IBTrACS databases and used as input in the HURRECON model (https://github.com/hurrecon-model/HurreconR) to generate wind data is provided in hurdat2-1851-2019-052520.txt. - The metadata file (Metadata_Meta-analysis_Litterfall-Mass.pdf) has the complete information on each variable included in the Litterfall_Mass.csv dataset, the data sources, and data processing information. - The metadata file (Metadata_Meta-analysis_Litterfall-Nutrients.pdf) has the complete information on each variable included in the Litterfall_Nutrients.csv dataset, the data sources, and data processing information.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

HFIR LEU High Density Silicide Dispersion Optimized Design Neutronics Analyses with PHAME

A high-fidelity neutronics model of the Oak Ridge National Laboratory High Flux Isotope Reactor (HFIR) with the low-enriched uranium (LEU) high-density silicide dispersion Optimized fuel design was updated and analyzed to generate reactor physics-based metrics to support follow-on thermal hydraulic and transient analyses of this design. The Python HFIR Analysis and Measurement Engine (PHAME) was also updated to enhance the automation capabilities of the framework developed and maintained to perform these reactor physics modeling and simulation efforts. The automated framework significantly increases the efficiency and reproducibility to design and thoroughly analyzes HFIR LEU core designs, changes, and uncertainties. Reactor physics metrics evaluated include but are not limited to fuel depletion, cycle length, fission rate density distributions, axial power peaking factors, kinetics data, reactivity coefficients, control element worths, heat deposition rates, and decay heat. These neutronics results provide essential input to follow-on steady state thermal, thermal hydraulic and reactor transient analyses, which are subject of other reports. The Optimized design operates at 95 MW to maintain HFIR’s current highly enriched uranium core performance level at 85 MW.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

Intermittency and Dissipative Structures Arising from Relativistic Magnetized Turbulence

Abstract Kinetic simulations of relativistic turbulence have significantly advanced our understanding of turbulent particle acceleration. Recent progress has highlighted the need for an updated acceleration theory that can account for particle acceleration within the plasma’s coherent structures. Here, we investigate how intermittency modeling connects statistical fluctuations in turbulence to regions of high-energy dissipation. This connection is established by employing a generalized She–Leveque model to characterize the exponents ζ p for the structure functions S p ∝ l ζ p . The fitting of the scaling exponents provides us with a measure of the codimension of the dissipative structures, for which we subsequently determine the filling fraction. We perform our analysis for a range of magnetizations σ and relative fluctuation amplitudes δ B 0 / B 0 . We find that increasing values of σ and δ B 0 / B 0 allow the turbulent cascade to break sheetlike structures into smaller regions of dissipation that resemble chains of flux ropes. However, as their dissipation measure increases, the dissipative regions become less volume filling. With this work, we aim to inform future turbulent acceleration theories that incorporate particle energization from interactions with coherent structures within relativistic turbulence.

Astronomy & Astrophysics↗
DOE OSTI2022

Investigation of the OECD/NEA PWR MOX/UO{sub 2} core transient benchmark using a coupled whole core pin-by-pin route in WIMS

This paper investigates the new CAMELOT coupling route developed in WIMS, which uses prepared tabulated cross-sections combined with the flux solver MERLIN module and the integrated sub-channel thermal hydraulics solver ARTHUR module to solve for the coupled core state. This coupling route is applied to the OECD/NEA pressurized water reactor (PWR) core transient benchmark exercise, involving 3-D neutronics and thermal-hydraulics analysis of a full-sized PWR core featuring MOX fuel at various core states, as well as simulating a control rod ejection transient. Analysis of the results for each part of the benchmark demonstrate that the high-resolution pin-by-pin analysis used in WIMS closely matches the other participants in the benchmark, with differences of between 20-200 pcm for k-effective, 1-2.5% RMS difference for assembly-averaged powers, 2-8% difference in rod worth calculations, < 200 ppm for calculations of critical boron concentration, zero difference in delayed neutron fraction calculation and a very close match for power response to the rod ejection transient. This is the first application of the WIMS-CAMELOT approach for the transient analysis of a full-sized reactor and as such further improvements such as reducing computational cost through parallelization, optimisation and memory reduction, are currently in active development. This paper demonstrates the accuracy and flexibility of the CAMELOT coupling route, providing a straightforward and consistent user image that can be easily used for a variety of modelling problems involving coupled neutronics and thermal-hydraulics. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Search for Lightly Ionizing Particles in SuperCDMS and simulation of neutron backgrounds

The Super Cryogenic Dark Matter Search (SuperCDMS) is a direct-detection dark matter search experiment that primarily aims to search for Weakly Interacting Massive Particles (WIMPs) using state-of-the-art solid-state detection technology. During its operation at the Soudan underground laboratory, germanium detectors were operated with high bias voltage mode known as the CDMS low ionization threshold experiment (CDMSlite) to achieve below-keV thresholds. CDMSlite, for being able to measure small energy depositions in detectors, also provides sensitivity to Lightly Ionizing Particles (LIPs) with very small fractional charges. This thesis will discuss an analysis to search LIPs with the data acquired in CDMSlite mode. An important component for LIPs search is the expected energy-deposition distributions for LIPs falling on the CDMSlite detector. In this thesis, a simulation framework to calculate the energy-deposition distributions is developed. This thesis presents first direct-detection limits on the intensity of cosmogenic LIPs with electric charges smaller than e /(3 × 10 5 ) as well as the strongest limits for charges ≤ e /160, with a minimum intensity of 1.36 × 10 -7 cm -2 s -1 sr -1 at charge e /160.In any rare-event search experiment, understanding background is crucial. Neutrons capable of mimicking dark matter signals are a major background for any dark matter search experiment. A simulation study to estimate the neutron background for an India based dark matter search experiment at Jaduguda Underground Science Laboratory (JUSL) is performed. The experiment at JUSL will be the first phase of a proposed Dark matter search at India-based Neutrino Observatory (DINO). It will be a direct detection experiment with primary aims to search for WIMPs as dark matter candidates. In this thesis, we discuss the methodology of estimating neutron flux at JUSL and report the results. The total neutron flux reaching the laboratory above 1 MeV energy threshold is found to be 5.76(±0.69) × 10 -6 cm -2 s -1 . The impact of neutron background on the sensitivity of the experiment to detect dark matter at JUSL is also discussed. The thesis is organized as follows. Chapter 1 provides a brief introduction to the Lightly Ionizing Particle (LIPs). The analysis to search LIPs in SuperCDMS is briefly outlined in the chapter. This chapter also discusses the importance of neutron background estimates in a dark matter search experiment, more specifically, in the context of a proposed India-based dark matter search experiment at Jaduguda Underground Science Laboratory. In Chapter 2, the SuperCDMS experiment is introduced. In Chapter 3, the framework developed to perform simulations for Lightly Ionizing Particles is presented. In Chapter 4, the LIPs search analysis with the CDMSlite data and the results are discussed. In Chapter 5, the simulation of neutron background and the feasibility of dark-matter search at JUSL is discussed. Finally, conclusions from all the results discussed in this thesis are presented in Chapter 6.

Banik, Samir↗
DOE OSTI2020

A coupled multipoint stress–multipoint flux mixed finite element method for the Biot system of poroelasticity

In this work, we present a mixed finite element method for a five-field formulation of the Biot system of poroelasticity that reduces to a cell-centered pressure–displacement system on simplicial and quadrilateral grids. A mixed stress–displacement–rotation formulation for elasticity with weak stress symmetry is coupled with a mixed velocity–pressure Darcy formulation. The spatial discretization is based on combining the multipoint stress mixed finite element (MSMFE) method for elasticity and the multipoint flux mixed finite element (MFMFE) method for Darcy flow. It uses the lowest order Brezzi–Douglas–Marini mixed finite element spaces for the poroelastic stress and Darcy velocity, piecewise constant displacement and pressure, and continuous piecewise linear or bilinear rotation. A vertex quadrature rule is applied to the velocity, stress, and stress–rotation bilinear forms, which block-diagonalizes the corresponding matrices and allows for local velocity, stress, and rotation elimination. This leads to a cell-centered positive-definite system for pressure and displacement at each time step. We perform error analysis for the semidiscrete and fully discrete formulations, establishing first order convergence for all variables in their natural norms. The numerical tests confirm the theoretical convergence rates and illustrate the locking-free property of the method.

42 ENGINEERING↗
DOE OSTI2026

Elevating SolTrace's Capabilities for the Next Generation of Concentrating Solar Analysis

SolTrace is an open-source Monte Carlo ray tracing software developed at NREL. SolTrace can characterize concentrating solar thermal (CST) collector optical performance and is CST technology agnostic. Shown in Fig. 1, SolTrace is a foundational tool in NREL's CST system and component modeling suite. SolTrace's generic surface elements can flexibly model novel collector and receiver designs to predict spatial and temporal flux distributions - critical to understand for CST component design, performance prediction, and system integration. Since its initial development, SolTrace has over 1,650 references on Google Scholar, over 9,800 downloads since 2017, and has served the CST research and development community as a benchmark of 3rd party verification. SolTrace provides users with many options for defining surface shape and boundaries. However, SolTrace provides limited documentation which can result in a steep learning curve for new users. Additionally, SolTrace lacks the computational performance required to evaluate optical performance of a CST system over the course of a year and/or iteratively over design parameters in a timely manner. To address this, we are working towards a new release of SolTrace that enables increased computational throughput by implementing ray tracing acceleration structures and enabling GPU parallelization. Additionally, we are working to improve SolTrace's usability, accessibility, and maintainability by (1) automating solar position time-dependent simulation processes, (2) creating general CST collector templates of grouped elements, (3) updating the user interface to better visualize model inputs and outputs, and (4) creating a user support network through forums, "how to" videos, and documentation.

14 SOLAR ENERGY↗
DOE OSTI2022

Numerical study of advanced solar receiver tubes based on a coupled thermo-mechanical analysis for concentrated solar power tower plant

The search for more sustainable energy to match the growing energy demand begins with finding more dispatchable resources such as solar energy. As one of the promising solar technologies, concentrated solar power (CSP) has a full capacity to store thermal energy for extended operation. Nevertheless, some key components in CSP systems usually face extreme environment, such as uneven solar flux, cyclic thermal expansion, structural degradation on the solar absorber tubes in a Concentrated Solar Power Tower (CSPT) Plant. In this study, we applied Multiphysics simulation to explore the benefits of introducing optimized fins for heat transfer enhancement and uniform temperature distribution, the goal is to improve the thermal efficiency of such advanced solar absorber tubes. The results of this study can supply design guidance for the manufacturing process of absorber tubes, and eventually can benefit the solar energy community for the next generation of molten salt based CSP system.

14 SOLAR ENERGY↗
DOE OSTI2025

Permeate fluxes from desalination of brines and produced waters: A reactive transport modeling study

The increasing interest in the use of membrane systems to desalinate inland brackish water, agricultural drainage, and industrially produced wastewater demands improved means of predicting desalination system performance under variable feedwater compositions. The interaction among water flow, solute transport, and chemical composition in these systems impacts permeate flux evolution. Here, an established multicomponent reactive transport simulator that accounts for these coupled processes is applied to compute osmotic pressure and permeate fluxes in reverse osmosis (RO) systems. The model is first validated by predicting permeate fluxes for a set of benchtop crossflow experiments subject to a range of feed flow rates and compositions, under fouling and non-fouling conditions. Results compare favorably with measured data that show that solutions with similar total dissolved solids concentrations but different compositions result in different permeate fluxes. The model is then applied to predict permeate fluxes from the desalination of produced waters using a commercial spiral wound RO module. For NaCl-dominant brines, at total dissolved salt concentrations (TDS) below about 70 g/L, permeate fluxes are inversely proportional to water mole fraction as the latter is a reasonable approximation of water activity (i.e. ideal mixing). In the case of Ca–Cl-, Na–CO3- and Na–SO4-dominant brines below about 70 g/L TDS, this relationship does not hold as well and tends to overpredict osmotic pressure and thus underpredict permeate fluxes. However, the opposite becomes true at higher TDS values for typical produced waters. The scaling potential of these waters is also computed by allowing the precipitation of minerals above their saturation limit on the RO membrane. This work demonstrates how reactive transport models developed for the analysis of waters from geological systems can be extended to improve process design, optimization, and control in desalination systems from produced waters and beyond.

Molins, Sergi↗
DOE OSTI2020

Advanced light source-Upgrade accumulator-ring gas bremsstrahlung production

Lawrence Berkeley National Laboratory is designing a new electron synchrotron for scientific research using synchrotron radiation. An upgrade of the Advanced Light Source (ALS) would position the nation to retain its international competitiveness in the soft X-ray sciences. The ALS is a 1.9 GeV storage ring operating at 500 mA of beam current. It is optimized to produce intense beams of soft X-rays, which offer spectroscopic contrast, nanometer-scale resolution, and broad temporal sensitivity. Upgrading the ALS storage ring from a triple-bend-achromat to a multibend-achromat (MBA) lattice design will provide a soft X-ray source that is up to 100–1,000 times brighter than today’s ALS and will provide a significantly higher fraction of coherent light in the soft X-ray region than is currently available at the ALS. Further, the upgraded ALS will offer the highest coherent flux of any existing or planned storage-ring facility, worldwide, up to a photon energy of 3.5 keV. The addition of an accumulator ring that enables on-axis, swap-out, recovery, and exchange of bunch trains is required to achieve this coherence and brightness. Analysis of the radiation generated and shielding required by the addition of this new accelerator is presented in this paper.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2026

Observed Land Surface Influence on Atmospheric Heat and Moisture Profiles During Interstorms

Land-atmospheric (L-A) feedbacks have historically been studied using models whose structure and parameterizations influence outcomes and insights. The representation of L-A feedbacks based on observations alone remains an ongoing challenge for understanding boundary layer development and precipitation. To address this gap, we use ground-based passive remote sensing and in-situ observations to present an analysis of the atmosphere during 103 interstorm soil moisture drydown events spanning nine warm seasons (2016–2024) in the U.S. Southern Great Plains region. By separating events based on local L-A coupling signals and characterizing the profiles of atmospheric heat and moisture to surface energy flux behavior, we investigate the physical mechanisms linking land surface processes to boundary layer development. We find that during interstorm drydowns, the atmospheric column follows a consistent pattern: moisture increases within the boundary layer, peaks near its top, and declines rapidly above, while warming occurs through the depth. Drydowns that shift toward evaporation produce stronger and deeper thermodynamic responses than cases dominated by sensible heating, which are weaker and shallower. Additionally, moisture is accumulated faster within the boundary layer during shorter drydowns, with longer drydowns representing slower, moisture-limited growth. Drydowns with wetter initial soil moisture will sustain stronger moistening within and above the boundary layer, accelerating buoyancy growth and convective potential toward the next storm. These results provide observational evidence linking surface flux evolution to boundary layer thermodynamics and offer a process-level benchmark for evaluating coupled L-A representations in models and demonstrating the influence of soil moisture on short-term weather forecasting skill.

Zhang, M. S. [Massachusetts Inst. of Technology (M↗
DOE OSTI2025

Selective semihydrogenation of acetylene in ethylene using defect-rich boron nitride catalyst from flux reconstruction

Efficient removal of trace acetylene from ethylene streams is essential for producing polymer-grade ethylene, yet achieving highly selective semihydrogenation without over-hydrogenation remains a long-standing challenge. A key barrier is the lack of a simple, low-cost catalyst that can activate hydrogen effectively while preventing ethylene from reacting further. Here we show that defect-rich boron nitride, prepared through a straightforward flux reconstruction method, serves as a highly selective and metal-free catalyst for acetylene semihydrogenation. The catalyst contains abundant open boron and nitrogen sites that enable efficient hydrogen activation and rapid release of ethylene, thereby avoiding over-hydrogenation. Experiments combined with isotope labeling and theoretical analysis reveal that these defects lower the energy barrier for hydrogen activation while accelerating product desorption. Our findings demonstrate a scalable strategy for defect engineering in boron nitride and highlight its potential as a robust, sustainable alternative to metal-based catalysts in industrial ethylene purification.

Wang, Tao [Oak Ridge National Laboratory (ORNL), O↗
DOE OSTI2023

Demonstrating Computational Equivalence Between Continuous and Discrete Adjoint Methods by Calculating Time-Dependent Adjoint Solutions with Neutron Diffusion Models

The continuous adjoint method and the discrete adjoint method are two alternative approaches used to calculate adjoint solutions for adjoint systems. The continuous adjoint method derives adjoint equations analytically from continuous forward equations and then solves the adjoint equations either analytically or numerically in a discretized form whereas the discrete adjoint method calculates the adjoint solutions directly from the discretized forward equations. With regard to the methodology development and calculation procedure, distinct differences are well recognized between the two methods. For certain reasons, both methods are exclusively preferred and commonly used by different computational communities, but limited studies clarify the connections between the two adjoint methods from either of the communities. Herein, this paper demonstrates the computational equivalence between the continuous and discrete adjoint methods by investigating time-dependent adjoint solutions to the two-group neutron diffusion model in nuclear reactor analysis problems using both methods. Adjoint solutions can be used to estimate system parameters for reactor safety analysis. Appropriate final state conditions for the adjoint systems are specified in both of the methods, and the conditions are clarified with proper physical explanations. With the help of an event-based case study on neutron diffusion models, the accuracy of the time-dependent adjoint fluxes obtained from both methods is verified, and the pros and cons of both adjoint methods are examined. More importantly, the computational equivalence of both methods is demonstrated when they are applied to multigroup neutron diffusion systems. The advantage of calculating time-dependent adjoint fluxes by directly solving time-dependent adjoint systems rather than taking steady-state approximations as in common practice is also demonstrated.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗