Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “scale coupling”

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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Simulating Hydrodynamics in Cosmology with CRK-HACC

Abstract We introduce CRK-HACC, an extension of the Hardware/Hybrid Accelerated Cosmology Code (HACC), to resolve gas hydrodynamics in large-scale structure formation simulations of the universe. The new framework couples the HACC gravitational N -body solver with a modern smoothed-particle hydrodynamics (SPH) approach called conservative reproducing kernel SPH (CRKSPH). CRKSPH utilizes smoothing functions that exactly interpolate linear fields while manifestly preserving conservation laws (momentum, mass, and energy). The CRKSPH method has been incorporated to accurately model baryonic effects in cosmology simulations—an important addition targeting the generation of precise synthetic sky predictions for upcoming observational surveys. CRK-HACC inherits the codesign strategies of the HACC solver and is built to run on modern GPU-accelerated supercomputers. In this work, we summarize the primary solver components and present a number of standard validation tests to demonstrate code accuracy, including idealized hydrodynamic and cosmological setups, as well as self-similarity measurements.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Breaking the Brønsted–Evans–Polanyi Relation with Dual-Metal Sites

Linear scaling relationships impose inherent limitations on catalyst activity; the Brønsted−Evans−Polanyi (BEP) relation, which correlates activation and reaction energies, is a prominent example. Here we report a dual-metal site catalyst (DMSC) on ceria that breaks the BEP relation for C−C coupling of methyl intermediates an elementary step in methane coupling to form ethane. The DMSC structure on CeO 2 (111) was discovered by density-functional theory (DFT) structural exploration and confirmed to be stable via ab initio thermodynamics and ab initio molecular dynamics. Homonuclear and heteronuclear DMSCs of Ni, Pd, Pt, Fe, Ru, Os, Co, Rh, and Ir (45 pairs in total) were examined for methyl affinity and methyl−methyl coupling activation energy. We found that many heteronuclear DMSCs break the BEP linear scaling due to a mixed low-affinity/high-affinity coadsorption of the two methyl groups, decoupling the step responsible for the activation energy (E a ) at the low-affinity site from the overall reaction energy (ΔE) determined by both sites. This mechanism of breaking the BEP relationship via the DMSCs offers a catalyst design principle for C−C coupling reactions.

