Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “RuN”

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 757 records · Page 42

Turbulence-parameter estimation for current-energy converters using surrogate model optimization

Surrogate models maximize information utility by building predictive models in place of computational or experimentally expensive model runs. Marine hydrokinetic current energy converters require large-domain simulations to estimate array efficiencies and environmental impacts. Meso-scale models typically represent turbines as actuator discs that act as momentum sinks and sources of turbulence and its dissipation. An OpenFOAM model was developed where actuator disc k-ε turbulence was characterized using an approach developed for flows through vegetative canopies. Turbine-wake data from laboratory flume experiments collected at two influent turbulence intensities were used to calibrate parameters in the turbulence-source terms in the k-ε equations. Additionally, parameter influences on longitudinal wake profiles were estimated using Gaussian process regression with subsequent optimization minimizing the objective function within 3.1% of those obtained using the full model representation, but for 74% of the computational cost (far fewer model runs). This framework facilitates more efficient parameterization of the turbulence-source equations using turbine-wake data.

16 TIDAL AND WAVE POWER↗

Multi-decade high-resolution regional hindcasts for wave energy resource characterization in U.S. coastal waters

Long-term, high-resolution, regional wave hindcast datasets were generated using unstructured-grid Simulating WAves Nearshore (SWAN) models for the U.S. coastal waters to support nearshore wave energy development in the U.S. including those bordering U.S. territorial islands. The model domains resolved the entire U.S. exclusive economic zones, with a spatial resolution of approximately 200 m nearshore. The regional SWAN models were driven by the global WAVEWATCH III® model outputs and run for a 42-year period from 1979 to 2020. Extensive model validations were performed using buoy observations and altimeter data. Regional resource characterization was performed based on hindcast data points at 2 km from shore and along the 100 m isobath. Aggregations of wave resource parameters were produced, and spatial and seasonal variations were analyzed for all the regions. Wave resource metrics recommended by international standards, including a 3-hour time series of six resource parameters, hourly frequency- and directionally resolved wave spectra at selected “virtual buoy” locations, and average-annual values of omni-directional wave power, significant wave height, and energy period are publicly disseminated through an Amazon Web Service and a Marine Energy Atlas web application tool to facilitate wave energy research and a wide range of coastal ocean applications.

16 TIDAL AND WAVE POWER↗

An Open-Source Virtual Testbed for a Real Net-Zero Energy Community

Net zero energy communities (NZECs) are critical to ensure sustainability and resilience of modernized power systems. System modeling helps overcome technical challenges in designing and operating NZECs. In this paper, we present the modeling work based on a real NZEC. Two sets of models are developed: higher-fidelity physics-based models considering the interaction between subsystems of the studied NZEC and capturing fast-dynamics; and lower-fidelity data-driven models requiring less resource to establish and/or run. All models are validated against measurements from this real NZEC. In addition, we create a simulation framework which streamlines the processes for simulation and thus allows using developed models to form a virtual testbed. To demonstrate the usage of the virtual testbed, a case study is conducted where a building-to-grid integration control is evaluated via simulation. The evaluation results suggest that the tested control significantly smooths the power draw of the studied community and doesn’t sacrifice the thermal comfort to a great extent.

Huang, Sen↗

MPLBM-UT: Multiphase LBM library for permeable media analysis

MPLBM-UT is a specialized lattice-Boltzmann library that makes running single- and two-phase flow simulations in porous media accessible to everyone. We provide a suite of tools to pre-process computational domains for simulation, to set up custom boundary conditions, to run simulations, to post-process simulation outputs, and to visualize simulation results and data. All of these tools are easily accessible to users through the mplbm_utils Python package included in and automatically installed with MPLBM-UT. The high-performance, highly parallel library Palabos is used as the solver backend. MPLBM-UT is easily deployed in a variety of systems, from laptops to supercomputer clusters. MPLBM-UT also features multiple examples and benchmark templates that allow for fast prototyping of different porous media problems. We also provide an interface for reading in different file types and downloading domains from the Digital Rocks Portal to perform simulations.

