Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “coupling method”

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 1,765 records · Page 98
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Crack opening calculation in phase-field modeling of fluid-filled fracture: A robust and efficient strain-based method

The phase-field method has become popular for the numerical modeling of fluid-filled fractures, thanks to its ability to represent complex fracture geometry without algorithms. However, the algorithm-free representation of fracture geometry poses a significant challenge in calculating the crack opening (aperture) of phase-field fracture, which governs the fracture permeability and hence the overall hydromechanical behavior. Although several approaches have been devised to compute the crack opening of phase-field fracture, they require a sophisticated algorithm for post-processing the phase-field values or an additional parameter sensitive to the element size and alignment. Here, we develop a novel method for calculating the crack opening of fluid-filled phase-field fracture, which enables one to obtain the crack opening without additional algorithms or parameters. Here we transform the displacement-jump-based kinematics of a fracture into a continuous strain-based version, insert it into a force balance equation on the fracture, and apply the phase-field approximation. Through this procedure, we obtain a simple equation for the crack opening which can be calculated with quantities at individual material points. We verify the proposed method with analytical and numerical solutions obtained based on discrete representations of fractures, demonstrating its capability to calculate the crack opening regardless of the element size or alignment.

58 GEOSCIENCES↗
DOE OSTI2024

Single-ancilla ground state preparation via Lindbladians

We design a quantum algorithm for ground state preparation in the early fault tolerant regime. As a Monte Carlo style quantum algorithm, our method features a Lindbladian where the target state is stationary. The construction of this Lindbladian is algorithmic and should not be seen as a specific approximation to some weakly coupled system-bath dynamics in nature. Our algorithm can be implemented using just one ancilla qubit and efficiently simulated on a quantum computer. It can prepare the ground state even when the initial state has zero overlap with the ground state, bypassing the most significant limitation of methods like quantum phase estimation. As a variant, we also propose a discrete-time algorithm, demonstrating even better efficiency and providing a near-optimal simulation cost depending on the desired evolution time and precision. Numerical simulations using Ising and Hubbard models demonstrate the efficacy and applicability of our method. Published by the American Physical Society 2024

Ding, Zhiyan (ORCID:000000018863403X)↗
DOE OSTI2025

Fast methods for multisite charge transfer. Processes II. Analytic nuclear gradients and nonadiabatic dynamics for cCASSCF(1,n) and cCASSCF(2n-1,n) wavefunctions

In this work we derive and implement analytic nuclear gradients and derivative couplings for a constrained complete active space self-consistent field with a small active space designed to model electron or hole transfer. Using a Lagrangian formalism, we are able to differentiate both the CASSCF energy and the constraint (which is required for smooth surfaces over a wide range of parameter space), and the resulting efficient algorithm can be immediately applied to nonadiabatic dynamics simulations of charge transfer processes. Here, we run initial surface-hopping simulations of a proton coupled electron transfer event for a phenoxyl–phenol system.

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

Coherent Magnon–Photon Coupling in the Magnetic Semiconductor CrSBr

Magnon-based hybrid quantum systems are promising candidates for quantum interconnects and quantum sensors, and they offer a rich platform for exploring nonlinear magnonics and cavity–photon interactions. Two-dimensional (2D) van der Waals magnets provide a compact, atomically flat geometry that can be easily integrated into existing quantum circuits, such as superconducting resonators and qubits. Among various 2D magnets, the magnetic semiconductor CrSBr is particularly unique due to its strong spin–exciton, spin–lattice, and magnon–exciton interactions. In this work, we demonstrate coherent coupling between antiferromagnetic (AFM) magnons in CrSBr and microwave photons in a niobium-(Nb)-based-on-chip resonator. We tuned the magnon–photon coupling strength by changing the number of CrSBr flakes integrated into the Nb microwave photon resonators. Furthermore, this work demonstrates the first step toward integrating layered van der Waals 2D magnets into superconducting microwave circuits, with full access for microwave and optical probing.

Electromagnetic radiation↗
DOE OSTI2022

Inline small-angle X-ray scattering-coupled chromatography under extreme hydrostatic pressure

