Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model scaling”

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 721 records · Page 40

The approximate second order coupled-cluster method based on a size-consistent Brillouin–Wigner partitioning

We present a variant of the approximate second order coupled-cluster method (CC2) with a two-parameter size-consistent Brillouin–Wigner (BW-s) partitioning instead of a Møller–Plesset (MP) partitioning for the unperturbed Hamiltonian, which we refer to as BWs-CC2. The computational complexity of this model scales identically to CC2 with molecular size. Conventional CC2 and its regularized BWs-CC2 variants, as well as conventional MP2 and two of its regularized BW-s2 variants, were assessed on a 535 element database spanning thermochemistry, non-covalent interactions, barrier heights, and isomerization energies. To ensure a well-defined model chemistry, the assessment was performed using internally stable spin-polarized Hartree–Fock (HF) orbitals in the finite aug-cc-pVQZ basis without counterpoise corrections. As a result of using stable orbitals, contrary to conventional wisdom, we find that CC2 substantially outperforms MP2 on molecules with significantly spin contaminated reference orbitals without a significant increase in error on systems with a spin-pure reference, showing the value of its single substitutions. While no single choice of regularization parameters can be optimal for all datasets, we find that BWs-CC2 generally outperforms both CC2 and BW-s2 with a single judicious parameter choice. Additional tests on dipole moments and bond lengths of diatomics provide further support for the utility of this choice. Furthermore, the main outliers and poorest performing cases are associated with large amounts of spin-contamination in the HF reference, which is indicative of systems with either strong correlation or extensive artificial symmetry breaking. Overall, these findings argue that the perception of the quality of the CC2 ground state should be reevaluated and that it can be further improved upon by the soundly based BWs-CC2 variant with the recommended parameter choice.

Correlation energy↗

Efficient wind farm layout optimization with the FLOWERS AEP model and analytic gradients

Wind farm layout optimization (WFLO) studies often aim to maximize the annual energy production (AEP) of a wind farm by choosing an arrangement of turbines that minimizes wake interactions. One way to reduce the cost of WFLO studies is by using more computationally efficient AEP models. The cost of standard AEP modeling approaches, based on the numerical integration of low-fidelity engineering wake models, scales poorly with the number of simulated discrete wind conditions. A second way to reduce cost when using a gradient-based algorithm is to supply exact gradient information instead of finite-difference estimates. However, analytical functions for the derivatives of AEP with respect to turbine positions are not always available in the conventional modeling approach. FLOWERS is a computationally inexpensive, analytical model for wind farm AEP that is specifically developed for WFLO applications. In this paper, we analyze the performance of the FLOWERS AEP model with analytic gradients in a layout optimization study compared with a reference optimization framework across three wind farm case studies. We find that the FLOWERS-based approach reduces computation time by a factor of 50–4000 and improves optimal AEP by about 0.3% with less than half of the variability in AEP across instances with randomized initial conditions. We also find the optimal layouts to be insensitive to model parameter tuning, making FLOWERS-based layout optimization a streamlined, user-friendly approach.

17 WIND ENERGY↗

An immersed interface method for microstructure-scale electrochemical battery models: numerical formulation and performance portable implementation

We present the numerical formulation, verification, and performance portable implementation of an immersed interface method for microstructure scale electrochemical modeling of batteries. The innovation in this approach is the resolution of chemical species and electrostatic potential discontinuities at active interfaces without the use of interface conforming unstructured grids. A unified formulation on Cartesian grids for all domains (electrodes and electrolyte) is used with interfacial flux conditions applied using volume fraction or “color” function gradients. We have developed one dimensional and two dimensional test cases with analytic solutions for electrochemical modeling using which we verified the consistency and accuracy of our scheme. Our solver is also validated against solutions from a macroscale model and an unstructured multi-subdomain solver for a full lithium ion cell. We then demonstrated the utility of our solver on an image-based complex battery electrode microstructure at high charging rate. Our technique also exhibits good scalability on distributed memory architectures using central processing units (CPU), with problem sizes up to 1.8 billion degrees of freedom and with 5400 ranks. Initial performance studies of our open-source performance portable solver showed about 70 times speed up using a graphics processing unit (GPU) compared to single compute core for a problem with 4 million cells.