97 MATHEMATICS AND COMPUTING↗

Dynamic compressed sensing for real-time tomographic reconstruction

Electron tomography has achieved higher resolution and quality at reduced doses with recent advances in compressed sensing. Compressed sensing (CS) exploits the inherent sparse signal structure to efficiently reconstruct three-dimensional (3D) volumes at the nanoscale from undersampled measurements. However, the process bottlenecks 3D reconstruction with computation times that run from hours to days. Here we demonstrate a framework for dynamic compressed sensing that produces a 3D specimen structure that updates in real-time as new specimen projections are collected. Researchers can begin interpreting 3D specimens as data is collected to facilitate high-throughput and interactive analysis. Using scanning transmission electron microscopy (STEM), we show that dynamic compressed sensing accelerates the convergence speed by ~3-fold while also reducing its error by 27% for a Au/SrTiO 3 nanoparticle specimen. Before a tomography experiment is completed, the 3D tomogram has interpretable structure within ~33% of completion and fine details are visible as early as ~66%. Upon completion of an experiment, a high-fidelity 3D visualization is produced without further delay. Additionally, we report reconstruction parameters that tune data fidelity can be manipulated throughout the computation without re-running the entire process.

47 OTHER INSTRUMENTATION↗

Process design tools and techno-economic analysis for capacitive deionization

Capacitive deionization (CDI) devices use cyclical electrosorption on porous electrode surfaces to achieve water desalination. Process modeling and design of CDI systems requires accurate treatment of the coupling among input electrical forcing, input flow rates, and system responses including salt removal dynamics, water recovery, energy storage, and dissipation. Techno-economic analyses of CDI further require a method to calculate and compare between a produced commodity (e.g. desalted water) versus capital and operational costs of the system. In this work, we demonstrate a new modeling and analysis tool for CDI developed as an installable Matlab program that allows direct numerical simulation of CDI dynamics and calculation of key performance and cost parameters. The program is provided for free and is used to run open-source Simulink models. The Simulink environment sends information to the program and allows for a drag and drop design space where users can connect CDI cells to relevant periphery blocks such as grid energy, battery, solar panel, waste disposal, and maintenance/labor cost streams. The program allows for simulation of arbitrary current forcing and arbitrary flow rate forcing of one or more CDI cells. We employ validated well-mixed reactor formulations together with a non-linear circuit model formulation that can accommodate a variety of electric double layer sub-models (e.g. for charge efficiency). The program includes a graphical user interface (GUI) to specify CDI plant parameters, specify operating conditions, run individual tests or parameter batch-mode simulations, and plot relevant results. The techno-economic models convert among dimensional streams of species (e.g. feed, desalted water, and brine), energy, and cost and enable a variety of economic estimates including levelized water costs.

42 ENGINEERING↗

Exploring the impacts of conformer selection methods on ion mobility collision cross section predictions

The prediction of structure dependent molecular properties, such as collision cross sections as measured using ion mobility spectrometry, are crucially dependent on the selection of the correct population of molecular conformers. Here, we report an in-depth evaluation of multiple conformation selection techniques, including simple averaging, Boltzmann weighting, lowest energy selection, low energy threshold reductions, and similarity reduction. Generating 50,000 conformers each for 18 molecules, we used the In Silico Chemical Library Engine (ISiCLE) to calculate the collision cross sections for the entire dataset. First, we employed Monte Carlo simulations to understand the variability between conformer structures as generated using simulated annealing. Then we employed Monte Carlo simulations to the aforementioned conformer selection techniques applied on the simulated molecular property—the ion mobility collision cross section. Based on our analyses, we found Boltzmann weighting to be a good tradeoff between precision and theoretical accuracy. Combining multiple techniques revealed that energy thresholds and root-mean-squared deviation-based similarity reductions can save considerable computational expense while maintaining property prediction accuracy. Molecular dynamic conformer generation tools like AMBER can continue to generate new lowest energy conformers even after tens of thousands of generations, decreasing precision between runs. This reduced precision can be ameliorated and theoretical accuracy increased by running density functional theory geometry optimization on carefully selected conformers.