Chen, Yiming [Vanderbilt Univ., Nashville, TN (Uni↗
DOE OSTI2021

A mesoscopic digital twin that bridges length and time scales for control of additively manufactured metal microstructures

We present our recent development of an integrated mesoscale digital twin (DT) framework for relating processing conditions, microstructures, and mechanical responses of additively manufactured (AM) metals. In particular, focusing on the laser powder bed fusion technique, we describe how individual modeling and simulation capabilities are coupled to investigate and control AM microstructural features at multiple length and time scales. We review our prior case studies that demonstrate the integrated modeling schemes, in which high-fidelity melt pool dynamics simulations provide accurate local thermal profiles and histories to subsequent AM microstructure simulations. We also report our new mechanical response modeling results for predicted AM microstructures. In addition, we illustrate how our DT framework has been validated through modeling–experiment integration, as well as how it has been practically utilized to guide and analyze AM experiments. Finally, we share our perspectives on future directions of further development of the DT framework for more efficient, accurate predictions and wider ranges of applications.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Localized basal area affects soil respiration temperature sensitivity in a coastal deciduous forest

Abstract. Soil respiration (Rs), the flow of CO2 from the soil surface to the atmosphere, is one of the largest carbon fluxes in the terrestrial biosphere. The spatial variability of Rs is both large and poorly understood, limiting our ability to robustly scale it in space. One factor in Rs spatial variability is the autotrophic contribution from plant roots, but it is uncertain how the presence of plants affects the magnitude and temperature sensitivity of Rs. This study used 1 year of Rs measurements to examine the effect of localized basal area on Rs in the growing and dormant seasons, as well as during moisture-limited times, in a temperate, coastal, deciduous forest in eastern Maryland, USA. In a linear mixed-effects model, tree basal area within a 5 m radius (BA5) exerted a significant positive effect on the temperature sensitivity of soil respiration. Soil moisture was the dominant control on Rs during the dry portions of the year, while soil moisture, temperature, and BA5 all exerted significant effects on Rs in wetter periods. Our results suggest that autotrophic respiration is more sensitive to temperature than heterotrophic respiration at these sites, although we did not measure these source fluxes directly, and that soil respiration is highly moisture sensitive, even in a record-rainfall year. The Rs flux magnitudes (0.46–15.0 µmol m-2 s-1) and variability (coefficient of variability 10 %–23 % across plots) observed in this study were comparable to values observed in similar forests. Six Rs observations would be required in order to estimate the mean across all study sites to within 50 %, and 518 would be required in order to estimate it to within 5 %, with 95 % confidence. A better understanding of the spatial interactions between plants and microbes, as well as the strength and speed of above- and belowground coupling, is necessary to link these processes with large-scale soil-to-atmosphere C fluxes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Algebraic discrete quantum harmonic oscillator with dynamic resolution scaling

We develop an algebraic formulation for the discrete quantum harmonic oscillator (DQHO) from the Hamiltonian for two, coupled QHOs and provide a physical picture for the Kravchuk function eigenstates of the oscillator. The familiar $\mathfrak{su}(2)$ structure of the coupled QHO Hamiltonian divides its spectrum into sets corresponding to the DQHO at different resolutions. In addition to energy ladder operators, the formulation allows for the introduction of resolution ladder operators connecting all DQHOs with different resolutions, thus enabling the dynamic scaling of the resolution of finite degree-of-freedom quantum simulations. The coherent state of the DQHO is constructed, and its expected position is proven to oscillate as a classical harmonic oscillator. The DQHO coherent state recovers that of the quantum harmonic oscillator at large resolution.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

High temporal frequency data from a four turbine, blade-resolved wind farm simulation with ExaWind

The data was generated with ExaWind (https://github.com/Exawind) which couples AMR-Wind (https://github.com/Exawind/amr-wind/), Nalu-Wind (https://github.com/Exawind/nalu-wind), TIOGA (https://github.com/Exawind/tioga), and OpenFAST (https://github.com/OpenFAST/openfast). This is a large-scale simulation of a blade-resolved wind farm using the ExaWind software stack. ExaWind couples together a background flow solver, AMR-Wind, and a near-body solver, Nalu-Wind, through an overset technique from the TIOGA application. Another application, OpenFAST, handles the structural dynamics of the turbine blades and towers, which informs the fluid-structure interaction of the wind turbines with the flow solvers. This particular simulation includes four blade-resolved wind turbines operating in a turbulent atmospheric boundary layer. The AMR-Wind solver uses 500 million cells and is being solved on 256 AMD GPUs of the Oakridge Leadership Computing Facility Frontier supercomputer. Each turbine is assigned its own Nalu-Wind solver with over 13 million elements per turbine and solved using 448 CPU cores, for a total of 1792 CPU cores. For each node, 56 cores contain Nalu-Wind, while 8 cores correspond to AMR-Wind operations on the GPUs. Consequently, ExaWind is entirely utilizing the CPUs and the GPUs of the nodes concurrently. The data used in the visualization is full flow field data output from the simulation. It is lossy-compressed to a specific accuracy using ZFP and written to disk every 16 time-steps to enable real-time flow visualization. The flow fields are sampled at a high temporal frequency to enable real-time, 24fps visualization. The flow fields are sampled every 12 simulation time steps (every 0.04132s).

17 WIND ENERGY↗
DOE OSTI2022

Assessing physics of ion temperature gradient turbulence via hierarchical reduced-model representations

In this work, the saturation physics of ion temperature gradient (ITG) turbulence is probed by studying how amplitudes and scalings with key parameters vary across a hierarchy of reduced models. The models derive from nonlinear fluid equations for toroidal ITG turbulence under approximations to the mode coupling interactions in wavenumber space and the representation of turbulent decorrelation. Mode coupling approximations include local-in-wavenumber treatments like the spectral density of flux in quasilinear theory, a truncation to three nonlinearly interacting waves, and the interactions in a cascade to high radial wavenumber mediated by a single zonal flow. Turbulent decorrelation treatments are based on the triplet correlation time with and without eddy damping. Model fidelity is assessed by the scalings and magnitudes of the squared amplitudes of unstable mode, stable mode, and zonal flow with respect to the flow-damping rate and temperature gradient. It is shown that all models reproduce fundamental scalings, provided they incorporate the coupling of unstable mode, stable mode, and zonal flow. Accurate amplitude prediction requires eddy damping in the triplet correlation time and proper representation of the zonal-flow drive by interactions associated with the radial wavenumber cascade.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Fragility of the nematic spin liquid induced by diagonal couplings in the square-lattice SU(3) model

Here, we present a large-scale density matrix renormalization group (DMRG) study of the spin-1 SU(3) bilinear-biquadratic model on the square lattice, which was suggested to host a nematic spin liquid state in recent DMRG calculations. We report that this spin liquid appears to strongly compete with a three-sublattice magnetic order. To further study the competition between the two states, and the reason of the emergent nematic spin liquid, we included an additional next-nearest-neighbor SU(3) symmetric interactions along one of the two plaquette diagonal directions. This allows to tune the square-lattice model to the triangular-lattice model. By computing spin correlation functions and various order parameters, we find that the three-sublattice order may develop at infinitesimal additional new coupling, at least within the precision of our study. Compared with the previous findings that the nematic spin liquid is stable in extended parameter regions with additional couplings that respect the lattice symmetries of the square lattice, we argue that here the diagonal couplings, which frustrate the bipartite-lattice structure, rapidly suppress the two-sublattice fluctuations and the three-sublattice order thus wins. This numerical result is consistent with the conjecture that the nematic spin liquid emerges from the competition between two- and three-sublattice fluctuations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Design and Operation of a Multi-Bed Catalytic Micro-Reactor for the Study of Co-Processing of Bio-Oils with VGO

An industry wide shift from fossil-based fuel to renewable fuel sources including biomass, municipal waste, and plastics will require new process monitoring methods to minimize transitional risks including off specification product formation and catalyst deactivation. This project aims to provide a machine learning based process monitoring tool composed of online, slipstream mass spectra for use in biomass refineries and co-processing in existing refineries allowing operators to monitor product qualities and adjust process conditions accordingly. In order to maximize the robustness of the tool, large volumes of data must be collected to fine tune model parameters which consists of both micro and pilot scale mass spectral data. Micro-scale data is collected with a multi-tube micro-reactor housing up to six catalysts in horizontal beds, coupled with a molecular beam mass spectrometer. A pyrolizer equipped with an auto-sampler streamlines the micro-scale data collection process. This type of pyrolizer/micro-reactor configuration does not exist on the market, and therefore had to be created for the purposes of this project. The design and commissioning of this reactor will be presented in detail. This reactor set-up is highly flexible and increases throughput of analysis. For catalyst testing, each bed can be individually selected simply by turning valves. For catalyst reduction and regeneration, simultaneous flow through all six beds is used. The reproducibility of the system was first assessed with whole biomass pyrolysis along with pyrolysis of calibration standards. Initial work on this system evaluated two FCC catalysts, equilibrium catalyst (E-cat), and a proprietary catalyst from Johnson Matthey specifically design for co-processing of bio-oil with vacuum gas oil (VGO). This work used model compounds and VGO which illuminated differences in products produced by the catalysts.

biomass↗
DOE OSTI2021

Quarkonium semiclassical transport in quark-gluon plasma: factorization and quantum correction

We study quarkonium transport in the quark-gluon plasma by using the potential nonrelativistic QCD (pNRQCD) effective field theory and the framework of open quantum systems. We argue that the coupling between quarkonium and the thermal bath is weak using separation of scales, so the initial density matrix of the total system factorizes and the time evolution of the subsystem is Markovian. We derive the semiclassical Boltzmann equation for quarkonium by applying a Wigner transform to the Lindblad equation and carrying out a semiclassical expansion. We resum relevant interactions to all orders in the coupling constant at leading power of the nonrelativistic and multipole expansions. The derivation is valid for both weakly coupled and strongly coupled quark-gluon plasmas. We find reaction rates in the transport equation factorize into a quarkonium dipole transition function and a chromoelectric gluon distribution function. For the differential reaction rate, the definition of the momentum dependent chromoelectric gluon distribution function involves staple-shaped Wilson lines. For the inclusive reaction rate, the Wilson lines collapse into a straight line along the real time axis and the distribution becomes momentum independent. The relation between the two Wilson lines is analogous to the relation between the Wilson lines appearing in the gluon parton distribution function (PDF) and the gluon transverse momentum dependent parton distribution function (TMDPDF). The centrality dependence of the quarkonium nuclear modification factor measured by experiments probes the momentum independent distribution while the transverse momentum dependence and measurements of the azimuthal angular anisotropy may be able to probe the momentum dependent one. We discuss one way to indirectly constrain the quarkonium in-medium real potential by using the factorization formula and lattice calculations. The leading quantum correction to the semiclassical transport equation of quarkonium is also worked out. The study can be easily generalized to quarkonium transport in cold nuclear matter, which is relevant for quarkonium production in eA collisions in the future Electron-Ion Collider.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Estimation of stress and stress-induced permeability change in a geological nuclear waste repository in a thermo-hydrologically coupled simulation

Geologic disposal is a promising solution for a safe permanent isolation of accumulated high-level nuclear waste from nuclear powerplants. The behavior of host rock is highly coupled thermally, hydromechanically, and chemically. Numerical simulations of such coupled phenomena for the extremely long term (>100,000 years) and large length scale (>1 km) of geologic disposal remain to be computationally challenging. In this study, a methodology has been developed to approximate stress and stress-induced permeability change in host rock using only thermo-hydrological (TH) variables. A coupled thermo-hydromechanical (THM) simulation is carried out using TOUGH-FLAC simulator to model THM behaviors of a generic nuclear waste repository, in order to evaluate the performance of the developed methodology, which is implemented in a coupled TH simulation. Results show that stress and permeability change estimated by the developed methodology in the TH-coupled simulation match those calculated in the THM-coupled simulation over the simulated timeframe of over 10,000 years. Details about errors in stress and permeability estimates accrued by the developed methodology are discussed here. The developed methodology will help incorporate stress-induced permeability change into existing TH simulators for the long-term radionuclide transport in geologic disposal.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2026

Importance of Spatially Continuous Urban Surface Properties in Urban‐Resolving Earth System Modeling

Accurate representation of urban properties and processes at higher resolutions in global modeling systems is essential for advancing our ability to capture the complexities of urban systems and informing effective resilience strategies. However, the prescription of coarse global-scale urban properties in most state-of-the-art Earth system models (ESMs) is limiting their potential for capturing urban signals as they advance toward kilometer-scale simulation capabilities. To bridge this gap in inadequate urban property representation and to advance urban-resolving Earth system modeling, this work integrates the newly-developed global 1 km-resolution facet-level urban surface property data set, U-Surf, into the land component of Community Earth System Model (CESM)—Community Terrestrial System Model (CTSM). The land-only CTSM simulations are validated against satellite measurements, ground-based urban weather stations, flux tower observations, and reanalysis data. Results demonstrate that the enhanced urban properties allow improved simulations of urban meteorology and surface energy fluxes compared to the default coarse-resolution categorical urban canopy parameters. Spatial scaling analysis reveals regime-dependent information loss during resolution aggregation, as well as substantial scale-dependent variations in urban surface energy flux representation. Furthermore, these findings have critical implications for coupled Earth system modeling when including the effect of land-atmosphere interaction. This work establishes a foundation for future urban-resolving kilometer-scale ESM development, which will enable systematic intra- and inter-city comparisons that inform urban adaptation strategies across diverse global urban environments.

Cheng, Yifan [University of Illinois Urbana-Champa↗
DOE OSTI2023

The ABC of RPV: classification of R-parity violating signatures at the LHC for small couplings

We perform a classification of all potential supersymmetric R-parity violating signatures at the LHC to address the question: are existing bounds on supersymmetric models robust, or are there still signatures not covered by existing searches, allowing LHC scale supersymmetry to be hiding? We analyze all possible scenarios with one dominant RPV trilinear coupling at a time, allowing for arbitrary LSPs and mass spectra. We consider direct production of the LSP, as well as production via gauge-cascades, and find 6 different experimental signatures for the $LL\overline{E}$-case, 6 for the $LQ\overline{D}$-case, and 5 for the $\overline{U}$ $\overline{D}$ $\overline{D}$-case; together these provide complete coverage of the RPV-MSSM landscape. This set of signatures is confronted with the existing searches by ATLAS and CMS. We find all signatures have been covered at the LHC, although not at the sensitivity level needed to probe the direct production of all LSP types. For the case of a dominant $LL\overline{E}$-operator, we use CheckMATE to quantify the current lower bounds on the supersymmetric masses and find the limits to be comparable to or better than the R-parity conserving case. Our treatment can be easily extended to scenarios with more than one non-zero RPV coupling.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Towards a multiscale approach for understanding irradiation induced swelling and creep in 316 stainless steels - A coupled cluster dynamics and crystal plasticity approach

Structural materials undergo mechanical degradation, in part due to irradiation-induced swelling and creep, under nuclear reactor conditions. While swelling results from the migration and clustering of irradiation-induced atomic-scale mobile defects, the interaction of mesoscale dislocations with these defects causes creep deformation. A coupled crystal plasticity (CP) and mean-field cluster dynamics (CD) approach is presented to investigate the effect of irradiation on the long-term mechanical behavior of 316 stainless steel, which are under consideration for use in nuclear reactors. The temporal evolution of Frenkel pairs and extended defect population, under a chosen irradiation flux and temperature, is predicted using the CD model. The impact of the irradiation defects on the stress state, and the resulting dislocation-mediated inelastic deformation, is modeled concurrently with the CP model. The inelastic deformation is irradiation flux dependent, and early-stage defect evolution determines the later-stage mechanical behavior in 316 stainless steel.

36 - MATERIALS SCIENCE↗
DOE OSTI2020

Zonal-Based Emission Source Term Model for Predicting Particulate Emission Factors in Wildfire Simulations

A physics/chemistry-based numerical model for predicting the emission of fine particles from wildfires is proposed. This model implements the fundamental mechanisms of soot formation in a combustion environment: soot nucleation, surface growth, agglomeration, oxidation, and particle fragmentation. These mechanisms occur on a scale too fine for the discretization of most wildfire models, which need to simulate landscape-scale dynamics. As a result this model implements a zonal approach, where the computed soot particle distribution is partitioned into process zones within a single resolved grid cell. These process zones include: an inception zone (for nucleation), a heating zone (for coagulation, surface growth, and fragmentation), a reaction zone (for oxidation), and a quenched zone (for atmospheric processes). Governing mechanisms are applied to the appropriate zones to predict total particle growth and emission. Additionally, the proposed model is implemented into HIGRAD/FIRETEC, a physics-based wildfire simulation code which couples interactions between fire, fuels, atmosphere, and topography on a landscape scale. Fire simulations among grasslands and conifer forests are performed and compared against experimental data for emission factors.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Modeling and Simulation of Fuel Dispersal During the Loss-of-Coolant Accident

This document is the compilation of the milestone portion to a larger end of project NEUP report. The executive summary of the modeling portion is provided below: In the event of cladding rupture during a postulated LOCA in a pressurized water reactor, fuel particles, along with fission gases, can be expelled into the reactor core from the fractured fuel rod, a phenomenon referred to as fuel dispersal. The initial stage of fuel dispersal is strongly influenced by the high-pressure ejection of fuel fragments, the size and geometry of the ruptured cladding, and the depressurization history of the fuel rod during the postulated LOCA transient. Depending on the location of the burst orifice relative to the quench front, the dispersal event represents an intricate three-phase flow and heat transfer phenomenon, where high-temperature fuel particles carried by the fission gases interact with the coolant within the narrow subchannels of the fuel assemblies, inducing localized phase change. Given the unique multiphysics nature of this phenomena, the current study develops a dedicated computational framework to predict the mass distribution and cooling of dispersing fuel particles, facilitating post-accident assessment and management of the fuel assemblies. Considering the scale of nuclear reactor applications, a continuum three-fluid model is proposed for simulating the transport of solids within the reactor core. With high-temperature fuel fragments within the liquid media, nucleation sites inducing phase changes are dispersed within the flow domain. Coupled with the fact that the transient dispersal event occurs on different time scales than other three-phase flow applications, this study derives a time-averaged three-fluid flow model without losing generality. The assumptions regarding the continuum treatment of the solid phase and the modeling of fuel dispersal behavior are incorporated to simplify the governing equations and derive applicable closure relations. The computational validation of the model was conducted using adiabatic experimental results obtained from ongoing research at Oregon State University, focusing on characterizing fuel dispersal behavior during simulated LOCA conditions. Settlement characteristics of the solids, quantified by the probability distribution of equivalent particles, closely matched the probability density functions reported in experimental studies. The transport of fuel particles within a scaled 5 × 5 lattice of a pressurized-water reactor rod bundle geometry was modeled through a two-fluid Eulerian framework. The required boundary conditions were evaluated from the fuel performance code BISON in a postulated large-break LOCA scenario. The modeling framework considered solid fuel particles as granular matter, interacting with the gaseous dry steam phase and fission gases through the governing interfacial momentum exchange between the participating fluids. The simulation results provided the volume fraction of the solids obtained at the bottom surface of the enclosing tank geometry. Postulated LOCA leading to fuel dispersal phenomena involves the strong coupling between fuel thermomechanics, cladding deformation, thermal-hydraulics, and fuel particle transport. Incorporation of such a strong coupling in numerical simulation is performed by coupling the multiphysics solvers. In the case of fuel dispersal, a strong coupled simulation can be performed by coupling the BISON code for fuel performance, the TRACE code for system-level thermal hydraulics, and fuel particle transport in Multiphysics Object-Oriented Simulation Environment (MOOSE). For such intricate infrastructure, the MOOSE Framework eases the data transfer between codes. The recent version of MOOSE has incorporated the Navier-Stokes module for the fluid flow. An exploratory exercise was done to gain familiarity with finite volume capabilities in the MOOSE framework to incorporate the Spalart-Allmaras (SA) turbulence model. New finite-volume and auxiliary kernels were introduced to assemble the SA transport equation, compute turbulent viscosity, and evaluate wall distance and diagnostic turbulence terms, fully integrated with existing Navier-Stokes modules. A turbulent lid-driven cavity at a Reynolds number of approximately 10,000 is used for verification. MOOSE shows the robust solver convergence and produces the turbulent features. But it underpredicts the velocity profile and turbulent quantities, emphasizing the need to develop improved SA near-wall treatments (e.g., low-Re corrections or wall functions) as a key direction for future work.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Coupled wave-current modeling for hydrodynamic load analysis of macroalgae cultivation farms

Large-scale cultivation of macroalgae is one of the most promising biofuel sources that could reduce our consumption of fossil fuels and would be particularly economically competitive when grown for the co-production of additional goods such as food and textiles. Current production of biofuel from algal biomass is limited by labor costs and lack of data on potential impacts of cultivation and harvesting on marine and coastal environments as well as a comprehensive characterization of the impact of environmental conditions on macroalgal feedstock and their bioproduct and biofuel yield and quality. Awarded by the United States (U.S.) Department of Energy Advanced Research Projects Agency-Energy (ARPA-E) Macroalgae Research Inspiring Novel Energy Resources (MARINER) program, which seeks to enable the U.S. as a global leader in the production of marine biomass, the Marine Biological Laboratory (MBL) has been developing a test system for tropical seaweed cultivation in the Gulf of Mexico and the Caribbean. An integral step in designing macroalgae cultivation farms is to adequately characterize the hydrodynamic climate at potential growth sites. In regions like the Gulf of Mexico and the Caribbean, which are highly exposed to tropical cyclones, it is critical to not only consider the day-to-day hydrodynamic conditions but to also assess the risk of extreme sea states to ensure the survivability of future macroalgae farms. Under these considerations, the Pacific Northwest National Laboratory (PNNL) is providing modeling support for MBL to inform the design and siting of the farm systems that have been proposed for their ARPA-E MARINER project. This report describes the development of a high-resolution coupled storm surge and wave model to simulate the hydrodynamics and wave climate at proposed macroalgae cultivation sites selected by MBL in Florida and Puerto Rico. Model results, including model validation, water level, current distributions, and sea states, are discussed for both selected sites in Florida and Puerto Rico coast. These simulations provide an accurate insight into the hydrodynamic conditions that a macroalgae farm is likely to experience during its operational lifetime, including current information to support fine-scale hydrodynamic load modeling, risk analysis and system design.

09 BIOMASS FUELS↗
DOE OSTI2020

Third-Order Møller–Plesset Theory Made More Useful? The Role of Density Functional Theory Orbitals

The practical utility of Møller–Plesset (MP) perturbation theory is severely constrained by the use of Hartree–Fock (HF) orbitals. It has recently been shown that the use of regularized orbital-optimized MP2 orbitals and scaling of MP3 energy could lead to a significant reduction in MP3 error. In this work, we examine whether density functional theory (DFT)-optimized orbitals can be similarly employed to improve the performance of MP theory at both the MP2 and MP3 levels. We find that the use of DFT orbitals leads to significantly improved performance for prediction of thermochemistry, barrier heights, noncovalent interactions, and dipole moments relative to the standard HF-based MP theory. Indeed, MP3 (with or without scaling) with DFT orbitals is found to surpass the accuracy of coupled-cluster singles and doubles (CCSD) for several data sets. We also found that the results are not particularly functional sensitive in most cases (although range-separated hybrid functionals with low delocalization error perform the best). As such, MP3 based on DFT orbitals thus appears to be an efficient, noniterative O(N 6 ) scaling wave-function approach for single-reference electronic structure computations. Scaled MP2 with DFT orbitals is also found to be quite accurate in many cases, although modern double hybrid functionals are likely to be considerably more accurate.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