25 ENERGY STORAGE↗

An Optimized Parameterization of Sub‐Grid Scale Advection for Convection Permitting Models

Convection‐permitting models (CPMs) explicitly resolve deep convection yet under‐resolve the organized lateral exchanges among drafts and their environment that control entrainment/detrainment, precipitation efficiency, and mesoscale structure. In this work, we introduce the Optimized Advection Scheme (OAS), which introduces a small rotation of the Cartesian frame of reference for the horizontal winds relative to other variables used in advection that induces cross‐gradient transport to mimic under‐resolved convective mixing. The rotation angle is selected to minimize the Kullback–Leibler divergence between the simulated and satellite observed precipitation intensity distributions, yielding a physically consistent perturbation that is computationally inexpensive and portable. Optimized Advection Scheme is implemented in WRF and evaluated over Amazon (April 2014). It shifts precipitation–precipitable‐water joint distributions toward lighter rain, reduces overly intense rates, and improves mesoscale convective system (MCS) lifetime and propagation. Mechanistically, the added cross‐gradient transport promotes convective detrainment and environmental mixing, which cools and moistens the mid‐troposphere, weakens downward momentum transport, alleviates excessive downwelling shortwave biases, and warms the surface temperature. The optimized rotation angle yields comparable improvements at 4‐km and 1‐km grid spacing, demonstrating resolution‐independent benefits across the CPM gray zone. By targeting the dynamical root of under‐mixed convective circulations, rather than tuning model microphysics or closures, OAS delivers robust, scale‐aware improvements in precipitation statistics, cloud vertical structure, and characteristics of MCS (MCSs), offering a practical pathway to more reliable CPM simulations for weather and climate applications.

CPM↗

Spatial Microsimulation and Activity Allocation in Python: An Update on the Likeness Toolkit

Understanding human security and social equity issues within human systems requires large-scale models of population dynamics that simulate high-fidelity representations of individuals and access to essential activities (work/school, social, errands, health). Likeness is a Python toolkit that provides these capabilities for Oak Ridge National Laboratory's (ORNL) UrbanPop spatial microsimulation project. In step with the initial development phase for Likeness (2021 - 2022), we built out several foundational examples of work/school and health service access. In this paper, we describe expansion and scaling of Likeness capabilities to metropolitan areas in the United States. We then provide an integrated demonstration of our methods based on a case study of Leon County, FL and perform validation exercises on 1) neighborhood demographic composition and 2) visits by demographic cohorts (gender/age) obtained from point of interest (POI) footfall data for essential services (grocery stores). Taking into account lessons learned from our case study, we scope improvements to our model as well as provide a roadmap of the anticipated Likeness development cycle into 2023 - 2024.

Tuccillo, Joe↗

Photosynthetic capacity is reduced by warming but unaffected by elevated CO2 in seedlings of five boreal tree species

Abstract Increasing atmospheric CO2 concentrations fuel global warming, with boreal regions warming at a faster rate than many other areas. Boreal forests are an important component of the global carbon cycle, yet we have little data on photosynthetic responses of boreal trees to elevated CO2 (EC) and warming. We grew seedlings of 5 widespread North American boreal tree species (from Betula, Larix, Picea, and Pinus) under current (410 ppm) or elevated (750 ppm) CO2 and either ambient (+0 °C) or increased (+4 °C or +8 °C) temperature, then measured photosynthetic traits over a range of leaf temperatures. Our results were generally consistent across species: photosynthetic capacity (maximum rates of Rubisco carboxylation, Vcmax, and electron transport, Jmax) was unaffected by EC but decreased under +8 °C warming. Accordingly, net photosynthesis measured at the growth CO2 concentration (Agrowth) was reduced under warming and increased under EC. The thermal optimum for Agrowth (ToptA) increased by ∼1.8 °C with EC but increased with warming in only two species. In contrast, the activation energies and thermal optima for Vcmax and Jmax, which are used to estimate photosynthesis in Earth System Models, were unaffected by growth environment. There were a few interactions between growth, CO2, and warming. These results suggest increased photosynthesis of widespread boreal tree species under EC may be offset by future reductions in photosynthetic capacity related to warming. We also show that the temperature sensitivities of parameters used to estimate global photosynthesis in large-scale models are generally unaffected by simulated climate change in these species.

