Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “partitioned methods”

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

Distributed PV Hosting Capacity Evaluation Considering Equitable PV Accommodation: Preprint

In distribution systems, distributed photovoltaic (DPV) hosting capacity is widely used by system operators to evaluate how much DPV the distribution networks can accommodate without violating operational constraints. With the proliferation of DPV deployment, if managed properly, it offers a great opportunity to mitigate the persistent energy inequity issue in electric grid. Otherwise, energy inequity can be exacerbated with unequal DPV access when hosting capacity improvement plans ignore the equity factor. To consider energy equity in distribution network planning for enhancing DPV accommodation, the hosting capacity evaluation should be improved to consider equity constraints endogenously. In this paper, a DPV hosting capacity evaluation method is proposed considering equitable DPV access among sub-regions of a distribution feeder partitioned by social-economic indicators such as the income levels. The results demonstrate that considering equity in DPV hosting capacity evaluation can provide valuable reference in the equitable energy transition. This work provides a valuable guidance for future distribution network upgrade planning for equitable DPV deployment.

distributed PV↗
DOE OSTI2023

A sharp interface Lagrangian-Eulerian method for flexible-body fluid-structure interaction

This paper introduces a sharp-interface approach to simulating fluid-structure interaction (FSI) involving flexible bodies described by general nonlinear material models and across a broad range of mass density ratios. This new flexible-body immersed Lagrangian-Eulerian (ILE) scheme extends our prior work on integrating partitioned and immersed approaches to rigid-body FSI. Our numerical approach incorporates the geometrical and domain solution flexibility of the immersed boundary (IB) method with an accuracy comparable to body-fitted approaches that sharply resolve flows and stresses up to the fluid-structure interface. Unlike many IB methods, our ILE formulation uses distinct momentum equations for the fluid and solid subregions with a Dirichlet-Neumann coupling strategy that connects fluid and solid subproblems through simple interface conditions. As in earlier work, we use approximate Lagrange multiplier forces to treat the kinematic interface conditions along the fluid-structure interface. This penalty approach simplifies the linear solvers needed by our formulation by introducing two representations of the fluid-structure interface, one that moves with the fluid and another that moves with the structure, that are connected by stiff springs. This approach also enables the use of multi-rate time stepping, which allows us to use different time step sizes for the fluid and structure subproblems. Our fluid solver relies on an immersed interface method (IIM) for discrete surfaces to impose stress jump conditions along complex interfaces while enabling the use of fast structured-grid solvers for the incompressible Navier-Stokes equations. The dynamics of the volumetric structural mesh are determined using a standard finite element approach to large-deformation nonlinear elasticity via a nearly incompressible solid mechanics formulation. This formulation also readily accommodates compressible structures with a constant total volume, and it can handle fully compressible solid structures for cases in which at least part of the solid boundary does not contact the incompressible fluid. Selected grid convergence studies demonstrate second-order convergence in volume conservation and in the pointwise discrepancies between corresponding positions of the two interface representations as well as between first and second-order convergence in the structural displacements. The time stepping scheme is also demonstrated to yield second-order convergence. To assess and validate the robustness and accuracy of the new algorithm, comparisons are made with computational and experimental FSI benchmarks. Test cases include both smooth and sharp geometries in various flow conditions. Furthermore, we also demonstrate the capabilities of this methodology by applying it to model the transport and capture of a geometrically realistic, deformable blood clot in an inferior vena cava filter.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

An Orthogonal Recursive Bisection (ORB) Based Time Advancement Algorithm for CFD-DEM Solvers

