Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “RANS”

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
NASA NTRSDate not supplied

New Capabilities and Improvements to the High-Order Glenn Flux Reconstruction Code

The Glenn Flux Reconstruction (GFR) code is a computational fluid dynamics (CFD) code under development at NASA Glenn Research Center. GFR is based on the high-order flux reconstruction (FR) method and provides a large-eddy simulation (LES) capability that is both accurate and efficient for complex aeropropulsion flows. Three significant new capabilities have been added to the code that improve its performance and functionality. First, a variety of explicit Runge-Kutta methods, including some with adaptive time stepping, were added to GFR with two methods offering a 33% improvement in time-to-solution. Second, GFR can now utilize fully unstructured, mixed-element meshes to more easily facilitate the grid generation process for complex geometries. Finally, a rotating reference frame capability has been added to GFR for solving rotating turbomachinery problems. A selection of results demonstrating these new capabilities are presented in this work. The Taylor-Green vortex problem is used to verify the new unstructured capability by showing similar accuracy and resolution for all element types. LES of the Turbulent Heat Flux Phase III (THX3) experiment with comparison to another high-order LES code and a popular Reynolds-averaged Navier-Stokes (RANS) code demonstrates the accuracy of the code for complex aeropropulsion flows. Finally, LES of a spacecraft cabin ventilation fan shows the ability of GFR to efficiently establish a fan performance map and identify operating points for further analysis at high orders of accuracy.

Direct Numerical Simulations↗
DOE OSTI2019

COMBUSTION AND EMISSION CHARACTERISTICS OF A LEAN AXIAL-STAGE COMBUSTOR

Gas turbine engines with high operational flexibility represent an important area of research. In support of this, the reduction of NOx emissions at high-temperature, full-load conditions are the main target of this research study. To overcome the challenge of decreasing the residence time of hot exhaust gas in the burner, the approach of axial-staged combustion has been used, which aims to burn at lean conditions. A main combustion stage with a conventional burner is followed by a jet-in-crossflow (JIC) stage downstream. This second combustion stage, the test section, burns well over a third of the total fuel amount, which is added by a jet inlet. A steady Computational Fluid Dynamics (CFD) model in Star-CCM+ was created to analyze the JIC stage. A structured mesh grid with 1.7 million cells and local refinement domains was used along with RANS equations, Realizable k-ε turbulence, Complex Chemistry model with Laminar Flame Concept (LFC) applied. The findings of the computational model include a detailed description of the flow field, focusing on vortex formation and penetration depth. Three vortex types and local pressure losses of up to 6% in the mixing domain were identified. Results were compared with similar data from literature and a very good match was yielded. Optical laser diagnostics should work well to capture the predicted flow penetration. The highest local reaction rates were observed inside a thin flame front, 4 to 6 diameters downstream of the jet. The jet showed a limited expansion into the available channel domain, resulting in low heat release and a methane conversion of only 52%. Providing an enhanced amount of oxygen was proven to increase the CH4 conversion significantly. Emissions were analyzed for an estimated species composition from the main burner. A CO production of 9.85‰ and the reduction of an incoming NOx amount by 15% were determined, outlining the potential of axial-staged combustion. Further numerical investigation along with experimental verification is in progress.

03 NATURAL GAS↗
DOE OSTI2019

MULTIDIMENSIONAL NUMERICAL SIMULATIONS OF REACTING FLOW IN A NON-PREMIXED ROTATING DETONATION ENGINE

