Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “3-D flow”

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 55 records · Page 3

UAE6 - Wind Tunnel Tests Data - UAE6 - Sequence K - Raw Data

Sequence K: Step AOA, Probes (P) This sequence was designed to quantify the 3-D blade static angle-of-attack response in the presence of rotational influences by varying the blade pitch angle. Sequence K used an upwind, rigid turbine with a 0° cone angle. The wind speeds ranged from 6 m/s to 20 m/s, and data were collected at yaw angles of 0° and 30°. The rotor rotated at 72 RPM. Blade and probe pressure measurements were collected. The teeter dampers were replaced with rigid links, and these two channels were flagged as not applicable by setting the measured values in the data file to -99999.99 Nm. The teeter link load cell was pre-tensioned to 40,000 N. The blade pitch angle ramped continuously at 0.18°/s over a wide range of increasing and decreasing pitch angles. A step sequence was also performed. The blade pitch was stepped 5°, the flow was allowed to stabilize, and the pitch angle was held for 8 seconds. Then the pitch angle step was repeated. Again, a wide range of pitch angles was obtained, both increasing and decreasing. The file lengths for this sequence varied from 96 seconds to 6 minutes, depending on the pitch angle range. Some short points were collected at 0° yaw and 3° pitch to verify the functionality of the instrumentation. The file name convention used the initial letter K, followed by two digits specifying wind speed, followed by two digits for yaw angle, followed by RU, RD, or ST, followed by the repetition digit. The angle-of-attack motion was differentiated by RU (ramp up), RD (ramp down), and ST (step down then step up). This sequence is related to sequences L and R.

17 WIND ENERGY↗

UAE6 - Wind Tunnel Tests Data - UAE6 - Sequence L - Raw Data

Sequence L: Step AOA, Parked (P) This sequence was designed to quantify the 3-D blade static angle-of-attack response in the absence of rotational influences by varying the blade pitch angle. This test sequence used an upwind, rigid turbine with a 0° cone angle. Wind speeds of 20 m/s and 30 m/s were used, and all data were collected at a yaw angle of 0°. The rotor was parked with the instrumented blade fixed at 0° azimuth, and the rotor lock was installed (see Appendix A). Blade and probe pressure measurements were collected. The teeter dampers were replaced with rigid links, and these two channels were flagged as not applicable by setting the measured values in the data file to -99999.99 Nm. The teeter link load cell was pre-tensioned to 40,000 N. The blade pitch angle ramped continuously at 0.18°/s over a wide range of increasing and decreasing pitch angles. A step sequence was also performed. The blade pitch was stepped 5°, the flow was allowed to stabilize, and the pitch angle was held for 8 seconds. Then the pitch angle step was repeated. Again, a wide range of pitch angles was obtained, both increasing and decreasing. The file lengths for this sequence varied from 5 to 10 minutes, depending on the pitch angle range. The file name convention used the initial letter L followed by two digits specifying wind speed, followed by 00 for yaw angle, followed by RU, RD, or ST, followed by the repetition digit. The angle-of-attack motion was differentiated by RU (ramp up), RD (ramp down), and ST (step down, then step up). This sequence is related to sequences K and R.

17 WIND ENERGY↗

UAE6 - Wind Tunnel Tests Data - UAE6 - Sequence R - Raw Data

