Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “TIME DEPENDENCE”

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,351 records · Page 75
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Size-Controlled Nanoparticles Embedded in a Mesoporous Architecture Leading to Efficient and Selective Hydrogenolysis of Polyolefins

A catalytic architecture, comprising a mesoporous silica shell surrounding platinum nanoparticles (NPs) supported on a solid silica sphere (mSiO2/Pt-X/SiO2; X is the mean NP diameter), catalyzes hydrogenolysis of melt-phase polyethylene (PE) into a narrow C23-centered distribution of hydrocarbons in high yield using very low Pt loadings (~10 -5 g Pt/g PE). During catalysis, a polymer chain enters a pore and contacts a Pt NP where the C–C bond cleavage occurs and then the smaller fragment exits the pore. mSiO 2 /Pt/SiO 2 resists sintering or leaching of Pt and provides high yields of liquids; however, many structural and chemical effects on catalysis are not yet resolved. Here, we report the effects of Pt NP size on activity and selectivity in PE hydrogenolysis. Time-dependent conversion and yields and a lumped kinetics model based on the competitive adsorption of long vs short chains reveal that the activity of catalytic material is highest with the smallest NPs, consistent with a structure-sensitive reaction. Remarkably, the three mSiO 2 /Pt-X/SiO 2 catalysts give equivalent selectivity. We propose that mesoscale pores in the catalytic architecture template the C 23 -centered distribution, whereas the active Pt sites influence the carbon–carbon bond cleavage rate. This conclusion provides a framework for catalyst design by separating the C–C bond cleavage activity at catalytic sites from selectivity for chain lengths of the products influenced by the structure of the catalytic architecture. The increased activity, selectivity, efficiency, and lifetime obtained using this architecture highlight the benefits of localized and confined environments for isolated catalytic particles under condensed-phase reaction conditions.

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

Time-Evolving Chirality Loss in Molecular Photodissociation Monitored by X-ray Circular Dichroism Spectroscopy

In this work, the ultrafast photoinduced chirality loss of 2-iodobutane is studied theoretically by time- and frequency-resolved X-ray circular dichroism (TRXCD) spectroscopy. Following an optical excitation, the iodine atom dissociates from the chiral center, which we capture by quantum non-adiabatic molecular dynamics simulations. At variable time delays after the pump, the resonant X-ray pulse selectively probes the iodine and carbon atom involved in the chiral dissociation through a selected core-to-valence transition. The TRXCD signal at the iodine L 1 edge accurately captures the timing of C–I photodissociation and thereby chirality loss, c.a 70 fs. The strong electric dipole–electric quadrupole (ED–EQ) response makes this signal particularly sensitive to vibronic coherence at the high X-ray regime. At the carbon K-edges, the signals monitor the molecular chirality of the 2-butyl radical photoproduct and the spin state of the iodine atom. The ED–EQ response is masked under the strong electric dipole-magnetic dipole response, making this signal intuitive for the electronic population. The evolution of the core electronic states and its chiral sensitivity is discussed. Overall, the element-specific TRXCD signal provides a detailed picture of molecular dynamics and offers a unique sensitive window into the time-dependent chirality of molecules.

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

Chemically Fueled Reinforcement of Polymer Hydrogels

Carbodiimide-fueled anhydride bond formation has been used to enhance the mechanical properties of permanently crosslinked polymer networks, giving materials that exhibit transitions from soft gels to covalently reinforced gels, eventually returning to the original soft gels. Temporary changes in mechanical properties result from a transient network of anhydride crosslinks, which eventually dissipate by hydrolysis. Over an order of magnitude increase in the storage modulus is possible through carbodiimide fueling. Here, the time-dependent mechanical properties can be modulated by the concentration of carbodiimide, temperature, and primary chain architecture. Because the materials remain rheological solids, new material functions such as temporally controlled adhesion and rewritable spatial patterns of mechanical properties have been realized.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Bottlebrush Block Copolymers at the Interface of Immiscible Liquids: Adsorption and Lateral Packing