As continuing discoveries highlight the surprising abundance and resilience of deep ocean and subsurface microbial life, the effects of extreme hydrostatic pressure on biological structure and function have attracted renewed interest. Biological small-angle X-ray scattering (BioSAXS) is a widely used method of obtaining structural information from biomolecules in solution under a wide range of solution conditions. Due to its ability to reduce radiation damage, remove aggregates, and separate monodisperse components from complex mixtures, size-exclusion chromatography-coupled SAXS (SEC-SAXS) is now the dominant form of BioSAXS at many synchrotron beamlines. While BioSAXS can currently be performed with some difficulty under pressure with non-flowing samples, it has not been clear how, or even if, continuously flowing SEC-SAXS, with its fragile media-packed columns, might work in an extreme high-pressure environment. Here we show, for the first time, that reproducible chromatographic separations coupled directly to high-pressure BioSAXS can be achieved at pressures up to at least 100 MPa and that pressure-induced changes in folding and oligomeric state and other properties can be observed. Further, the apparatus described here functions at a range of temperatures (0°C–50°C), expanding opportunities for understanding biomolecular rules of life in deep ocean and subsurface environments.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Probing the Potential Energy Profile of the I + (H 2 O) 3 → HI + (H 2 O) 2 OH Forward and Reverse Reactions: High Level CCSD(T) Studies with Spin-Orbit Coupling Included

Three different pathways for the atomic iodine plus water trimer reaction I + (H 2 O) 3 → HI + (H 2 O) 2 OH were preliminarily examined by the DFT-MPW1K method. Related to previous predictions for the F/Cl/Br + (H 2 O) 3 reactions, three pathways for the I + (H 2 O) 3 reaction are linked in terms of geometry and energetics. To legitimize the results, the “gold standard” CCSD(T) method was employed to investigate the lowest-lying pathway with the correlation-consistent polarized valence basis set up to cc-pVQZ(-PP). According to the CCSD(T)/cc-pVQZ(-PP)//CCSD(T)/cc-pVTZ(-PP) results, the I + (H 2 O) 3 → HI + (H 2 O) 2 OH reaction is predicted to be endothermic by 47.0 kcal mol -1 . The submerged transition state is predicted to lie 43.7 kcal mol -1 above the separated reactants. The I···(H 2 O) 3 entrance complex lies below the separated reactants by 4.1 kcal mol -1 , and spin-orbit coupling has a significant impact on this dissociation energy. The HI···(H 2 O) 2 OH exit complex is bound by 4.3 kcal mol -1 in relation to the separated products. Compared with simpler I + (H 2 O) 2 and I + H 2 O reactions, the I + (H 2 O) 3 reaction is energetically between them in general. It is speculated that the reaction between the iodine atom and the larger water clusters may be energetically analogous to the I + (H 2 O) 3 reaction. The iodine reaction I + (H 2 O) 3 is connected with the analogous valence isoelectronic bromine/chlorine reactions Br/Cl + (H 2 O) 3 but much different from the F + (H 2 O) 3 reaction. Significant difference with other halogen systems, especially for barrier heights, are seen for the iodine systems.

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

High precision tests of QCD without scale or scheme ambiguities: The 40th anniversary of the Brodsky–Lepage–Mackenzie method