Nielson, Felicity F.↗

Potential Adoption and Benefits of Co-Optimized Multimode Engines and Fuels for U.S. Light-Duty Vehicles

Exploring a diverse portfolio of technologies for decarbonization is crucial to understanding the potential impacts of different technological solutions and their associated environmental implications. Using high-octane, high-sensitivity biofuel blends in co-optimized multimode engines can increase engine efficiency and reduce vehicle emissions. Here, the multimode engine research focuses on the benefits of light-duty vehicle engines, which can operate in multiple modes depending on the vehicle's load. Low-temperature combustion can improve efficiency and reduce emissions (such as those from oxides of nitrogen and particulate matter) during low-load operation, while spark ignition performance is maintained in high-load operation. These advanced engines can be optimized to run on blends of biobased fuels. This analysis models scenarios for potential market adoption of co-optimized multimode vehicles fueled by three different bioblendstocks: ethanol, isopropanol, and isobutanol. An integrated modeling approach is used to forecast the energy and environmental impacts of the deployment of co-optimized multimode vehicles and fuels in the light-duty sector over the 2020-to-2050 time horizon. The multidisciplinary approach combines vehicle sales modeling, system dynamics modeling of the biorefining industry, and life cycle assessment to estimate the emissions and energy benefits. The models consider market forces such as consumer preferences for vehicle attributes, biofuel supply and demand dynamics subject to biorefinery capacity build-out and bioresource constraints, and forecasted changes to the U.S. bulk energy system over time. Market adoption of co-optimized vehicles is evaluated across a wide parameter space for incremental vehicle cost and engine efficiency improvement. This analysis reveals that the deployment of co-optimized multimode fuels and vehicles results in up to a 5% reduction in annual sector-wide life cycle greenhouse gas (GHG) emissions by 2050, relative to a business-as-usual scenario, but is also indicates environmental trade-offs, such as higher life cycle water-use. Emission benefits could potentially increase beyond 2050, as the new technologies penetrate the market and gain a foothold. Results also show that, under certain circumstances, vehicles with engines co-optimized for use with high-octane, high-sensitivity biofuel blends can be cost-competitive with conventional gasoline, while reducing GHG emissions. Our modeling results indicate that co-optimized multimode fuels and engines can be strategically leveraged in tandem with electrification to decarbonize the light-duty sector. Co-optimized vehicles could play a role in the early years of the time horizon, while electric vehicles (EVs) could become more competitive in the later years, highlighting the complementary benefits of these technologies for GHG reductions.

Oke, Doris↗

Simple but tricky: Investigations of terephthalic acid purity obtained from mixed PET waste

In this study, we report for the first time, the basic depolymerization of mixed waste-polyethylene terephthalate (PET) by hydrolysis and subsequent terephthalic acid monomer recovery at high purity using benign reaction conditions. Several conditions were tested for depolymerization such as PET chips size, concentration of aqueous sodium hydroxide (20 or 30%), organic co-solvent (ethylene glycol or ethanol), temperature at which the reaction was run and duration of the heating. More importantly, several batches of PET were utilized as starting materials including a commercial PET, chopped PET obtained from clean bottles whose caps and labels were removed and the purity of the product from each condition was evaluated via nuclear magnetic resonance (1HNMR and 13CNMR), differential scanning calorimetry (DSC) and powder Xray diffraction (XRD). Unsurprisingly, the conversion of PET is dependent on the particle size varying from 100% conversion for fine powder to 73% conversion for (300 µm, mesh 6 – 20 or mesh 14-20). Ethanol appears to be more efficient as a co-solvent than ethylene glycol, with higher PET depolymerization conversions (94% versus 75-80%), shorter reaction times (2h versus 6h) and lower temperatures (80 °C versus 110 °C). The terephthalic acid (TPA) recovered appeared to have only subtle differences among the batches, most notably a pink color when the reaction was run in ethanol/base. The DSC of the compounds produced in ethylene glycol water appear to display a melting point (280-288 °C) while the samples prepared in ethanol as well as a commercial sample did not. Overall, the purity of the various TPA batches is comparable, and similar to commercial TPA, demonstrating the utility of the method to depolymerize realistic waste streams. The method is simple, demonstrated on multigram scale (15-30g) and allows for the complete removal of waste other than PET unaffected by alkaline conditions.