Amphiphilic bottlebrush block copolymers (BBCPs), having a hydrophilic bottlebrush polymer (BP) linked covalently to a hydrophobic BP, were found to segregate to liquid-liquid interfaces to minimize the free energy of the system. The key parameter influencing the outcome of the experiments is the ratio between the degree of polymerization of the backbone (NBB) and that of the side-chain brushes (NSC). Specifically, a spherical, star-like configuration results when NBB < NSC, while a cylindrical, bottlebrush-like shape is preferred when NBB > NSC. Dynamic interfacial tension (γ) and fluorescence recovery after photobleaching (FRAP) measurements show that the BBCP configuration influences the areal density and in-plane diffusion at the fluid interface. The characteristic relaxation times associated with BBCP adsorption (τA) and reorganization (τR) were determined by fitting time-dependent interfacial tension measurements to a sum of two exponential relaxation functions. Both τA and τR initially increased with NBB up to 92 repeat units, due to the larger hydrodynamic radius in solution and slower in-plane diffusivity, attributed to a shorter cross-sectional diameter of the side-chains near the block junction. This trend reversed at NBB = 190, with shorter τA and τR attributed to increased segregation strength and exposure of the bare water/toluene interface due to tilting and/or wiggling of the backbone chains, respectively. The adsorption energy barrier decreased with higher NBB, due to a reduced BBCP packing density at the fluid interface. In conclusion, this study provides fundamental insights into macromolecular assembly at fluid interfaces, as it pertains to unique bottlebrush block architectures.

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

Overdestabilization vs Overstabilization in the Theoretical Analysis of f-Orbital Covalency

The complex nature of the f-orbital electronic structures and their interaction with the chemical environment pose significant computational challenges. Advanced computational techniques that variationally include scalar relativities and spin–orbit coupling directly at the molecular orbital level have been developed to address this complexity. Among these, variational relativistic multiconfigurational multireference methods stand out for their high accuracy and systematic improvement in studies of f-block complexes. Additionally, these advanced methods offer the potential for calibrating low-scaling electronic structure methods such as density functional theory. However, studies on the Cl K-edge X-ray absorption spectra of the [Ce(III)Cl 6 ] 3– and [Ce(IV)Cl 6 ] 2– complexes show that time-dependent density functional theory with approximate exchange–correlation kernels can lead to inaccuracies, resulting in an overstabilization of 4f orbitals and incorrect assessments of covalency. In contrast, approaches utilizing small active space wave function methods may understate the stability of these orbitals. The results herein demonstrate the need for large active space, multireference, and variational relativistic methods in studying f-block complexes.

chemical calculations↗
DOE OSTI2024

Mixed Valence {Ni 2+ Ni 1+ } Clusters as Models of Acetyl Coenzyme A Synthase Intermediates

Acetyl coenzyme A synthase (ACS) catalyzes the formation and deconstruction of the key biological metabolite, acetyl coenzyme A (acetyl-CoA). The active site of ACS features a {NiNi} cluster bridged to a [Fe4S4] n+ cubane known as the A-cluster. The mechanism by which the A-cluster functions is debated, with few model complexes able to replicate the oxidation states, coordination features, or reactivity proposed in the catalytic cycle. In this work, we isolate the first bimetallic models of two hypothesized intermediates on the paramagnetic pathway of the ACS function. The heteroligated {Ni 2+ Ni 1+ } cluster, [K(12-crown-4) 2 ][1], effectively replicates the coordination number and oxidation state of the proposed “A red ” state of the A-cluster. Addition of carbon monoxide to [1] - allows for isolation of a dinuclear {Ni 2+ Ni 1+ (CO)} complex, [K(12-crown-2) n ][2] (n = 1–2), which bears similarity to the “A NiFeC ” enzyme intermediate. Structural and electronic properties of each cluster are elucidated by X-ray diffraction, nuclear magnetic resonance, cyclic voltammetry, and UV/vis and electron paramagnetic resonance spectroscopies, which are supplemented by density functional theory (DFT) calculations. Calculations indicate that the pseudo-T-shaped geometry of the three-coordinate nickel in [1] – is more stable than the Y-conformation by 22 kcal mol –1 , and that binding of CO to Ni 1+ is barrierless and exergonic by 6 kcal mol –1 . UV/vis absorption spectroscopy on [2] - in conjunction with time-dependent DFT calculations indicates that the square-planar nickel site is involved in electron transfer to the CO π*-orbital. Further, we demonstrate that [2] - promotes thioester synthesis in a reaction analogous to the production of acetyl coenzyme A by ACS.

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