Plant Sciences↗

EAGLES Liquid Cloud Testbed Large Eddy Simulation Library (v3)

This library consists of large eddy simulation (LES) model output using the PINACLES codebase coupled to the Hebrew University Fast Spectral Bin Microphysics scheme representing shallow, liquid phase clouds from four global liquid cloud testbed regions. This data is particularly suitable for examining interactions among cloud, aerosol, and turbulence in warm, boundary-layer clouds and for developing, assessing and refining parameterizations for coarser-scale models.

Kaul, Colleen M↗

Inference offers a metric to constrain dynamical models of neutrino flavor transformation

The multimessenger astrophysics of compact objects presents a vast range of environments where neutrino flavor transformation may occur and may be important for nucleosynthesis, dynamics, and a detected neutrino signal. Development of efficient techniques for surveying flavor evolution solution spaces in these diverse environments, which augment and complement existing sophisticated computational tools, could leverage progress in this field. To this end we continue our exploration of statistical data assimilation (SDA) to identify solutions to a small-scale model of neutrino flavor transformation. SDA is a machine learning formula wherein a dynamical model is assumed to generate any measured quantities. Specifically, we use an optimization formulation of SDA wherein a cost function is extremized via the variational method. Regions of state space in which the extremization identifies the global minimum of the cost function will correspond to parameter regimes in which a model solution can exist. Our example study seeks to infer the flavor transformation histories of two monoenergetic neutrino beams coherently interacting with each other and with a matter background. We require that the solution be consistent with measured neutrino flavor fluxes at the point of detection, and with constraints placed upon the flavor content at various locations along their trajectories, such as the point of emission, and the locations of the Mikheyev-Smirnov-Wolfenstein resonances. We show how the procedure efficiently identifies solution regimes and rules out regimes where solutions are infeasible. Overall, results in this work intimate the promise of this “variational annealing” methodology to efficiently probe an array of fundamental questions that traditional numerical simulation codes render difficult to access.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Multiscale modeling and nonlinear model predictive control for flue gas desulfurization

The primary source of sulfur dioxide (SO 2 ) emissions is flue gas from fossil fuels-based power plants. SO 2 emissions are known to not only cause health issues, but also have an adverse effect on the environment in various ways. Several Flue Gas Desulfurization (FGD) technologies have been incorporated in power plants. The most popular technology is Wet FGD, where a limestone slurry is used to absorb SO 2 from the flue gas. A detailed droplet scale model describing the instantaneous and finite rate chemistry is developed. The ill-posed Differential-Algebraic Equation (DAE) droplet model is reformulated to a well-posed index-1 DAE through index reduction. The droplet model is integrated with the bulk phase by incorporating gas-liquid mass transfer, and an oxidation reactor model to simulate the dynamic operation of the counter-current spray scrubber. As a result, the model has a well-conditioned Jacobian and overcomes the modeling challenges of previous works and enables numerical solution without requiring carefully selected initialization or specialized solution procedures. The model is successfully validated using power plant measurements, and nonlinear model predictive control (NMPC) studies are demonstrated to optimize recycle stream flowrates to minimize pumping costs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Modeling Nanoscale Ohmics in Carbon Supports of Fuel Cell Cathodes

Here, reducing platinum (Pt) loading in polymer electrolyte fuel cells (PEFCs) while meeting performace requirements is critical to their widespread deployment. However, significant polarization losses manifest at higher current densities in cathodes with lower Pt content. The morphology of the carbon supports in PEFC cathodes affects the location of Pt deposition into the micropits or onto the surface of the carbon support, translating into different kinetic and transport resistances. In this work, we derive an agglomerate scale model that differentiates the sink terms for Pt on the surface and in the pits of carbon supports. We develop an approach to assess the impact of nanoscale ohmic resistance to Pt in the micropits arising from weakly ionic solution in the carbon support on PEFC performance. Effectiveness factors relating the actual reaction rate to the maximum reaction rate (had all the Pt been exposed) are derived and embedded into a one-dimensional catalyst layer model. Parameters in the catalyst layer model are tuned based on experimental local oxygen transport resistances. Subsequently, we estimate bounds for the micropore resistances based on geometric and physical arguments. Lastly, polarization curves are simulated to assess the effect of the micropore resistance in fully-humidified and oxygen-rich environments.