Sequence R: Step AOA, No Probes (P) This sequence was designed to quantify the effect of the five-hole probes on the 3-D blade static angle-of-attack response in the presence of rotational influences by repeating Sequence K without five-hole probes. This test sequence used an upwind, rigid turbine with a 0° cone angle. The wind speeds ranged from 6 m/s to 20 m/s, and data were collected at yaw angles of 0° and 30°. The rotor rotated at 72 RPM. Blade pressure measurements were collected. The five-hole probes were removed and the plugs were installed. Plastic tape 0.03-mm-thick was used to smooth the interface between the plugs and the blade. The teeter dampers were replaced with rigid links, and these two channels were flagged as not applicable by setting the measured values in the data file to –99999.99 Nm. The teeter link load cell was pre-tensioned to 40,000 N. During post-processing, the probe channels were set to read –99999.99. The blade pitch angle ramped continuously at 0.18°/s over a wide range of increasing and decreasing pitch angles. A step sequence was also performed. The blade pitch was stepped 5°; the flow was allowed to stabilize; and the pitch angle was held for 5 seconds. Then the pitch angle step was repeated. Again, a wide range of pitch angles was obtained, both increasing and decreasing. The file lengths for this sequence varied from 96 seconds to 6 minutes, depending on the pitch angle range. Some short points were collected at 0° yaw and 3° pitch to ascertain the functionality of the instrumentation and repeatability over time. The file name convention used the initial letter R, followed by two digits specifying wind speed, followed by two digits for yaw angle, followed by RU, RD, or ST, followed by the repetition digit. The angle of attack motion was differentiated by RU (ramp up), RD (ramp down), and ST (step down, then step up). This sequence is related to Sequences K and L.

17 WIND ENERGY↗

Investigation of thermal hydraulic behavior of the High Temperature Test Facility's lower plenum via large eddy simulation

A high-fidelity computational fluid dynamics (CFD) analysis was performed using the Large Eddy Simulation (LES) model for the lower plenum of the High–Temperature Test Facility (HTTF), a ¼ scale test facility of the modular high temperature gas-cooled reactor (MHTGR) managed by Oregon State University. In most next–generation nuclear reactors, thermal stress due to thermal striping is one of the risks to be curiously considered. This is also true for HTGRs, especially since the exhaust helium gas temperature is high. In order to evaluate these risks and performance, organizations in the United States led by the OECD NEA are conducting a thermal hydraulic code benchmark for HTGR, and the test facility used for this benchmark is HTTF. HTTF can perform experiments in both normal and accident situations and provide high-quality experimental data. However, it is difficult to provide sufficient data for benchmarking through experiments, and there is a problem with the reliability of CFD analysis results based on Reynolds–averaged Navier–Stokes to analyze thermal hydraulic behavior without verification. To solve this problem, high-fidelity 3-D CFD analysis was performed using the LES model for HTTF. It was also verified that the LES model can properly simulate this jet mixing phenomenon via a unit cell test that provides experimental information. As a result of CFD analysis, the lower the dependency of the sub-grid scale model, the closer to the actual analysis result. In the case of unit cell test CFD analysis and HTTF CFD analysis, the volume-averaged sub-grid scale model dependency was calculated to be 13.0% and 9.16%, respectively. As a result of HTTF analysis, quantitative data of the fluid inside the HTTF lower plenum was provided in this paper. As a result of qualitative analysis, the temperature was highest at the center of the lower plenum, while the temperature fluctuation was highest near the edge of the lower plenum wall. The power spectral density of temperature was analyzed via fast Fourier transform (FFT) for specific points on the center and side of the lower plenum. FFT results did not reveal specific frequency-dominant temperature fluctuations in the center part. It was confirmed that the temperature power spectral density (PSD) at the top increased from the center to the wake. The vortex was visualized using the well-known scalar Q-criterion, and as a result, the closer to the outlet duct, the greater the influence of the mainstream, so that the inflow jet vortex was dissipated and mixed at the top of the lower plenum. Additionally, FFT analysis was performed on the support structure near the corner of the lower plenum with large temperature fluctuations, and as a result, it was confirmed that the temperature fluctuation of the flow did not have a significant effect near the corner wall. In addition, the vortices generated from the lower plenum to the outlet duct were identified in this paper. It is considered that the quantitative and qualitative results presented in this paper will serve as reference data for the benchmark.

97 MATHEMATICS AND COMPUTING↗

Infrasonic directivity of monopole, dipole and bipole ground-surface reflected sources

