Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Re”

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,189 records · Page 66
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Spectroscopic, Magnetic, and Computational Investigations on a Series of Rhenium(III) Cyclopentadienide β-diketiminate Halide and Pseudohalide Complexes

Here the low-valent rhenium(I) salt Na[Re(η 5 -Cp)(BDI)] (BDI = N,N'-bis(2,6-diisopropylphenyl)-3,5-dimethyl-β-diketiminate) was shown to react with various halide and chalcogenolate reagents, leading to the isolation of a series of rhenium(III)-halide and -pseudohalide complexes: Re(X)(η 5 -Cp)(BDI) (1-X, X = F, Cl, Br, I; 2, X = OTf) and Re(ER)(η 5 -Cp)(BDI) (3-ER, ER = SBn, SeBn, TePh). The 1 H NMR spectra of these complexes displayed sharp resonances shifted several ppm from typical diamagnetic regions, as well as distinct chemical shift trends down both the halide (with the exception of F) and chalcogenolate series, with both Cp and BDI backbone peaks moving downfield with increasing atomic number. Subsequent magnetic susceptibility measurements of the rhenium(III) halides 1-Cl, 1-Br, and 1-I indicated that these complexes display temperature-independent paramagnetism (TIP), with χ TIP values of 7.41(44) × 10 -4 to 1.50(51) × 10 -3 cm 3 mol -1 . Multireference complete active space self-consistent field (CASSCF) computations incorporating spin-orbit state mixing revealed small energy separations (1.7-3.0 kcal/mol) between thermally isolated ground states and the first excited spin-orbit states for the rhenium halides, confirming that low-lying excited states are responsible for the observed TIP behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Luminescent Properties of DOBDC Containing MOFs: The Role of Free Hydroxyls

A novel metal–organic framework (MOF), Mn-DOBDC, has been synthesized in an effort to investigate the role of both the metal center and presence of free linker hydroxyls on the luminescent properties of DOBDC (2,5-dihydroxyterephthalic acid) containing MOFs. Co-MOF-74, RE-DOBDC (RE–Eu and Tb), and Mn-DOBDC have been synthesized and analyzed by powder X-ray diffraction (PXRD) and the fluorescent properties probed by UV–Vis spectroscopy and density functional theory (DFT). Mn-DOBDC has been synthesized by a new method involving a concurrent facile reflux synthesis and slow crystallization, resulting in yellow single crystals in monoclinic space group C2/c. Mn-DOBDC was further analyzed by single-crystal X-ray diffraction (SCXRD), scanning electron microscopy–energy-dispersive spectroscopy (SEM-EDS), and photoluminescent emission. Results indicate that the luminescent properties of the DOBDC linker are transferred to the three-dimensional structures of both the RE-DOBDC and Mn-DOBDC, which contain free hydroxyls on the linker. In Co-MOF-74 however, luminescence is quenched in the solid state due to binding of the phenolic hydroxyls within the MOF structure. Finally, Mn-DOBDC exhibits a ligand-based tunable emission that can be controlled in solution by the use of different solvents.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Operando XAS and DFT Uncover Structure-Performance Relationships in Re/TiO 2 for Selective CO 2 Hydrogenation to Methanol

The conversion of CO 2 into value-added chemicals, such as methanol, offers a promising pathway toward a renewable energy future. However, a precise kinetic control and a highly selective catalyst are necessary to overcome the thermodynamic preference for CO 2 hydrogenation to methane. Rhenium-based catalysts, particularly Re/TiO 2 , demonstrate high activity and selectivity for methanol under high-pressure conditions. For example, at 100 bar and 200 °C, a methanol selectivity of 97−99% was obtained. Catalysts with 1 wt % Re and 5 wt % Re/ TiO 2 were used to study the effect of cluster sizes. At 250 °C, the 1 wt % catalyst achieves 97% selectivity at 23% conversion, whereas 5 wt % Re/TiO 2 achieves 74% selectivity at 40% conversion, corresponding to a drop in space-time yield from 65 to 16 g CH 3 OH ·g Re −1 ·h −1 , respectively. X-ray absorption spectroscopy provided insights into the structure of the active sites, while density functional theory calculations revealed the effects of cluster size on the energy barriers for H 2 activation, CH 3 OH dissociation, and CH 3 OH desorption, all of which directly influence conversion and selectivity. These results underscore the importance of balancing cluster size for optimal catalyst performance and provide insights into the design of efficient and selective catalysts for renewable methanol production.

XAS↗
DOE OSTI2020

Electron Redistribution within the Nitrogenase Active Site FeMo-cofactor During Reductive Elimination of H2 to Achieve N=N Triple-Bond Activation

Nitrogen fixation by nitrogenase begins with accumulation of four reducing equivalents at the active-site FeMo-cofactor (FeMo-co), generating a state (denoted E4(4H)) with two [Fe-H-Fe] bridging hydrides. Recently, photolytic reductive elimination (re) of the E4(4H) hydrides showed that enzymatic re of E4(4H) hydride yields an H2-bound complex (E4(H2,2H)), in a process corresponding to a formal 2-electron reduction of the metal-ion core of FeMo-co. The resulting electron-density redistribution from Fe-H bonds to the metal ions themselves enables N2 to bind with concomitant H2 release, a process illuminated here by QM/MM molecular dynamics simulations. What is the nature of this redistribution? Although E4(H2,2H) hasn’t been trapped, cryogenic photolysis of E4(4H) provides a means to address this question. Photolysis of E4(4H) causes hydride-re with release of H2, generating doubly-reduced FeMo-co (denoted E4(2H)*), the extreme limit of the electron-density redistribution upon formation of E4(H2,2H). Here we examine the doubly-reduced FeMo-co core of the E4(2H)* limiting-state by 1H, 57Fe, and 95Mo ENDOR to illuminate the partial electron-density redistribution upon E4(H2,2H) formation during catalysis, complementing these results with corresponding DFT computations. Inferences from the E4(2H)* ENDOR results as extended by DFT computations include: (i) the Mo-site participates negligibly, and overall it is unlikely that Mo changes valency throughout the catalytic cycle; (ii) two distinctive E4(4H) 57Fe signals are suggested as associated with ‘anchors’ of one bridging hydride, two others with anchors of the second, with NBO-analysis identifying one anchor of each hydride as a major recipient of electrons released upon breaking Fe-H bonds.

Lukoyanov, Dmitriy↗
DOE OSTI2022

Correlating Thermodynamic and Kinetic Hydricities of Rhenium Hydrides

The kinetics of hydride transfer from Re( R bpy)(CO) 3 H (bpy = 4,4'-R-2,2'-bipyridine; R = OMe, t Bu, Me, H, Br, COOMe, CF 3 ) to CO 2 and seven different cationic N-heterocycles were determined. Additionally, the thermodynamic hydricities of complexes of the type Re( R bpy)(CO) 3 H were established primarily using computational methods. Linear free-energy relationships (LFERs) derived by correlating thermodynamic and kinetic hydricities indicate that, in general, the rate of hydride transfer increases as the thermodynamic driving force for the reaction increases. Kinetic isotope effects range from inverse for hydride transfer reactions with a small driving force to normal for reactions with a large driving force. Hammett analysis indicates that hydride transfer reactions with greater thermodynamic driving force are less sensitive to changes in the electronic properties of the metal hydride, presumably because there is less buildup of charge in the increasingly early transition state. Bronsted α values were obtained for a range of hydride transfer reactions and along with DFT calculations suggest the reactions are concerted, which enables the use of Marcus theory to analyze hydride transfer reactions involving transition metal hydrides. It is notable, however, that even slight perturbations in the steric properties of the Re hydride or the hydride acceptor result in large deviations in the predicted rate of hydride transfer based on thermodynamic driving forces. This indicates that thermodynamic considerations alone cannot be used to predict the rate of hydride transfer, which has implications for catalyst design.

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

Aqueous electrocatalytic CO 2 reduction using metal complexes dispersed in polymer ion gels

We use fac-[Re(bpy)(CO) 3 Cl] ([Re–Cl]) dispersed in polymer ion gel (PIG) ([Re]–PIG) to carry out electrocatalytic CO 2 reduction in water. Electrolysis at -0.68 V vs. RHE in a CO 2 -saturated KOH and K 2 CO 3 solution produces CO with over 90% Faradaic efficiency. The PIG electrode is readily combined with water oxidation catalysts to generate a full electrochemical cell. Additionally, we provide evidence that the PIG electrode can be implemented with other molecular catalysts.

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

Controlled site coverage of strong metal–support interaction (SMSI) on Pd NP catalysts

Here, strong metal–support interaction catalysts have been shown to improve desired product selectivity at the cost of fractional rates due to active site coverage. The goal of this study was to determine if the active site coverage of metallic nanoparticles could be controlled to lower levels than have been previously reported in SMSI catalysts with the aim of improving the rate while maintaining high selectivity. 2Pd–XTi/SiO 2 (2 wt% Pd, X wt% Ti) strong metal–support interaction (SMSI) catalysts with Ti loadings between 0–1.0 wt% were synthesized to control Pd nanoparticle coverage. Calcination at 450 °C and reduction at 550 °C were sufficient for forming ~2 nm sized Pd particles in all catalysts. Increasing the Ti loading from 0.1 to 1.0 wt% increased the surface coverage from 40 to 85% at a fixed reduction temperature of 550 °C. The IR spectra of the SMSI catalysts were similar with a high fraction of linear bonded CO which was much higher than that of Pd nanoparticles of similar size. The SMSI overlayer could be removed by oxidation at 350 °C and re-reduction at 200 °C. EXAFS of the oxidized catalysts indicates that nearly full oxidation of the metallic nanoparticle was required to remove the SMSI overlayer. Oxidation temperatures from 30 to 300 °C partially oxidized the Pd nanoparticles and subsequent re-reduction at 200 °C partially decreases the SMSI coverage. The fractional surface coverage was determined by measuring the rate of propylene hydrogenation with and without the SMSI overlayer. Increasing the reduction temperature from 200 to 550 °C increased the SMSI coverage from 0 to 85% depending on the Ti loading and temperature. In conclusion, after reduction at 550 °C and oxidation at 350 °C, the range of coverages varied between ~10% with 0.1 wt% Ti after re-reduction at 300 °C and ~85% with 1 wt% Ti after reduction at 550 °C.

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

Advanced thermal/environmental barrier coatings of high-entropy rare earth disilicates tuned by strong anharmonicity of Eu 2 Si 2 O 7

Advancing thermal/environmental barrier coating (TEBC) materials with integrated thermal-mechanical functions is paramount for safeguarding SiC-based ceramic matrix composites (CMCs) in high-efficiency gas turbines. Herein, we employ a synergistic approach, combining density functional theory (DFT) methods and combinatorial chemistry techniques, to design high-performance and low-cost RE 2 Si 2 O 7 (RE = rare earth elements) TEBC materials tailored for enhanced compatibility with SiC-based CMCs. Expanding on phase stability of alloying pure RE 2 Si 2 O 7 , the investigation extends to the mechanical and thermal properties of solid solution systems, including Er 1/2 Y 3/4 Yb 3/4 Si 2 O 7 , Gd 1/4 Er 1/4 Y 3/4 Yb 3/4 Si 2 O 7 , and Eu 1/4 Er 1/4 Y 3/4 Yb 3/4 Si 2 O 7 . The solid solution systems exhibit a major reduction in lattice thermal conductivity relative to their pure counterparts, achieving ultralow values of 0.25 to 0.39 W m −1 K −1 at 1500 K. Furthermore, the coefficients of thermal expansion (CTE) of these solid solutions are precisely tuned within the desired range for SiC (4.4 to 5.5 × 10 −6 K −1 ), while maintaining good mechanical properties. Here, in particular, the addition of Eu 2 Si 2 O 7 demonstrates to be an important variable to the tuning of CTE and lattice thermal conductivity by leveraging its strong anharmonicity, presenting a pioneering avenue for fine-tuning material properties. In summary, this research not only identifies promising TEBC materials with superior thermal properties, but also introduces a valuable computational material design methodology for the rapid discovery of complex materials for harsh environments.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Beyond real: alternative unitary cluster Jastrow models for molecular electronic structure calculations on near-term quantum computers

Near-term quantum devices require wavefunction ansätze that are expressive while also of shallow circuit depth in order to both accurately and efficiently simulate molecular electronic structure. While the unitary coupled cluster ansatz (e.g., UCCSD) has become a standard, the high gate count associated with the implementation of this limits its feasibility on noisy intermediate-scale quantum (NISQ) hardware. k -Fold unitary cluster Jastrow (uCJ) ansätze mitigate this challenge by providing O( kN 2 ) circuit scaling and favorable linear depth circuit implementation. Previous work has focused on the real orbitalrotation (Re-uCJ) variant of uCJ, which allows an exact (Trotter-free) implementation. Here we extend and generalize the k -fold uCJ framework by introducing two new variants, Im-uCJ and g-uCJ, which incorporate imaginary and fully complex orbital rotation operators, respectively. Similar to Re-uCJ, both of the new variants achieve quadratic gate-count scaling. Our results focus on the simplest k = 1 model, and show that the uCJ models frequently maintain energy errors within chemical accuracy (∼1 kcal mol −1 ). Both g-uCJ and Im-uCJ are more expressive in terms of capturing electron correlation and are also more accurate than the earlier Re-uCJ ansatz. We further show that Im-uCJ and g-uCJ circuits can also be implemented exactly, without any Trotter decomposition. Numerical tests using k = 1 on H 2 , H 3 + , Be 2 , C 2 H 4 , C 2 H 6 and C 6 H 6 in various basis sets confirm the practical feasibility of these shallow Jastrow-based ansätze for applications on near-term quantum hardware.

Tkachenko, Nikolay V. [University of California, B↗
DOE OSTI2020

Enhancement in hard magnetic properties of (Nd, Pr)–Fe–B melt-spun ribbons

The coercivity of RE 2 Fe 14 B-type permanent magnets is strongly influenced by the microstructural features such as grain boundary (GB) phases as well as grain sizes. Here, we have combined micromagnetic simulations and experiments to elucidate the role of excess RE (Nd/Pr) in determining the resulting hard magnetic properties of Nd–Pr–Fe–B melt-spun ribbons. The intrinsic coercivity (H c ) at room temperature significantly enhanced from 9.7 kOe to 15.3 kOe with the increase in the Nd/Pr content. Furthermore, the effect of non-magnetic grain refining refractory carbide (TiC) on both the microstructure and magnetic hardening was studied. The addition of TiC showed a very high coercivity H c of up to 19.0 kOe at room temperature. Micromagnetic simulation indicates that the coercivity enhancement is mainly due to the reduction of inter-grain magnetic interaction, which is due to the RE-rich nonmagnetic grain boundary (GB) phase and/or TiC distributed at the GB. This work provides useful information on the roles of non-magnetic grain boundary phases for improving the coercivity of Nd–Pr–Fe–B magnets. Combined with experimental and modeling results, we have discussed the mechanism responsible for the enhancements in coercivity and the suitability of the alloys for high-performance permanent magnet development.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Flow and heat transfer in the wake of a triangular arrangement of spheres

This research work seeks to investigate the influence of spacing and heat transfer on the wake behavior of a triangular arrangement of spheres. Four experimental configurations have been investigated at three Reynolds numbers, Re 1 = 350, Re 2 = 700, and Re 3 = 1050. Two isothermal cases were investigated with spacing between the spheres of zero and one sphere diameter, and two cases were investigated with an applied heat flux at the same spacing conditions. The time resolved particle image velocimetry results revealed various flow phenomena including flow separations, von Karman vortex shedding, and Kelvin–Helmholtz instabilities. The turbulent statistics reveal the effect of proximity and heat transfer on the time averaged values of the wake size, turbulent strengths, and Reynolds shear stress in the wake of each sphere, namely, the laminarization effects from the addition of heat and the suppression of the lead sphere wake from the proximity of the trailing spheres. These results are complemented by the application of proper orthogonal decomposition (POD) to the flow fields, which extracts the coherent structures from the flow. Here the modes that describe the coherent structures are extracted and described in detail, which provide further insight into effects of the experimental conditions on the temporal behavior of the flow. Many of the low order modes are found to be associated in pairs, corresponding to asymmetric structures or advection of a given structure downstream. The capability of POD to produce reduced order models of the flow is then utilized to facilitate vortex identification analysis. A turbulent kinetic energy based mode truncation criteria, which has been found to enhance vortex identification capability, is applied to select the POD modes and temporal coefficients to be used in the reduced order modeling. The reconstructed velocity fields are then analyzed with vortex identification algorithms to extract the vortex cores and boundaries. The combination of these approaches allows the study of the effect of proximity and heat transfer on the vortex characteristics, such as size, strength, and distribution.

42 ENGINEERING↗
DOE OSTI2022

Particle deceleration for collective QED signatures

Frequency upshifts have been proposed as a first experimental signature of collective effects in quantum electrodynamic cascade generated electron–positron pair plasmas. Since the high effective masses of generated pairs will reduce any frequency change, stopped pairs at a minimal Lorentz factor in the lab frame were thought to be the dominant contribution to the laser upshift. However, we demonstrate that only considering stopped particles unduly neglects the contributions of particles re-accelerated in the laser propagation direction. Re-accelerated particles should, on a per particle basis, affect the laser more strongly and over a much longer timescale. To maximize particle contributions to the laser upshift, we consider a Laguerre–Gaussian (LG) mode laser beam to better reflect generated pairs. Furthermore, the LG mode does not have an advantage in particle deceleration and re-acceleration when compared against a Gaussian beam, but the LG mode can maintain particle contributions for a longer duration, allowing for more pair density accumulation. Deceleration with a structured beam to keep pairs within the laser should create a larger upshift, thereby lowering the demands on the driving laser.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

The impact of collisionality on the runaway electron avalanche during a tokamak disruption

The exponential growth (avalanching) of runaway electrons (REs) during a tokamak disruption continues to be a large uncertainty in RE modeling. The present work investigates the impact of tokamak geometry on the efficiency of the avalanche mechanism across a broad range of disruption scenarios. It is found that the parameter ν *,crit ⁠, describing the collisionality at the critical energy to run away, delineates how toroidal geometry impacts RE formation. In particular, utilizing a reduced but self-consistent description of plasma power balance, it is shown that for a high-density deuterium-dominated plasma, ν *,crit is robustly less than one, resulting in a substantial decrease in the efficiency of the RE avalanche compared to predictions from slab geometry. In contrast, for plasmas containing a substantial quantity of neon or argon, ν *,crit ≳1, no reduction of the avalanche is observed due to toroidal geometry. This sharp contrast in the impact of low- vs high-Z material results primarily from the relatively strong radiative cooling from high-Z impurities, enabling the plasma to be radiatively pinned at low temperatures and thus large electric fields, even for modest quantities of high-Z material.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Resistive hose modes in tokamak runaway electron beams II

Resistive hose instabilities of runaway electron (RE) beams immersed in resistive background plasmas are examined with analytic and numerical calculations. The RE beam-plasma equilibria considered are characteristic of the situation observed post-thermal quench in a tokamak disruption. An analytic linear dispersion relation is presented for the case of a uniform RE current density profile with a sharp boundary in cylindrical geometry. Initial value linear calculations for a more general profile in toroidal geometry find that reducing aspect ratio increases the resistive hose mode growth rates with fixed safety factor profile. Nonlinear calculations in cylindrical and toroidal geometry find that the resistive hose instability-driven fluctuations relax the gradient of the current density profile. In toroidal geometry, changes to the magnetic topology are observed as a result of the resistive hose activity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Phase stability and coercivity in La 2 Fe 14 B magnet

Critical rare-earth free La 2 Fe 14 B (2:14:1) has the potential to be a gap permanent magnet. However, La 2 Fe 14 B decomposes into La, α-Fe, and La 4 Fe 4 phases below 1067 K. The phase stability and coercivity have been studied in La 2 Fe 14 B magnet using first principles DFT (density functional theory) calculation and micromagnetic simulation. For a perfect La 2 Fe 14 B cube (edge length of 256 nm) without any structural defects and soft magnetic secondary phases, the coercivity (8.5 kOe) is reduced to ~40% of its magnetocrystalline anisotropy field ( H A = 20 kOe). Further, the coercivity sharply reduces to 3.2 kOe upon forming a thin layer (2 nm) of α-Fe on the surface of the La 2 Fe 14 B cube particle. The DFT calculations indicate that a partial replacement of La by other rare-earth (RE) elements can enhance the structural stability of 2:14:1. The gains in cohesive energy are 0.75, 0.10, and 0.33 eV per formula unit in (La 0.5 RE 0.5 ) 2 Fe 14 B with RE = Ce, Pr, and Nd, respectively. Stabilizing the 2:14:1 structure and mitigating the formation of soft magnetic structural defects or impurity phases such as α-Fe is necessary to develop La 2 Fe 14 B based magnet, which can be moderately achieved via partial substitution of La by other rare earth elements such as Ce, Pr, and Nd.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Harnessing interpretable and unsupervised machine learning to address big data from modern X-ray diffraction

The information content of crystalline materials becomes astronomical when collective electronic behavior and their fluctuations are taken into account. In the past decade, improvements in source brightness and detector technology at modern X-ray facilities have allowed a dramatically increased fraction of this information to be captured. Now, the primary challenge is to understand and discover scientific principles from big datasets when a comprehensive analysis is beyond human reach. We report the development of an unsupervised machine learning approach, X-ray diffraction (XRD) temperature clustering (X-TEC), that can automatically extract charge density wave order parameters and detect intraunit cell ordering and its fluctuations from a series of high-volume X-ray diffraction measurements taken at multiple temperatures. We benchmark X-TEC with diffraction data on a quasi-skutterudite family of materials, (Ca x Sr 1–x ) 3 Rh 4 Sn 13 , where a quantum critical point is observed as a function of Ca concentration. We apply X-TEC to XRD data on the pyrochlore metal, Cd 2 Re 2 O 7 , to investigate its two much-debated structural phase transitions and uncover the Goldstone mode accompanying them. We demonstrate how unprecedented atomic-scale knowledge can be gained when human researchers connect the X-TEC results to physical principles. Specifically, we extract from the X-TEC–revealed selection rules that the Cd and Re displacements are approximately equal in amplitude but out of phase. This discovery reveals a previously unknown involvement of 5d 2 Re, supporting the idea of an electronic origin to the structural order. Our approach can radically transform XRD experiments by allowing in operando data analysis and enabling researchers to refine experiments by discovering interesting regions of phase space on the fly.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Fluid inertia controls mineral precipitation and clogging in pore to network-scale flows

Mineral precipitation caused by fluid mixing presents complex control and predictability challenges in a variety of natural and engineering processes, including carbon mineralization, geothermal energy, and microfluidics. Precipitation dynamics, particularly under the influence of fluid flow, remain poorly understood. Combining microfluidic experiments and three-dimensional reactive transport simulations, we demonstrate that fluid inertia controls mineral precipitation and clogging at flow intersections, even in laminar flows. We observe distinct precipitation regimes as a function of Reynolds number (Re). At low Reynolds numbers (Re < 10), precipitates form a thin, dense layer along the mixing interface, which shuts precipitation off, while at high Reynolds numbers (Re > 50), strong three-dimensional flows significantly enhance precipitation over the entire intersection, resulting in rapid clogging. When injection rates from two inlets are uneven, flow symmetry-breaking leads to unexpected flow bifurcation phenomena, which result in enhanced concurrent precipitation in both downstream channels. Finally, we extend our findings to rough channel networks and demonstrate that the identified inertial effects on precipitation at the intersection scale are also present and even more dramatic at the network scale. This study sheds light on the fundamental mechanisms underlying mixing-induced mineral precipitation and provides a framework for designing and optimizing processes involving mineral precipitation.

Science & Technology - Other Topics↗
DOE OSTI2023

Mesh and model requirements for capturing deep-stall aerodynamics in low-Mach-number flows

Here, this work presents a comprehensive computational fluid dynamics investigation of the effects of grid resolution and turbulence-model choice for capturing the unsteady three-dimensional aerodynamic performance of NACA 0012 and 0021 airfoils, with specific focus on the deep-stall regime. At high angles of attack (α), wind turbine blades routinely experience vortex-induced vibrations, which can cause significant structural damages. Accurate predictions of post-stall aerodynamics can identify the frequencies at which such vibrations maybe triggered. In this context, the NACA 0012 airfoil simulations are conducted at a chord-based Reynolds number, Re c =2×10 6 , with the k-ω Shear-Stress Transport Reynolds-Averaged Navier-Stokes (RANS) and Improved Delayed Detached Eddy Simulation (IDDES) hybrid RANS-Large Eddy Simulation turbulence models. The effect of mesh resolution both in the wall-normal and spanwise directions is investigated. Only the IDDES model with a minimum spanwise resolution of 24 cells per chord length correctly predicts the aerodynamic forces. Spectral analysis shows the peak primary shedding frequency at α=30°, which signifies the end of the stall region. In the post-stall regime, both lift and drag frequencies drop asymptotically with increasing α. The Strouhal number, based on normalised chord length, remains nearly constant in this region. Based on this study, NACA 0021 airfoil runs are performed with IDDES for Re c =2.7×10 5 and 2.0×10 6 on the finest wall-normal mesh and three spanwise grids. Simulations conducted on the finer spanwise grids demonstrate grid independence and show good agreement with experiments. The effect of varying Rec on the airfoil frequency statistics is investigated. Additionally, comparison studies are presented to investigate the impact of airfoil thickness on the frequency content at Re c =2.0×10 6 . The results from the study provide guidance on the choice of mesh resolution with the IDDES model to accurately capture aerodynamic quantities for complex industrial applications.

17 WIND ENERGY↗