25 ENERGY STORAGE↗

Assessing the Role of Interfacial and Metal Sites in Pt/TiO2-Catalyzed Acetic Acid Hydrodeoxygenation

Thermochemical conversion of biomass to produce drop-in quality biofuels typically involves hydrodeoxygenation (HDO) steps following catalytic fast pyrolysis (CFP) to remove excess oxygen and create a more stable bio-oil product. HDO involves co-feeding the CFP vapor-phase product and H2 gas over a bi-functional catalyst. Reducible metal oxide-supported noble-metal catalysts (e.g., Pt/TiO2) are promising materials for HDO chemistry, with recent work aiming to elucidate the role of various Pt/TiO2 actives sites (i.e., Pt-metal, TiO2-support, Pt-TiO2-interfacial sites) in the competing desired deoxygenation and undesired decarboxylation/decarbonylation reactions for important classes of CFP vapor model compounds. Carboxylic acids are one important class of bio-derived compounds that has not been studied extensively for HDO on Pt/TiO2, particularly at the atomic level. Past experimental work evaluating Pt/C and Pt/TiO2 catalysts for acetic acid HDO (AA-HDO) demonstrated selectivity toward C-C and C-O bond-dissociation products, respectively. This work utilizes atomic-scale modeling to discern the role of Pt-metal and Pt-TiO2-interfacial sites in promoting key C-C bond-breaking, C-O bond-breaking, and (de)hydrogenation steps in AA-HDO. Using Pt(111) and Pt-TiO2-interface surface models, adsorption and reaction energetics calculated by density functional theory provide fundamental insights into the role of interfacial sites and oxygen vacancies in promoting desired deoxygenation pathways over undesired decarboxylation/decarbonylation pathways.

BIOMASS FUELS,INORGANIC, ORGANIC, PHYSICAL, AND AN↗

Revisiting a Drag Partition Model For Canopy-Like Roughness Elements

Turbulent flows over a large surface area (S) covered by n obstacles experience an overall drag due to the presence of the ground and the protruding obstacles into the flow. The drag partition between the roughness obstacles and the ground is analyzed using an analytical model proposed by Raupach (Boundary-Layer Meteorol 60:375-395, 1992) and is hereafter referred to as R92. The R92 is based on the premise that the wake behind an isolated roughness element can be described by a shelter area A and a shelter volume V. The individual sizes of A and V without any interference from other obstacles can be determined from scaling analysis for the spread of wakes. To upscale from an individual roughness element to n/S elements where wakes may interact, R92 adopted a background stress re-normalizing instead of reducing A or V with each element addition. This work demonstrates that R92’s approach results in a linear background stress reduction in A and V only when the ratio of n/S is small, due to a low probability of wake interactions. This probabilistic nature suggests that up-scaling from individual to multiple roughness elements can be re-formulated using stochastic averaging methods proposed here. The two approaches are shown to recover R92 under plausible conditions. An alternative scaling for the shelter volume is also proposed here using thermodynamic arguments of work and dissipation though the final outcome remains similar to R92. Comparisons between R92 and available data spanning more than two decades after R92 on blocks and vegetation-like roughness elements confirm the practical utility of R92. The agreement between R92 and this updated databases of experiments and simulations confirm the potential use of R92 in large-scale models provided that the relevant parameters accommodate certain features of the roughness element type (cube versus vegetation-like) and, to a lesser extent, their configuration throughout S. Last, a comparison between R92 and models based on first-order closure principles with constant mixing length suggests that R92 can outperform such models when evaluated across a wide range of roughness densities.

canopy turbulence↗