The time integration of the granular phase in coupled computational fluid dynamics (CFD) – discrete element method (DEM) simulations presents a unique computational challenge brought about by the large variations in particle collisional time scales. Particles in the dilute regions of the computational domain can be advanced with large time steps while dense regions require much smaller time increments. However, the time step size in most solvers is globally set as the limit for accuracy and stability imposed by the collisions and is typically orders of magnitude less than that required away from collisions. This work addresses this precise issue and provides a strategy to avoid the use of a global conservative small time step size for the entire set of particles.A novel time stepping algorithm for CFD-DEM solvers using a partitioning approach using orthogonal recursive bisection (ORB) that allows for variable time steps among particles is described and its computational performance is compared against baseline explicit methods, typically used in several CFD-DEM solvers. ORB has advantages of being relatively quick and easy to update incrementally and has the required heuristic behavior (i.e., it will split the region in half with a cluster on each side) when groups of particles are well separated (clustered). The algorithm presented in this work uses a local time stepping approach to resolve collisional time scales for subsets of particles that are present at the leaves of the ORB, thereby resulting in substantial reduction of computational cost. The parallel implementation of this method where a ``knapsack” algorithm is used in tandem with ORB for effective load-balancing is also presented, where a best possible partitioning is obtained based on number of particles and local time-stepping costs. The algorithm is tested against benchmark problems with varying particle distributions that include fluidized bed and riser flow scenarios. Preliminary results indicate that the approach is 2-3X faster than traditional explicit methods for problems that involve both dense and dilute regions, while maintaining the same level of accuracy.

adaptive timestepping↗
DOE OSTI2020

A one-loop test of the near-AdS 2 /near-CFT 1 correspondence

We analyze quantum fluctuations around black hole solutions to the Jackiw-Teitelboim model. We use harmonic analysis on Euclidean AdS 2 to show that the logarithmic corrections to the partition function are determined entirely by quadratic holomorphic differentials, even when conformal symmetry is broken and harmonic modes are no longer true zero modes. Our quantum-corrected partition function agrees precisely with the SYK result. We argue that our effective quantum field theory methods and results generalize to other theories of two-dimensional dilaton gravity.

2D Gravity↗
DOE OSTI2023

A Robust Parallel Distributed State Estimation for Large Scale Distribution Systems

The growing need and interest in real-time monitoring of large distribution networks motivated by the rapid population of renewable sources, EVs and etc. demand a computationally efficient state estimation framework. Furthermore, this paper presents an improved computational framework for implementing a robust state estimator using a multi-core processor. The main contribution of the paper is the proposed computational framework along with two partitioning strategies which enable fast and robust state estimation for large scale radial and/or meshed distribution systems. Formulation of the proposed method and its implementation are described in detail. Performance of the estimator is tested by simulations first using a small 84-bus radial distribution system. Then the method’s scalability is demonstrated by simulations on two very large scale distribution networks one configured radially and the other meshed each containing over 12,500 buses.

42 ENGINEERING↗
DOE OSTI2020

Impact of the adhesion-force lever-arm “a” on the Rock ‘n’ Roll resuspension model and how to compute it from contact mechanics

The paper proposes a methodology to fill a gap in the context of mechanistic particle resuspension models related to the adhesion and lift force lever arms ( a ): to date, only one experimental observation exists from a centrifuge experiment with a specific particle/surface pair which found r P / a = 100 and assumed that the contact points were symmetric with respect to the particle center. However, the Reeks and Hall’s (2001) model is highly sensitive to this parameter, which therefore cannot be used as-is for different particle diameters or materials. This paper proposes a methodology to compute the adhesion and lift force lever arms ( a A and a L , respectively) based on Atomic Force Microcopy (AFM) surface roughness measurements and contact mechanics arguments. In practice, a spherical, perfectly smooth particle is brought in contact (computationally) with the measured surface roughness and is then rotated by the point of contact to determine the two or more asperities over which the particle would sit on the surface. A method to compute the asperity radius of curvature is devised to compute the asperity deformation and the correct partitioning of the total deformation between the particle and the asperity is presented. For the glass surface under consideration, the asperity deformation was negligible, i.e., it was below the van der Waals radius, which is the limit of applicability of contact mechanics. The Reeks and Hall’s (2001) model is then extended to account for the joint probability of a A and a L . Significant differences in particle resuspended after a day in common atmospheric conditions are found between the approach that uses the lever-arm distributions or their mean values. Also, important differences are found when comparing the lever-arm distribution and r P / a = 100. Hence, one should be cautious of using r P / a = 100 for surface/particle pairs differing from those in the original study. Finally, limits of applicability of the theories used are presented together with considerations regarding how to extend the method to non-spherical particles and different materials.

42 ENGINEERING↗
DOE OSTI2020

Elemental characterization of crystalline silicotitanate following Hanford tank waste processing