Infrasound (acoustic waves below 20 Hz) can be used to detect, locate and quantify activity in the atmosphere such as volcanic eruptions and anthropogenic explosions. Attempts to quantify volcanic eruption parameters such as exit velocity, plume height and mass flow rate using infrasound data depend strongly on assumptions of the acoustic source type. Infrasonic sources may produce omnidirectional or directional wavefields, while propagation effects, such as interaction with topography, can induce further wavefield directivity that is measured by field instrumentation. Limited sampling of these wavefields can hinder our ability to infer the underlying source, and thus our understanding of the eruption characteristics. Equivalent sources are often used to represent acoustic source mechanisms and resultant wavefields. In this study, we review equivalent acoustic sources as they pertain to infrasonic scale and wavelengths commonly encountered in very local (⁠<5 km range) geophysical field deployments. We highlight the equivalent infrasonic bipole source that can be induced by ground-reflection of an elevated monopole; we are not aware of any prior infrasound studies that use the bipole source concept. We use analytical and numerical methods to explore source directivity of monopole, dipole and bipole ground-reflected sources at infrasonic frequencies as well as the additional directivity complications introduced by interactions with topography. We illustrate that for typical volcano-infrasound wavelengths, increasing height above the ground as well as increasing source frequency leads to increased wavefield directivity. Numerical modelling using a simple omnidirectional monopole source embedded in topography further illustrates that both horizontal and vertical infrasound directionality can be induced by topography at the distance scales appropriate for local volcano infrasound monitoring. Information summarized in this analytical and numerical exploration of infrasound directivity may be used to help guide future volcano-infrasound field deployments intended to estimate source parameters or quantify wavefield directivity. Analytic solutions for simple whole-space or half-space atmospheres provide useful formulations for planning or initially analysing geophysical field-scale experimental data; however, especially at very local distances from the source (⁠<5 km), 3-D simulations are necessary to account for complex topography commonly encountered in volcano-infrasound applications.

Infrasound↗

Biofouling of filtration membranes in wastewater reuse: In situ visualization with confocal laser scanning microscopy

The filtration membranes utilized in advanced treatment of municipal wastewater are prone to biological fouling. Permeability loss and water recovery limitations due to fouling necessitate additional pretreatment and membrane cleaning, which raise the cost of water reuse. Better fouling mitigation calls for a deeper understanding of biofilm behavior on water reuse membranes. This study uses confocal microscopy to monitor biofouling under realistic conditions. To simulate water reuse, a miniature flow cell with microfiltration and ultrafiltration membranes was operated under typical filtration conditions using secondary effluent from California's East Bay Municipal Utility District. The growing biofilms were stained for extracellular polymeric substances and live and dead cells, and they were monitored in situ under filtration conditions with a confocal laser scanning microscope to produce 2-D and 3-D images. In situ image series revealed biofilm growth and removal processes including internal and external fouling, extracellular polymeric substance production, wrinkling, delamination, and film re-deposition. Additionally, changes in biofilm morphology between in situ and ex situ images highlighted limitations of ex situ imaging. Here, these results provide insight into the physical and biological mechanisms of biofouling in hopes of informing the development of improved techniques for biofouling mitigation.

ultrafiltration↗

Scale interactions and anisotropy in Rayleigh–Taylor turbulence

We study energy scale transfer in Rayleigh–Taylor (RT) flows by coarse graining in physical space without Fourier transforms, allowing scale analysis along the vertical direction. Two processes are responsible for kinetic energy flux across scales: baropycnal work Λ, due to large-scale pressure gradients acting on small scales of density and velocity; and deformation work Π, due to multiscale velocity. Our coarse-graining analysis shows how these fluxes exhibit self-similar evolution that is quadratic-in-time, similar to the RT mixing layer. We find that Λ is a conduit for potential energy, transferring energy non-locally from the largest scales to smaller scales in the inertial range where Π takes over. In three dimensions, Π continues a persistent cascade to smaller scales, whereas in two dimensions Π rechannels the energy back to larger scales despite the lack of vorticity conservation in two-dimensional (2-D) variable density flows. This gives rise to a positive feedback loop in 2-D RT (absent in three dimensions) in which mixing layer growth and the associated potential energy release are enhanced relative to 3-D RT, explaining the oft-observed larger α values in 2-D simulations. Despite higher bulk kinetic energy levels in two dimensions, small inertial scales are weaker than in three dimensions. Moreover, the net upscale cascade in two dimensions tends to isotropize the large-scale flow, in stark contrast to three dimensions. Furthermore, our findings indicate the absence of net upscale energy transfer in three-dimensional RT as is often claimed; growth of large-scale bubbles and spikes is not due to ‘mergers’ but solely due to baropycnal work Λ.

