Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Static analysis”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Theory of Ion-Mediated Segmental Localization, Activated Structural Relaxation, and the Glass Transition in Polymerized Ionic Liquids

Here we combine polymer integral equation theory for structural correlations with a force-level theory for activated segmental relaxation to construct a microscopic theory for how ions modify the onset of segmental localization, the polymer alpha relaxation time, and glass transition temperature, T g , in dry polymerized ionic liquids (PolyILs). Our most important findings concerning the behavior of T g are (i) a weak dependence on anion–cation Coulomb attraction strength due to strong decoupling of ion and segmental dynamics, (ii) strong reduction with increasing mobile ion size due to both plasticization effects and weakening of ion-modified monomer caging, and (iii) increase with backbone rigidity and intrachain dynamic cooperativity. We quantitatively estimate T g variations as ion size and polymer static and dynamic rigidity are widely varied. The experimental data of 19 PolyILs are organized into a specific pattern consistent with the theoretical predictions with glass transition temperatures varying over a very wide range (~250 K) as observed experimentally. Based on the theoretical results, suggestions are made concerning the search for high conductivity Li- or Na-based PolyILs based on further lowering of T g via new polymer synthesis. Calculations of the alpha time as a function of packing fraction reveal the dynamic fragility increases for larger mobile ions, a more flexible backbone, and/or a smaller degree of intrachain dynamic cooperativity, as a consequence of the increasing importance of collective elasticity in determining the activation barrier. An analysis of transient segmental localization and a crossover to activated dynamics reveals similar, but not identical, trends as predicted at the laboratory vitrification timescale.

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

Hyper-fidelity depletion coupled with discrete pebble motion in pebble bed reactors

Pebble bed reactors have raised new interest during the past decade due to their attractive characteristics. Therefore, accurate simulations must be performed to better understand these systems and ensure optimal and safe designs. Most current methods use lower fidelity approaches with representative unit-cells or macro-zones with uniform fluxes, which have accuracy and flexibility limitations. A novel hyper-fidelity method for pebble bed reactors depletion is presented and internally couples Serpent 2 and a pseudo-motion routine. Pseudo-motion is applied handling vertical shifts of compositions in a static pebble bed, random reinsertion of used pebbles, insertion of fresh pebbles and used pebbles discarding. Associated with individual depletion to correctly determine the flux spectrum and composition in each pebble, this hyper-fidelity approach paves the way towards more accurate depletion calculation in pebble bed reactors. Using this method, a demonstration is completed on a small-scale reactor. In this application, the core reaches equilibrium, and the following data is extracted: core-wise parameters evolution, pebble-wise spatial and statistical distribution. Discarded pebbles are analyzed, and relevant information is shown. This work proves the feasibility of hyper fidelity depletion with Serpent 2, and the range of use for this method: reactor design and analysis for equilibrium and slow transients, lower fidelity methods validation and feeding fuel performance, thermal-hydraulics, or waste management models. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2025

Data-driven modeling of dynamic occupant thermostat override behavior for demand response applications

Buildings consume nearly 40% of global energy and produce similar emissions. Whiletechnological advances address efficiency, occupant behavior causes energy use variations up to 300% between identical buildings. This gap between predicted and actual building performance impacts building design, operations, and grid demand management programs. Through analyses of smart thermostat data from 1,400 single-occupant homes, the researchdemonstrates that occupants respond to 8°F thermostat setpoint changes within a median of 15 minutes, while 2°F changes trigger responses within a median of 30 minutes. This highlights an understudied temporal relationship between thermostat setbacks and response time of occupant behaviors. Models of such behavior dynamics are required to incorporate occupant impacts into building performance simulation. A key contribution of this dissertation is the Thermal Frustration Theory (TFT), which positsthat thermal discomfort driven behaviors are caused by the time-accumulation of discomfort, not simply a temperature deviation threshold or a delay from an initiating event. Using a dataset of 634 thermostats, each with 25+ manual setpoint changes, a comparative analysis of TFT and comfort zone and a delayed response theories demonstrated that personalized TFT models better predict when manual setpoint change occur. This was measured by the area under the curve statistical measure (AUC); all three models perform similarly by a Matthews Correlation Coefficient measure. Higher AUC performance is especially important for modeling occupant behavior in demand response programs where false negatives of rare occupant interactions could adversely affect grid stability. EnergyPlus based simulations were conducted with TFT-derived occupant models, demonstrating the ability to identify parameters of known TFT models from only data observable with smart thermostats, even under the presence of noise from routine overrides. Overall, the dissertation highlights that thermostat interactions are neither static,instantaneous, nor driven solely by the environment. Instead, temporal accumulation of discomfort and routine-based behavior play important roles. The methodology and results offer a pathway towards more accurate modeling of human-building interactions for policy assessment, building design, and demand response programs.