To expedite vitrification of low activity waste at the Hanford Site, a Tank-Side Cesium Removal system is being evaluated. This method utilizes an inorganic ion exchange (IX) media, crystalline silicotitanate (CST), to remove 99.9% of the 137Cs from tank waste supernate. Laboratory scale testing of the ion exchange process has been conducted by Pacific Northwest National Laboratory to evaluate Cs removal with CST in supernates collected from Hanford tanks 241-AP-107 and 241-AW-102 under prototypic plant operating conditions. In an effort to understand the selectivity and partitioning of other analytes toward CST, specifically Resource Conservation and Recovery Act (RCRA) hazardous metals, selected +2 cations, and transuranic elements, a method for the complete dissolution of CST and Cs removal from the digestate was applied to the spent CST post-processing Hanford tank waste. The CST was shown to adsorb a relevant fraction of Ca, Cd, Pb, Sr, and potentially Fe by IX and/or other mechanisms. The spent CST also adsorbed actinides, including U, Pu, Np, and Am. This work discusses digestion of CST and removal of Cs from the digestate allowing effective analysis and provides insight as to what analytes may be competing for IX sites on the CST media.

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

Multireference Embedding and Fragmentation Methods for Classical and Quantum Computers: From Model Systems to Realistic Applications

One of the primary challenges in quantum chemistry is the accurate modeling of strong electron correlation. While multireference methods effectively capture such correlation, their steep scaling with system size prohibits their application to large molecules and extended materials. Quantum embedding offers a promising solution by partitioning complex systems into manageable subsystems. In this Review, we highlight recent advances in multireference density matrix embedding and localized active space self-consistent field approaches for complex molecules and extended materials. We discuss both classical implementations and the emerging potential of these methods on quantum computers. Here, by extending classical embedding concepts to the quantum landscape, these algorithms have the potential to expand the reach of multireference methods in quantum chemistry and materials.

Algorithms↗
DOE OSTI2025

Robust Multi-fidelity Bayesian Optimization with Deep Kernel and Partition

Multi-fidelity Bayesian optimization (MFBO) is a powerful approach that utilizes lowfidelity, cost-effective sources to expedite the exploration and exploitation of a high-fidelity objective function. Existing MFBO methods with theoretical foundations either lack justification for performance improvements over single-fidelity optimization or rely on strong assumptions about the relationships between fidelity sources to construct surrogate models and direct queries to low-fidelity sources. To mitigate the dependency on cross-fidelity assumptions while maintaining the advantages of low-fidelity queries, we introduce a random sampling and partition-based MFBO framework with deep kernel learning. This framework is robust to cross-fidelity model misspecification and explicitly illustrates the benefits of low-fidelity queries. Our results demonstrate that the proposed algorithm effectively manages complex cross-fidelity relationships and efficiently optimizes the target fidelity function.