Dynamic density functional theory of polymers with salt in electric fields

Here we present a dynamic density functional theory for modeling the effects of applied electric fields on the local structure of polymers with added salt (polymer electrolytes). Time-dependent equations for the local electrostatic potential and volume fractions of polymer, cation, and anion of added salt are developed using the principles of linear irreversible thermodynamics. For such a development, a field theoretic description of the free energy of polymer melts doped with salts is used, which captures the effects of local variations in the dielectric function. Connections of the dynamic density functional theory with experiments are established by relating the three phenomenological Onsager’s transport coefficients of the theory to the mutual diffusion of electrolyte, ionic conductivity, and transference number of one of the ions. The theory is connected with a statistical mechanical model developed by Bearman and Kirkwood [J. Chem. Phys. 28, 136 (1958)] after relating the three transport coefficients to friction coefficients. The steady-state limit of the dynamic density functional theory is used to understand the effects of dielectric inhomogeneity on the phase separation in polymer electrolytes. The theory developed here provides not only a way to connect with experiments but also to develop multi-scale models for studying connections between local structure and ion transport in polymer electrolytes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Quasi-mode evolution in a stochastic magnetic field

Abstract We present a multi-scale model of quasi-mode evolution in a stochastic magnetic field. The similarity between a quasi-mode and a ballooning mode enables us to address the challenges arising from the disparate geometries in the theories of ballooning modes in the presence of resonant magnetic perturbations. We obtain useful insights into our understanding of ballooning mode dynamics in a stochastic background. To maintain quasi-neutrality at all scales, the beat between the quasi-mode and the stochastic magnetic field drives microturbulence, which drives the turbulent background that promotes mixing and damps the quasi-mode. As a result of the broad mode structure of the quasi-mode, the turbulent viscosity and the turbulent diffusivity produced by the microturbulence are larger than those in our related study on resistive interchange modes. The stochastic magnetic field can also enhance the effective plasma inertia and reduce the effective drive, thereby slowing the mode growth. A nontrivial correlation between the microturbulence and the magnetic perturbations is shown to develop. This could account for the reduction in the Jensen–Shannon complexity of pedestal turbulence in the Resonant Magnetic Perturbation Edge-Localized Mode suppression phase observed in recent experiments. Directions for future experimental and theoretical studies are suggested.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Flexible User-Defined Domain Decomposition in Kilometer-Scale E3SM Land Model Simulation

The Energy Exascale Earth System Model (E3SM) Land Model (ELM) has been extended to kilometer-scale (km-ELM) resolutions, enabling high-fidelity simulations of terrestrial processes at 1 km x 1 km grid spacing. In ELM, domain decomposition partitions the computational domain across processors, ensuring efficient parallel execution. Currently, round-robin decomposition is applied, providing a straightforward way to distribute computational workload. As ELM continues evolving at the kilometer-scale (km-scale), particularly with integrating lateral flow modeling, decomposition strategies must also account for the increased workload and data movement. This paper introduces a flexible user-defined domain decomposition framework, allowing users to customize domain partitioning based on application requirements. The impact of different decomposition strategies is evaluated across various applications concerning computation, communication, and I/O. Results demonstrate that while 1D partitioning yields superior I/O performance, k-nearest neighbors (KNN) clustering effectively reduces inter-process communication overhead. This study lays the groundwork for scalable partitioning in large-scale land surface simulations, enhancing next-generation Earth system modeling.

Wang, Dali [ORNL] (ORCID:0000000168065108)↗

A unified framework of stabilized finite element method for solving the Boltzmann transport equation

This paper presents a unified framework of stabilized finite element method for solving the Boltzmann transport equation. Unlike the traditional Petrov-Galerkin finite element method which modifies the test function to construct the stabilization term, we derive the stabilization methods from the standard Galerkin weak form with Sub-grid scale model. The basic idea of this method is to decompose the unknowns into its numerical solution and residual, with an approximation for the residual and embeds it in the Galerkin weak form to yield a stabilized variational formula. Different approximations of the residual lead to different stabilization methods, all the frequently used stabilized methods, including the Streamline Upwinding Petrov-Galerkin (SUPG) method, Galerkin/Least-Square (GLS) method, and Algebraic Sub-Grid Scale (ASGS) method can be obtained from this framework. The similarities and differences of the different approximations are compared in this paper. The numerical results show that the behaviors of the different methods area similar with the same stabilization parameter, and all these stabilized techniques can obtain a correct and stable solution. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