Sharma, Kunind [Northeastern University] (ORCID:00↗
DOE OSTI2022

Air classification of forest residue for tissue and ash separation efficiency

The goal of this Case Study was to evaluate the performance of air classification of logging residues toward meeting conversion CMAs for carbon and ash contents, as compared to the static status quo Base Case system in which the residues are first dried and then ground in a hammer mill with a 6 mm screen and fines < 1.18 mm are removed. Also considered were moisture and ash impacts on throughput and Overall Operating Effectiveness (OOE), as well as delivered feedstock cost and minimum fuel selling price (MFSP). Laboratory data on the impacts of fan speed and moisture content on the separation efficiency of soil ash, needles and bark from white wood were received from FCIC Subtask 5.2. Average throughput and energy consumption data were obtained from the Bioenergy Feedstock National User Facility (BFNUF) for the same air classifier. These data were utilized to develop the necessary response surface equations to perform throughput analysis using discrete event simulation. Because the Base Case status quo system utilizes drying prior to grinding, we modeled the Case Study with drying prior to air classification and subsequent grinding of the separated white wood to isolate the individual quality and cost impacts of air classification relative to the Base Case system. Key takeaways from this Case Study are that air classification improves the quality of the final material, but increases the production cost, especially when lights are disposed; this becomes a tradeoff between increased conversion yield and the extra cost. Removing material should be done as early in the process as possible. Each operation that occurs prior to removing the material increases the cost of the disposed material and leads to wasted energy expenditures. As processes are included or modified, the impact on monetary and energy cost should be included in the decision process. Identifying alternative uses and the associated value for material that is separated from the feedstock stream will have a significant impact on the delivered cost of the material. Although there was an increase in production cost by adding air classification, there was a resulting benefit to the MFSP when meeting or exceeding all of the conversion CMAs; this was an important finding that should be explored further with FCIC Subtask 8.3 and potentially FCIC Task 6.

09 BIOMASS FUELS↗
DOE OSTI2021

Internal Microstructure Dictates Interactions of Polymer-grafted Nanoparticles in Solution

Understanding the effects of polymer brush architecture on particle interactions in solution is requisite to enable the development of functional materials based on self-assembled polymer-grafted nanoparticles (GNPs). Static and dynamic light scattering of polystyrene-grafted silica particle solutions in toluene reveals that the pair interaction potential, inferred from the second virial coefficient, A 2 , is strongly affected by the grafting density, σ, and degree of polymerization, N, of tethered chains. In the limit of intermediate σ (~0.3 to 0.6 nm -2 ) and high N, A 2 is positive and increases with N. This confirms the good solvent conditions and can be qualitatively rationalized on the basis of a pair interaction potential derived for grafted (brush) particles. In contrast, for high σ > 0.6 nm -2 and low N, A 2 displays an unexpected reversal to negative values, thus indicating poor solvent conditions. These findings are rationalized by means of a simple analysis based on a coarse-grained brush potential, which balances the attractive core-core interactions and the excluded volume interactions imparted by the polymer grafts. The results suggest that the steric crowding of polymer ligands in dense GNP systems may fundamentally alter the interactions between brush particles in solution and highlight the crucial role of architecture (internal microstructure) on the behavior of hybrid materials. The effect of grafting density also illustrates the opportunity to tailor the physical properties of hybrid materials by altering geometry (or architecture) rather than a variation of the chemical composition.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Greenhouse gas emissions from the global transportation of crude oil: Current status and mitigation potential

Global crude-oil transportation contributes a significant portion of greenhouse gas (GHG) emissions in the marine transportation sector. In this work, we first compile a detailed country-level global crude-oil transportation network in 2018 and estimate that the direct and well-to-hull GHG emissions related to crude transportation were 97 and 109 million metric tons, respectively. Combining with the country-specific crude recovery GHG intensities, the consumption-based well-to-country-gate crude-oil GHG intensities are derived for individual countries, ranging from 2.99 to 27.32 g CO 2 eq/MJ, with a global crude-volume-weighted average of 8.67 g CO 2 eq/MJ. We then project the global crude transportation emissions at the regional level in 2050 under a static (no change) scenario (based on current ship energy efficiency) and a sustainable-development (SD) scenario (based on the International Energy Agency's projections of ship energy efficiency and penetration of alternative marine fuels). Results show that the global well-to-hull GHG emissions related to crude transportation would be 82 and 59 million metric tons in 2050 in the static and SD scenarios, respectively. To further evaluate the impact of potential fuel-switching on decarbonizing the crude oil transportation sector, we estimate the GHG emissions for 20 fuel/production options in 2050. We find that, in comparison to the static scenario, ~50% reduction in global well-to-hull GHG emissions from crude transportation could be achieved under the SD scenario if green ammonia further replaces conventional ammonia. The methodology developed here can be applied to other commodities to estimate the emissions associated with their global marine transportation and to evaluate the potential emission mitigation options.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Mesoscopic structure of mixed type domain walls in multiaxial ferroelectrics

The structure of a 180° uncharged rotational domain wall in a multiaxial ferroelectric film is studied in the framework of an analytical Landau-Ginzburg-Devonshire (LGD) approach. Finite element modeling (FEM) is used to solve numerically the system of the coupled nonlinear Euler-Lagrange (EL) second-order differential equations for two components of polarization. We show that the structure of the domain wall and corresponding metastable or stable phase of the film are controlled by a single parameter—the dimensionless ferroelectric anisotropy factor μ. We fitted the static profile of a solitary domain wall, calculated by FEM, with kinklike functions for polarization components, and extracted the five μ-dependent parameters from the fitting to FEM curves. The surprisingly high accuracy of the fitting results for two polarization components in the entire μ range allows us to conclude that the analytical functions, which are trial functions in the direct variational method, can be treated as a high-accuracy variational solution of the static EL equations. We further derive exact two-component analytical solutions of the static EL equations for a polydomain 180° domain structure in a multiaxial ferroelectric film. Using these, we derive analytical expressions for the system free energy and analyze its dependence on the film thickness and boundary conditions at the film surfaces. The single-domain state is ground for zero polarization derivative at the surfaces, while the polydomain states minimize the system energy for zero polarization at the surfaces. Counterintuitively, the energy of the polydomain states split into two levels, “0” and “1,” for zero polarization at the surfaces, and each of the levels contains a large number of close-energy sublevels, whose structure is characterized by a different number and type of domain walls. The analytical solutions can become a useful tool for Bayesian analysis of high-resolution scanning transmission electron microscopy images in ferroelectric films.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Challenges and Solutions for Real-Time Phasor Modeling of Large-scale Distribution Network with High PV Penetration

The conversion process of a practical large-scale feeder data from a quasi-static time series (OpenDSS) model into a real-time phasor model (ePHASORSIM in Opal-RT) is discussed in the paper. The process is implemented using an open source Python software. Previous reported implementations for the conversion process lead to several errors when applied to a larger-scale system such as the one considered here. Hence in this work, we describe the common issues in this conversion and propose a customized solution to enhance the efficiency of the conversion and reduce the complexity in the process. A quantitative validation of the enhanced conversion process is presented in this work using an actual high PV penetration feeder model that consists of 2230 buses, and using actual load and PV profile data. After a detailed analysis, this customized conversion software will be made available as an open source tool and is expected to be helpful for researchers who want to pursue a similar conversion. Solutions to various observed issues such as identifying the lines due to islanded network, representation of full impedance model of transformer/lines as sequential models, complexity in the representation of single phase buses/lines as three phase buses/lines to make it compatible with the simulator platform are discussed. Comparison of power flow and time series simulation results obtained from both OpenDSS and ePHASORsim models show very low errors, validating the accuracy of the proposed conversion process.

14 SOLAR ENERGY↗
DOE OSTI2026

Protein Adhesion on Semi-Fluorinated Polystyrene Surfaces in Static and Dynamic Measurements

Reducing protein adhesion is a critical strategy in fouling-resistant material innovation, with broad applications spanning biomedical and healthcare devices, biosensors, industrial and environmental systems, and other important technological domains. Here, in this study, we elucidated protein adhesion behavior on polystyrene-based thin films by neutron reflectometry (NR) and quartz crystal microbalance with dissipation (QCM-D), using both lysozyme and bovine serum albumin (BSA) as model proteins. To this end, semifluorinated polystyrene thin films with gradient wettability and surface energy were fabricated through dry processing using plasma oxidation and gas-phase deposition. Although it is believed that a fully fluorinated alkyl chain offers extremely low surface energy, thus rejecting foulants, and has been used in many fouling-resistant surface designs, enhanced protein–surface interactions were observed consistently in NR and QCM-D results, due to the combined effects of surface morphology and chemistry. On the contrary, depositing shorter fluorinated silane onto a hydrophilic PS surface contributed to a more homogeneous nanoscale fluorine coating, resulting in less initial protein adsorption and improved surface recovery. Comparative analysis of proteins with different sizes on the nanopatterned semifluorinated surface revealed the influence of molecular characteristics on surface interactions. Lysozyme, being smaller and more compact, showed faster adsorption kinetics and higher surface coverage but largely reversible binding, whereas BSA, with its larger and more flexible structure, formed broader and more stable interfacial layers. This study fills the gap in understanding protein adhesion within the range of hydrophobicity (water contact angle ∼90°), as current strategies often associate with extreme hydrophilic and superhydrophobic surfaces due to hydration or low-surface-energy rejection mechanisms, respectively. It also provides in-depth insights into current combinatorial fouling-resistant surface design.

Yuan, Yue [Oak Ridge National Laboratory (ORNL), O↗
DOE OSTI2020

Contact Angle Measurements Using Sessile Drop and Micro-CT Data from Six Sandstones

Numerous sessile drop and micro-computed tomography (micro-CT) studies have been conducted to quantify Geologic Carbon Storage formation wettability by measuring static contact angles (θ); however, the influence of pore geometry remains unknown. In this work, six sandstones (Bandera Brown, Berea, Bentheimer, Mt. Simon, Navajo, and Nugget) are used to measure θ using the two aforementioned experimental methods at identical testing conditions (45°C and 12.41 MPa). The range of θ measured at in situ conditions (micro-CT) exceeds those the range at ex situ (sessile drop method) conditions for all sandstones. However, when droplets with more representative in situ diameters are analyzed, θ averages show ex situ θ exceed those of in situ θ. Pore geometry does influence local θ, but the size of ex situ droplets relative to pore size appears to influence θ. This is an important to consider for future sessile drop studies used for analysis of CO 2 behavior in carbon storage reservoirs.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Backbone Stitching in Bottlebrush Copolymer Mesodomains and the Impact of Side Chain Crystallization

We synthesized bottlebrush statistical copolymers (BSCPs) having poly­(ethylene oxide) (PEO) and poly­(dimethylsiloxane) (PDMS) side chains attached to a polynorbornene backbone. Small-angle X-ray scattering analysis showed that for densely grafted BSCPs, the scattering length density gradually transitions between the PEO and PDMS domains. For loosely grafted BSCPs, the polymer backbone formed a distinct mesodomain, with a lower electron and mass density than both the PEO and PDMS domains. The bottlebrush backbone essentially “stitches” the PEO and PDMS side chains, looping back and forth from the PEO to PDMS domains with the backbone segments oriented normal to the domain interfaces. Self-consistent field theory (SCFT) calculations validated the stitching of the backbone driven by the microphase separation of PEO and PDMS, along with a strong segmental order of the side chains in the melt. The reduced birefringence upon PEO crystallization suggests the disruption of the strong segmental order by the crystallization. Both the static intrinsic and the form birefringences of the BSCPs decreased upon PEO crystallization. Solid-state NMR confirmed the rigidity of PEO crystallites and the bottlebrush backbone. Self-assembly of BSCPs containing polyhedral oligomeric silsesquioxane (POSS) pendent groups was also evaluated by X-ray scattering, showing the formation of lamellar microdomains that inhibited POSS crystallization.

Hu, Mingqiu↗
DOE OSTI2021

Specular-reflection contributions to static and dynamic radiation forces on circular cylinders

Interest in the response of highly reflecting objects in water to modulated acoustical radiation forces makes it appropriate to consider contributions to such forces from perfectly reflecting objects to provide insight into radiation forces. The acoustic illumination can have wavelengths much smaller than the object's size, and objects of interest may have complicated shapes. Here, the specular contribution to the oscillating radiation force on an infinite circular cylinder at normal incidence is considered for double-sideband-suppressed carrier-modulated acoustic illumination. The oscillatory magnitude of the specular force decreases monotonically with increasing modulation frequency, and the phase of the oscillating force depends on the relative phase of the sidebands. The phase dependence on the modulation frequency can be reduced with the appropriate selection of a sideband relative-phase parameter. That is a consequence of the significance of rays that are incident on the cylinder having small impact parameters that are nearly backscattered. For one choice of a relative sideband phase, a prior partial wave series (PWS) solution is available, which supports the specular analysis when the PWS is evaluated for a rigid cylinder. The importance of specular contributions for aluminum cylinders in water is noted. A specular analysis for an analogous spherical reflector is also summarized.

Marston, Philip L. (ORCID:0000000306778537)↗
DOE OSTI2021

Implementation and Verification of Cable Bending Stiffness in MoorDyn

The relatively large motions experienced by floating wind turbines and wave energy converters pose a challenge for power cables, whose internal components provide significant bending resistance and are sensitive to deformation. The behavior and associated design considerations of power cables in these highly dynamic applications make coupled analysis relevant for design. Bending stiffness capabilities have recently been added to the lumped-mass mooring dynamics model MoorDyn to enable simulation of dynamic power cables. MoorDyn is a common modeling choice for floating wind energy simulation (often coupled with OpenFAST) and floating wave energy converter simulation (often coupled with WEC-Sim) but the model’s previous line elasticity formulation only considered axial stiffness. To properly capture the dynamics of power cables, a bending stiffness model has been added that approximates cable curvature based on the difference in tangent vectors of adjacent elements. The resulting bending moment is realized by applying forces on adjacent nodes, enabling cable modeling while leaving the underlying lumped-mass formulation unchanged. In this paper, the new bending stiffness implementation is verified in static conditions against analytical solutions and then in a dynamic power cable scenario in comparison with the commercial simulator OrcaFlex. The dynamic scenario uses prescribed motions and includes wave loadings on the cable. Results indicate correct implementation of bending stiffness and show close agreement with OrcaFlex.

50 EE - Wind and Water Power Program - Water (EE-4↗
DOE OSTI2024

Machine Learning for Predictive Performance Analysis in Charged Particle Beam Tools

Imaging methods driven by probes, electrons, and ions have played a dominant role in modern science and engineering. Opportunities for machine vision and AI that focus on consumer problems like driving and feature recognition, are now presenting themselves for automating aspects of the scientific processes. This proposal aims to enable and drive discovery in ultra-low energy implantation by taking advantage of faster processing, flexible control and detection methods, and architecture-agnostic workflows that will result in higher efficiency and shorter scientific development cycles. Custom microscope control, collection and analysis hardware will provide a framework for conducting novel in situ experiments revealing unprecedented insight into surface dynamics at the nanoscale. Ion implantation is a key capability for the semiconductor industry. As devices shrink, novel materials enter the manufacturing line, and quantum technologies transition to being more mainstream. Traditional implantation methods fall short in terms of energy, ion species, and positional precision. Here we demonstrate 1 keV focused ion beam Au implantation into Si and validate the results via atom probe tomography. We show the Au implant depth at 1 keV is 0.8 nm and that identical results for low energy ion implants can be achieved by either lowering the column voltage, or decelerating ions using bias – while maintaining a sub-micron beam focus. We compare our experimental results to static calculations using SRIM and dynamic calculations using binary collision approximation codes TRIDYN and IMSIL. A large discrepancy between the static and dynamic simulation is found that is due to lattice enrichment with high stopping power Au and surface sputtering. Additionally, we demonstrate how model details are particularly important to the simulation of these low-energy heavy-ion implantations. Finally, we discuss how our results pave a way to much lower implantation energies, while maintaining high spatial resolution.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Photoexcited energy relaxation and vibronic couplings in π-conjugated carbon nanorings

Conjugated carbon nanorings exhibit unique photophysical properties that, combined with their tunable sizes and conformations, make them suitable for a variety of practical applications. These properties are intimately associated to their strained, bent and sterically hindered cyclic structures. Herein we perform a comparative analysis of the photoinduced dynamics in carbon nanorings composed of nine phenyl units([9]CPP) and nine naphthyl units ([9]CN) respectively. The sterically demanding naphthyl units lead to large dihedral angles between neighboring units. Nevertheless, the ultrafast electronic and vibrational energy relaxation and redistribution is found to be similar for both systems. We observe that vibronic couplings, introduced by nonadiabatic energy transfer between electronic excited states, ensure the intramolecular vibrational energy redistribution through specific vibrational modes. The comparative impact of the internal conversion process on the exciton spatial localization and intra-ring migration indicates that naphthyl units in [9]CN achieve more efficient but less dynamical self-trapping compared to that of phenyl units in [9]CPP. That is, during the photoinduced process, the exciton in [9]CN is more static and localized than the exciton in [9]CPP. Furthermore, the internal conversion processes take place through a specific set of middle- to high-frequency normal modes, which directly influence the spatial exciton redistribution during the internal conversion, self-trapping and intra-ring migration.

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

Novel experimental probes of QCD in SIDIS and e + e - annihilation (Final Technical Report)

The research addressed with this award seeks to advance our understanding of the structure and dynamics underlying the properties of visible matter. In our current understanding the nucleons (protons and neutrons) are not fundamental but are comprised of quarks and gluons, which are collectively called partons. A static quark picture fails to explain the properties of the nucleons, such as their mass and intrinsic spin, which are thought to emerge dynamically from the quark-gluon interactions via the strong force. In this work, novel observables employing correlations of particles produced in the scattering of high energy electrons off protons at the CLAS12 experiment at Jefferson Lab were analyzed to probe quark-gluon interactions. The ultimate goal of this line of inquiry is to be able to describe the properties of protons and neutrons from first principles, similar to how studying the hydrogen atom has led to the formulation of the theory of Quantum Electrodynamics. Because quarks cannot be observed directly but only as part of more complex composite particles, a smaller, complimentary part of this work was the analysis of particle production in electron-positron annihilation at the Belle II experiment to understand the production of particles from initial quarks. This takes advantage of the fact that in e + e - annihilation the initial quark dynamics is known, unlike in the scattering of nucleons. Several new applications using Machine Learning algorithms for event tagging and reconstruction to support this program were developed as part of this award. In addition, we had a significant role in the development of the physics program for the future Electron-Ion Collider, which is a new collider to be build in the US within the next decade.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Quantum-classical hybrid algorithm for the simulation of all-electron correlation

While chemical systems containing hundreds to thousands of electrons remain beyond the reach of quantum devices, hybrid quantum-classical algorithms present a promising pathway toward a quantum advantage. Hybrid algorithms treat the exponentially scaling part of the calculation-the static correlation-on the quantum computer and the non-exponentially scaling part-the dynamic correlation-on the classical computer. While a variety of algorithms have been proposed, the dependence of many methods on the total wave function limits the development of easy-to-use classical post-processing implementations. Here, we present a novel combination of quantum and classical algorithms, which computes the all-electron energy of a strongly correlated molecular system on the classical computer from the 2-electron reduced density matrix (2-RDM) evaluated on the quantum device. Significantly, we circumvent the wave function in the all-electron calculations by using density matrix methods that only require input of the statically correlated 2-RDM. Although the algorithm is completely general, we test it with two classical density matrix methods, the anti-Hermitian contracted Schrödinger equation (ACSE) and multiconfiguration pair-density functional theories, using the recently developed quantum ACSE method for simulating the statically correlated 2-RDM. Furthermore, we obtain experimental accuracy for the relative energies of all three benzyne isomers and thereby demonstrate the ability of the developed algorithm to achieve chemically relevant and accurate results on noisy intermediate-scale quantum devices.

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

Updraft dynamics and microphysics: on the added value of the cumulus thermal reference frame in simulations of aerosol–deep convection interactions

One fundamental question about atmospheric moist convection processes that remains debated is whether, or under which conditions, a relevant variability in background aerosol concentrations may have a significant dynamical impact on convective clouds and their associated precipitation. Furthermore, current climate models must parameterize both the microphysical and the cumulus convection processes, but this is usually implemented separately, whereas in nature there is a strong coupling between them. As a first step to improve our understanding of these two problems, we investigate how aerosol concentrations modify key properties of updrafts in eight large-eddy-permitting regional simulations of a case study of scattered convection over Houston, Texas, in which convection is explicitly simulated and microphysical processes are parameterized. Dynamical and liquid-phase microphysical responses are investigated using the following two different reference frames: static cloudy updraft grid cells versus tracked cumulus thermals. In both frameworks, we observe the expected microphysical responses to higher aerosol concentrations, such as higher cloud number concentrations and lower rain number concentrations. In terms of the dynamical responses, both frameworks indicate weak impacts of varying aerosol concentrations relative to the noise between simulations over the observationally derived range of aerosol variability for this case study. On the other hand, results suggest that thermals are more selective than cloudy updraft grid cells in terms of sampling the most active convective air masses. For instance, vertical velocity from thermals is significantly higher at upper levels than when sampled from cloudy updraft grid points, and several microphysical variables have higher average values in the cumulus thermal framework than in the cloudy updraft framework. In addition, the thermal analysis is seen to add rich quantitative information about the rates and covariability of microphysical processes spatially and throughout tracked thermal lifecycles, which can serve as a stronger foundation for improving subgrid-scale parameterizations.

54 ENVIRONMENTAL SCIENCES↗