Metal–Ligand Covalency in the Valence Excited States of Metal Dithiolenes Revealed by S 1s3p Resonant Inelastic X-ray Scattering

Metallo dithiolene complexes with biological and catalytic relevance are well-known for having strong metal–ligand covalency, which dictates their valence electronic structures. We present the resonant sulfur Kβ (1s3p) X-ray emission spectroscopy (XES) for a series of Ni and Cu bis(dithiolene) complexes to reveal the ligand sulfur contributions to both the occupied and unoccupied valence orbitals. While S K-edge X-ray absorption spectroscopy played a critical role in identifying the covalency of the unoccupied orbitals of metal dithiolenes, the present focus on XES explores the occupied density of states. For a series of [Cu(mnt) 2 ] n– and [Ni(mnt) 2 ] n– anions and dianions, a comparison of the nonresonant and resonant S Kβ XES spectra highlights the dramatic improvement in spectral resolution and corresponding ability to differentiate subtle changes in occupied electronic structure across the series. Furthermore, the use of resonant inelastic X-ray scattering (RIXS) probes the valence excited states and the core–valence couplings of the complexes. By employing a theoretical approach based on time-dependent density functional theory to interpret the RIXS spectra, we reveal how metal–ligand covalency influences the excited state energies and covalencies. We identify the low energy excited states as having the same symmetry as the nominal “ligand field” or “d–d” states that typically dominate the photophysics of 3d metal complexes but with significant metal– ligand charge transfer character dictated by their covalency. These results suggest that strong metal–ligand covalency can be used to influence the charge-transfer photochemistry of first row transition metal complexes.

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

Structural Evidence of Interanionic Hydrogen Bonding in Phosphoric Acid Solutions

Interanionic hydrogen bonding (IAHB) is a noncovalent interaction between like-charged ions that challenges conventional electrostatic understanding. This study provides direct structural evidence of IAHB in concentrated aqueous phosphoric acid (PA) solutions, which exhibit >60% dissociation under these conditions. Oxygen K-edge X-ray absorption fine structure spectroscopy, combined with electron affinity time-dependent density functional theory calculations, reveals the formation of stable, cyclic phosphate-phosphate IAHB dimers at PA concentrations ≥7 M. Extended X-ray absorption fine structure data show distinct long-range ordering consistent with these dimers, and near-edge X-ray absorption fine structure spectra confirm a concentration-dependent transition from monomeric to dimeric species. Energy decomposition analysis through density functional theory shows that the formation of solution-phase IAHB is energetically favored and is attributed to polarization of and the charge transfer between the two fragments driven by the surrounding solvent molecules, in addition to the permanent electrostatics. These findings offers crucial structural insights into the H-bonded networks in concentrated PA, highlighting the critical role of solvent in facilitating anion–anion association.

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

Helical Complex Ladder Polymer with Amplification of Asymmetry

Here, we report a Cu–salen helical complex ladder polymer (HLP-Cu) built from covalently fused Cu(Salen) units. Point chirality from enantiopure diamines directs the formation of a rigid, π-conjugated helix. Compared with a planar analogue (LLP-Cu) and the Cu(Salen) monomer, HLP-Cu exhibits up to 5-fold stronger circular dichroism and a maximum dissymmetry factor of 7.4 × 10 –3 . Titration experiments suggest a cooperative amplification of asymmetry arising from extended conjugation and helical geometry during helix formation. Time-dependent density functional theory calculations are employed to support this hypothesis, showing that increasing conjugation length and helical order enhance the rotational strength of contributing excited states via elevated magnetic transition moments and improved alignment of electric and magnetic transition moments. Our findings provide a new polymer platform for designing functional chiral materials with tunable chiroptical and spin-dependent properties.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Plasmon-Assisted Electrochemical Epoxidation using Water as an Oxidant