A key issue in making precise predictions in QCD is the uncertainty in setting the renormalization scale μ r and thus determining the correct values of the QCD running coupling α s (μ r ) at each order in the perturbative expansion of a QCD observable. It has often been conventional to simply set the renormalization scale to the typical scale of the process Q and vary it in the range μ r $\in$ [Q/2, 2Q] in order to estimate the theoretical error. This is the practice of Conventional Scale Setting (CSS). The resulting CSS prediction will however depend on the theorist’s choice of renormalization scheme and the resulting pQCD series will diverge factorially. It will also disagree with renormalization scale setting used in QED and electroweak theory thus precluding grand unification. A solution to the renormalization scale-setting problem is offered by the Principle of Maximum Conformality (PMC), which provides a systematic way to eliminate the renormalization scale-and-scheme dependence in perturbative calculations. The PMC method has rigorous theoretical foundations, it satisfies Renormalization Group Invariance (RGI) and preserves all self-consistency conditions derived from the renormalization group. The PMC cancels the renormalon growth, reduces to the Gell-Mann–Low scheme in the N c → 0 Abelian limit and leads to scale- and scheme-invariant results. The PMC has now been successfully applied to many high-energy processes. In this article we summarize recent developments and results in solving the renormalization scale and scheme ambiguities in perturbative QCD. In particular, we present a recently developed method the PMC ∞ and its applications, comparing the results with CSS. The method preserves the property of renormalizable SU(N)/U(1) gauge theories defined as Intrinsic Conformality (iCF). This property underlies the scale invariance of physical observables and leads to a remarkably efficient method to solve the conventional renormalization scale ambiguity at every order in pQCD. This new method reflects the underlying conformal properties displayed by pQCD at NNLO, eliminates the scheme dependence of pQCD predictions and is consistent with the general properties of the PMC. A new method to identify conformal and β-terms, which can be applied either to numerical or to theoretical calculations is also shown. We present results for the thrust and C-parameter distributions in e + e - annihilation showing errors and comparison with the CSS. We also show results for a recent innovative comparison between the CSS and the PMC ∞ applied to the thrust distribution investigating both the QCD conformal window and the QED N c → 0 limit. In order to determine the thrust distribution along the entire renormalization group flow from the highest energies to zero energy, we consider the number of flavors near the upper boundary of the conformal window. In this flavor-number regime the theory develops a perturbative infrared interacting fixed point. These results show that PMC ∞ leads to higher precision and introduces new interesting features in the PMC. In fact, this method preserves with continuity the position of the peak, showing perfect agreement with the experimental data already at NNLO. We also show a detailed comparison of the PMC ∞ with the other PMC approaches: the multi-scale-setting approach (PMCm) and the single-scale-setting approach (PMCs) by comparing their predictions for three important fully integrated quantities R e+e- , R$_Τ$ and Γ (H→$b\bar{b}$) up to the four-loop accuracy.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Development of a capacitance measurement for pulsed magnetic fields

Capacitance measurements are crucial for probing the electrical properties of materials. In this study, we develop and implement a capacitance measurement technique optimized for pulsed magnetic fields. Our approach employs an auto-balancing bridge method, leveraging a high-bandwidth transimpedance amplifier to mitigate parasitic contributions from coaxial cables. This technique enables precise capacitance measurements in rapidly changing magnetic fields, as demonstrated in experiments on the magnetoelectric material NiCo 2 TeO 6 . The results reveal strong magnetoelectric coupling, including a pronounced hysteresis in capacitance that coincides with magnetization measurements and an enhanced energy dissipation peak at high sweep rates. Compared to traditional LCR meter measurements in DC fields, our method exhibits excellent agreement while providing additional insight into field-induced phase transitions. This work establishes a robust methodology for capacitance measurements under extreme conditions and opens new opportunities for studying multiferroic and correlated electron systems under high magnetic fields.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Separation and determination of ultratrace rhenium quantities in molybdenum matrix

The production of very pure molybdenum, free of Re even at ultratrace levels, is extremely important in meeting target purity specifications during the production of 99 Mo for medical applications. Here we propose two methods for separating ultratrace levels of Re from a bulk Mo matrix using either solvent extraction or extraction chromatography, combined with inductively coupled plasma quadrupole mass spectrometry (ICP-MS) detection. Re is extracted (D > 70) from solutions in sulfuric acid as a complex with tri-n-butyl phosphate. Scrubbing and washing steps result in total decontamination factors for Mo up to 1,700 for solvent extraction and 3,500 for extraction chromatography. Re is stripped into a 3 M NH 4 OH matrix and analyzed by ICP-MS. Detection (L D ) and quantification (L Q ) limits were optimized by matrix-matching to allow detection of Re in strip solutions at levels of L D = 0.1 ng Re/g-Mo and L Q = 0.3 ng Re/g-Mo in Mo powder samples. By employing these two extraction methods, excellent recoveries of Re from bulk Mo are achieved, and the L Q is improved by a factor of 5 x 10 4 .

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

Ultrafast charge transfer coupled to quantum proton motion at molecule/metal oxide interface