Cosimbescu, Lelia↗

Mineral Scale Formation during Crossflow Reverse Osmosis at Constant Flux and Constant Transmembrane Pressure Conditions

Mineral scale formation on membrane surfaces is a significant challenge in reverse osmosis water purification. Laboratory fouling experiments are typically run such that the transmembrane pressure (TMP) is fixed, and the permeate flux decreases over time as scales accumulate on the membrane surface. However, this change in flux means that the hydrodynamic conditions at the membrane surface are continuously changing, which could affect crystallization and foulant deposition processes. Operating under constant permeate flux conditions, in contrast, is advantageous because it keeps the hydrodynamic conditions relatively consistent, making it possible to compare how membrane properties (e.g., surface chemistry) affect fouling propensity. Industrial reverse osmosis operations are not run strictly in either constant TMP or constant flux mode; while they may start at a constant TMP, feed pressure may be periodically adjusted to maintain permeate water production within a specified range. Here, the scarcity of constant permeate flux reverse osmosis scaling experiments reported in the literature frustrates efforts to compare membrane fouling processes under constant TMP and constant flux conditions. For the first time, the evolution of the fouling layer resistance was compared as a function of cumulative permeate volume per membrane area during constant TMP and constant flux reverse osmosis filtrations. Scaling experiments were conducted by challenging commercial reverse osmosis membranes with a model feed solution nearly saturated with calcium sulfate dihydrate (gypsum). At low fluxes, the increase in fouling layer resistance was quantitatively similar for the two operational modes. In contrast, at high fluxes, the fouling layer resistance increased more rapidly in constant flux filtration than in constant TMP filtration. The mechanism of scale formation in constant TMP and constant flux operation was self-limiting and self-reinforcing, respectively.

36 MATERIALS SCIENCE↗

A Deep Learning-Driven Sampling Technique to Explore the Phase Space of an RNA Stem-Loop

The folding and unfolding of RNA stem-loops are critical biological processes; however, their computational studies are often hampered by the ruggedness of their folding landscape, necessitating long simulation times at the atomistic scale. Here, we adapted DeepDriveMD (DDMD), an advanced deep learning-driven sampling technique originally developed for protein folding, to address the challenges of RNA stem-loop folding. Although tempering- and order parameter-based techniques are commonly used for similar rare-event problems, the computational costs or the need for a priori knowledge about the system often present a challenge in their effective use. DDMD overcomes these challenges by adaptively learning from an ensemble of running MD simulations using generic contact maps as the raw input. DeepDriveMD enables on-the-fly learning of a low-dimensional latent representation and guides the simulation toward the undersampled regions while optimizing the resources to explore the relevant parts of the phase space. We showed that DDMD estimates the free energy landscape of the RNA stem-loop reasonably well at room temperature. Our simulation framework runs at a constant temperature without external biasing potential, hence preserving the information on transition rates, with a computational cost much lower than that of the simulations performed with external biasing potentials. Here, we also introduced a reweighting strategy for obtaining unbiased free energy surfaces and presented a qualitative analysis of the latent space. This analysis showed that the latent space captures the relevant slow degrees of freedom for the RNA folding problem of interest. Finally, throughout the manuscript, we outlined how different parameters are selected and optimized to adapt DDMD for this system. We believe this compendium of decision-making processes will help new users adapt this technique for the rare-event sampling problems of their interest.