42 ENGINEERING↗

Computational Study on the Effect of Thermal Boundary Conditions and Axial Aspect Ratio on Catalytic Oxidative Coupling of Methane

Packed-bed reactor models with a large reactor-diameter-to-particle-diameter ratio (N > 10) are commonly used in catalysis research. Their use is underpinned by strong simplifying assumptions including unidirectional plug flow, negligible radial gradients, small axial diffusion, and isothermal operation. However, the complex exothermic–endothermic nature of reactions like oxidative coupling of methane (OCM) coupled with heat transfer from the catalyst surface and convective heat losses from the reactor wall provides an interesting case study to evaluate packed-bed assumptions and understand axial and radial temperature variations that may impact conversion and selectivity. This computational study investigates the effect of thermal boundary conditions on catalyst performance in laboratory OCM reactors (N = ~13) through computational fluid dynamics (CFD) simulation. Further, the 1-D model typically used for modeling plug flow reactors assumes negligible radial temperature and species concentration gradients. This assumption is evaluated for OCM reactors over a range of geometric axial aspect ratios, defined as the reactor length divided by the diameter, comparing the 1-D with 3-D CFD results. 1-D model predictions of catalyst performance start to deviate from the 3-D model predictions as the axial aspect ratio decreases, specifically as diameter increases at constant length. In the work, the 1-D model assumption is only found to be valid for axial aspect ratios higher than 1.6, a result considered to be relevant for all exothermic catalytic reactors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Convective Momentum Transport and Its Impact on the Madden‐Julian Oscillation in E3SM‐MMF

Abstract Convective momentum transport (CMT) is the process of vertical redistribution of horizontal momentum by small‐scale turbulent flows from moist convection. Traditional general circulation models (GCMs) and their multiscale modeling framework (MMF) versions poorly represent CMT due to insufficient information of subgrid‐scale flows at each GCM grid. Here the explicit scalar momentum transport (ESMT) scheme for representing CMT is implemented in the Energy Exascale Earth System Model‐Multiscale Modeling Framework (E3SM‐MMF) with embedded 2‐D cloud‐resolving models (CRMs), and verified against E3SM‐MMF simulations with 3‐D CRMs and observations. The results show that representing CMT by ESMT helps reduce climatological mean precipitation model bias over the western Pacific and the ITCZ regions, which is attributed to the weakened mean easterlies over the Pacific. Also, CMT from simulations with 2‐D and 3‐D CRMs impose a similar impact on Kelvin waves by reducing their variability and slowing down their phase speed, but opposite impacts on the Madden‐Julian Oscillation (MJO) variability. The ESMT scheme readily captures the climatological mean spatial patterns of the zonal and meridional components of CMT and their variability across multiple time scales, but shows some differences in estimating its magnitude. CMT mainly affects the MJO by decelerating its winds in the free troposphere, but accelerates its near‐surface winds. This study serves as a prototype for implementing CMT scheme in the MMF simulations, highlighting its crucial role in reducing model bias in mean state and spatiotemporal variability.

54 ENVIRONMENTAL SCIENCES↗

An experimentally informed design process for future inertial confinement fusion facilities