Understanding how the nuclear quantum effects (NQEs) in the hydrogen bond (H-bond) network influence the photoexcited charge transfer at semiconductor/molecule interface is a challenging problem. By combining two kinds of emerging molecular dynamics methods at the ab initio level, the path integral–based molecular dynamics and time-dependent nonadiabatic molecular dynamics, and choosing CH 3 OH/TiO 2 as a prototypical system to study, we find that the quantum proton motion in the H-bond network is strongly coupled with the ultrafast photoexcited charge dynamics at the interface. The hole trapping ability of the adsorbed methanol molecule is notably enhanced by the NQEs, and thus, it behaves as a hole scavenger on titanium dioxide. The critical role of the H-bond network is confirmed by in situ scanning tunneling microscope measurements with ultraviolet light illumination. It is concluded the quantum proton motion in the H-bond network plays a critical role in influencing the energy conversion efficiency based on photoexcitation.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Machine Learning-Based Identification of the Interface Regions for Coupling Local and Nonlocal Models

Local-nonlocal coupling approaches provide a means to combine the computational efficiency of local models and the accuracy of nonlocal models. However, the coupling process can be challenging, requiring expertise to identify the interface between local and nonlocal regions. Here, this study introduces a machine learning-based approach to automatically detect the regions in which the local and nonlocal models should be used. The method uses loading functions evaluated at grid points to decide the model selection at those points. Training of the networks is based on datasets provided by classes of loading functions for which reference coupling configurations are computed using accurate coupled solutions, where accuracy is measured in terms of the relative error between the solution to the coupling approach and the solution to the nonlocal model. We study two approaches that vary in data structure. The first, the full-domain input data approach, uses the entire load vector and outputs a complete label vector, performing a global classification. The second, a window-based approach, processes loads into windows and addresses the problem as a node-wise classification where each window's central point is classified individually. The classification problems are solved via deep learning algorithms based on convolutional neural networks. The performance of these approaches is studied on one-dimensional numerical examples using F1-scores and accuracy metrics. Notably, the windowing approach achieves an accuracy of 0.96 and an F1-score of 0.97, highlighting its potential to automate coupling processes effectively and enhance computational efficiency in material science applications.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Domain Knowledge Guided Bayesian Optimization For Autonomous Alignment Of Complex Scientific Instruments

Bayesian Optimization (BO) is a powerful tool for optimizing complex non-linear systems. However, its performance degrades in high-dimensional problems with tightly coupled parameters and highly asymmetric objective landscapes, where rewards are sparse. In such needle-in-a-haystack scenarios, even advanced methods like trust-region BO (TurBO) often lead to unsatisfactory results. We propose a domain knowledge guided Bayesian Optimization approach, which leverages physical insight to fundamentally simplify the search problem by transforming coordinates to decouple input features and align the active subspaces with the primary search axes. We demonstrate this approach's efficacy on a challenging 12-dimensional, 6-crystal Split-and-Delay optical system, where conventional approaches, including standard BO, TuRBO and multi-objective BO, consistently led to unsatisfactory results. When combined with an reverse annealing exploration strategy, this approach reliably converges to the global optimum. The coordinate transformation itself is the key to this success, significantly accelerating the search by aligning input co-ordinate axes with the problem's active subspaces. As increasingly complex scientific instruments, from large telescopes to new spectrometers at X-ray Free Electron Lasers are deployed, the demand for robust high-dimensional optimization grows. Our results demonstrate a generalizable paradigm: leveraging physical insight to transform high-dimensional, coupled optimization problems into simpler representations can enable rapid and robust automated tuning for consistent high performance while still retaining current optimization algorithms.

FOS: Computer and information sciences↗
DOE OSTI2021

Network-Cognizant Time-Coupled Aggregate Flexibility of Distribution Systems Under Uncertainties

Increasing integration of distributed energy resources (DERs) within distribution feeders provides unprecedented flexibility at the distribution-transmission interconnection. To exploit this flexibility and to use the capacity potential of aggregate DERs, feasible substation power injection trajectories need to be efficiently characterized. This paper provides an ellipsoidal inner approximation of the set of feasible power injection trajectories at the substation such that for any point in the set, there exists a feasible disaggregation strategy of DERs for any load uncertainty realization. The problem is formulated as one of finding the robust maximum volume ellipsoid inside the flexibility region under uncertainty. Though the problem is NP-hard even in the deterministic case, this paper derives novel approximations of the resulting adaptive robust optimization problem based on optimal second-stage policies. The proposed approach yields less conservative flexibility characterization than existing flexibility region approximation formulations. The efficacy of the proposed method is demonstrated on a realistic distribution feeder.