Gupta, Ayush↗

A Theory of Ultrafast Charge Transfer Relaxation with Non-Innocent Solvent Molecules

In this work we revisit the photodynamics of tetracyanoethylene-hexamethylbenzene (TCNE- HMB), the molecular complex studied by Hochstrasser et al. [J. Chem. Phys. 100, 4797–4810, 1994] that has long challenged the applicability of Marcus’s theory of elec- tron transfer for predicting photochemical reactions. Using a novel black-box electronic structure algorithm (time-dependent density functional theory with one double, TD- DFT-1D) to efficiently run molecular dynamics that can treat charge recombination, we run ab initio surface hopping molecular dynamics and confirm that, for a polar solvent, charge recombination rates can be incredibly fast (indeed faster than the sol- vent relaxation time); for non-polar solvents, the rate is much slower. We demonstrate that, although Marcus theory cannot be directly applied, these nonequilibrium (and sometimes incredibly fast) photoexcited dynamics can be effectively explained within a two-state model without any evidence of a transition through a conical intersection. Most importantly, for this paradigmatic model system, we are able to identify two nuclear coordinates of interest (rather than the single coordinate predicted by Marcus or a full set of internal quantum modes studied by Bixon and Jortner): the solvent relaxation in the first shell (that strongly modulates the energies of the charge trans- fer state and differentiates time scales for relaxation) and a nuclear displacement in the TCNE-HMB complex arising from a handful of vibrations that induces non-Born Oppenheimer motion and eventually facilitates an abrupt electronic transition to the ground state. Altogether, these findings suggest a tractable generalization of Marcus theory for future simulations of photochemistry with non-innocent solvent environ- ments in the spirit of a Hamiltonian suggested by Stuchebrukhov (J. Chem. Phys. 107, 3821, 1997).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Effect of Flexibility in Molecular Simulations of Carbon Dioxide Adsorption and Diffusion in a Cuprous Triazolate Framework

Using fixed atom force fields to model gas adsorption in flexible metal organic frameworks (MOFs) is known to pose difficulties in accurately reproducing and predicting experimental results. This paper studies the difference in accuracy between flexible and fixed atom force fields in reproducing CO 2 adsorption measurements in MAF-2 ([Cu(etz)]∞ (MAF-2, Hetz) 3,5-diethyl-1,2,4-triazole), an NbO-type triazolate scaffold with a bcu cavity system and attached ethyl groups. The flexible force field used to run the hybrid molecular dynamics and grand canonical Monte Carlo calculations were generated using the QuickFF software incorporating van der Waals parameters from the Universal Force Field (UFF) and density derived electrostatic and chemical (DDEC) charges. The fixed atom force field used was composed of UFF van der Waals parameters together with DDEC charges. The calculations were run at 298 K and at pressures of 0.1, 0.3, 0.5, 0.8, and 1 bar. It was observed that for this MOF the rigid force field overpredicted gas adsorption, whereas the flexible force field values closely matched experimental data. In the flexible structure, the freely moving ethyl groups of MAF-2 hindered adsorption, reducing the interaction energy between CO 2 and the N atoms of the triazolate framework as well as reducing the size of the largest cavity diameter. Here, the combination of these factors led to improved prediction of adsorption values with the flexible forcefield as compared to the rigid forcefield, demonstrating the need for modeling MOFs in a way more indicative of their behavior.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

High-Resolution Comonomer Sequencing of Blocky Brominated Syndiotactic Polystyrene Copolymers Using 13C NMR Spectroscopy and Computer Simulations