The achievement of ignition in the laboratory has renewed interest in defining the requirements for a future high-gain inertial confinement fusion (ICF) facility. Our best chance of predicting future ICF performance is with 3-D radiation hydrodynamic simulations that have been benchmarked against experimental data, but their high computational cost is prohibitive for use in practical design studies. We introduce a hierarchical approach where 3-D simulations are tuned to match experimental measurements and used to train 3-D degradation models in 1-D simulations allowing for accurate predictions over the entire OMEGA direct-drive database. A genetic algorithm was used in combination with the trained 1-D simulations to search for optimal direct-drive implosion designs at driver energies ranging from 20 kJ to 10 MJ. As the fidelity of 3-D codes improves, this approach will provide a viable experimentally informed tool for defining the next ICF facility.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Multiphysics Analysis of Li Cooled Divertor Substrate During Loss of Coolant Accident (LOCA)

In the ongoing study of potential designs for liquid metal (LM) plasma-facing components (PFCs), so-called “slow” and “fast” Li flow divertor concepts are under investigation. In the previous studies on design and analysis of the slow Li flow divertor and comparison with the fast Li flow divertor, the magnetohydrodynamics (MHD)/heat transfer effects of the Li flowing inside the substrate as a second coolant were comprehensively investigated under the normal steady-state operation conditions. Here, in the present study, the multiphysics analysis is extended to the unsteady abnormal divertor scenario where the Li layer on top of the substrate does not provide full coverage or even totally disappears for a certain period of time, so that the substrate becomes directly exposed to the incident high plasma heat flux. Such an unwanted event may happen regardless of the concept of the divertor and is worth detailed investigations, typically referred to as a loss of coolant accident (LOCA). To address this situation, a simplified scoping analysis is conducted first in 2-D, and then an integrated 3-D modeling is performed using a time-dependent multiphysics model in COMSOL Multiphysics that integrates LM MHD, heat transfer, and solid mechanics. The main goal is to evaluate conditions under which the major material limits, such as the maximum allowable temperature, stress, and displacement of the substrate, can still be met. It was shown that the maximum time over which the substrate of RAFM steel can retain structural integrity during the LOCA is around 0.2 ∼ 0.3 s. Any divertor concept that utilizes RAFM steel as a substrate material and liquid Li as a second coolant should take such a permitted time into consideration.

divertor↗

Physical properties, internal structure, and the three‐dimensional petrography of CI chondrites

physical properties and the nature of their breccation, we investigated nine samples of the Ivuna and Orgueil CI chondrites ranging in size from 1 mm to 4 cm in approximate diameter. The combined mass of unique material investigated in this work is 113 g. For our investigations, we use ideal gas pycnometry, 3-D laser scanning, x-ray computed microtomography (μCT), and accompanying digital data extraction techniques. We found that the bulk density of the samples ranged from 1.61 to 2.10 g cm −3 . Larger samples tend to have a lower bulk density. Grain density (ranging from 2.44 to 2.55 g cm −3 ) is significantly less variable than the bulk density in our samples and the quantity of porosity (ranging from 14.6% to 33.8%) is the dominant factor in determining the bulk density of CI chondrite material. Our μCT results show that the visible porosity across all sizes of our CI chondrite samples is in the form of cracks, but these cracks can account for less than two-thirds of the porosity in the CI chondrites. Other porosity is not visible, even at μCT resolutions of 2.7 μm voxel edge −1 and we conclude that it is sub-micron in nature. It is not clear if the cracks seen in our samples are indigenous to the chondrites or are a result of terrestrial processes. We also find that the CI chondrites are excellent examples of the fractal-like nature of brecciation, where clasts can be observed at all scales we imaged. The breccias are composed of sub-equant-shaped and sub-rounded-textured clasts like melt-free impact breccias on other solar system bodies. From our μCT volume and digital data extraction, we determine that the Ivuna CI chondrite breccia is organized: the mostly sub-equant clasts within our ~2 cm chunk of Ivuna have a mean diameter of 1.33 mm and their aligned longest axes define a lineation structure. We speculate that the lineation was imparted after fragmentation of the clasts by slight shear on the parent asteroid which could be the result of seismic-related granular flow or mild non-axial impact-related compaction. These data will help to place returned asteroidal material from asteroids 162173 Ryugu and 101955 Bennu and the CI chondrites into a mutual geological context.

CI chondrite↗

Acoustic Codes in 2D Spherical Coordinate