Over the last two decades, detonation based propulsion has received a great deal of attention as a potential means to achieve significant improvement in the performance of air-breathing and rocket engines. Detonative combustion mode is particularly interesting due to the resulting pressure gain from reactants to products, faster heat release, decreased entropy generation, more available work and higher thrust compared to conventional deflagrative combustion. Rotating detonation engine (RDE) is one such novel combustor concept. Realistic RDE configurations utilize separate fuel and air injection schemes, hence are not perfectly premixed. Moreover, RDE performance is governed by a large number of design parameters and operating conditions. In this context, computational fluid dynamics (CFD) has the potential to enhance the understanding of RDE combustion and aid future development/optimization of this technology. In the present work, a CFD model was developed to simulate a representative non-premixed RDE combustor. Unsteady Reynolds-Averaged Navier-Stokes (RANS) simulations were performed for the full combustor geometry (including the separate fuel and air injection ports), with hydrogen as fuel and air as the oxidizer. Adaptive mesh refinement (AMR) was incorporated to achieve a trade-off between model accuracy and computational expense. A finite-rate chemistry model along with a 10-species detailed kinetic mechanism was employed to describe the H2-Air combustion chemistry. Two operating conditions were simulated, corresponding to the same global equivalence ratio of unity but different fuel and air mass flow rates. For both conditions, the capability of the model to capture the essential detonation wave dynamics was assessed. A validation study was performed against experimental data available on detonation wave frequency/height, reactant fill height, oblique shock angle, axial pressure distribution in the channel, and fuel/air plenum pressure. The CFD model predicted the sensitivity of these wave characteristics to the operating conditions with good accuracy, both qualitatively and quantitatively. The present CFD model offers a potential capability to perform rapid design space exploration and/or performance optimization studies for realistic full-scale RDE configurations.

Pal, Pinaki↗
DOE OSTI2019

Development of WEC Design Loads: A Comparison of Numerical and Experimental Approaches: Preprint

The development of accurate design loads is a critical part of the design of a wave energy converter (WEC). In this paper, we evaluate the impact of different approaches to determining extreme wave loading on the Triton WEC using a combination of mid-fidelity and high-fidelity numerical modeling tools, complemented by scaled physical model tests. The mid-fidelity approach used is a time-domain model based on linearized potential flow hydrodynamics while the high-fidelity modeling tool is an unsteady RANS CFD model. A 1:30 scale physical model of the Triton WEC was tested at the Oregon State University Large Wave Flume. A comparison will be presented between the design loads predicted by the mid-fidelity model, the CFD model, and the physical model tests, and suggestions for best practices will be offered.

16 TIDAL AND WAVE POWER↗
DOE OSTI2019

An Evaluation of Advanced Tools for Distributed Wind Turbine Performance Estimation: Preprint

We evaluate various classes of distributed wind turbine performance tools across two sites in the USA. The tools class ranges from the simple mass conservation model to coupled Reynolds Averaged Navier-Stokes (RANS) model, all initiated by the WIND Toolkit dataset. The resource estimation at the site is evaluated against measured data at the mast location. Taking a sample 100 kW wind turbine and constant losses, we evaluate Openwind, Continuum, and WindNinja tools and document Annual Energy Production (AEP) and time series statistics associated with the performance estimation of the wind turbine. Using a methodology that is consistent and unbiased across the three advance options currently available in the industry, we elaborate results at the two sample locations and discuss the probable sources of discrepancy in the AEP estimates. The discussion includes additional values these tools may bring into the energy assessment process to enhance the owners' confidence over the distributed wind power systems.

17 WIND ENERGY↗
DOE OSTI2020

Multiphase Reacting Flow Simulations and Optimization of Commercial-Scale Aerobic Bioreactors