This work demonstrates the first high-resolution comonomer sequencing of Blocky brominated syndiotactic polystyrene (sPS-co-sPS-Br) copolymers based on pentad assign- ments of the quaternary carbon region of the nuclear magnetic resonance spectrum. Copolymers containing p-bromostyrene (Br-Sty) units were prepared in matched sets using postpolymerization bromination methods carried out in the heterogeneous gel state (Blocky) and homogeneous solution state (Random). Quantitative information from the quaternary carbon spectra, heteronuclear multiple bond correlation spectroscopy, electronic structure calculations, and simulated statistically random copolymers was correlated to confirm the carbon resonance assignments for all 20 possible pentad comonomer sequences. Using the experimental pentad sequence prevalences, a computer code was developed to simulate chains with microstructures typical of each sample as a means to visually represent the copolymer blockiness with quantitative precision. Based on the microstructure and distribution of run lengths in these chains, the simulations revealed that the Blocky copolymers contain a high degree of blockiness. By comparing the run lengths in the simulated chains to the average number of styrene units in a crystalline segment of sPS (found by small-angle X-ray scattering), copolymer crystallizability was predicted. For the simulated Blocky B-21% (21 mol % Br-Sty) chain, the probability of randomly selecting a styrene unit in a crystallizable block was 25.8%, while that in the simulated Random R-18% was zero, in excellent agreement with the experimental crystallization behavior measured by differential scanning calorimetry. Additionally, these predictions confirmed that the simulated chains accurately represent the ensemble of chains in their respective copolymer samples. Furthermore, each simulated Blocky chain contained one or more long sPS blocks that paralleled the measured 38-40 styrene units spanning a crystalline segment within the sPS/CCl4 gel. This finding affirmed that the long sPS segments originated from the precise lamellar structure within the heterogeneous gel morphology (i.e., block length is correlated with lamellar thickness). Overall, the ability to tailor the copolymer microstructure through control of the semicrystalline gel morphology opens the door to synthesizing ordered copolymers by postpolymerization functionalization processes with unprecedented levels of compositional control.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dynamic Life-Cycle Analysis of Fast Pyrolysis Biorefineries: Impacts of Feedstock Moisture Content and Particle Size

The commercialization of biorefineries has been met with significant obstacles due to the technical difficulties of handling solid biomass feedstocks, which are highly variable in physical properties and chemical compositions. An understanding of the principles of unit operations along the supply chain is key to the success of the biofuel industry. This study applies a dynamic life-cycle analysis (DLCA) methodology by developing quantitative relationships between the inputs and outputs for each unit operation based on scientific understanding of the causal effects. Using the DLCA, we assessed system sustainability of drop-in fuel production from fast pyrolysis of pine residues followed by hydro-processing. Life-cycle greenhouse gas (GHG) emissions were calculated for 4641 runs involving two key feedstock parameters: moisture content after field drying and particle size of the feed to pyrolyzer. Moisture contents varied between 25% and 35%, while particle sizes varied between 0.5 and 5 mm in these runs. The life-cycle GHG emissions of these cases vary from 22 to 40 g/MJ. Close examination of each unit operation's contribution to the total GHG emissions reveals that low moisture content after field drying produces savings in energy consumption during feedstock transportation and on-site drying. However, although small particle size leads to overall higher fuel yield, it also requires a significant amount of energy for feedstock size reduction and fuel production. The energy penalty outweighs the benefit of increased fuel yields, especially when fine particles (<1 mm) are used for pyrolysis; thus, small particle size leads to increased GHG emissions overall. The results highlight the trade-offs between the energy demand for preprocessing and the conversion yields, which can be addressed with DLCA.

09 BIOMASS FUELS↗

Autonomous Nanoparticle Synthesis Guided by In Situ Multiscale Structural Characterization