Olefin epoxidation, an important industrial reaction, often uses hazardous oxidants, causing challenges in waste disposal. Here, we demonstrate the use of water as an oxidant by a plasmon-assisted electrochemical strategy. The electrocatalyst is comprised of a hybrid of a water oxidation catalyst, manganese oxide, and plasmonic gold nanoparticles. Visible-light irradiation of the electrocatalyst enhanced the epoxidation of 4-styrenesulfonate 5-fold as compared to dark conditions at the same temperature. From electrochemical analyses conducted under plasmon excitation conditions, complemented by real-time time-dependent density functional tight binding simulations, it is found that the plasmonic boost of the electrochemical styrene epoxidation is due to energetic holes generated by the excitation of localized surface plasmon resonances of gold nanoparticles. These photogenerated holes activate adsorbed water for oxidation and enhance the binding of 4-styrenesulfonate at interfacial sites. Here, this work demonstrates a proof of concept and establishes the mechanistic basis for plasmon-assisted activation of water as an O atom source for electrochemical epoxidations.

Gold↗
DOE OSTI2020

Enhanced Menshutkin S N 2 Reactivity in Mesoporous Silica: The Influence of Surface Catalysis and Confinement

A significant enhancement in the Menshutkin S N 2 reaction between 1-methylimidazole (MeIm) and methyl thiocyanate (MeSCN) is observed when the reaction is confined in the nanoscale silica pores of MCM41 and SBA15. The experiments in the silica pores are conducted without surrounding bulk reaction mixture. The influences of temperature, pore radius, and surface chemistry on the kinetics of the confined reaction are analyzed with time-dependent infrared spectroscopy, molecular dynamics simulations, and ab initio calculations. The rate constant of the pseudo-first order reaction increases with decreasing pore size, and the activation energy is found to decrease by 5.6 kJ/mol in the smallest pore studied (2.8 nm) relative to the bulk reaction. The rate constant dependence on pore size is accurately described by a two-state model in which molecules within the 4.6 Å interfacial layer experience a 2.4-fold rate constant increase relative to those reacting at the bulk rate further away from the interface. The removal of polar silanol groups from the silica surface via passivation with trimethylsilyl chloride results in bulk-like kinetics despite a reduction in the pore diameter, demonstrating the role of silanols as catalytic sites. Here, electronic structure calculations of the energy profile on a model silica surface confirm that silanol groups, particularly those of the vicinal type, can reduce the activation energy and reaction endothermicity through the donation of hydrogen bonds to the reactant, transition state, and product complexes.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Evaluation of Waveform Profiles for Traveling Wave Ion Mobility Separations in Structures for Lossless Ion Manipulations

Structures for lossless ion manipulations (SLIM) have recently enabled a powerful implementation of traveling wave ion mobility spectrometry (TWIMS) for ultra-high resolution separations. However, optimization of experimental parameters remains challenging. Most TWIMS separations have utilized square-shaped waveforms applied by time-dependent voltage stepping across a repeating sets of electrodes; however, alternative waveforms may provide further improvements to resolution. In this work, we characterize five waveforms (including square and sine) in terms of their transmission efficiency, IMS resolution, and resolving power, and explore the effects of TW amplitude and speed on the performance of each. We found, consistent with previous work, separations were generally improved with higher TW amplitudes, moderately improved by lower speeds (limited by ion ‘surfing’ with the waves), and found decreases in signal intensity at the extremes of operating conditions. The triangle and asymmetric ‘ramp forward’ shaped profiles were found to provide modestly greater resolution and resolving power, an observation we tentatively attribute to their relatively uniform fields and minimal low-field regions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Acoustic Energy Release During the Laboratory Seismic Cycle: Insights on Laboratory Earthquake Precursors and Prediction

Abstract Machine learning can predict the timing and magnitude of laboratory earthquakes using statistics of acoustic emissions. The evolution of acoustic energy is critical for lab earthquake prediction; however, the connections between acoustic energy and fault zone processes leading to failure are poorly understood. Here, we document in detail the temporal evolution of acoustic energy during the laboratory seismic cycle. We report on friction experiments for a range of shearing velocities, normal stresses, and granular particle sizes. Acoustic emission data are recorded continuously throughout shear using broadband piezo‐ceramic sensors. The coseismic acoustic energy release scales directly with stress drop and is consistent with concepts of frictional contact mechanics and time‐dependent fault healing. Experiments conducted with larger grains (10.5 μm) show that the temporal evolution of acoustic energy scales directly with fault slip rate. In particular, the acoustic energy is low when the fault is locked and increases to a maximum during coseismic failure. Data from traditional slide‐hold‐slide friction tests confirm that acoustic energy release is closely linked to fault slip rate. Furthermore, variations in the true contact area of fault zone particles play a key role in the generation of acoustic energy. Our data show that acoustic radiation is related primarily to breaking/sliding of frictional contact junctions, which suggests that machine learning‐based laboratory earthquake prediction derives from frictional weakening processes that begin very early in the seismic cycle and well before macroscopic failure.