Aerobic bioreactors are widely used in the synthesis of value-added products in pharmaceutical and biotechnology industries, wherein microbial bioreactions aid in the conversion of sugars to macromolecules. These reactors are also being actively investigated as a cost-effective pathway for the production of low-value commodities such as biofuels and animal feed. Reactor scale-up is one of the important challenges when designing these systems. Gas and liquid phase transport, mass-transfer, and mixing physics at large length scales can significantly affect microbial conversion rates, particularly when the microbial reaction requires a narrow set of conditions. These phenomena are difficult to study in small-scale bench-top reactors that are typically well-mixed. Predictive computational fluid dynamics (CFD) based simulations can therefore aid in the design and optimization of these reactors. This work presents multiphase Euler-Euler CFD simulations of commercial-scale (~ 500 m3) aerobic bioreactors. Our mathematical model treats the gas and liquid as interpenetrating phases. This approach reduces the computational complexity of tracking individual gas bubbles that are several orders of magnitude smaller than reactor dimensions. We solve the Reynolds averaged Navier-Stokes (RANS) multiphase equations that account for phase and chemical species transport, interphase mass and momentum transfer and uses a phenomenological model for oxygen uptake by microbes. We use a customized solver derived from open-source CFD toolbox, OpenFOAM, to perform these simulations, which has been validated against small-scale reactors in our previous work. This work examines the performance of three different reactor designs, viz. bubble column reactor, airlift reactor with an internal draft tube, and a stirred-tank reactor with Rushton impellers. Reactor oxygen mass-transfer coefficient, gas hold-up, and oxygen distribution are critically analyzed among reactors, and sensitivity studies pertaining to gas flow rates and reactor geometry will be presented.

09 BIOMASS FUELS↗
DOE OSTI2019

A Lifeline-Based Approach for Work Requesting and Parallel Particle Advection

Particle advection, a fundamental building block for many flow visualization algorithms, is very difficult to parallelize efficiently. That said, work requesting is a promising technique to improve parallel performance for particle advection. With this work, we introduce a new work requesting-based method which uses the Lifeline scheduling method. To evaluate the impact of this new algorithm, we ran 92 experiments, running at concurrencies as high as 8192 cores, data sets as large as 17 billion cells, and as many as 16 million particles, comparing against other work requesting scheduling methods. Overall, our results show that Lifeline has significantly less idle time than other approaches, since it reduces the number of failed attempts to request work.

Pugmire, Dave↗
DOE OSTI2020

Measurement of the EMC Effect of the Helium-3 Nucleus at Jefferson Lab

The MARATHON (MeAsurement of the Fn 2 /Fp 2 , d/u RAtios and A=3 EMC Effect in Deep Inelastic Scattering Off the Tritium and Helium MirrOr Nuclei) experiment ran in the Hall A Facility of the Thomas Jefferson National Accelerator Facility (JLab). The experiment used a 10.59 GeV electron beam from the Continuous Electron Beam Accelerator Facility (CEBAF) of JLab to measure the ratio of the F2 structure functions of Helium-3 and Deuterium. Data were taken in the kinematic range of 0.195 lte x lte 0.825, 2.5 (GeV/c)2 lte Q2 lte 13 (GeV/c)2, and 3.5 (GeV/c2)^2 lte W^2 lte 13 (GeV/c2)^2. This is the first measurement of the 3He EMC effect that is purely in the Deep Inelastic Scattering region. The results of this experiment constitute an important, complimentary addition, to the limited data presently available for the EMC effect of light nuclei. The MARATHON results are compared to the world data for light, medium, and heavy nuclei. This measurement is considered an essential component to uncovering the nature of the EMC effect.

Hague, Tyler↗
DOE OSTI2020

Measurement of the Ratio of the Neutron to Proton Structure Functions, and the Three-Nucleon EMC Effect in Deep Inelastic Electron Scattering Off Tritium and Helium-3 Mirror Nuclei