Autonomous synthesis platforms promise rapid exploration of vast parameter spaces; yet, integrating in situ structural characterization in closed-loop synthesis optimization remains challenging. We demonstrate a realization of such a closed-loop platform coupled with a droplet-flow microreactor, in situ X-ray scattering methods (SAXS/WAXS), and Gaussian process optimization to synthesize citrate-reduced Au nanoparticles with targeted characteristics. The system efficiently explored ∼19,000 synthesis recipes through 365 experiments, achieving precise control over size (4–60 nm) and polydispersity (σ < 0.11) across large citrate/gold ratios, exceeding traditional synthesis boundaries (1–10). Beyond confirming classical Turkevich–Frens trends, partial-dependence analysis revealed strong nonlinear coupling among precursor, citrate, and pH effects. Combining quantitative SAXS/WAXS analysis with electron microscopy characterization, we uncovered that crystallite size (d c ) and particle size (d) follow d c = 0.18d + β, where synthesis chemistry controls the intercept β while maintaining a universal slope. This parallel-band structure enables independent tuning of crystallite domain size at fixed particle diameter through a combination of chloride, gold precursor, citrate, and pH contributions (cross-validated Spearman ρ = 0.7 ± 0.1). High-resolution electron microscopy shows multiple lattice-fringe orientations within single particles, directly confirming polycrystalline domains and the ability to tune d c at the fixed d. The platform’s validation includes indistinguishable static versus flowing measurements, stable droplet transport at 100 °C, and <5% run-to-run variation, establishing a robust framework for mapping and controlling multiscale nanoparticle structure across expansive chemical spaces. In conclusion, the developed closed-loop platform can be applied to a borad range of nanosyntheis processes.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Microphysical Sensitivity of Superparameterized Precipitation Extremes in the Contiguous United States Due to Feedbacks on Large-Scale Circulation

Superparameterized (SP) global climate models have been shown to better simulate various features of precipitation relative to conventional models, including its diurnal cycle as well as its extremes. While various studies have focused on the effect of differing microphysics parameterizations on precipitation within limited-area cloud-resolving models, we examine here the effect on contiguous U.S. (CONUS) extremes in a global SP model. We vary the number of predicted moments for hydrometeor distributions, the character of the rimed ice species, and the representation of raindrop self-collection and breakup. Using a likelihood ratio test and accounting for the effects of multiple hypothesis testing, we find that there are some regional differences, particularly during spring and summer in the Southwest and the Midwest, in both the current climate and a warmer climate with uniformly increased sea surface temperatures. These differences are most statistically significant and widespread when the number of moments is changed. To determine whether these results are due to (fast) local effects of the different microphysics or the (slower) ensuing feedback on the large-scale atmospheric circulation, we run a series of short, 5-day simulations initialized from reanalysis data. We find that the differences largely disappear in these runs and therefore infer that the different parameterizations impact precipitation extremes indirectly via the large-scale circulation. Finally, we compare the present-day results with hourly rain gauge data and find that SP underestimates extremes relative to observations regardless of which microphysics scheme is used given a fixed model configuration and resolution.

54 ENVIRONMENTAL SCIENCES↗

The Leading Modes of Asian Summer Monsoon Variability as Pulses of Atmospheric Energy Flow

Monsoon rainfall anomalies are often organized into dynamical modes that are more predictable than rainfall at specific locations. Here a linear response function (LRF) is constructed to extract the neutral modes (NM) of the dynamical atmospheric system and their optimal forcing that govern precipitation variability, providing a dynamical underpinning for the empirically derived, recurrent rainfall patterns. The leading mode (NM1) features a shift in the Atlantic ITCZ, a drying over the Indian Ocean and Indian subcontinent, and a “South-Flood-North Drought (SFND)” pattern over eastern China, a pattern inferable from changes in the atmospheric energy flow forced by the corresponding optimal forcing and further enhanced by the water vapor and cloud feedbacks. Moreover, the slow evolving component of the NM1 in response to CO2 forcing can be skillfully predicted by running the ocean heat uptake through the LRF, implicating a deterministic ocean dynamic source for the NM1-like trend under climate warming.

monsoon mode of variability, EBM, climate change, ↗