Finite-difference methods are widely used to simulate infrasound propagation in the atmosphere. Flexibility of finite-difference scheme allows implementation of highly heterogeneous media for sound propagation as well as complex source models for sound generation. While full 3-D finite-difference methods have been utilized for local infrasound propagation with pronounced topography, 2-D modeling approach has been preferred for regional and global propagation as full 3-D methods generally require enormous computational resources. Infrasound propagation is often simulated with a second-order finite difference scheme. This lowest-order finite-difference scheme is robust and straightforward to implement complex boundary conditions, but the solution includes large error with numerical dispersion and dissipation. This large numerical error may make the second-order finite-difference not optimal for long range infrasound propagation modeling as the numerical dispersion degrades the accuracy of the solution unacceptably. Here, we developed a high-order finite-difference solver for long-range infrasound simulation. The high-order scheme is particularly popular for linear wave modeling in aeroacoustics owing to its low-dispersive and low-dissipative behavior. We develop and evaluate a high-order finite difference scheme in 2-D axisymmetric coordinates. The axisymmetry allows to approximate 3-D spherical sound propagation and amplitude attenuation by a 2-D method. AC2Dr is developed to simulate infrasound propagation in realistic atmosphere, but can be used for linear acoustic waves in general materials with background flow. AC2Dr in an axisymmetric coordinates allows for spherical radiation of acoustic waves from compact sources.

Sjogreen, Bjorn↗

A Process with Decoupled Absorber Kinetics and Solvent Regeneration through Membrane Dewatering and In-Column Heat Transfer (Final Report)

This report summarizes the work conducted on project DE-FE0031604 where University of Kentucky Center for Applied Energy (UK CAER) has validated its intensified CO 2 capture process through substantial enhancements to the kinetics of the absorption process and energy reductions by absorber temperature profile modification, dewatering and heat integration technologies for achieving significant capital and operating cost reductions. To address DOE’s objective of improving post-combustion CO 2 capture technology and reducing associated cost, UK CAER employed an intensified process which combined three key aspects targeted at overcoming inherent limitations or barriers in the conventional CO 2 capture and desorption process. The process designed to be independent of the type of solvent used, included (1) the use of 3-D printed two-channel structured packing material to control the temperature profile and increase the CO 2 absorption rate in the absorber, (2) a zeolite membrane dewatering unit for dewatering of the carbon-rich solvent to decouple solvent concentration needs for CO 2 absorption and desorption, and (3) a rich-split feed with two-phase flow heat transfer prior to the stripper that provided a secondary point of vapor generation to provide energy savings in steam extraction and solvent regeneration. The project was executed over two budget periods. This involved testing of individual process components which included the advanced heat transfer packing and the dewatering membrane on UK CAER’s 30 liter per minute (L/min) CO 2 (3” Column) capture bench unit with simulated flue gas in the first budget period. Subsequent scaled-up testing of these components together with the split-feed configuration were also tested in UK CAER’s 0.1 MWth CO 2 capture unit with coal-derived flue gas in the second-budget period. Long term studies were done during this period to assess process and solvent performance over extended duration. Project partners Lawrence Livermore National Laboratory (LLNL) and Media & Process Technology (MPT) led the development efforts for the advanced packing material and dewatering membrane respectively. Data from the long term testing was used as input for an environmental, health and safety (EH&S) assessment for the process and scaled technology performed by ALL4 LLC. Trimeric Corporation also completed a techno-economic analysis (TEA) for the UK CAER technology which was compared to the DOE reference Case B12B. Tests on the 3” column capture unit showed that the advanced heat transfer packing could be used to lower the bulge temperature in the absorber, and this was also proven in the scaled testing in the 0.1 MWth CO 2 capture unit. The bulge temperature could be lowered by >10 °C, changing the temperature profile in the absorber, and showed potential to enhance absorption with the ability to tailor the profile to provide conditions suitable for a solvent’s properties and kinetics. Conditions for short term evaluation of a 19” zeolite dewatering membrane on the 3” column capture unit yielded desirable fluxes and sustained rejection rates of >80%. However, for the scaled testing of six membrane modules consisting of 21 parts of 31-inch-long membrane tubes in each bundle (surface area 0.3 m 2 ), over a more extended duration, similar rejection rates could not be achieved. With the split-feed of the rich stream to the stripper, improved heat recovery minimized waste heat exiting the top of the stripper. The stripper exhaust temperatures could be reduced by >10 °C; reducing the amount of water vaporization contributed to lowering the regeneration energy by ~ 15%. The energy benefit could be sustained from the long term monitoring of the solvent performance. The solvent properties were not significantly impacted over the long-term operations. The benefits of the UK CAER process demonstrated experimentally were mostly validated from the TEA comparing a commercial scale application of the technology to DOE reference Case B12B. The cost of CO 2 capture for the UK CAER technology was estimated to be ~$34.97/tonne of CO 2 captured; a reduction of 23% compared to Case B12B. The increase in cost of electricity was also shown to be 16% lower than that of Case B12B. The total parasitic demand was also shown to be 11% lower. The key drivers for the benefits are a result of the process intensification approaches employed in the UK CAER technology for enhanced solvent performance, effective heat recovery and improved energy performance. The EH&S assessment did not find any major environmental concerns or barriers to the full scale implementation of the technology.