The proton and neutron structure functions Fp 2 and Fn 2, respectively are fundamental to understanding many studies in nuclear physics. They provide important information about quark distributions. For example, the ratio Fn 2/Fp 2 is one of the best measurements to find the ratio of d quark over u quark distribution inside the proton. While the calculations of structure functions and quark distributions are non-perturbative, they can be determined by the parameterization of experimental data. The understanding of Fn 2/Fp 2 and d/u as x -> 1 has a large influence on global fits and parameterization, and can be used to distinguish the non-perturbative models which give different predictions. However, Fn 2/Fp 2 measured using deuteron and hydrogen targets has large nuclear uncertainties at large x, because the nuclear effects in the deuteron become significant at large x. The MARATHON experiment, which ran in spring 2018 using the upgraded 11 GeV Jefferson Lab electron beam, employs a novel method. It performed deep inelastic scattering off tritium and helium-3 mirror nuclei to measure Fn 2/Fp 2 over the range x = 0.17 to x = 0.82. Since tritium and helium-3 are mirror nuclei, theoretical uncertainties largely cancel out in the ratio. The extracted Fn 2/Fp 2 has much smaller uncertainties compared with previous experiments at large x. The MARATHON experiment also provided results on the EMC effect for tritium and helium-3 nuclei. The results are considered essential for understanding the EMC effect. This thesis describes the MARATHON experiment, and presents results for Fn 2/Fp 2, and the EMC effect for tritium and helium-3.

Liu, Hanjie↗
DOE OSTI2020

Deep Exclusive p 0 Electroproduction Measured in Hall A at Jefferson Lab with the Upgraded CEBAF

Our understanding of the proton's structure in terms of its fundamental degrees of freedom, quarks and gluons, is incomplete. Generalized Parton Distributions (GPDs) can enhance our understanding of the partonic structure of the proton as it includes the correlation between the momentum and transverse position of partons. Exclusive electroproduction of the pion (ep→e'p'π 0 ) is one of the potential channels to access the GPDs. GPDs are defined at the Bjorken limit where the total amplitude can be factorized in a hard scattering of the photon off the quark and a soft part described by the GPDs. At sufficiently high Q 2 , a factorization scheme for exclusive meson production processes exists but only for longitudinally polarized virtual photons. As long as the longitudinally polarized virtual photons dominate the cross-section, we can have factorization and extract the GPDs. But the existing data for π 0 production at low Q 2 and low x B show the dominance of transverse part of the cross-section over the longitudinal one. Some transversity GPD models have found that an effective way to factorize the contribution from transversely polarized virtual photons to explain the dominance of σ T over σ L , shown by existing data. The experiment E12-06-114 ran in experimental Hall A during 2014-2016. The primary goal of this experiment was to measure the DVCS cross-section but the experimental configuration also allows the measurement of deep exclusive π 0 production. Taking advantage of the upgraded CEBAF at Jefferson Lab, we explored the high Bjorken-x (0.60) and high Q 2 regime for the first time. The extracted cross-section terms for π 0 production helps us to test the validity of QCD factorization in our kinematic regime. Our π 0 production results can be used to test and improve the existing models of proton structure. Cross-section results will be presented at two different Q 2 values for fixed x B of 0.60.

Karki, Bishnu↗
DOE OSTI2020

AGC-4 Experiment Irradiation Monitoring Data Qualification Final Report

