Search NASA⌕ Search

SEARCH · Search NASA

Results for “Complexation”

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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Redox-active dinuclear oxorhenium(V) pyrazolate complexes

Four new structurally similar dinuclear oxorhenium(V) complexes, [{Re(O)X(PPh 3 )} 2 (μ-O)(μ-4-x'-pz) 2 ], where pz = pyrazolate anion, X = X' = Cl (1) and Br (4), X = Cl, X' = Br (2), and X = Br, X' = Cl (3), have been synthesized and characterized. Little variation in spectroscopic features – 1 H NMR, IR, UV–Vis – exists among the four complexes. All complexes possess a bent Re-O-Re core as well as distorted octahedral coordination geometry around the rhenium centers. Finally, a reversible one-electron electrochemical process is observed at approximately 0.84 V vs. Fc + /Fc in all four complexes; however, changing the terminal halide from chloride to bromide slightly destabilizes the oxidized Re(VI) center.

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

Molecular structures of dissolved and colloidal As V –Fe III complexes and their roles in the mobilization of As V under strongly acidic conditions

The effect of high concentration of iron (Fe III ) on the speciation and mobility of arsenic (As) under strongly acidic conditions remains unclear. This work studied the redistribution and speciation of As V and Fe III at Fe/As molar ratio of 1–14 and pH 1.5–2.0 in the dissolved, colloidal, and solid phases. Results showed that the elevated Fe III induced the decomposition of the precipitated poorly crystalline ferric arsenate by forming dissolved (< 3 kDa) and colloidal (3 kDa–0.1 µm) As–Fe complexes. The fraction of particulate As (> 0.1 µm) decreased from 70–90% to less than 20% when the Fe/As molar ratio increased from 1 to 14. The particle size of the bulk samples decreased significantly with the increase of Fe III concentration. The FTIR results suggested that AsV in dissolved/colloidal As–Fe complexes dominantly occurred as HAsO 4 2– species. The EXAFS results indicated that each HAsO 4 2– coordinated with approximately two Fe atoms in dissolved/colloidal As–Fe complexes at Fe/As ≥ 2. Furthermore, the findings suggest that high aqueous Fe III concentration can promote the mobility of As by forming dissolved/colloidal Fe–As complexes in acidic waters, potentially accelerating As transport from source to downstream in acid mine drainage systems.

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

Framework for predicting the fractionation of complex liquid feeds via polymer membranes

The separation of complex liquid hydrocarbon mixtures was recently demonstrated using the glassy polymer SBAD-1, showing that small molecule fractionation is possible by such organic membrane materials. Here, in this work, we develop a framework that will enable workable predictions of permeate flux and composition in complex hydrocarbon liquids through intrinsically porous glassy polymers. The predictions are made by incorporating experimentally-derived unary sorption and diffusion parameters in a Maxwell-Stefan framework coupled with multicomponent sorption models and various distinct diffusion phenomena. Across the range of sorption and diffusion phenomena considered, both the conventional Flory-Huggins model and the proposed Langmuir + Flory-Huggins sorption model combined with a simple average guest diffusivity or a more complex free-volume theory-based transport resulted in the lowest prediction error for three chosen multicomponent separations. The proposed Maxwell-Stefan framework simply requires pure component transport parameters to allow a fast approximation of the separation of multicomponent liquid hydrocarbon feeds that can potentially be extended to more complex feeds such as crude oil fractions.

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

A chiral alkali metal capped Ni 4 cubane complex: Synthesis, structure, magnetic and catalytic bromination studies

A new chiral NaO 2 Ni 4 O 2 N 2 cluster has been developed with a distorted capped cubane-type structure. The structure was determined by single crystal X-ray diffraction. The tetradentate Schiff-base ligand was synthesized by the condensation of o-vanillin and $\textit{D}$-(-)–threo-2-amino-1-(4-nitrophenyl)-1,3-propanediol. X-ray structure reveals that the complex possesses a Ni 4 cubane cluster where sodium ion is anchored on the face of the cubane structure. Continuous shape measurement studies revealed that all the Ni centres are in distorted octahedral geometry and sodium ion adopts the spherical square pyramidal geometry. Further structural exploration reveals that the complex shows local distortions of the coordination geometries at each metal site, and these distortions govern the magnetic anisotropy associated with each Ni II ion. For the present case, an elongation ($d_{str}$ > 0) of the octahedral coordination sphere results in a local easy-plane anisotropy, with a positive $D_{Ni(II)}$ parameter. Temperature dependant magnetization measurements showed a predominant ferromagnetic interaction amongst the four Ni centres. This complex is also catalytically active towards bromination of aromatic phenolic compounds. Initial catalytic activity was verified with the conversion of phenol red to bromophenol blue. Different phenolic substances show good conversion to their corresponding bromo derivatives. The complex showed good catalytic activity towards bromination of aromatic phenolic compounds, offering up to 99% conversion with 100% selectively to their corresponding bromo derivatives.

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

Structure and Mechanism of Acetylation by the N-Terminal Dual Enzyme NatA/Naa50 Complex

NatA co-translationally acetylates the N termini of over 40% of eukaryotic proteins and can associate with another catalytic subunit, Naa50, to form a ternary NatA/Naa50 dual enzyme complex (also called NatE). The molecular basis of association between Naa50 and NatA and the mechanism for how their association affects their catalytic activities in yeast and human are poorly understood. Here, we determined the X-ray crystal structure of yeast NatA/Naa50 as a scaffold to understand coregulation of NatA/Naa50 activity in both yeast and human. We find that Naa50 makes evolutionarily conserved contacts to both the Naa10 and Naa15 subunits of NatA. These interactions promote catalytic crosstalk within the human complex, but do so to a lesser extent in the yeast complex, where Naa50 activity is compromised. These studies have implications for understanding the role of the NatA/Naa50 complex in modulating the majority of the N-terminal acetylome in diverse species.

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

Community Data Mining Approach for Surface Complexation Database Development

This paper presents a comprehensive data-to-model workflow, including a findable, accessible, interoperable, reusable (FAIR) community sorption database (newly developed LLNL Surface Complexation/Ion Exchange (L-SCIE) database) along with a data fitting workflow to efficiently optimize surface complexation reaction constants with multiple surface complexation model (SCM) constructs. This workflow serves as a universal framework to mine, compile, and analyze large numbers of published sorption data as well as to estimate reaction constants for parameterizing reactive transport models. Here the framework includes (1) data digitization from published papers, (2) data unification including unit conversions, and (3) data-model integration and reaction constant estimation using geochemical software PHREEQC coupled with the universal parameter estimation code PEST. We demonstrate our approach using an analysis of U(VI) sorption to quartz based on a first L-SCIE implementation, concluding that a multisite SCM construct with carbonate surface species yielded the best fit to community data. Surface complexation reaction constants extracted from this approach captured all available sorption data available in the literature and provided insight into previously published reaction constants and surface complexation model constructs. The L-SCIE sorption database presented herein allows for automating this approach across a wide range of metals and minerals and implementing novel machine learning approaches to reactive transport in the future.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2020

Iron(II) Complexes of an Anionic Bis(ylide)diphenylborate Ligand

Double deprotonation of the salt [Ph 2 B(PMe 3 ) 2 ][OTf] provides access to a bis(ylide)diphenylborate ligand that is readily transferred in situ to iron(II). Depending on the reaction stoichiometry, both the “ate” complex [Ph 2 B(Me 2 PCH 2 ) 2 Fe(μ-Cl) 2 Li(THF) 2 ] and the homoleptic complex [Ph 2 B(Me 2 PCH 2 ) 2 ] 2 Fe can be prepared from FeCl 2 (THF) 1.5 . Further reaction of with FeCl 2 (THF) 1.5 produces the chloride-bridged dimer [Ph 2 B(Me 2 PCH 2 ) 2 Fe(μ-Cl) 2 Fe(CH 2 PMe 2 ) 2 BPh 2 ]. Attempts to reduce or alkylate provide as the only isolable product, likely a consequence of the low steric hindrance of the bis(ylide)diphenylborate ligand. On the other hand, reaction of with the strong field ligand CN t Bu provides the six-coordinate, diamagnetic complex [Ph 2 B(Me 2 PCH 2 ) 2 Fe(CN t Bu) 4 ][Cl]. Here, the electronic structure calculations for the bis(ylide)diphenylborate ligand and homoleptic complex suggest that the C(ylide) atoms are strong σ-donors with little π-bonding character. These initial results suggest the potential for this bis(ylide)diphenylborate ligand in coordination chemistry.

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

Restructuring of Benzimidazole-Based Copper Complexes during Electrochemical CO 2 Reduction

Copper-based molecular catalysts have been widely used for the electrochemical reduction of carbon dioxide (ERCO 2 ). However, establishing a correlation between structure and catalytic performance has been challenging due to the dynamic nature of these materials under ERCO 2 working conditions. Here, we report on the restructuring of the copper complexes [Cu(bzimpy)Cl 2 ] and [Cu(pyrben) 2 (NO 3 )]NO 3 during ERCO 2 in an acidic medium, evidencing its effects on the observed activity, selectivity, and stability. In situ X-ray absorption spectroscopy (XAS) suggested that, at sufficiently negative potentials, the Cu 2+ centers of both complexes are reduced to undercoordinated Cu 0 species. Pb underpotential deposition of post-ERCO 2 samples indicated that such Cu 0 species also present a dominant contribution of the (100) domain, which is in line with the high Faradaic efficiency toward C 2 products (∼45%) shown by both complexes at −1.24 V vs RHE. By combining online electrochemical mass spectrometry, in situ XAS, electron paramagnetic resonance (EPR), and quantitative product analysis, we found that [Cu(bzimpy)Cl 2 ] exhibits a partially reversible restructuring, regenerating the main complex moiety, [Cu(bzimpy)] 2+ , and its geometry.

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

Olefin Coupling Catalyzed by (Pybox)Os Complexes via Osmacyclopentane Intermediates: Comparison with Isoelectronic (Phebox)Ir

(Pybox)Os is found to catalyze alkene hydrovinylation, effecting the dimerization of ethylene, tail-to-tail coupling of propene and 1-butene, and cross-coupling of ethylene with higher α-olefins. This reactivity contrasts with the previously reported dehydrogenative coupling of ethylene to give butadiene catalyzed by the isoelectronic fragment (Phebox)Ir. The reaction mechanism was investigated through computational and experimental means. Both the Os- and Ir-catalyzed reactions proceed through a [2 + 2 + 1] cyclization of the corresponding bis-olefin complex to yield an experimentally observed metallacyclopentane intermediate. In both cases, the metallacyclopentane undergoes β- H elimination, via a dechelated κ 2 -pincer-ligated intermediate, to yield a σ−π-but-3-enyl hydride complex or derivative. Both the greater reactivity and the distinct chemoselectivity of the Os system relative to the Ir system are attributable to C−H reductive elimination by the σ−π-but-3-enyl hydride having a barrier for Os much lower than that for Ir. This lower barrier to C−H elimination for Os is unexpected given that the thermodynamic driving force for elimination is much less for Os than for Ir. Computational studies of model complexes were conducted, comparing (Pybox)Os(L)(CH 3 )(H) with the isoelectronic (Phebox)Ir(L)(CH 3 )(H). The results indicate that the more facile kinetics with Os relative to Ir may be general for C−H elimination from six-coordinate d 6 complexes of the two metals, as well as for the microscopic reverse, i.e., C−H addition to the corresponding four-coordinate d 8 species.

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

Structural chemistry of penta- and hexanitrato thorium($\tiny{IV}$) complexes isolated using N–H donors

In this study, a series of fifteen tetravalent thorium phases were prepared. The compounds were isolated from acidic aqueous nitrate solutions using protonated nitrogen heterocycles of varying hydrogen-bond donation strength. Structural analysis via single crystal X-ray diffraction showed that the structures are built from pentanitrato, [Th(NO 3 ) 5 (H 2 O) 2 ] 1- , and hexanitrato, [Th(NO 3 ) 6 ] 2- , molecular units, with the latter being far more prevalent in the solid state. The vibrational properties of the compounds were examined using Raman and IR spectroscopy; the spectra are dominated by stretches characteristic of nitrate and the organic ions. The relative energetics of nitrate complexation was examined using electronic structure theory. These results confirmed that there are clear thermodynamic sinks for the penta- and hexanitrato structural units that were observed experimentally. Additionally, electrostatic surface potentials (ESPs) were calculated in an effort to better understand the counterion stabilization of the complexes. The ESP surfaces showed that the position of the water and nitrate molecules and the coordination geometry of the metal complex had a clear effect on the polarizability of the two structural motifs. Despite limited speciation of the Th–nitrate structural units, the compounds exhibit rich supramolecular chemistry resulting from hydrogen bonding of the Th complexes with the organic N–H donors and π–π stacking interactions from the protonated N-heterocycles.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Highly-destabilized ligand field excited states of iron carbene complexes and their relation to charge transfer state lifetimes

Lifetimes of photoexcited charge transfer (CT) states in transition metal chromophores are influenced by low-lying ligand field (LF) excited states, especially for 3d metal complexes. To manipulate interactions between LF and CT states, it is important to be able to control LF excited state energies using tunable synthetic variables. In this report, we use Fe 2p3d L 3 -edge resonant inelastic X-ray scattering (RIXS) to measure LF excited state energies of three homoleptic iron chromophores coordinated by strong-field N-heterocyclic carbenes (NHCs). We investigate the effect of oxidation state and ligand scaffold on LF energies and covalency parameters. A cyclometalated bis(NHC) ligand affords both high LF excited state energies (and thus high 10 Dq) as well as high metal–ligand covalency compared to other iron complexes with very strong-field ligands. However, for the set of complexes investigated, we do not observe meaningful correlation between the LF excited state energies and the CT excited state lifetimes. These results illustrate that targeting long-lived CT excited states necessitates control of multiple molecular excited state properties, with destabilization of the LF excited state energies proving necessary, but insufficient, to control the CT excited state lifetime in Fe carbene complexes.

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

Structural characterization of an activin class ternary receptor complex reveals a third paradigm for receptor specificity

TGFβ family ligands, which include the TGFβs, BMPs, and activins, signal by forming a ternary complex with type I and type II receptors. For TGFβs and BMPs, structures of ternary complexes have revealed differences in receptor assembly. However, structural information for how activins assemble a ternary receptor complex is lacking. In this paper, we report the structure of an activin class member, GDF11, in complex with the type II receptor ActRIIB and the type I receptor Alk5. The structure reveals that receptor positioning is similar to the BMP class, with no interreceptor contacts; however, the type I receptor interactions are shifted toward the ligand fingertips and away from the dimer interface. Mutational analysis shows that ligand type I specificity is derived from differences in the fingertips of the ligands that interact with an extended loop specific to Alk4 and Alk5. The study also reveals differences for how TGFβ and GDF11 bind to the same type I receptor, Alk5. For GDF11, additional contacts at the fingertip region substitute for the interreceptor interactions that are seen for TGFβ, indicating that Alk5 binding to GDF11 is more dependent on direct contacts. In support, we show that a single residue of Alk5 (Phe 84 ), when mutated, abolishes GDF11 signaling, but has little impact on TGFβ signaling. The structure of GDF11/ActRIIB/Alk5 shows that, across the TGFβ family, different mechanisms regulate type I receptor binding and specificity, providing a molecular explanation for how the activin class accommodates low-affinity type I interactions without the requirement of cooperative receptor interactions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Circuit complexity near critical points

Here, we consider the Bose–Hubbard model in two and three spatial dimensions and numerically compute the quantum circuit complexity of the ground state in the Mott insulator and superfluid phases using a mean field approximation with additional quadratic fluctuations. After mapping to a qubit system, the result is given by the complexity associated with a Bogoliubov transformation applied to the reference state taken to be the mean field ground state. In particular, the complexity has peaks at the O(2) critical points where the system can be described by a relativistic quantum field theory. Given that we use a Gaussian approximation, near criticality the numerical results agree with a free field theory calculation. To go beyond the Gaussian approximation we use general scaling arguments that imply that, as we approach the critical point t → t c , there is a non-analytic behavior in the complexity c 2 (t) of the form |c 2 (t) – c 2 (t c )| ~ |t – t c | νd , up to possible logarithmic corrections. Here d is the number of spatial dimensions and ν is the usual critical exponent for the correlation length ξ ~ |t – t c | –ν . As a check, for d = 2 this agrees with the numerical computation if we use the Gaussian critical exponent $v=\frac{1}{2}$. Finally, using AdS/CFT methods, we study higher dimensional examples and confirm this scaling argument with non-Gaussian exponent ν for strongly interacting theories that have a gravity dual.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Delocalized excitons in natural light-harvesting complexes

Natural organisms such as photosynthetic bacteria, algae, and plants employ complex molecular machinery to convert solar energy into biochemical fuel. An important common feature shared by most of these photosynthetic organisms is that they capture photons in the form of excitons typically delocalized over a few to tens of pigment molecules embedded in protein environments of light harvesting complexes (LHCs). Delocalized excitons created in such LHCs remain well protected despite being swayed by environmental fluctuations and are delivered successfully to their destinations over 100 nanometer distances in about 100 ps times. Decades of experimental and theoretical investigation have produced a large body of information offering insight into major structural, energetic, and dynamical features contributing to LHCs’ extraordinary capability to harness photons using delocalized excitons. The objective of this review is (i) to provide a comprehensive account of major theoretical, computational, and spectroscopic advances that have contributed to this body of knowledge, and (ii) to clarify the issues concerning the role of delocalized excitons in achieving efficient energy transport mechanisms. The focus of this review is on three representative systems: the Fenna-Matthews-Olson complex of green sulfur bacteria, the light harvesting 2 complex of purple bacteria, and phycobiliproteins of cryptophyte algae. Although we offer a more in-depth and detailed description of theoretical and computational aspects, major experimental results and their implications are also assessed in the context of achieving excellent light-harvesting functionality. Future theoretical and experimental challenges to be addressed in gaining a better understanding and utilization of delocalized excitons are also discussed.

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

Mesoscale to Microscale Simulations over Complex Terrain with the Immersed Boundary Method in the Weather Research and Forecasting Model

Improvements to the Weather Research and Forecasting (WRF) Model are made to enable multiscale simulations over highly complex terrain with dynamically downscaled boundary conditions from the mesoscale to the microscale. Over steep terrain, the WRF Model develops numerical errors that are due to grid deformation of the terrain-following coordinates. An alternative coordinate system, the immersed boundary method (IBM), has been implemented into WRF, allowing for simulations over highly complex terrain; however, the new coordinate system precluded nesting within mesoscale simulations using WRF’s native terrain-following coordinates. Here, the immersed boundary method and WRF’s grid-nesting framework are modified to seamlessly work together. This improved framework for the first time allows for large-eddy simulation over complex (urban) terrain with IBM to be nested within a typical mesoscale WRF simulation. Simulations of the Joint Urban 2003 field campaign in Oklahoma City, Oklahoma, are performed using a multiscale five-domain nested configuration, spanning horizontal grid resolutions from 6 km to 2 m. These are compared with microscale-only simulations with idealized lateral boundary conditions and with observations of wind speed/direction and SF 6 concentrations from a controlled release from intensive observation period 3. The multiscale simulation, which is configured independent of local observations, shows similar model skill predicting wind speed/direction and improved skill predicting SF 6 concentrations when compared with the idealized simulations, which require use of observations to set mean flow conditions. Use of this improved multiscale framework shows promise for enabling large-eddy simulation over highly complex terrain with dynamically downscaled boundary conditions from mesoscale models.

42 ENGINEERING↗
DOE OSTI2018

Thermodynamic parameters for the complexation of technetium(IV) with EDTA

Technetium-99 is a high yield (~6% fission yield) fission product and long-lived (2.13×10 5 year half-life) component of nuclear waste that will be disposed of in a geological repository. Some 99 Tc has been released into the environment due to nuclear fuel and weapon production activities at sites such as Hanford, WA. Strongly complexing ligands such as ethylenediamine- N , N , N ', N '-tetraacetic acid (EDTA) are known to increase Tc(IV) solubility and mobility in environmental systems and an accurate quantification of the complexation of Tc(IV) with EDTA is important for predicting its behavior in a geological repository. A liquid–liquid extraction system utilizing 0.2 M TOPO in dodecane was used to measure the stability constants of Tc(IV)-EDTA in 0.50 m NaNO 3 at variable temperatures (14.0±0.1, 25.0±0.1, and 32.0±0.1°C). The acid dependence of the apparent stability constants in the p C H range of 2.00–2.70 indicated the formation of TcO(EDTA) 2- (log β 101 =17.9±0.3, 25.0±0.1°C) and a protonated complex TcO(H)(EDTA) - (log β 111 =20.5±0.1, 25.0±0.1°C). The associated thermodynamic parameters Δ r G 101 =-101.7±0.4 kJ·mol -1 , Δ r H 101 =-47±9 kJ·mol -1 , Δ r S 101 =179±36 J·mol -1 ·K -1 , Δ r G 111 =-117.2±0.3 kJ·mol -1 , Δ r H 111 =-23±5 kJ·mol -1 , and Δ r S 111 =315±63 J·mol -1 ·K -1 (0.50 m NaNO 3 , 25.0±0.1°C) were determined by van’t Hoff analysis. The formation of each Tc(IV)–EDTA complex is exothermic and present favorable entropy terms.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2019

Molecular Studies of the Protein Complexes Involving Cis -Prenyltransferase in Guayule ( Parthenium argentatum ), an Alternative Rubber-Producing Plant

Guayule ( Parthenium argentatum ) is a perennial shrub in the Asteraceae family and synthesizes a high quality, hypoallergenic cis -1,4-polyisoprene (or natural rubber; NR). Despite its potential to be an alternative NR supplier, the enzymes for cis -polyisoprene biosynthesis have not been comprehensively studied in guayule. Recently, implications of the protein complex involving cis -prenyltransferases (CPTs) and CPT-Binding Proteins (CBPs) in NR biosynthesis were shown in lettuce and dandelion, but such protein complexes have yet to be examined in guayule. In this work, we identified four guayule genes - three PaCPTs ( PaCPT1-3 ) and one PaCBP , whose protein products organize PaCPT/PaCBP complexes. Co-expression of both PaCBP and each of the PaCPTs could complemented the dolichol (a short cis -polyisoprene)-deficient yeast, whereas the individual expressions could not. Microsomes from the PaCPT/PaCBP -expressing yeast efficiently incorporated 14 C-isopentenyl diphosphate into dehydrodolichyl diphosphates; however, NR with high molecular weight could not be synthesized in in vitro assays. Furthermore, co-immunoprecipitation and split-ubiquitin yeast 2-hybrid assays using PaCPTs and PaCBP confirmed the formation of protein complexes. Of the three PaCPTs , guayule transcriptomics analysis indicated that the PaCPT3 is predominantly expressed in stem and induced by cold-stress, suggesting its involvement in NR biosynthesis. The comprehensive analyses of these PaCPTs and PaCBP here provide the foundational knowledge to generate a high NR-yielding guayule.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Efficacy of the Cell Perturbation Method in Large-Eddy Simulations of Boundary Layer Flow over Complex Terrain

A challenge to simulating turbulent flow in multiscale atmospheric applications is the efficient generation of resolved turbulence motions over an area of interest. One approach is to apply small perturbations to flow variables near the inflow planes of turbulence-resolving simulation domains nested within larger mesoscale domains. While this approach has been examined in numerous idealized and simple terrain cases, its efficacy in complex terrain environments has not yet been fully explored. Here, we examine the benefits of the stochastic cell perturbation method (CPM) over real complex terrain using data from the 2017 Perdigão field campaign, conducted in an approximately 2-km wide valley situated between two nearly parallel ridges. Following a typical configuration for multiscale simulation using nested domains within the Weather Research and Forecasting (WRF) model to downscale from the mesoscale to a large-eddy simulation (LES), we apply the CPM on a domain with horizontal grid spacing of 150 m. At this resolution, spurious coherent structures are often observed under unstable atmospheric conditions with moderate mean wind speeds. Results from such an intermediate resolution grid are often nested down for finer, more detailed LES, where these spurious structures adversely affect the development of turbulence on the subsequent finer grid nest. We therefore examine the impacts of the CPM on the representation of turbulence within the nested LES domain under moderate mean flow conditions in three different stability regimes: weakly convective, strongly convective, and weakly stable. In addition, two different resolutions of the underlying terrain are used to explore the role of the complex topography itself in generating turbulent structures. We demonstrate that the CPM improves the representation of turbulence within the LES domain, relative to the use of high-resolution complex terrain alone. During the convective conditions, the CPM improves the rate at which smaller-scales of turbulence form, while also accelerating the attenuation of the spurious numerically generated roll structures near the inflow boundary. During stable conditions, the coarse mesh spacing of the intermediate LES domain used herein was insufficient to maintain resolved turbulence using CPM as the flow develops downstream, highlighting the need for yet higher resolution under even weakly stable conditions, and the importance of accurate representation of flow on intermediate LES grids.

17 WIND ENERGY↗