20 FOSSIL-FUELED POWER PLANTS↗

CFD Simulations to Support Pronghorn Modeling of a Molten Salt Fast Reactor

To assist with the deployment of next generation molten salt reactors (MSRs), the DOE NEAMS program is developing advanced simulation capabilities. As part of this work, the high-fidelity CFD tool, Nek5000 was used to inform model development for the engineering scale tool, Pronghorn. The particular case chosen was the EVOL molten salt fast reactor concept. This concept uses an open core design, which is well known to be sensitive to the particular geometry and can have complex flow behavior that varies across a range of conditions. This makes it an idea candidate for demonstration of the Hi2Lo concept of using higher fidelity models to benchmark lower fidelity models. Simulations of the EVOL design were performed using both LES (high-fidelity) and RANS (moderate-fidelity) in Nek5000. The LES results indicated discrepancies in 2-D axisymmetric RANS, indicating that the core is better modeled in RANS with a 3-D wedge representation of a portion of the core, which can accurately account for the effect of the inlet channels. Finally, results from the 3-D RANS were used to modify an existing turbulence model in Pronghorn. Improvement in the calibrated Pronghorn model was demonstrated. Future work is suggested to focus on expanding the LES calculation to provide a more direct point of comparison for both the RANS and Pronghorn models as well as incorporating more of the relevant MSR physics (such as delayed neutron precursor tracking) into the model.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

UAE6 - Wind Tunnel Tests Data - UAE6 - Sequence N - Raw Data

Sequence N: Sin AOA, Rotating (P) This sequence was designed to quantify the blade 3-D unsteady aerodynamic response in the presence of rotational influences by varying blade pitch angle. Test sequence N used an upwind, rigid turbine with a 0° cone angle. The wind speed was 15 m/s, and the yaw angle was 0° throughout the sequence. The rotor rotated at 72 RPM. Blade and probe pressure measurements were collected. The teeter dampers were replaced with rigid links, and these two channels were flagged as not applicable by setting the measured values in the data file to –99999.99 Nm. The teeter link load cell was pre-tensioned to 40,000 N. Sinusoidal angle-of-attack histories were designed to emulate 2-D S809 dynamic data obtained at the Ohio State University (OSU) wind tunnel (Reuss Ramsey et al. 1995) and the 1P yawed flow angle-of-attack variations the turbine could be expected to encounter in routine operation. The reduced frequency (K), mean angle of attack (αm), and oscillation amplitude (αω) were computed for each of the five primary span locations. The file lengths varied according to the time required to complete 40 oscillations at the specified frequency. The file name convention used the sequence designation N, followed by two digits for the span location (30, 47, 63, 80, or 95), followed by two digits representing the condition specified in the test matrix, followed by the repetition digit. Due to time constraints, most of the test conditions corresponding to the 63% span location were eliminated. Test points with a static pitch angle of 3° were collected throughout the sequence to ascertain instrumentation fidelity.