The Graphite Technology Development Program ran a series of six experiments to quantify the effects of irradiation on nuclear grade graphite. This report focuses on the fourth experiment, Advanced Graphite Creep 4 (AGC 4). The Advanced Reactor Development (ARD) Technology Development Office (TDO) Program for research and development activities require documentation of qualified monitoring data to design and license the first high-temperature reactor nuclear plant. Qualified data meets the requirements for use as described in the experiment planning and quality assurance documents. Failed data do not meet the requirements and provide no useable information. Trend data may not meet all requirements but still provide some useable information. Use of Trend data requires assessment of how any deficiencies affect a particular use of the data. AGC-4 began with Advanced Test Reactor (ATR) Cycle 157D on May 30, 2015. After irradiating the graphite for two cycles, the capsule was removed from the reactor after ATR Cycle 158A, which ended on January 2, 2016, due to interference with another experiment. Irradiation was resumed with Cycle 162A on October 7, 2017 after the interfering experiment was removed from the reactor. AGC-4 irradiation completed by the end of Cycle 166B on January 10, 2020. Between Cycles 162A and 166B, AGC-4 capsule was removed from the reactor core during two Powered Axial Locator Mechanism (PALM) Cycles 163A and 165A to avoid overheating. All thermocouples (TCs) have functioned throughout the AGC-4 experiment. A total of 9,256,184 out of all 11,167,213 temperature records (or 82.9%) are Qualified for use by the ARD Program and 1,911,029 missing values are Failed. Argon, helium, and total gas flow data were within expected ranges, except only eight out-of-range values occurred during outages. A total of 21,240,627 out of all 22,288,5103 gas flow rates (or 92.8%) are Qualified for use by the ARD Program and 1,644,476 are Failed records mostly due to missing values. Discharge gas line moisture values were consistently low during full ATR power, except for 162B, when moisture content was increasing to more than 200 ppmv by the end of this cycle. During outages, moisture content reached as high as 700 ppmv. 1,089,941 out of a total of 1,107,331 moisture values (or 98.4%) are within the measurement range of the instrument and are Qualified for use by the ARD Program and 17,390 missing moisture values are Failed. Graphite creep specimens were subjected to one of three loads: 393, 491, or 589 lbf. For a brief period during Cycle 157D between 12:19 on June 2, 2015, and 08:23 on June 11, 2015, the load cells were wired incorrectly, resulting in missing stack load data. Missing stack loads were estimated from measured ram pressures using regression equations developed from the existing data from Cycle 157D. Estimated stack loads during this period are considered to be an accurate representation of actual load applied to the stacks. These loads deviate slightly from the planned loads. This deviation does not prevent the data from being Qualified for use but must be taken into account when analyzing the effect of load on creep. 6,095,728 out of a total of 6,403,985 moisture values (or 95.2%) are within the measurement range of the instrument and are Qualified for use by the ARD Program and 308,257 missing load values are Failed. Stack displacement increased consistently throughout the eight cycles, with total displacement reached highest value of 2.4 in by the end of irradiation. During ATR outages, a set of pneumatic rams raised the stacks of graphite creep specimens to ensure the specimens were not stuck within the test train. This stack raising was performed seven times throughout irradiation. All stacks were raised successfully each time. 4,744,974 out of a total of 6,094,513 displacement values (or 77.9%) are within the measurement range of the instrument and are Qualified for use by the ARD Program and 1,349,539 displacement values are Failed mostly due to missing values. Analyses were conducted on correlations between TCs to look for trends and step changes that might indicate instrument degradation or failure. Correlation analysis was used to identify instances when TCs form short circuits, referred to as virtual junctions, which result in TCs reporting temperatures from some location in the capsule other than the location where they were intended to read. No evidence of virtual junctions was found. Analyses were also conducted on control charts of temperature differences of between two TCs, which are expected to behave consistently throughout the entire irradiation period. Upward or downward trend over time indicates at least one TC in the pair was drifted. Examining control charts for all 66 possible pairs out of twelve TCs installed in the AGC-4 capsules reveal no clear drift failures occurred, except unstable behavior of two TCs in Zone 3, TC-7 and TC-8, over irradiation time. In conclus

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

SST ..kappa..-..omega.. Simulations of the Atmospheric Boundary Layer Including the Coriolis Effect

For large-scale structures in the atmospheric boundary layer (ABL), the Coriolis effect and near-wall behavior can have a meaningful impact. For example, both the Coriolis effect and blade boundary layer impact how muchpower wind farms produce. RANS simulations of the ABL typically use the ..kappa..-epsilon turbulence model, which has been developed to accurately capture the Coriolis effect but typically does not perform well near walls. The SST ..kappa..-..omega.. turbulence model performs well near walls but does not accurately capture the Coriolis effect. This work modifies SST ..kappa..-..omega.. to accurately model the Coriolis effect. We discuss the similarities and differences in how to modify ..kappa..-epsilon and SST ..kappa..-..omega.. for the Coriolis effect. Finally, we compare ..kappa..-epsilon and SST ..kappa..-..omega.. simulations of the ABL, including the Coriolis effect, with large eddy simulations and measurements.