aggregates↗
DOE OSTI2022

A reaction–diffusion model for grayscale digital light processing 3D printing

We report that digital light processing (DLP) 3D printing is an additive manufacturing process that utilizes light patterns to photopolymerize a liquid resin into a solid. Due to the accuracy of modern digital micromirror devices (DMD) and recent advances in resin chemistry, it is now possible to create functionally graded structures using different light intensity values, also known as grayscale DLP (g-DLP). Different intensities of light lead to differences in the polymer crosslinking density after curing, which ultimately produces a part with gradients of material properties. However, g-DLP is a complicated process. First, the DLP printing is a highly coupled chemical and physical process that involves light propagation, chemical reactions, species diffusion, heat transfer, volume shrinkage, and changes in mechanical behaviors of the curing resin. Second, in g-DLP, light gradients create strong in plane gradients of chemical species concentrations in the curing liquid resin due to the strong dependence of light intensity on the rate of monomer crosslinking. Furthermore, light gradients through the depth create concentration gradients due to the degree of cure dependent light absorption and the use of photoabsorbers. These complex physical features of the printing process must be understood in order to properly control printing parameters such as light exposure time, printing speed, and grayscale variations to achieve accurate mechanical properties. In this paper, a photopolymerization reaction–diffusion model is developed and used in conjunction with experiments to investigate the coupled effects of light propagation, chemical reaction rates, and species diffusion during g-DLP 3D printing. The model is implemented numerically utilizing the finite difference method and simulation results are compared to experimental findings of simple printed structures. The agreement between experimental and model predictions of simple quantities of interest, such as geometric feature sizes, shows that the model can capture the overcure due to free-radical and other species diffusion during printing when grayscale patterns are employed. This model lays the groundwork for future extensions that can incorporate more complex coupled physics such as heat transfer, volume shrinkage, and material property evolution, which are critically important in utilizing g-DLP 3D printing for the fabrication of high-performance parts which excellent geometric and material property tolerances.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Standard Analytical Methods for Pyrolysis Bio-Oils

