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.
MAESTROeX is a massively parallel, finite-volume C++/F90 solver for low Mach number astrophysical flows. Our code utilizes a low Mach number equation set allowing for more efficient, long-time integration of highly subsonic flows compared to compressible approaches. The recommended range of applicability is for flows where the Mach number does not exceed ~0.1.
The ndrop_mfrsr value-added product (VAP) provides an estimate of the cloud droplet number concentration of overcast water clouds retrieved from cloud optical depth from the multi-filter rotating shadowband radiometer (MFRSR) instrument and liquid water path (LWP) retrieved from the microwave radiometer (MWR). When cloud layer information is available from vertically pointing lidar and radars in the Active Remote Sensing of Clouds (ARSCL) product, the VAP also provides estimates of the adiabatic LWP and an adiabatic parameter (beta) that indicates how divergent the LWP is from the adiabatic case. quality control (QC) flags (qc_drop_number_conc), an uncertainty estimate (drop_number_conc_toterr), and a cloud layer type flag (cloud_base_type) are useful indicators of the quality and accuracy of any given value of the retrieval. Examples of these major input and output variables are given in sample plots in section 6.0.
The Los Alamos QRNG is a hardware-based high-performance Random Number Generator capable of generating 200 Mbit/s or more of true random numbers. The device harvests entropy from fluctuations in an optical source that arise from quantum mechanical properties of light. These quantum effects are irreducibly random; the resulting numbers are unpredictable and beyond the influence of any adversary. Qrypt, Inc., launched in 2017, has amassed multiple quantum entropy sources to create high-quality random keys at scale. The company is engaging with Los Alamos through license and a Cooperative Research and Development Agreement to facilitate the transition of QRNG technology and deploy the technology into the marketplace.
The Los Alamos Quantum Random Number Generator (QRNG) is a hardware-based, high-performance Random Number Generator capable of generating 200 Mbit/s or more of true random numbers. Like flipping a coin, it is very much random and essential for information security like encrypting data on the internet, checking email, or purchasing something from an online vendor. The device harvests entropy from fluctuations in an optical source that arise from quantum mechanical properties of light. Qrypt, Inc., a company launched in 2017, began making strategic investments and developing partnerships to advance cutting-edge quantum hardware solutions. One of those key investments was licensing QRNG from Los Alamos and subsequently collaborating with advanced quantum materials and technology researcher Dr. Raymond Newell to create high-quality random keys at scale.
This study evaluates reactor physics aspects of the proposed High Flux Isotope Reactor (HFIR) permanent beryllium reflector number 5 design. HFIR is a pressurized, light water–cooled, and beryllium reflected research reactor that operates at 85 MW for cycle lengths of approximately 24 to 26 days. A new permanent reflector design is highly desirable for increased versatility in irradiation experiments, to arrange the vertical experiment facilities (VXFs) to minimize their impact on neutron scattering if loaded with neutronabsorbing experiments, to enhance reflector thermal-structural-hydraulic performance, and to simplify the complex fabrication process. The current number 4 and concept number 5 reflectors have 22 and 28 VXFs, respectively, and are modeled and analyzed in this study with the MCNP5 and SCALE 6.1.3 nuclear simulation codes.
This is the final technical report for DOE DE-SC0021034 award entitled “Application of aerosol-ice nucleating particle closure to establish the leading parameters governing ice crystal number concentration under commonly observed mixed-phase cloud conditions”. The overarching objective of this project was to gain a predictive understanding of the processes that lead from aerosols to ice-nucleating particles (INPs) to ice crystal number concentrations (ICNCs) in stratiform mixed-phase clouds where supercooled water droplets and ice crystals coexist, and to inform climate model parameterizations. During this four-year project we addressed the three main research goals of this objective: 1) Study the aerosol-INP link; 2) Develop an aerosol-cloud column model (AC-1D); and 3) Study the INP-ICNC link. Those activities significantly advanced our predictive understanding of how to represent INPs in could models and its implications on ice crystal number concentrations. The results have been published in five peer-reviewed publications with three more in review and preparation. It allowed for collaboration resulting in additional fourteen peer-reviewed publications that are concerned with the objective of this research project. This research has been communicated to the scientific community at national and international meetings and conferences by numerous talks, seminars, and posters including fifteen invited presentations. This project allowed training of four graduate students and one postdoctoral research associate resulting in one M.S. thesis and advancement to Ph.D. candidacy.
In this letter, we reassess the underlying physics of the number sum rule for dihadron fragmentation functions. We will argue that, currently, there are no settled constraints on what constitutes a valid number density interpretation for multihadron fragmentation functions. Imposing overly restrictive criteria might lead to misinterpretating the data. Most importantly, and on the basis of phenomenological analyses, the slightly varying definitions used in previous work are not excluded from possessing legitimate number density interpretations (up to the usual issues with ultraviolet divergences and renormalization), so long as they are paired with appropriate factorization theorems. We advocate for further theoretical analyses to be challenged with experimental data, available at JLab or at the future EIC.
Atmospheric new particle formation (NPF) is a crucial process driving aerosol number concentrations in the atmosphere; it can significantly impact the evolution of atmospheric aerosol and cloud processes. This study analyses at least 1 year of asynchronous particle number size distributions from six different locations in India. We also analyze the frequency of NPF and its contribution to cloud condensation nuclei (CCN) concentrations. We found that the NPF frequency has a considerable seasonal variability. At the measurement sites analyzed in this study, NPF frequently occurs in March–May (pre-monsoon, about 21% of the days) and is the least common in October–November (post-monsoon, about 7% of the days). Considering the NPF events in all locations, the particle formation rate (J SDS ) varied by more than 2 orders of magnitude (0.001–0.6 cm –3 s –1 ) and the growth rate between the smallest detectable size and 25nm (GR SDS-25nm ) by about 3 orders of magnitude (0.2–17.2nm h –1 ). We found that J SDS was higher by nearly 1 order of magnitude during NPF events in urban areas than mountain sites. GR SDS did not show a systematic difference. Our results showed that NPF events could significantly modulate the shape of particle number size distributions and CCN concentrations in India. The contribution of a given NPF event to CCN concentrations was the highest in urban locations (4.3 × 10 3 cm –3 per event and 1.2 × 10 3 cm –3 per event for 50 and 100nm, respectively) as compared to mountain background sites (2.7 × 10 3 cm –3 per event and 1.0 × 10 3 cm –3 per event, respectively). We emphasize that the physical and chemical pathways responsible for NPF and factors that control its contribution to CCN production require in situ field observations using recent advances in aerosol and its precursor gaseous measurement techniques.
Accurate quantification of the vertical distribution of cloud condensation nuclei (CCN) number concentrations is critical for improving our understanding of aerosol–cloud interactions. Ground-based Raman lidars operated by the Atmospheric Radiation Measurement (ARM) program, together with surface CCN measurements, are used to retrieve vertically resolved CCN number concentrations (Retrieved Number concentration of CCN, RNCCN). These retrievals rely on several assumptions, including that aerosol composition is vertically homogeneous. To assess this assumption, we developed and tested a framework to infer the dominant aerosol classes/types at different altitudes. This was done by applying a k-Nearest-Neighbors (kNN) algorithm to lidar ratio and linear depolarization ratio measurements from Raman lidar. We evaluated the framework using aircraft aerosol and CCN measurements from the ARM Holistic Interactions of Shallow Clouds, Aerosols, and Land Ecosystems (HI-SCALE) field campaign. The results show that RNCCN performance degrades as vertical aerosol complexity increases, i.e., RNCCN agrees with the aircraft CCN in vertically homogeneous conditions, but closure decreases in layered aerosol structures. To generalize beyond individual examples, we introduce a metric (heterogeneity index) that quantifies the vertical complexity by assessing the variation of inferred aerosol classes/types. Case-level statistics show a tendency for RNCCN and aircraft differences to increase with this metric. By detecting retrievals that are likely compromised by aerosol vertical heterogeneity, the proposed framework improves the interpretability and effective use of RNCCN used for long-term evaluation of models and aerosol–cloud interactions.
Abstract. Optimizing turbine layout is a challenging problem that has been extensively researched in the literature. However, optimizing the number of turbines within a given boundary has not been studied as extensively and is a difficult problem because it introduces discrete design variables and a discontinuous design space. An essential step in performing wind power plant layout optimization is to define the objective function, or value, that is used to express what is valuable to a wind power plant developer, such as annual energy production, cost of energy, or profit. In this paper, we demonstrate the importance of selecting the appropriate objective function when optimizing a wind power plant in a land-constrained site. We optimized several different wind power plants with different wind resources and boundary sizes. Results show that the optimal number of turbines varies drastically depending on the objective function. For a simple, one-dimensional, land-based scenario, we found that a wind power plant optimized for minimal cost of energy produced just 72 % of the profit compared to the wind power plant optimized for maximum profit, which corresponded to a loss of about USD 2 million each year. This paper also compares the performance of several different optimization algorithms, including a novel repeated-sweep algorithm that we developed. We found that the performance of each algorithm depended on the number of design variables in the problem as well as the objective function.
Octane sensitivity (OS), defined as the research octane number (RON) minus the motor octane number (MON) of a fuel, has gained interest among researchers due to its effect on knocking conditions in internal combustion engines. Compounds with a high OS enable higher efficiencies, especially within advanced compression ignition engines. RON/MON must be experimentally tested to determine OS, requiring time, funding, and specialized equipment. Thus, predictive models trained with existing experimental data and molecular descriptors (via quantitative structure-property relationships (QSPRs)) would allow for the preemptive screening of compounds prior to performing these experiments. Here, the present work proposes two methods for predicting the OS of a given compound: using artificial neural networks (ANNs) trained with QSPR descriptors to predict RON and MON individually to compute OS (derived octane sensitivity (dOS)), and using ANNs trained with QSPR descriptors to directly predict OS. Twenty-five ANNs were trained for both RON and MON and their test sets achieved an overall 6.4% and 5.2% error, respectively. Twenty-five additional ANNs were trained for both dOS and OS; dOS calculations were found to have 15.3% error while predicting OS directly resulted in 9.9% error. A chemical analysis of the top QSPR descriptors for RON/MON and OS is conducted, highlighting desirable structural features for high-performing molecules and offering insight into the inner mathematical workings of ANNs; such chemical interpretations study the interconnections between structural features, descriptors, and fuel performance showing that connectivity, structural diversity, and atomic hybridization consistently drive fuel performance.
The purpose of this agreement between the National Renewable Energy Laboratory (NREL) and Gevo, under the DOE Small Business Vouchers Pilot, is to develop a predictive octane blending model for isobutanol and blendstocks for oxygenated blending (BOBs). Validated models of ethanol octane blending effects are currently used in refinery blending models. While it is known that isobutanol increases octane when blended into BOBs, the effect is non-linear, and dependent on hydrocarbon blendstock properties. Because no predictive model currently exists for isobutanol and BOBs, blenders must perform expensive and time-consuming tests to determine the octane effect on each batch of finished fuel. NREL will make basic measurements of research octane number (RON) and motor octane number (MON) of finished gasoline-isobutanol blends prepared in a broad range of BOBs. The key properties and detailed composition of BOBs utilized will be measured by NREL. Using these measured fuel properties, NREL will develop a predictive octane blending model for isobutanol in BOBs that can be used by terminals to further reduce barriers to market penetration of isobutanol.
We extend previous studies of the conformal 0+1d kinetic non-equilibrium attractor in relaxation time approximation by enforcing number conservation through the introduction of a dynamical fugacity (chemical potential). We derive two coupled integral equations for the effective temperature and fugacity which are then solved numerically to obtain the exact solution. The resulting solutions exhibit convergence to a unique nonequilibrium attractor when the scaled moments of the distribution function are plotted as a function of the rescaled time w = τ /τeq. This occurs even though the system is out of chemical equilibrium at late times. In addition, compared to the case where number conservation was not imposed, we find that the moments converge to their respective attractors more quickly, particularly for moments with m = 0. Finally, we compare the resulting attractor moments with predictions from different hydrodynamic frameworks.
Current research octane number (RON) and motor octane number (MON) gasoline performance characterization techniques use dated, complex engine testing methodology and limit researchers’ ability to easily characterize small volumes of experimental fuels. A novel methodology is presented that correlates measured ignition delay (ID) time to RON in an Advanced Fuel Ignition Delay Analyzer (AFIDA) constant-volume combustion chamber device at a single pressure/temperature condition, with an r2 of 0.99 and standard error (SE) of 1.0. The correlation of the slope of the ID time between two additional temperature points to octane sensitivity (S) produces an r2 of 0.97 and SE of 0.69; however, fuels with S>12 are indistinguishable. These results are based on methodology calibration using 31 primary and toluene reference fuels containing 0%-40% ethanol with RON values ranging from 85 to 113. Validation of these methods using a 102-sample fuel matrix spanning an array of base fuels and additive chemistry designed to test the robust applicability of the method, along with pump gasoline and high-octane surrogate blend samples, demonstrates an r2 of 0.94 and SE of 1.3 for the RON correlation over all samples, whereas the equivalent S correlation produces an r2 of 0.78 and SE of 1.2 by excluding two additives, 3-pentanone and diisobutylene, which displayed poor S correlation results. This novel AFIDA analysis method can be performed in 1 h and with 40 mL of fuel, offering significant improvements in time and volume requirements over traditional techniques.
The objective of this paper is to propose and validate a model for predicting the transitional film Reynolds number between the droplet and jet modes of falling film flows. The proposed model relies on thermodynamic availability with no fitting to mode transition data. The model directly predicts the numerical values of the transitional film Reynolds numbers at different operating conditions. The results of the model are validated against 43 experimental data points from the literature, which cover 10 working fluids and a wide range of tube geometries. The accuracy of the proposed model in predicting the validation data is -on average- more than 20 % better than that of the widely used Re = a (Ga) b and Re = a (Ga) b (S/ξ) c models. The model also achieves low computation time (less than 30 seconds for the 43 validation data points). The model presented in this paper has the potential of integration with heat transfer, mass transfer, and pressure drop models for optimizing the design and operation of falling film heat and mass exchangers.
In this work, 3D simulations of oil jets impinging on a flat, heated wall are presented. The numerical setup uses the Volume of Fluid (VoF) method to model the two-phase flow. A careful grid definition across the liquid film, along with the use of the Conjugate Heat Transfer (CHT) approach allowed local heat transfer to be solved with fine resolution at the wall. Variations of liquid flow rate, liquid temperature and surface temperature allow to cover a wide range of local Reynolds and Prandtl numbers (226 < Re < 2850, 77 < Pr < 161). Resulting surface-averaged heat transfer compares very well with experimental measurements conducted in a previous study. In-depth analysis of the flow has identified expected features from the literature. In particular, the impact of jet axial velocity profiles on the heat transfer distribution in the stagnation zone was clearly stated. The increase in heat transfer when warming the liquid film was also reproduced and explained by a decrease in oil viscosity and an increase in film velocity. All those effects were taken into account in correlations for stagnation and local values of Nusselt number. A grid sensitivity study was also conducted, showing that if the grid solving the thermal boundary layer in the stagnation zone can be coarsened without impacting local and surface-averaged predictions of heat transfer, a minimum resolution (2 to 3 cells) within the thermal boundary layer is however required for an accurate prediction of heat transfer.
Designing alternative fuels for advanced compression ignition engines necessitates a predictive model for cetane number (CN). In this study, the physics-informed graph neural networks are introduced for a reliable CN prediction by considering molecular features pertinent to the physical properties of molecules that affect CN. The reliability of measured data is another key factor to consider for improving the predictive model. Various experimental instruments for measuring CN exist, including standard and non-standard methods. In this regard, a systematic data quality analysis was carried out for the total 630 CNs collected from literature and new measurements in this study using Advanced Fuel Ignition Delay Analyzer (AFIDA). The results from this data curation process were reflected in the model by imposing lower sample weights on the data coming from less reliable measurement techniques. This approach effectively maximized the prediction accuracy while incorporating data from all available sources. Using the sample weights decreased the mean absolute error (MAE) up to 0.8 CN units. The accuracy was also improved by introducing the CN-related physical properties (the number of hydrogen bond donors and acceptors); the test set MAE is 5.74 and 7.01 for the model with and without such properties, respectively. Investigating molecular structural effects on CN was also carried out to gain chemical insights into factors used to design new fuel candidates. The dimensionality reduction analysis of feature vectors showed a clear clustering in terms of functional groups and CN and the structural effect derived from the model was consistent with the physicochemical insights. Finally, this physics-informed model and data curation would be helpful for accurate CN prediction and inform rational fuel design.
The stochastic theory of neutron transport is extended to describe the cumulative distribution of fission numbers and deposited fission energy in a multiplying assembly. Solutions for the probability distributions are obtained using analytical approximations and Monte Carlo simulation in lumped geometry and in symmetric homogeneous and heterogeneous spheres. The results show the development of a power-law tail in the steady state fission number and deposited energy distributions when the medium is critical, independent of the fission neutron multiplicity distribution and domain heterogeneity. In contrast, the asymptotic decay is faster than exponential in subcritical media due to rapid chain extinction and in supercritical media due to the increasing probability of chain divergence. Here, a formal asymptotic analysis of the problem in lumped geometry with an arbitrary fission neutron multiplicity confirms the existence of power-law tails at critical.