atmospheric boundary layer↗
DOE OSTI2022

A Strong-QCD Regime Measurement of the Proton’s Spin Structure

The theory of the strong force, Quantum Chromodynamics (QCD) remains one of the most important ways to understand the fundamental properties of ordinary matter. However, at low momentum transfer Q 2 , in the regime where the strong force becomes extremely strong, our understanding of QCD for ordinary nucleons becomes hazy. Several cutting edge theories such as Chiral Perturbation Theory (χPT) and Lattice QCD have provided valuable predictions in this regime, but Lattice QCD has not yet extended predictions of many important quantities to this kinematic region, and Chiral Perturbation Theory has faced several important disagreements with experimental data for the neutron over the last several decades. It is therefore of extreme importance to have a benchmark of experimental data in the low energy regime for the proton’s behavior, as a test of leading theories for the behavior of QCD in this regime. The E08-027 (g2p) experiment ran at Jefferson Lab in 2012 with the goal of collecting this valuable data, and though I was still completing my undergraduate studies at the time, I became involved in the analysis in 2015 and built on the previous work to complete it and analyze the exciting results. This experiment achieved a high precision measurement of the spin structure functions g1 and g2 for the proton, quantities which describe the internal spin structure of the proton. These measurements were taken in the valuable low Q 2 region described above, and used to extract several moments of these spin structure functions which can be directly compared to the cutting-edge predictions of Chiral Perturbation Theory. Though the experiment’s timeline was such that I didn’t have a chance to work directly on the experimental setup, I had the opportunity to acquire hands-on experience working on a polarized target at UNH which is very similar to the crucial polarized target used in the g2p experiment. Full details of the g2p experiment and my experimental work at the University of New Hampshire are presented in this thesis, as well as a detailed description of the analysis process and the exciting benchmark results, which serve as a direct test of all current and future theories of QCD in the low-Q 2 regime.

Ruth, David↗
DOE OSTI2023

Transient Efficiency, Flexibility, and Reliability Optimization of Coal-Fired Power Plants - Final Program Review

This program developed an advanced model-based monitoring and model-predictive control algorithms for a coal fired power plant (CFPP), and deployed these algorithms in a real-time platform to demonstrate performance benefits for transient flexibility and plant operation efficiency. More specifically, the objectives were successfully achieved through a combination of (i) developing a high-fidelity transient plant model in Apros, which was used as a high-fidelity plant simulation between $100-50\% TMCR$ where TMCR denotes the turbine maximum continuous rating, i.e., baseload, (ii) developing a very fast physics-based reduced-order model (ROM) of the plant, which ran more than $100\times$ faster than real-time, enabling its use as real-time embedded model for model-based estimation (MBE) and model predictive control (MPC) (iii) implementing a real-time MBE based on ROM using a robust unscented Kalman filter (UKF) to continuously tune the ROM to match the measurements from high-fidelity Apros plant model despite significant plant-model mismatch, and thus, obtain a Digital Twin of the plant (iv) designing and implementing a real-time MPC with dual objectives of transient plant load tracking with high ramp rates and minimizing coal consumption, i.e., improving plant efficiency in the baseload-partload operation range of $100-50\% TMCR$. Each key element above was developed and tested individually, and has been reported in corresponding Topical Reports. Finally, all the individual elements were integrated in an overall closed-loop system, that was successfully tested in desktop Simulink test harness simulations with ROM or high-fidelity model as the plant. Thereafter, the Simulink implementation was used to auto-generate C-code and deploy as real-time Docker microservice containers in Linux, and validate that they can run in real-time in the hardware-in-the loop (HIL) setup and produce the same results as in Simulink. The results of the integrated simulation tests in Simulink as well as the real-time HIL deployment are documented in this final report, showing good load tracking for load ramps at $3-4\%/min$ ramp rates, and achieving up to $5.5\%$ reduction in coal relative to baseline operation at $50\% TMCR$ load. The desktop and HIL simulations show successful performance of the overall model based estimation and control solution and achieve the key objectives of the program for flexible, efficient and reliable operation of sub-critical coal fired power plants.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2022