58 GEOSCIENCES↗
DOE OSTI2020

Amorphization and Plasticity of Olivine During Low-Temperature Micropillar Deformation Experiments

Experimentally quantifying the viscoplastic rheology of olivine at the high stresses and low temperatures of the shallow lithosphere is challenging due to olivine's propensity to deform by brittle mechanisms at these conditions. In this study, we use microscale uniaxial compression tests to investigate the rheology of an olivine single crystal at room pressure and temperature. Pillars with nominal diameters of 1.25 μm were prepared using a focused ion beam milling technique and were subjected to sustained axial stresses of several gigapascal. The majority of the pillars failed after dwell times ranging from several seconds to a few hours. However, several pillars exhibited clear evidence of plastic deformation without failure after 4–8 hr under load. The corresponding creep strain rates are estimated to be on the order of 10-6 to 10 -7 s -1 . The uniaxial stresses required to achieve this deformation (4.1–4.4 GPa) are in excellent agreement with complementary data obtained using nanoindentation techniques. Scanning transmission electron microscopy observations indicate that deformation occurred along amorphous shear bands within the deformed pillars. Electron energy loss spectroscopy measurements revealed that the bands are enriched in Fe and depleted in Mg. We propose that inhomogeneities in the cation distribution in olivine concentrate stress and promote the amorphization of the Fe-rich regions. Finally, the time dependence of catastrophic failure events suggests that the amorphous bands must grow to some critical length scale to generate an unstable defect, such as a shear crack.

58 GEOSCIENCES↗
DOE OSTI2019

Mathematical Reconstruction of Land Carbon Models From Their Numerical Output: Computing Soil Radiocarbon From C Dynamics

Radiocarbon ( 14 C) is a powerful tracer of the global carbon cycle that is commonly used to assess carbon cycling rates in various Earth system reservoirs and as a benchmark to assess model performance. Therefore, it has been recommended that Earth System Models (ESMs) participating in the Coupled Model Intercomparison Project Phase 6 report predicted radiocarbon values for relevant carbon pools. However, a detailed representation of radiocarbon dynamics may be an impractical burden on model developers. Here, we present an alternative approach to compute radiocarbon values from the numerical output of an ESM that does not explicitly represent these dynamics. The approach requires computed 12 C stocks and fluxes among all carbon pools for a particular simulation of the model. From this output, a time-dependent linear compartmental system is computed with its respective state-transition matrix. Using transient atmospheric 14 C values as inputs, the state-transition matrix is then applied to compute radiocarbon values for each pool, the average value for the entire system, and component fluxes. We demonstrate the approach with ELMv1-ECA, the land component of an ESM model that explicitly represents 12 C, and 14 C in 7 soil pools and 10 vertical layers. Results from our proposed method are highly accurate (relative error <0.01%) compared with the ELMv1-ECA 12 C and 14 C predictions, demonstrating the potential to use this approach in CMIP6 and other model simulations that do not explicitly represent 14 C.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Adaptive Multirate Mass Transfer (aMMT) Model: A New Approach to Upscale Regional‐Scale Transport Under Transient Flow Conditions

