Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Coefficient of performance”

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,261 records · Page 70

PERFORMANCE CHARACTERIZATION OF THE DIRECT CONTACT MEMBRANE DISTILLATION MODULE OPERATED WITH HIGHLY CONCENTRATED SOLUTION

The experimental study was conducted to characterize the performance of the flat sheet Direct Contact Membrane Distillation. The length degradation was studied in a short and long module with and without spacers. The experimental study is validated using the predicted flux for a wide range of operating conditions using deionized water. Three different membranes are used in the parametric study for feed Reynolds number varied between 500 and 2000, and inlet feed temperature varied between 45°C and 80°C. PTFE membrane was shown to be more effective than the PVDF due to its higher permeability coefficient and higher hydrophobicity. Fouling study was conducted on the selected PTFE membrane with various feed Reynolds numbers between 500 and 2000 and feed inlet temperature between 50°C and 80°C. Membranes used in the fouling study were characterized using scanning electron microscopy (SEM) and X-ray photoelectron spectroscopy (XPS) for the morphology and chemical composition of the fouling layer. 36 hours of operation with 10 wt% NaCl (100,000 ppm) at feed inlet temperatures 80°C and 60°C and fixed Reynolds 2000 demonstrated a drop in flux by 31% and 27%, respectively. Low flow rate and high temperature are favorable conditions for fouling formation and reduced permeate flux for highly concentrated solutions. The spacers improved the flux by 30.8% and 27.6% at the operating Reynolds numbers 2000 the short and long channel modules. This study helps to design and optimize DCMD modules for industrial-size desalination and water treatment applications.

Kim, Deliya↗

Underlying mechanism of structural transformation between GaSb and GaAs response to intense electronic excitation

Ion irradiation of semiconductors has emerged as a promising approach for fabricating self-organized nanosystems with high atomic precision, despite often being accompanied by undesirable phenomena. Exploring the mechanisms underlying structural transformations is crucial for assessing nanostructure array types under complex irradiation environments. By quantitatively calculating the thermodynamically driven processes and analyzing the impact of intrinsic structural parameters, distinct structural transformations in response to intense electronic excitation are systematically investigated in gallium antimonide (GaSb) and gallium arsenide (GaAs) systems. In high-energy regimes, the nanofibers layer of GaSb exhibits intriguing structural discrepancy, characterized by partial nanofibers with coherent boundaries, interspersed nanopores accompanied by antisite defects and Ga precipitates, distinguishing to a series of discontinuous latent tracks that emerged within cylindrical trajectories in GaAs. Furthermore, significant diffusion behaviors of the nanohillocks are discovered in GaAs, with higher average roughness than GaSb, driven by the gradient stress distribution influenced by the free-surface effects. The deposition energy for melting phase formation, Gibbs free energy, and Ga diffusion coefficients contribute to the distinctive structural features, evidencing relatively stable morphological configurations and higher irradiation resistance in GaAs. Consequently, special optoelectronic properties associated with structural discrepancies facilitate the design and optimization of material functionalities by irradiation technologies.

36 MATERIALS SCIENCE↗

Replicating HCCI-like autoignition behavior: What gasoline surrogate fidelity is needed?