Measurement of H(e,e’) and D(e,e’) Cross-Sections for Bloom-Gilman Quark-Hadron Duality Study

The E12-10-002 experiment which ran in 2018 in Hall C at the Thomas Jefferson National Accelerator Facility measured the inclusive electron scattering cross-sections for hydrogen and deuterium. This was the first experiment conducted with the new Super High-Momentum Spectrometer and accompanying instrumentation. These measurements extend the previous Hall C measurements up to Q2 = 16 GeV2. This dissertation presents the analyses of the new detectors and shows the preliminary results for F2 structure function and the study of quark-hadron duality in the four nucleon-resonance regions.

Sun, Abel↗
DOE OSTI2021

Neutron Structure Functions at Large Bjorken X

The E12-10-002 (F2) experiment ran in Hall C at JLab to extract the F2 structure functions of proton and neutron from inclusive H(e, e0) and D(e, e0) reactions in the high Bjorken x region. Those extractions put constraints on the Parton Distribution Functions PDFs (especially at large Bjorken x), and facilitate the study of Quark Hadron Duality. The F2 experiment measurements covers a large kinematic range in x from 0.2 to 1.0, and in Q2 from 4 to 16 GeV2 . While the proton F2 structure function has been studied extensively through inelastic electron-proton scattering, much less is known about neutron structure due to the unavailability of high density, free neutron targets. The BONuS12 experiment was proposed to measure the neutron F2 on a nearly free neutron within a weakly bound deuteron target via the spectator tagging method. Tagging the slow backward moving spectator protons minimizes both off-shell and final-state interaction effects, and the measured proton momentum is used to correct for the initial-state momentum of the struck neutron. The recoil detector is used to detect spectator protons with momenta 70 < ps < 150 MeV/c. This dissertation outlines the design, construction and testing of the recoil detector, Radial Time Projection Chamber (RTPC) and discusses the results of the H(e, e0) and D(e, e0) cross sections and the F2 structure functions obtained from the Hall C F2 data analysis.

Nadeeshani, Sooriyaarachchilage↗
DOE OSTI2023

Gas-Liquid Flow Modeling for Renewable Fuels Production