17 WIND ENERGY↗

UAE6 - Wind Tunnel Tests Data - UAE6 - Sequence O - Raw Data

Sequence O: Sin AOA, Parked (P) This sequence was designed to quantify the blade 3-D unsteady aerodynamic response in the absence of rotational influences by varying blade pitch angle. This test sequence used an upwind, rigid turbine with a 0° cone angle. The wind speed ranged from 18.9 m/s to 39.3 m/s, and the yaw angle was 0° throughout the sequence. The rotor was parked with the instrumented blade fixed at 0° azimuth, and the rotor lock was installed (see Appendix A). Blade and probe pressure measurements were collected. The teeter dampers were replaced with rigid links, and these two channels were flagged as not applicable by setting the measured values in the data file to -99999.99 Nm. The teeter link load cell was pre-tensioned to 40,000 N. Sinusoidal angle-of-attack histories were designed to emulate 2-D S809 dynamic data obtained at the OSU wind tunnel (Reuss Ramsey et al. 1995) as well as 1P yawed flow angle-of-attack variations the turbine could be expected to encounter in routine operation. The reduced frequency (K), mean angle of attack (αm), and oscillation amplitude (αω) were computed for each of the five primary span locations. The file lengths varied according to the time required to complete 40 oscillations at the specified frequency. The file name convention used the sequence designation O, followed by two digits for the span location (30, 47, 80, or 95), followed by two digits representing the condition specified in the test matrix, followed by the repetition digit. Due to time constraints, the test conditions corresponding to the 63% span location were eliminated. Test points with a static pitch angle of 3° were collected throughout the sequence to ascertain instrumentation fidelity.

17 WIND ENERGY↗

Capturing plume behavior in complex terrain: an overview of the Nevada National Security Site Meteorological Experiment (METEX21)

METEX21 was an atmospheric tracer release experiment executed at the Department of Energy’s Nevada National Security Site (NNSS) in the southwestern U.S to study terrain-induced wind and thermodynamic conditions that influence local-scale (<5-km) plume transport under varying atmospheric forcing conditions. Meteorological observations were collected using 10-m tall meteorological towers, 2-m tall tripods with 3-d sonic anemometers, a 3-m tall eddy covariance flux tower, Doppler profiling lidars, Doppler scanning lidars, weather-balloon launched radiosondes, and a tethered balloon equipped with wind, temperature, and aerosol sensors at heights up to 800 m. A smoke tracer was released along three transects in the horizontal and vertical directions and observed with video cameras, aerosol sensors and lidars (via aerosol backscatter). The observations showed evidence of large-scale/synoptic transience as well as local-scale upslope and downslope flows, along-axis valley flows, recirculation eddies on leeward slopes, and periods of strong shear and veer aloft. The release days were classified as either synoptically-driven or locally-driven, and a single case day is presented in detail for each. Synoptically-forced days show relatively narrow smoke plumes traveling down the valley from north to south (with the predominant wind direction), with little deviation in transport direction regardless of the elevation or ground locations of the smoke releases, except near the presence of leeside recirculation eddies. Locally-forced days exhibit a wider range of plume behavior due to the combination of thermally-induced valley and slope flows, which are often flowing in different cardinal directions, and wind shear found aloft at higher altitudes and elevations. We saw evidence of smoke lofting on top of the mesas due to strong upslope flows on these days. A major finding of this experiment was the effectiveness of scanning lidars to measure 2-dimensional plume transport out to a 2–3 km distance; much farther than could be visibly observed. METEX21 was the first of three planned tracer experiments at NNSS, and future experiments will incorporate multiple tracers to improve individual plume identification so that finer resolution flow details can be attained from these measurements, as well as deploy a larger suite of meteorological instrumentation, including more temperature profiling data.

54 ENVIRONMENTAL SCIENCES↗