This work seeks to characterize the fidelity needed in a gasoline surrogate with the intent to replicate the complex autoignition behavior exhibited within advanced combustion engines, and specifically Homogeneous Charge Compression Ignition (HCCI). A low-temperature gasoline combustion (LGTC) engine operating in HCCI mode and a rapid compression machine (RCM) are utilized to experimentally quantify fuel reactivity, through autoignition and preliminary heat release characteristics. Fuels considered include a research grade E10 U.S. gasoline (RD5-87), three multi-component surrogates (PACE-1, PACE-8, PACE-20), and a binary surrogate (PRF88.4). Each fuel was studied at lean/HCCI-like conditions covering a wide range of temperatures and pressures that are representative of naturally aspirated to high boost engine operation. Detailed chemical kinetic modeling is also undertaken using a recently updated gasoline surrogate kinetic model to simulate the RCM experiments and to provide chemical insight into surrogate-to-surrogate differences. The LGTC engine experiments demonstrate nearly identical reactivity between PACE-20 and RD5-87 across studied conditions, while faster phasing is seen for both PACE-1 and PACE-8 due to their stronger intermediate- and low-temperature heat release (ITHR/LTHR) at naturally aspirated and boosted conditions, respectively. The RCM experiments reveal typical low-temperature, negative temperature coefficient (NTC) and intermediate-temperature autoignition behaviors at all pressure conditions for RD5-87, which are qualitatively reproduced by all surrogates. Quantitative discrepancies in both autoignition and preliminary heat release are observed for all surrogates, while their ability to replicate RD5-87 autoignition behavior follows the order of PACE-20 > PACE-1 > PACE-8 > PRF88.4. Excellent mapping is obtained between the LGTC engine and the RCM, where the engine pressure-time trajectories can be characterized by the regimes represented by the RCM autoignition isopleths. The kinetic model performs commendably when simulating both autoignition and preliminary heat release of PACE-20, while typically overpredicting ignition delay times for PACE-1, PACE-8 and PRF88.4 at high-pressure and low-temperature/NTC conditions. Sensitivity and rate of production (ROP) analyses highlight surrogate-to-surrogate differences in the governing chemical kinetics where n-pentane initiates rapid OH branching at a faster rate and an earlier timing for PACE-20 than iso-pentane does for PACE-1 and PACE-8, making it computationally more reactive than the other surrogates. The current study highlights the need to include non-standardized properties, such as the lean/HCCI-like autoignition characteristics, in addition to ASTM properties (e.g., RON, MON) as metrics of fuel reactivity and targets to be matched when formulating high-fidelity surrogates that fully capture gasoline advanced combustion behavior such as HCCI-like autoignition.

42 ENGINEERING↗

Neutron elastic and inelastic scattering differential cross sections on carbon