BioRT‐HBV 1.0: A Biogeochemical Reactive Transport Model at the Watershed Scale

Abstract Reactive Transport Models (RTMs) are essential tools for understanding and predicting intertwined ecohydrological and biogeochemical processes on land and in rivers. While traditional RTMs have focused primarily on subsurface processes, recent watershed‐scale RTMs have integrated ecohydrological and biogeochemical interactions between surface and subsurface. These emergent, watershed‐scale RTMs are often spatially explicit and require extensive data, computational power, and computational expertise. There is however a pressing need to create parsimonious models that require minimal data and are accessible to scientists with limited computational background. To that end, we have developed BioRT‐HBV 1.0, a watershed‐scale, hydro‐biogeochemical RTM that builds upon the widely used, bucket‐type HBV model known for its simplicity and minimal data requirements. BioRT‐HBV uses the conceptual structure and hydrology output of HBV to simulate processes including advective solute transport and biogeochemical reactions that depend on reaction thermodynamics and kinetics. These reactions include, for example, chemical weathering, soil respiration, and nutrient transformation. The model uses time series of weather (air temperature, precipitation, and potential evapotranspiration) and initial biogeochemical conditions of subsurface water, soils, and rocks as input, and output times series of reaction rates and solute concentrations in subsurface waters and rivers. This paper presents the model structure and governing equations and demonstrates its utility with examples simulating carbon and nitrogen processes in a headwater catchment. As shown in the examples, BioRT‐HBV can be used to illuminate the dynamics of biogeochemical reactions in the invisible, arduous‐to‐measure subsurface, and their influence on the observed stream or river chemistry and solute export. With its parsimonious structure and easy‐to‐use graphical user interface, BioRT‐HBV can be a useful research tool for users without in‐depth computational training. It can additionally serve as an educational tool that promotes pollination of ideas across disciplines and foster a diverse, equal, and inclusive user community.

Sadayappan, Kayalvizhi↗

High-Voltage DC MULTI-terminal SIMulation (HVDC MULTISIM): Technical Program Summary

Expansion of the power grid in the USA is essential to fulfilling the rising energy demand of the country. In particular, the transmission grid, a network of electrical energy corridors that enable the flow of large amounts of energy from the point of generation to the point of consumption, needs significant expansion to cater to this growth of electrification. The current infrastructure is also dated, and any new transmission corridor should be based on a vision of building a modern infrastructure that is futureproof. High voltage DC (HVDC) transmission technology falls in this category: it is the most economical way to build large transmission lines and relies on sophisticated electronics and flexible controls, rather than just passive components like switchgear and transformers. This project is aimed at establishing the technical and economic feasibility of a network of HVDC lines that are interconnected to form a multi-terminal DC network. HVDC transmission is not a new technology and United States has a number of such lines; however, these are point to point transmissions and do not form a DC grid. Modeling and analyzing DC grids is challenging since there are no established modeling tools and the traditional ways of modeling AC grids fall short of providing the fidelity required the fast dynamics of the DC grid. This project aims to build a software-in-the-loop simulator for a multi-terminal HVDC grid, called MULTISIM and, establish new control and protection algorithms to operate such a connected DC system. NLR, a key partner, established the core framework of the MULTISIM simulator and validated its functionality. The GE Vernova team, created a full-scale model of a four-terminal HVDC network using PSCAD (Power Systems Computer Aided Design) software and established the baseline performance of the system during normal and fault operation. The project was started on 10/12024 with a kick-off meeting held on 12/11/2024 and ran through two quarters till termination. All the tasks, milestones and deliverables during this period were met. Several interim reports detailing the various tasks were submitted. The following sections provide details on the completed tasks till the project pause and termination after Q2.

24 POWER TRANSMISSION AND DISTRIBUTION↗