Aerobic/anaerobic and gas fermentation pathways have emerged as promising new technologies for the generation of renewable fuels/chemicals from biomass derived sugars, and mixtures of greenhouse/energy rich gas streams (CO2/CH4/H2/CO) via microbial action. Example pathways include sugars-to-ethanol conversion, biomethanation (CO2/H2 to CH4), biogas upgrading, CO fermentation and wet-waste conversion. Gas and liquid phase transport, mass-transfer, and mixing physics at large length scales can significantly affect microbial conversion rates, particularly when the microbial reaction requires a narrow set of conditions. These phenomena are difficult to study in small-scale bench-top reactors that are typically well-mixed. Predictive computational fluid dynamics (CFD) based simulations can therefore aid in the scale-up, design and optimization of these reactors. This work presents multiphase Euler-Euler CFD simulations of at-scale (~500 m3) bioreactors. Our mathematical model treats the gas and liquid as interpenetrating phases. This approach reduces the computational complexity of tracking individual gas bubbles that are several orders of magnitude smaller than reactor dimensions. We solve the Reynolds averaged Navier-Stokes (RANS) multiphase equations that account for phase and chemical species transport, interphase mass and momentum transfer and uses a phenomenological model for gas uptake by microbes. We use a customized solver derived from open-source CFD toolbox, OpenFOAM [1], to perform these simulations, which has been validated against small-scale reactors in our previous work [2]. There is currently a knowledge-gap regarding bubble-size distributions when using gas mixtures with vastly different properties, which can have a significant impact overall mass-transfer. For example, hydrogen bubbles are more buoyant compared to other relatively heavier gases (CO2/CH4/CO), resulting in a large distribution of residence times and bubble sizes. This work therefore develops a deeper understanding of bubble dynamics and interphase mass transfer in such heterogenous gas mixtures through well-resolved computational models. We use a population balance model (PBM) for bubble-size-distribution modeling that is validated against small-scale experiments in our solver with an uncertainty quantification study for bubble coalescence and break-up model parameters. Results pertaining to multiple simulations of gas-fermentation reactors are presented where gas mixtures with varying compositions of CO2/CH4/CO/H2 are imposed at the sparger boundaries. The spatio-temporal variations in bubble-size distribution and mass transfer coefficient are analyzed for varying superficial velocities and gas-compositions for varying sizes of bubble-column and airlift reactors. This work will also examine the performance of different reactor designs, viz. bubble column reactor, airlift reactor with an internal draft tube, and a stirred-tank reactor with Rushton impellers. Reactor mass-transfer coefficient, gas hold-up, and dissolved gas distribution are critically analyzed among reactors, and sensitivity studies pertaining to gas flow rates and reactor geometry will be presented. [1] Weller, H., Tabor, G., Jasak, H. and Fureby, C., A tensorial approach to computational continuum mechanics using object-oriented techniques, Computers in physics, 12, 6, 620--631, 1998. [2] Rahimi, M., Sitaraman, H., Humbird, D. and Stickel, J., Computational fluid dynamics study of full-scale aerobic bioreactors: Evaluation of gas-liquid mass transfer, oxygen uptake, and dynamic oxygen distribution, Chemical Engineering Research and Design, 139: 283-295.

BIOMASS FUELS↗
DOE OSTI2023

The Importance of Addressing Disagreements Between Nominal and Effective Treatments During Bat Mortality Minimization Validation Studies

As the wind energy industry grows, so too does our need for effective and low-cost bat mortality minimization solutions. Despite knowledge gaps in our understanding of what drives bats to collide with spinning turbines, minimization solutions have shown success during validation studies. Curtailment has been consistently effective, reducing bat mortality from 33%-79%, depending on the curtailment scenario and species present. Several validation studies have demonstrated deterrent solutions, specifically ultrasonic deterrents, to be effective, but overall, more variable compared to curtailment solutions. Further, a study combining curtailment with ultrasonic deterrents produced encouraging results, finding that adding ultrasonic deterrents to turbines that were designated to curtail, significantly reduced mortality rates compared to curtailment only turbines. Despite positive results, the strength of inference achieved with validation studies is limited by methodological constraints associated with appropriately assigning fatalities to treatments. Because of the temporal separation between when a bat collides with a turbine and when it is discovered during a ground-based carcass survey, we must assign carcasses to a treatment that ran during the prior night(s). This process may introduce errors at multiple stages. First, for studies that rotate treatments among turbines, mortality surveyors must be confident that carcasses are 'fresh' such that mortalities are correctly assigned to the treatment from the previous night(s). Second, we must reconcile any misalignments between how we assign treatments and how we implement treatments (e.g., turbines assigned with a deterrent treatment may be implemented as another treatment when deterrent devices are not operating as designed). Finally, for curtailment solutions, it is critical to recognize that measured effects are a function of not just the treatment as implemented but the proportion of the night the treatment is realized (e.g., if wind speeds are greater than the curtailment treatment cut-in speed for the entire night, we cannot expect there to me any differences in mortality reduction relative to control turbines). Using a dataset collected between June and October 2017, that rotated 3 treatments (Deterrent Only, 5 m/s Curtail only, and 5 m/s Curtail & Deterrent) and 1 control condition across 16 turbines each night, we explored the importance and implications of accounting for potential errors in assigning fatalities to treatments. We present preliminary results comparing the mortality associated with error-corrected treatments and control conditions highlighting how the measured effect of a treatment (mortality) greatly depends on site specific implementation.

bats↗