There has been significant recent interest in the production of renewable fuels and chemicals from biomass and waste feedstocks. Pyrolysis pathways produce a liquid bio-oil product, which must be processed further, or upgraded, to yield fuel or chemical products. Bio-oils are very complex and often unstable samples, and research and development on upgrading processes needs reliable analytical information. In particular, chemical characterization techniques are needed to quantify both functional groups and individual compounds present in bio-oils. Reliable analytics are also needed to enable the bioenergy industry, as industrial facilities often have different analytical needs and capabilities than research facilities. In this presentation, we will discuss the development of a suite of standard analytical methods for pyrolysis bio-oils. Analytical methods to be discussed include: Determination of Carbon, Hydrogen, Nitrogen, and Oxygen in bio-oils; Accelerated Aging of Fast Pyrolysis Bio-oil using Carbonyl Titration; Determination of Water Content in Bio-oils by Volumetric Karl Fischer Titration; Determination of Carbon Functional Groups; Elemental Analysis of Bio-oils by Inductively Coupled Plasma Optical Emission Spectroscopy (ICP-OES) - Na, K, Mg, Ca, S, P, and Fe; Determination of Phenolic Groups in Bio-oils using Revised Folin-Ciocalteu Methods: Single Cuvette and Plate Reader; Corrosivity of Bio-oils: Screening Test using Metal Leaching; Determination of Biogenic Content by 14C Measurement using Liquid Scintillation Counter. These new analytical methods are publicly available as Laboratory Analytical Procedures (https://www.nrel.gov/bioenergy/bio-oil-analysis.html), along with previously developed standard methods: GC-MS, Acid Titration, Carbonyl Titration, and 31P NMR. Additionally, the development of diffusion ordered NMR for characterization of bio-oil molecular weight will be discussed. Collectively, this suite of analytical methods represents the most comprehensive set of standard methods available for pyrolysis bio-oils. These standard methods are commonly used by the bioenergy community, and provide reliable information that enables research, scaleup, and industrial processing of biomass to produce renewable fuels and chemicals.

analytical↗
DOE OSTI2022

Charge ordering in Ir dimers in the ground state of Ba 5 AlIr 2 O 11

It has been well established experimentally that the interplay of electronic correlations and spin-orbit interactions in Ir 4+ and Ir 5+ oxides results in insulating J eff = 1/2 and J eff = 0 ground states, respectively. However, in compounds where the structural dimerization of iridium ions is favorable, the direct Ir d-d hybridization can be significant and takes a key role. Here, in this study, we investigate the effects of direct Ir d-d hybridization in comparison with electronic correlations and spin-orbit coupling in Ba 5 AlIr 2 O 11 , a compound with Ir dimers. Using a combination of ab initio many-body wave-function quantum chemistry calculations and resonant inelastic x-ray scattering experiments, we elucidate the electronic structure of Ba 5 AlIr 2 O 11 . We find excellent agreement between the calculated and the measured spin-orbit excitations. Contrary to expectations, the analysis of the many-body wave function shows that the two Ir (Ir 4+ and Ir 5+ ) ions in the Ir 2 O 9 dimer unit in this compound preserve their local J eff character close to 1/2 and 0, respectively. The local point group symmetry at each of the Ir ions plays an important role, significantly limiting the direct d-d hybridization. Our results emphasize that minute details in the local crystal field environment can lead to dramatic differences in the electronic states in iridates and 5d oxides in general.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

DeepM&Mnet: Inferring the electroconvection multiphysics fields based on operator approximation by neural networks

Electroconvection is a multiphysics problem involving coupling of the flow field with the electric field as well as the cation and anion concentration fields. Here, we use electroconvection as a benchmark problem to put forward a new data assimilation framework, the DeepM&Mnet, for simulating multiphysics and multiscale problems at speeds much faster than standard numerical methods using pre-trained neural networks. We first pre-train DeepONets that can predict independently each field, given general inputs from the rest of the fields of the coupled system. DeepONets can approximate nonlinear operators and are composed of two sub-networks, a branch net for the input fields and a trunk net for the locations of the output field. DeepONets, which are extremely fast, are used as building blocks in the DeepM&Mnet and form constraints for the multiphysics solution along with some sparse available measurements of any of the fields. We demonstrate the new methodology and document the accuracy of each individual DeepONet, and subsequently we present two different DeepM&Mnet architectures that infer accurately and efficiently 2D electroconvection fields for unseen electric potentials. Furthermore, the DeepM&Mnet framework is general and can be applied for building any complex multiphysics and multiscale models based on very few measurements using pre-trained DeepONets in a “plug-and-play” mode.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Distributed ADMM Using Private Blockchain for Power Flow Optimization in Distribution Network With Coupled and Mixed-Integer Constraints

The optimization problem for scheduling distributed energy resources (DERs) and battery energy storage systems (BESS) integrated with the power grid is important to minimize energy consumption from conventional sources in response to demand. Conventionally this optimization problem is solved in a centralized manner, limiting the size of the problem that can be solved and creating a high communication overhead because all the data is transferred to the central controller. These limitations are addressed by the proposed distributed consensus-based alternating direction method of multiplier (DC-ADMM) optimization algorithm, which decomposes the optimization problem into subproblems with private cost function and constraints. The distribution feeder is partitioned into low coupling subnetworks/regions, which solves the private subproblem locally and exchanges information with the neighboring regions to reach consensus. The relaxation strategy is employed for mixed-integer and coupled constraints introduced in the optimal power flow (OPF) problem by stationary and transportable BESS because DC-ADMM convergence is only guaranteed for strict convex problems. The information exchange and synchronization between subnetworks/regions are vital for distributed optimization. In this work, both of these aspects are addressed by the blockchain. The smart contract deployed on the blockchain network acts as a mediator for secure data exchange and synchronization in distributed computation. The blockchain-based distributed optimization problem’s effectiveness is tested for a 0.5-MW laboratory microgrid for one hour ahead and day-ahead for the IEEE 123-bus and EPRI J1 test feeders, and results are compared with a centralized solution.

25 ENERGY STORAGE↗