Elastic and inelastic neutron scattering angular distributions were measured on natural carbon samples to confirm existing experimental data and evaluations in the fast neutron region and to guide improvements in resonance parameters, where needed. Sixty-four (n,n') differential cross section measurements were performed at 45 incident neutron energies between 0.5 and 8.0 MeV. Experimental angle-integrated elastic scattering cross sections are consistent with ENDF/B-VIII.0 values with the exception of the region from 3.2 to 4.0 MeV where our results are ~3% higher. In the 3.4 to 3.6 MeV region our differential cross sections are slightly lower at forward angles and somewhat higher at backward angles than the ENDF/B-VIII.0 calculations; however, our results are consistent with measured data from other research groups. The first- through fourth-order elastic scattering Legendre coefficients from fits to these experimental data are consistent with ENDF/B-VIII.0 values across this range of incident neutron energies. Inelastic scattering cross sections were measured at 12 incident neutron energies between 5.6 and 7.8 MeV. Angle-integrated cross sections agree well with the 1978 Perey and ENDF/B-VIII.0 values; however, the angular-distribution shapes deviate somewhat from the ENDF/B-VIII.0 calculations. These data will allow refinement of the resonance parameter description in this difficult energy region.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Mixed Precision Fermi-Operator Expansion on Tensor Cores from a Machine Learning Perspective

Here we present a second-order recursive Fermi-operator expansion scheme using mixed precision floating point operations to perform electronic structure calculations using tensor core units. A performance of over 100 teraFLOPs is achieved for half-precision floating point operations on Nvidia’s A100 tensor core units. The second-order recursive Fermi-operator scheme is formulated in terms of a generalized, differentiable deep neural network structure, which solves the quantum mechanical electronic structure problem. We demonstrate how this network can be accelerated by optimizing the weight and bias values to substantially reduce the number of layers required for convergence. We also show how this machine learning approach can be used to optimize the coefficients of the recursive Fermi-operator expansion to accurately represent the fractional occupation numbers of the electronic states at finite temperatures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Motorette Thermal Performance Testing and Modeling for an Electric Motor with Additively Manufactured Hollow Conductors with Integrated Heat Pipes

This paper discusses the design, build and test of a motorette to characterize the thermal performance of an additively manufactured coil integrated with heat pipes. The motorette is powered using a variable frequency AC power supply and cooled through two independent cooling systems. The first cooling system cools the stator core using forced convection of air flowing over a finned heat sink mounted on the stator outer diameter. The second cooling system cools the heat pipes using forced convection of 50-50 water ethylene glycol (WEG) mixture flowing through condenser chambers. A 3D thermal FEA model of the setup is built and heat transfer coefficients (HTC) of convective boundaries are computed using coolant flow rates from test data and empirical equations. Temperature at locations on the coil and heat sink are recorded and compared to thermal model predicted values. A maximum temperature error of 22.7% occurs at 180ARMS and 800 Hz operating point for the AM dual coil. The thermal model provides higher temperatures compared to test results and hence it is on the conservative side but in general, there is good correlation between test and model results. Lastly, opportunities for improvement to both test data measurement and 3D thermal FEA modeling are discussed.

additive manufacturing↗

Numerical modeling of impurity powder injection in W7-X

In this work, we present numerical simulation results of powder injection in W7-X using the EMC3-EIRENE and DIS codes. First, we model powder injection experiments performed in W7-X with the Probe Mounted Powder Injector. The simulation results qualitatively agree with visible imaging measurements. Secondly, we perform predictive simulations to guide the installation of an Impurity Powder Dropper in W7-X, allowing to choose in between several available non-vertical ports to maximize the amount of powder penetrating into the plasma, as well as the verticality of the port, to minimize sticking of the powders in the in-vessel stainless steel guiding tube. Port AEM41 is selected as the best candidate for IPD installation. The robustness of the simulation results has been verified for different plasma densities, powder materials and sizes, powder friction coefficient and changes in the plasma flow.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

QRF4P-NRT: Probabilistic Post-Processing of Near-Real-Time Satellite Precipitation Estimates Using Quantile Regression Forests

Accurate and reliable near-real-time satellite precipitation estimation is of great importance for operational large-scale flood forecasting and drought monitoring. The state-of-the-art precipitation post-processing model is based on a deterministic approach to construct relationships between satellites estimates and ground observations. We propose a probabilistic postprocessor, the Probabilistic Post-Processing of Near-Real-Time Satellite Precipitation Estimates using Quantile Regression Forests (QRF4P-NRT), based on quantile modeling, yielding both deterministic and probabilistic predictions. The experimental design incorporates different solutions of near-real-time predictors to further improve the model performance. Using the Integrated Multi-satellitE Retrievals Early Run for Global Precipitation Measurement Mission (IMERG-E) product as an example, we illustrate that the proposed method significantly improves the overall quality of the raw IMERG-E and is also superior to the bias-corrected product (IMERG Final Run, IMERG-F) at daily scale in a complex mountain basin. Evaluations of the corrected IMERG-E, raw IMERG-E, and IMERG-F using ground observation show that the corrected IMERG-E improves correlation coefficients (0.7), mean error (-0.14 mm/day) and root mean square error (3.3 mm/day) relative to the raw IMERG-E (0.31, -0.72 and 5.5 mm/day) and IMERG-F (0.34, -0.09 and 6.0 mm/day). The error decomposition further confirms that the QRF4P-NRT improves on the various deficiencies of the raw IMERG-E product. The ensemble assessment also demonstrates that the quantile outputs provide reliable prediction spread and sharp prediction intervals. The promising results indicate the great potential of the proposed method for probabilistic post-processing for near-real-time satellite precipitation estimates, and for further applications such as hydrological ensemble forecasting.

54 ENVIRONMENTAL SCIENCES↗

An LES study of secondary motion and wall shear stresses in a pipe bend

In this study, large-eddy simulations (LES) of a single-phase, turbulent flow in a 90 ° pipe bend are performed at three Reynolds numbers (5300, 27 000, and 45 000) to investigate the correlation between secondary flow motion and wall shear stresses, which is suspected to be a potential mechanism responsible for material erosion. The isothermal flows are validated against available experimental and numerical data first. The snapshot proper orthogonal decomposition (POD) is applied for the medium and high Reynolds number flows to identify the secondary flow motions and the oscillation of the Dean vortices that are found to cause swirl-switching. Distinguished frequencies of the POD time coefficients at Strouhal numbers of 0.25 and 0.28 are identified for Reynolds numbers at 27 000 and 45 000, respectively. Moreover, shear stress on the pipe wall and the associated power spectral density are obtained and shown to have the same oscillating frequency as the swirl-switching.

42 ENGINEERING↗

Quantitative assessment of selenium diffusion and passivation in CdSeTe solar cells probed by spatially resolved cathodoluminescence

The introduction of selenium in CdSeTe/CdTe solar cells has led to improved device performances attributed to the passivation of bulk defects. In this work, high-resolution cathodoluminescence experiments are performed on a series of CdSeTe/CdTe thin films with different Se concentrations to quantify the mechanisms and the passivation role of Se. We demonstrate a universal dependence between the Se concentration and the radiative efficiency and a ten-fold enhancement of the luminescence between CdTe and CdSe 0.4 Te 0.6 . Raw luminescence maps are converted into maps of the Se concentration, revealing its graded profile within the stack. We demonstrate the diffusion of Se along CdTe grain boundaries induced by the cadmium chloride annealing treatment and determine the diffusion coefficients, which are more than eight times higher at grain boundaries than in grain interiors. These results provide microscopic insights into the distribution of Se and its impact on the passivation of CdSeTe/CdTe solar cells.

14 SOLAR ENERGY↗

MIC-DP: A Scalable Correlation-Aware Differential Privacy Framework for High-Dimensional Data

Conventional differential privacy (DP) assumes record independence, limiting effectiveness on real-world datasets with temporal, spatial, or structural correlations. These dependencies undermine privacy guarantees and degrade utility in domains like healthcare, IoT, and smart city analytics. We propose Maximum Information Correlated Differential Privacy (MIC-DP), a novel framework that dynamically calibrates noise based on statistical dependencies. MIC-DP uses the Maximum Information Coefficient (MIC) to capture both linear and nonlinear correlations without explicit modeling, enabling adaptive sensitivity adjustment and improved privacy–utility trade-offs. Evaluations on healthcare (MIMIC), demographic (ACI), and synthetic datasets show that MIC-DP reduces mean absolute error (MAE) by up to 5.2% under strict privacy budgets (ϵ≤1), with aggregate utility improvements reaching 18% across datasets and evaluation metrics. MIC-DP provides formal (ϵ,δ)-privacy guarantees, scales efficiently with feature count, and supports deployment in moderate-scale, privacy-sensitive applications. Its tunable performance and runtime efficiency make MIC-DP suitable for privacy-sensitive applications where low-latency analytics and strong privacy guarantees must coexist. These results demonstrate MIC-DP’s effectiveness as a correlation-aware solution for practical DP.

Yang, Wenjun [Univ. of Washington, Tacoma, WA (Uni↗

Vertical and Lateral Etch Survey of Ferroelectric AlN/Al 1-x Sc x N in Aqueous KOH Solutions

Due to their favorable electromechanical properties, such as high sound velocity, low dielectric permittivity and high electromechanical coupling, Aluminum Nitride (AlN) and Aluminum Scandium Nitride (Al 1-x Sc x N) thin films have achieved widespread application in radio frequency (RF) acoustic devices. The resistance to etching at high scandium alloying, however, has inhibited the realization of devices able to exploit the highest electromechanical coupling coefficients. In this work, we investigated the vertical and lateral etch rates of sputtered AlN and Al 1-x Sc x N with Sc concentration x ranging from 0 to 0.42 in aqueous potassium hydroxide (KOH). Etch rates and the sidewall angles were reported at different temperatures and KOH concentrations. We found that the trends of the etch rate were unanimous: while the vertical etch rate decreases with increasing Sc alloying, the lateral etch rate exhibits a V-shaped transition with a minimum etch rate at x = 0.125. By performing an etch on an 800 nm thick Al 0.875 Sc 0.125 N film with 10 wt% KOH at 65 °C for 20 min, a vertical sidewall was formed by exploiting the ratio of the {$10\overline{11}$ planes and {$1\overline{1}00$ planes etch rates. This method does not require preliminary processing and is potentially beneficial for the fabrication of lamb wave resonators (LWRs) or other microelectromechanical systems (MEMS) structures, laser mirrors and Ultraviolet Light-Emitting Diodes (UV-LEDs). It was demonstrated that the sidewall angle tracks the trajectory that follows the {$\overline{1}2\overline{12}$ of the hexagonal crystal structure when different c/a ratios were considered for elevated Sc alloying levels, which may be used as a convenient tool for structure/composition analysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Progress on Computational Modeling of Water-Based NSTF (FY2021)

This report documents the FY21 progress and achievements made in the computational analyses of the water-based NSTF. Both system-level and high-fidelity Computational Fluid Dynamics (CFD) analyses were performed to gain a complete understanding of the complex flow and heat transfer phenomena in natural convection systems. The progress on the waterbased NSTF experimental testing is summarized in a companion report (ANL-ART-230). As a continuation of progress from previous years, in FY21 the RELAP5 model of the NSTF was first updated to include a representative heated cavity that incorporated convection and conduction means of heat transfer. Two different cavity models were developed, one that models the radiation heat exchange among all surfaces in one enclosure (Cavity Model 1), and another that employs multiple radiation enclosures with each considering the surfaces at one axial level (Cavity Model 2). This updated RELAP5 deck was benchmarked with single-phase test data to tune form loss coefficients of the elbows and tees, along with heat transfer coefficients within the cavity and off the external insulation panels. This tuned model was then found to accurately predict fluid temperatures within 2.5% and both system and riser flow rates within 6% during single-phase, steady-state simulations. For transient simulations, there is a slight over prediction of the flow rate, fluid temperatures, and heater temperatures during the heat-up period, which indicate that the transient heat loss may be underestimated. Following, the capabilities of RELAP5-3D were examined, which feature a conduction enclosure model that is not present in the current version of RELAP5-MOD3.3 used for all previous NSTF simulations. As an initial comparison of the two codes, single-phase simulations were run that used similar input decks. The results indicate that the single-phase flow predicted by the two codes was very similar, however RELAP5-3D appears to predict a higher system flow during the two-phase transient when compared to the flow rates predicted by RELAP5-MOD3.3. Furthermore, RELAP5-3D run did not capture the density wave oscillations or geysering phenomena, which are large oscillation flow instabilities observed in the experimental facility. To further the prediction capabilities for capturing two-phase instabilities, Options 55, 57, 58, 61, which are implemented by default in RELAP5-MOD3.3 but not in RELAP5-3D, were investigated. When enabling Option 61, RELAP5-3D agrees well with MOD3.3 in predicting the total system flow, while other options have negligible effect on the predicted total system flow by RELAP5-3D. However, none of these options enables RELAP5-3D to produce the same flow instabilities. For future steps, efforts are currently underway to investigate RELAP5-3D’s different behavior in predicting low-pressure two-phase instabilities. STAR-CCM+ was used to perform CFD analysis of the water-based test facility and provide simulation results for comparisons to experimental data. Efforts first focused on describing the characteristics of flow within the heated cavity. A study was performed on whether possible leakage flow between the two sides of the cavity could substantially impact and explain the results. An additional second assumption was examined which defined adiabatic side walls, so physical insulation structures were added to the models to allow conduction across this outer insulation. The results suggest that potential porosity of the insulation could be a major factor in the temperature distribution within the cavity and should be considered in any studies where it is desirable to know the cavity wall temperatures to higher accuracy.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

The performance of ensemble-based free energy protocols in computing binding affinities to ROS1 kinase

Optimization of binding affinities for compounds to their target protein is a primary objective in drug discovery. Herein we report on a collaborative study that evaluates a set of compounds binding to ROS1 kinase. We use ESMACS (enhanced sampling of molecular dynamics with approximation of continuum solvent) and TIES (thermodynamic integration with enhanced sampling) protocols to rank the binding free energies. The predicted binding free energies from ESMACS simulations show good correlations with experimental data for subsets of the compounds. Consistent binding free energy differences are generated for TIES and ESMACS. Although an unexplained overestimation exists, we obtain excellent statistical rankings across the set of compounds from the TIES protocol, with a Pearson correlation coefficient of 0.90 between calculated and experimental activities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Neural network based fast prediction of β N limits in HL-2M

Artificial neural networks (NNs) are trained, based on the numerical database, to predict the no-wall and ideal-wall β N limits, due to onset of the n = 1 ( n is the toroidal mode number) ideal external kink instability, for the HL-2M tokamak. The database is constructed by toroidal computations utilizing both the equilibrium code CHEASE (Lütjens et al 1992 Comput. Phys. Commun. 69 287) and the stability code MARS-F (Liu et al 2000 Phys. Plasmas 7 3681). The stability results show that (1) the plasma elongation generally enhances both β N limits, for either positive or negative triangularity plasmas; (2) the effect is more pronounced for positive triangularity plasmas; (3) the computed no-wall β N limit linearly scales with the plasma internal inductance, with the proportionality coefficient ranging between 1 and 5 for HL-2M; (4) the no-wall limit substantially decreases with increasing pressure peaking factor. Furthermore, both the NN model and the convolutional neural network (CNN) model are trained and tested, producing consistent results. The trained NNs predict both the no-wall and ideal-wall limits with as high as 95% accuracy, compared to those directly computed by the stability code. Additional test cases, produced by the Tokamak Simulation Code (Jardin et al 1993 Nucl. Fusion 33 371), also show reasonable performance of the trained NNs, with the relative error being within 10%. The constructed database provides effective references for the future HL-2M operations. The trained NNs can be used as a real-time monitor for disruption prevention in the HL-2M experiments, or serve as part of the integrated modeling tools for ideal kink stability analysis.

Physics↗

Automatically generated model for light alkene combustion

Light alkenes are common combustion intermediates for a variety of fuels. Therefore, understanding their oxidation and pyrolysis chemistry is key to building detailed mechanisms for heavier fuels. This work was focused on the development and evaluation of a detailed kinetic mechanism suitable for the combustion of light alkenes up to C 4 without the use of tuned parameters, instead, the parameter values come from first principles or direct measurements. The generated mechanism accurately estimates the laminar burning velocity (S u ) and ignition delay time (IDT) of light alkenes available in the literature, which represent fundamental combustion properties at a wide range of conditions. Because each parameter is thought to have a physically realistic value, not tuned to these measurements, the new model could be used as a sub-mechanism in models for other applications. The reaction network was generated with the open-source Reaction Mechanism Generator (RMG) software. Sensitivity analyses were performed under wide ranges of temperatures and pressures, allowing for the identification of the most impactful species and reactions. Based on these, a comprehensive thermochemistry database, including calculations on 550 molecules performed in this work at the CBS-QB3 level of theory, and a kinetic library, including theoretically-derived reaction rates retrieved from the literature, were built and used in the mechanism generation. The developed mechanism was compared against several existing detailed kinetic mechanisms for ethene, propene, 1-butene, 2-butene, and isobutene. Here, the newly generated model is the most accurate among the ones analyzed, in terms of fractional bias and normalized mean square error. Hence, this new model was used to analyze the chemistry of alkene combustion. Key rate coefficients were compared, to identify the cause of deviations between the models and possible areas for further improvements.

09 BIOMASS FUELS↗

An Algorithm to Retrieve Total Precipitable Water Vapor in the Atmosphere from FengYun 3D Medium Resolution Spectral Imager 2 (FY-3D MERSI-2) Data

The atmosphere has substantial effects on optical remote sensing imagery of the Earth’s surface from space. These effects come through the functioning of atmospheric particles on the radiometric transfer from the Earth’s surface through the atmosphere to the sensor in space. Precipitable water vapor (PWV), CO 2 , ozone, and aerosol in the atmosphere are very important among the particles through their functioning. This study presented an algorithm to retrieve total PWV from the Chinese second-generation polar-orbiting meteorological satellite FengYun 3D Medium Resolution Spectral Imager 2 (FY-3D MERSI-2) data, which have three near-infrared (NIR) water vapor absorbing channels, i.e., channel 16, 17, and 18. The algorithm was improved from the radiance ratio technique initially developed for Moderate-Resolution Imaging Spectroradiometer (MODIS) data. MODTRAN 5 was used to simulate the process of radiant transfer from the ground surfaces to the sensor at various atmospheric conditions for estimation of the coefficients of ratio technique, which was achieved through statistical regression analysis between the simulated radiance and transmittance values for FY-3D MERSI-2 NIR channels. The algorithm was then constructed as a linear combination of the three-water vapor absorbing channels of FY-3D MERSI-2. Measurements from two ground-based reference datasets were used to validate the algorithm: the sun photometer measurements of Aerosol Robotic Network (AERONET) and the microwave radiometer measurements of Energy’s Atmospheric Radiation Measurement Program (ARMP). The validation results showed that the algorithm performs very well when compared with the ground-based reference datasets. The estimated PWV values come with root mean square error (RMSE) of 0.28 g/cm 2 for the ARMP and 0.26 g/cm 2 for the AERONET datasets, with bias of 0.072 g/cm 2 and 0.096 g/cm 2 for the two reference datasets, respectively. The accuracy of the proposed algorithm revealed a better consistency with ground-based reference datasets. Thus, the proposed algorithm could be used as an alternative to retrieve PWV from FY-3D MERSI-2 data for various remote sensing applications such as agricultural monitoring, climate change, hydrologic cycle, and so on at various regional and global scales.

54 ENVIRONMENTAL SCIENCES↗

Assessment of Core Physics Characteristics of Extended Enrichment and Higher Burnup LWR Fuels using the Polaris/PARCS Two-Step Approach. Vol. I: PWR Fuel

Nuclear fuel with extended enrichment (235U enrichment within 5-8 wt%) is one of the evolutionary changes that have been pursued in recent years by commercial light water reactor operators and fuel vendors to improve the fuel cycle economy and operation performance of a nuclear plant. This work assesses the performance of the Polaris/PARCS two-step approach in core physics modeling of the pressurized water reactor cores with extended enrichment fuel, referred to as “LEU+” in this report. A representative LEU+ core with a 24-month fuel cycle developed by Southern Nuclear Company (SNC) was modeled using this two-step approach. A representative LEU core with an 18-month fuel cycle was also modeled to provide a reference for the LEU+ core. As expected, significantly more burnable poison absorbers were used in the LEU+ core to accomodate its higher fuel enrichment. Nine different fuel assembly types were modeled using Polaris for each core to generate the assembly cross sections, which were then processed by GenPMAX to prepare the cross-section data for PARCS. The average specific powers of each fuel batch in each core were derived from VERA results and higher specific powers in fresh assemblies were found in the LEU core due to its less total uranium loading included in the VERA LEU model, given that the total core power was assumed to be the same for both cores. PARCS models were developed to simulate the steady-state operations of both cores. PARCS results on the LEU+ core were first compared with the VERA results for verification purpose; good agreements were seen in soluble boron and burnup distribution results, indicating that the Polaris/PARCS modeling and simulation were correctly implemented. Core physics parameters calculated by PARCS, at zero power physics tests, beginning of cycle (BOC), and end of cycle conditions (EOC), were compared between the LEU+ core and the LEU core, including soluble boron concentration, burnup distributions, assembly and pin power peaking factors, fuel temperature reactivity coefficients, moderator temperature and density reactivity coefficients, control rod worth, and shut down margin. The main differences in PARCS results between the LEU+ and the LEU cores are summarized below: 1)The critical boron concentrations were found to be much higher in the LEU+ core than in the LEU core (1582 vs. 1335 ppm for peak values). 2)Higher assembly radial power peaking factors (1.4 vs. 1.3 for peak values), 2D pin peaking factors (1.53 vs. 1.42 for peak values), and 3D pin peaking factors (1.89 vs. 1.81 for peak values) were found in the LEU+ core than in the LEU core. 3)Significantly higher reactivity coefficients of moderator temperature (and density) were found in LEU+ than in LEU.4)Significantly lower control rod worth at EOC were found in LEU+ than LEU for all but one control banks.5)Significantly lower shut down margins were found in the LEU+ core than in the LEU core. The Polaris/GenPMAX/PARCS code suite was found to be capable of modeling the LEU+ PWR core for steady-state operations and no unexpected results in core physics parameters were observed, in spite of that a) several bugs in PARCS were identified and workarounds were used; b) several features were found lacking in the current version of PARCS that would be useful for core modeling. A list of requests for bug fixes and feature upgrades for PARCS originated from this work were transmitted to the code developers. The assessments on the performance of the Polaris/PARCS two-step approach in core modeling for boiling water reactor with LEU+ fuel is ongoing.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