Abstract The long‐term evaluation of regional‐scale groundwater quality needs efficient upscaling methods for transient flow. Upscaling techniques, such as the Multirate Mass Transfer (MRMT) method with constant upscaling parameters, have been used for transport with steady‐state flow, yet the upscaling parameters (i.e., rate coefficients) may be time dependent. This study proposed and validated an adaptive MRMT (aMMT) method by allowing the mass transfer coefficients in MRMT to change with the flow field. Advective‐dispersive contaminant transport simulated in a 3‐D heterogeneous medium was used as a reference solution. Equivalent transport under homogeneous flow conditions was evaluated by applying the MRMT and aMMT models for upscaling. The relationship between mass transfer coefficients and flow rates was fitted under steady‐state flow driven by various hydraulic gradients. A power law relationship was obtained, which was then used to update the mass transfer coefficients in each stress period under transient flow conditions in the aMMT method. Results indicated that for advection‐dominated transport, both the MRMT and aMMT methods can upscale the anomalous transport dynamics affected by subgrid heterogeneity under transient flow conditions. Whereas for diffusion‐dominated systems, the MRMT model failed to capture the tails of tracer breakthrough curves after the boundary condition changed, but the results from the aMMT model were significantly improved. However, if the overall flow direction changed, both MRMT and aMMT failed to represent the breakthrough curve tail generated by the heterogeneous system. The results point toward a promising path for upscaling transport in complex aquifers with transient flow.

Guo, Zhilin↗
DOE OSTI2021

Covariation of Slab Tracers, Volatiles, and Oxidation During Subduction Initiation

Abstract Subduction‐related lavas have higher Fe 3+ /∑Fe than midocean ridge basalts (MORB). Hypotheses for this offset include imprint from subducting slabs and differentiation in thickened crust. These ideas are readily tested through examination of the time‐dependent evolution of slab‐derived signatures, thickening crust of the overriding plate, and evolving redox during subduction initiation. Here, we present Fe 3+ /ΣFe and volatile element abundances of volcanic glasses recovered from International Ocean Discovery Program (IODP) Expedition 352 to the Izu‐Bonin‐Mariana (IBM) forearc. The samples include forearc basalts (FAB) that are stratigraphically overlain by low‐ and high‐silica boninite lavas. The FAB glasses have 0.18–0.85 wt% H 2 O, 75–233 ppm CO 2 , S contents controlled by saturation with a sulfide phase (602–1,386 ppm), Ba/La from 3.9‐10, and Fe 3+ /ΣFe ratios from 0.136 to 0.177. These compositions are similar to MORB and suggest that decompression melting of dry and reduced mantle dominates the earliest stages of subduction initiation. Low‐ and high‐silica boninite glasses have 1.51–3.19 wt% H 2 O, CO 2 below detection, S contents below those required for sulfide saturation (5–235 ppm), Ba/La from 11 to 29, and Fe 3+ /∑Fe from 0.181 to 0.225. The compositions are broadly similar to modern arc lavas in the IBM arc. These data demonstrate that the establishment of fluid‐fluxed melting of the mantle, which occurs in just 0.6–1.2 my after subduction initiation, is synchronous with the production of oxidized, mantle‐derived magmas. The coherence of high Fe 3+ /∑Fe and Ba/La ratios with high H 2 O contents in Expedition 352 glasses and the modern IBM arc rocks strongly links the production of oxidized arc magmas to signatures of slab dehydration.

Brounce, Maryjo↗
DOE OSTI2022

The FastEddy® Resident-GPU Accelerated Large-Eddy Simulation Framework: Moist Dynamics Extension, Validation and Sensitivities of Modeling Non-Precipitating Shallow Cumulus Clouds

Herein we describe the moist dynamics formulation implemented within the graphics processing unit-resident large-eddy simulation FastEddy® model, which includes a simple saturation adjustment scheme for condensation and evaporation processes. Two LES model intercomparison exercises for non-precipitating shallow cumulus clouds are simulated in order to validate this model extension, including a static forcing and a time-dependent forcing case. Overall, we find our dynamical, thermodynamical and microphysical quantities, along with turbulence variability and fluxes, to be commensurate with the corresponding model intercomparison results. In addition, sensitivities to specific model settings are investigated. Among these settings, it is shown that boundary layer and cloud layer structure and characteristics are sensitive to use of higher-order advection schemes impacting the vertical distribution of cloud content and associated turbulence statistics. Increasing the timescale of the saturation scheme leads to enhanced liquid water presence and decreases vertical velocity variance within the cloud deck. In some cases, these sensitivities agree with the model-to-model variability reported in the intercomparison exercises, highlighting the important role of specific model implementation choices in the context of shallow cumulus convection simulations. These analyses and findings also provide the basis for future extensions and applications of FastEddy® for modeling moist convection and precipitation scenarios.

54 ENVIRONMENTAL SCIENCES↗