Zhang, Fengxue [University of Chicago, Illinois, U↗
DOE OSTI2022

Path Integrals for Nonadiabatic Dynamics: Multistate Ring Polymer Molecular Dynamics

This review focuses on a recent class of path-integral-based methods for the simulation of nonadiabatic dynamics in the condensed phase using only classical molecular dynamics trajectories in an extended phase space. Specifically, a semiclassical mapping protocol is used to derive an exact, continuous, Cartesian variable path-integral representation for the canonical partition function of a system in which multiple electronic states are coupled to nuclear degrees of freedom. Building on this exact statistical foundation, multistate ring polymer molecular dynamics methods are developed for the approximate calculation of real-time thermal correlation functions. As a result, the remarkable promise of these multistate ring polymer methods, their successful applications, and their limitations are discussed in detail.

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

Deep learning for characterizing the self-assembly of three-dimensional colloidal systems

Creating a systematic framework to characterize the structural states of colloidal self-assembly systems is crucial for unraveling the fundamental understanding of these systems’ stochastic and non-linear behavior. The most accurate characterization methods create high-dimensional neighborhood graphs that may not provide useful information about structures unless these are well-defined reference crystalline structures. Dimensionality reduction methods are thus required to translate the neighborhood graphs into a low-dimensional space that can be easily interpreted and used to characterize non-reference structures. We investigate a framework for colloidal system state characterization that employs deep learning methods to reduce the dimensionality of neighborhood graphs. The framework next uses agglomerative hierarchical clustering techniques to partition the low-dimensional space and assign physically meaningful classifications to the resulting partitions. In this work, we first demonstrate the proposed colloidal self-assembly state characterization framework on a three-dimensional in silico system of 500 multi-flavored colloids that self-assemble under isothermal conditions. We next investigate the generalizability of the characterization framework by applying the framework to several independent self-assembly trajectories, including a three-dimensional in silico system of 2052 colloidal particles that undergo evaporation-induced self-assembly.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Efficient shallow Ritz method for 1D diffusion problems

This paper studies the shallow Ritz method for solving the one-dimensional diffusion problem. It is shown that the shallow Ritz method improves the order of approximation dramatically for non-smooth problems. To realize this optimal or nearly optimal order of the shallow Ritz approximation, we develop a damped block Newton (dBN) method that alternates between updates of the linear and non-linear parameters. Per each iteration, the linear and the non-linear parameters are updated by exact inversion and one step of a modified, damped Newton method applied to a reduced non-linear system, respectively. The computational cost of each dBN iteration is $\mathcal{O}$(n). Starting with the non-linear parameters as a uniform partition of the interval, numerical experiments show that the dBN is capable of efficiently moving mesh points to nearly optimal locations. In conclusion, to improve the efficiency of the dBN further, we propose an adaptive damped block Newton (AdBN) method by combining the dBN with the adaptive neuron enhancement (ANE) method [28].

Diffusion problems↗
DOE OSTI2020

Economic Extraction and Recovery of REEs and Production of Clean Value-Added Products from Low-Rank Coal Fly Ash

The University of North Dakota (UND) Energy & Environmental Research Center (EERC) teamed with Pacific Northwest National Laboratory (PNNL) , the North Dakota Industrial Commission Lignite Research Program, and commercial partners Basin Electric Power Cooperative, Great River Energy, and Southern Company to execute a Cooperative Agreement funded by the U.S. Department of Energy’s (DOE’s) National Energy Technology Laboratory (NETL) focused on identifying unique pathways and pretreatments to extract rare-earth elements (REEs) from low-rank coal (LRC) ash in a more economical and environmentally benign manner than current methods. The EERC drew upon the unique and deep knowledge of the team regarding LRC geochemistry, ash transformation mechanisms during combustion, and fly ash chemistry in the approach to understand the portioning of REEs in LRC ash. A comprehensive characterization effort was undertaken to fully establish the form, associations, and partitioning of the REEs and other elements/minerals of interest in fly ash, as well as the ash chemistry, mineralogy, and morphology. Using this information, the team tested methods of REE extraction to take advantage of the unique and advantageous properties of LRC ash. Multiple complementary approaches were evaluated that combined both novel and proven methods to develop an extractive process for REEs and other high-value metals from LRC ash that provide technological and economic progress on the current state of the art. Methodologies were also evaluated for the potential to synergistically remove toxic metals while extracting other high-value metals from the ash, generating a clean, safe ash material as a value-added product. This would eliminate waste and offer significant additional revenue potential. Based on laboratory-scale testing, an initial high-level technical and economic analysis to estimate operating expenses and product revenues to guide future phases of technology development was implemented. By focusing on minimizing process steps and consumptive use of extraction solutions, a tunable, economically viable process for REE extraction from LRC ash was developed. The tunable process has the ability to be adjusted to adapt to the differing ash chemical and physical properties that are characteristic of LRC ash. This is key to deployment as a commercial process. In demonstrating the process, the EERC was able to illustrate mixed REE concentration greater than the goal of 2 wt%. High-level economic analysis of the process with LRC ash of moderate REE levels demonstrated the potential of the process, but recovery of just the REEs was insufficient to cover operational costs. With the additional extraction of other high-value metals found in LRC ash, the process economics approach would break even, and it is anticipated that with process optimization, the developed process could potentially be economically sound. The same economic analysis performed utilizing LRC ash with high levels of REEs showed the recovery of the REEs alone would make the process economically viable even without further optimization.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2020

Optimized Lie–Trotter–Suzuki decompositions for two and three non-commuting terms

Lie–Trotter–Suzuki decompositions are an efficient way to approximate operator exponentials exp ( t H ) when H is a sum of n (non-commuting) terms which, individually, can be exponentiated easily. They are employed in time-evolution algorithms for tensor network states, digital quantum simulation protocols, path integral methods like quantum Monte Carlo, and splitting methods for symplectic integrators in classical Hamiltonian systems. Here, we provide optimized decompositions up to order t 6 . The leading error term is expanded in nested commutators (Hall bases) and we minimize the 1-norm of the coefficients. For n = 2 terms, several of the optima we find are close to those in McLachlan (1995). Generally, our results substantially improve over unoptimized decompositions by Forest, Ruth, Yoshida, and Suzuki. We explain why these decompositions are sufficient to efficiently simulate any one- or two-dimensional lattice model with finite-range interactions. This follows by solving a partitioning problem for the interaction graph.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Comparison of measurement techniques and sorption of radium-226 in low and high salinity aqueous samples

Human activities have the potential to redistribute radium (Ra) in the marine environment in a manner that may necessitate monitoring or management of subsequent human or environmental exposures. There is therefore a need to identify accurate and accessible techniques for Ra measurement in high salinity samples and to describe the distribution of Ra in estuarine and marine environments, but most efforts in these areas have focused on low salinity matrices. In addition, rapid and reliable measurements are crucial for time-sensitive samples such as short-lived isotopes or emergency situations. The objective of this study is to describe the limits of detection, cost, and relative ease for measurement of Ra in both low and high salinity aqueous samples via three analytical methods: liquid scintillation counting (LSC), high purity germanium (HPGe) gamma spectrometry, and inductively coupled plasma mass spectrometry (ICP-MS). To contextualize these measurements for real-world scenarios, the partitioning of 226 Ra to substrates relevant to the marine environment was also characterized. Although HPGe detection with solid phase extraction had the lowest limit of detection for low salinity samples (0.27 Bq L −1 ), poor 226 Ra recovery for high salinity samples and high materials costs make this method prohibitive for many users. Limits of detection for high salinity samples were lower for LSC (1.28 Bq L −1 ) than for ICP-MS without dilution (11.4 Bq L −1 ), but significant and unexpected degradation of the high salinity LSC standards was observed after six months. Furthermore, our preferred measurement method for high salinity Ra samples is ICP-MS with sample dilution as necessary to reduce matrix effects.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2022

Facile alignment estimation in carbon nanotube films using image processing

Whether a macroscopic assembly of carbon nanotubes can exhibit the one-dimensional properties expected from individual nanotubes critically depends on how well the nanotubes are aligned inside the assembly. Therefore, a simple and accurate method for assessing the degree of alignment is desired for the rapid characterization of carbon nanotube films and fibers. Here, we present an end-to-end solution for determining the global and local spatial orientation of carbon nanotubes in films within a short amount of time using a fast, precise, and economical approach based on an image processing method applied to scanning electron microscopy images. Further, we first use Laplacian edge enhancement filtering for improving the appearance of edge regions, which is followed by image partitioning into multiple blocks to capture the nanoscale orientation characteristics and total variation-based image decomposition of these image blocks. We then perform a 2D-fast Fourier transform on the image decomposed textural components of these edge-enhanced image blocks to determine the orientation distribution, which is utilized to estimate the 2D nematic order parameter. To show the effectiveness of our method, we corroborated our results against results obtained with other state-of-the-art image processing and experimental techniques.

2D-FFT↗
DOE OSTI2021

Local indistinguishability and incompleteness of entangled orthogonal bases: Method to generate two-element locally indistinguishable ensembles

Highlights: • Local indistinguishability of orthogonal quantum states. • Bipartite and multipartite quantum ensembles. • Unextendible entangled bases and uncompletable entangled bases. • Two-element locally indistinguishable quantum ensembles. • Multiparty unextendible entangled bases with unextendibility in all partitions. We relate the phenomenon of local indistinguishability of orthogonal states with the properties of unextendibility and uncompletability of entangled bases for bipartite and multipartite quantum systems. We prove that all two-qubit unextendible entangled bases are of size three and they cannot be perfectly distinguished by separable measurements. We identify a method of constructing two-element orthogonal ensembles, based on the concept of unextendible entangled bases, that can potentially lead to information sharing applications. Two-element ensembles form the fundamental unit of ensembles, and yet does not offer locally indistinguishable ensembles for pure state elements. Going over to mixed states does open this possibility, but can be difficult to identify. The method provided using unextendible entangled bases can be used for their systematic generation. In multipartite systems, we find a class of unextendible entangled bases for which the unextendibility property remains conserved across all bipartitions. We also identify nonlocal operations, local implementation of which require entangled resource states from a higher-dimensional quantum system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