Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Distributed computation”

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,621 records · Page 90
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Adaptive Computing (AC) (Open Source) [SWR-24-106]

The Adaptive Computing (AC) software stack supports goal-based computing, for which a simulation workload is created on the fly, adapting to the results of calculations. Application-specific code defines an objective, which may be to solve an optimization problem or to train a surrogate model with minimal uncertainty. Then, the AC driver decides where in the design parameter space to run simulations to best achieve that objective. This process is iterative and online; as new data is returned from simulations, the AC driver chooses new simulations to run. The AC driver can strategically run simulations on distributed hardware resources (including high performance computing machines, cloud resources, and edge devices) to maximize throughput and obey resource constraints.

Griffin, Kevin [National Laboratory of the Rockies↗
DOE OSTI2018

Electric Power Distribution System Model Simplification Using Segment Substitution

Quasi-static time-series (QSTS) simulation is used to simulate the behavior of distribution systems over long periods of time (typically hours to years). The technique involves repeatedly solving the load-flow problem for a distribution system model and is useful for distributed energy resource (DER) planning. When a QSTS simulation has a small time step and a long duration, the computational burden of the simulation can be a barrier to integration into utility workflows. One way to relieve the computational burden is to simplify the system model. The segment substitution method of simplifying distribution system models introduced in this paper offers model bus reduction of up to 98% with a simplification error as low as 0.2% (0.002 pu voltage). In contrast to existing methods of distribution system model simplification, which rely on topological inspection and linearization, the segment substitution method uses black-box segment data and an assumed simplified topology.

Reiman, Andrew P. (ORCID:0000000302042960)↗
DOE OSTI2021

Standardizing Performance Metrics for Building-Level Electrical Distribution Systems

Building-level electrical distribution systems comprise a myriad of current-carrying equipment, conversion devices, and protection devices that deliver power from the utility or local distributed energy resources to end-use building loads. Electric power has traditionally been generated, transmitted, and distributed in alternating current (AC). However, the last decade has seen a significant increase in the integration of native direct current (DC) equipment that has elevated the importance of DC distribution systems. Numerous studies have comparatively examined the performance of various electrical distribution systems in buildings but have failed to achieve uniform conclusions, primarily because of a lack of consistent and analogous performance evaluation methods. This paper aims to fill this gap by providing a standard set of metrics and measurement boundaries to consistently evaluate the performance of AC, DC, or hybrid AC/DC electrical distribution systems. The efficacy of the proposed approach is evaluated on a representative medium-sized commercial office building model with AC distribution and an equivalent hybrid AC/DC and DC distribution model, wherein the AC distribution model is concluded to be the most efficient. The simulation results show variation in computed metrics with different selected boundaries that verify the effectiveness of the proposed approach in ensuring consistent computation of the performance of building-level electrical distribution systems. This paper provides an initial set of guidelines for building energy system stakeholders to adopt appropriate solutions, thus leading to more efficient energy systems.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2018

Asymptotic Efficiency of New Distribution-Free Tests of Symmetry for Generalized Skew Alternatives

The Bahadur efficiency of new nonparametric tests of symmetry recently proposed by Nikitin and Ahsanullah is calculated. In contrast to this result, where only location alternatives were discussed, in the present paper generalized skew alternatives are of interest. It is shown that the new tests are highly efficient for a large class of skew alternatives. The problem of most favorable alternatives is also studied.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2012

Three-dimensional Reconstruction and Homogenization of Heterogeneous Materials Using Statistical Correlation Functions and FEM

In this study, a previously developed reconstruction methodology is extended to three-dimensional reconstruction of a three-phase microstructure, based on two-point correlation functions and two-point cluster functions. The reconstruction process has been implemented based on hybrid stochastic methodology for simulating the virtual microstructure. While different phases of the heterogeneous medium are represented by different cells, growth of these cells is controlled by optimizing parameters such as rotation, shrinkage, translation, distribution and growth rates of the cells. Based on the reconstructed microstructure, finite element method (FEM) was used to compute the effective elastic modulus and effective thermal conductivity. A statistical approach, based on two-point correlation functions, was also used to directly estimate the effective properties of the developed microstructures. Good agreement between the predicted results from FEM analysis and statistical methods was found confirming the efficiency of the statistical methods for prediction of thermo-mechanical properties of three-phase composites.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Nuclide inventory characterization of EBR-II MOX fuel test pins

Characterization of 3 mixed oxide (MOX) fuel pins from the Experimental Breeder Reactor II (EBR-II) SPA-2/-2B experiment program was performed in support of the Advanced Reactor Experiments for Sodium Fast Reactor Fuels (ARES) joint project between Idaho National Laboratory (INL) and the Japan Atomic Energy Agency (JAEA). The characterization efforts were performed to enable transient testing in the Transient Reactor Test (TREAT) facility at INL and will include post irradiation examination (PIE) and other analyses. MOX fuel pins UW02009, 02011, and 06024 were historically irradiated from 1989 to 1994 for a total of 822.42 effective full power days (EFPDs) and peak burnup of approximately 14.3 at.% (atomic percent of heavy metal atoms, ∼134.4 GWd/t). They had since been placed in storage until a future use was identified. Characterization of the fission products and activation of each fuel pin was necessary to enable removal from storage and subsequent irradiation testing or PIE activities. The initial fuel pin geometry and composition were established using historical documentation. Irradiation history details were utilized to compute irradiated and decayed fuel pin masses, atomic densities, and activities. A pellet-wise axial neutron flux was also computed to assist in evaluation of the axial distribution of {sup 235}U and {sup 239}Pu, per pellet, within each of the test fuel pins. The results computed herein supported experiment design for advanced high-burnup MOX fuel designs under slow transient overpower (MOXTOP) conditions, which was Phase I of the ARES project. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2014

Computer modeling of electrical and thermal performance during bipolar pulsed radiofrequency for pain relief

Purpose: Pulsed RF (PRF) is a nonablative technique for treating neuropathic pain. Bipolar PRF application is currently aimed at creating a “strip lesion” to connect the electrode tips; however, the electrical and thermal performance during bipolar PRF is currently unknown. The objective of this paper was to study the temperature and electric field distributions during bipolar PRF. Methods: The authors developed computer models to study temperature and electric field distributions during bipolar PRF and to assess the possible ablative thermal effect caused by the accumulated temperature spikes, along with any possible electroporation effects caused by the electrical field. The authors also modeled the bipolar ablative mode, known as bipolar Continuous Radiofrequency (CRF), in order to compare both techniques. Results: There were important differences between CRF and PRF in terms of electrical and thermal performance. In bipolar CRF: (1) the initial temperature of the tissue impacts on temperature progress and hence on the thermal lesion dimension; and (2) at 37 °C, 6-min of bipolar CRF creates a strip thermal lesion between the electrodes when these are separated by a distance of up to 20 mm. In bipolar PRF: (1) an interelectrode distance shorter than 5 mm produces thermal damage (i.e., ablative effect) in the intervening tissue after 6 min of bipolar RF; and (2) the possible electroporation effect (electric fields higher than 150 kV m{sup −1}) would be exclusively circumscribed to a very small zone of tissue around the electrode tip. Conclusions: The results suggest that (1) the clinical parameters considered to be suitable for bipolar CRF should not necessarily be considered valid for bipolar PRF, and vice versa; and (2) the ablative effect of the CRF mode is mainly due to its much greater level of delivered energy than is the case in PRF, and therefore at same applied energy levels, CRF, and PRF are expected to result in same outcomes in terms of thermal damage zone dimension.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Fast resonance decays in nuclear collisions

Here, we present a new method to calculate the final decay spectrum of direct resonance decays directly from hydrodynamic fields on a freeze-out surface. The method is based on identifying components of the final particle spectrum that transform in an irreducible way under rotations in the fluid-restframe. Corresponding distribution functions can be pre-computed including all resonance decays. Just a few of easily tabulated scalar functions then determine the Lorentz invariant decay spectrum from each space-time point, and simple integrals of these scalar functions over the freeze-out surface determine the final decay products. This by-passes numerically costly event-by-event calculations of the intermediate resonances. The method is of considerable practical use for making realistic data to model comparisons of the identified particle yields and flow harmonics, and for studying the viscous corrections to the freeze-out distribution function.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2016

Investigating Evaporation in Gas Diffusion Layers for Fuel Cells with X-ray Computed Tomography

Understanding evaporation in porous media and the associated water distribution for a given saturation is critical for optimizing many different technologies including polymer-electrolyte fuel cells. In these devices, heat and mass-transport are coupled due to the two-phase transport of water and operating temperatures from subzero to 80°C. Especially critical is understanding phase change in the mixed wettability, carbon gas-diffusion layers (GDLs). While previous works have measured evaporation rates empirically for a given saturation, there remains a need to explore the mechanisms governing evaporation, which are tied directly to the internal water distribution. Here in this paper, liquid-water evaporation rates in GDLs are measured in situ using synchrotron X-ray computed tomography (CT). X-ray CT allows visualizing the evaporating water-front location and interfacial water/air surface area, thereby enabling true surface-area based evaporation rates. It is found that the overall specific evaporation rate is essentially constant as a function of saturation and that the water/air interfacial area scales almost linearly with saturation. To isolate transport and kinetic contributions to the overall evaporation rate, we systematically varied gas flow rate and composition. A three-dimensional mathematical model with direct meshes of liquid-water evaporation fronts from the X-ray CT studies allowed for the determination that the evaporation is transport limited. The overall results provide insight into evaporation phenomena in porous media.

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

Human Dosimetry and Preliminary Tumor Distribution of (superscript)18F-Fluoropaclitaxel in Healthy Volunteers and Newly Diagnosed Breast Cancer Patients Using PET/CT

{sup 18}F-fluoropaclitaxel is a radiolabeled form of paclitaxel, a widely used chemotherapy agent. Preclinical data suggest that {sup 18}F-fluoropaclitaxel may be a reasonable surrogate for measuring the uptake of paclitaxel. As a substrate of P-glycoprotein, a drug efflux pump associated with multidrug resistance, {sup 18}F-fluoropaclitaxel may also be useful in identifying multidrug resistance and predicting tumor response for drugs other than paclitaxel. After informed consent was obtained, 3 healthy volunteers and 3 patients with untreated breast cancer (neoadjuvant chemotherapy candidates, tumor size > 2 cm) received an intravenous infusion of {sup 18}F-fluoropaclitaxel and then underwent PET/CT. Healthy volunteers underwent serial whole-body imaging over an approximately 3-h interval, and organ {sup 18}F residence times were determined from the time-activity curves uncorrected for decay to determine dosimetry. Radiation dose estimates were calculated using OLINDA/EXM software. For breast cancer patients, dynamic imaging of the primary tumor was performed for 60 min, followed by static whole-body scans at 1 and 2 h after injection. Dosimetry calculations showed that the gallbladder received the highest dose (229.50 {mu}Gy/MBq [0.849 rad/mCi]), followed by the small and large intestines (161.26 {mu}Gy/MBq [0.597 rad/mCi] and 184.59 {mu}Gy/MBq [0.683 rad/mCi]). The resultant effective dose was 28.79 {mu}Gy/MBq (0.107 rem/mCi). At approximately 1 h after injection, an average of 42% of the decay-corrected activity was in the gastrointestinal system, with a mean of 0.01% in the tumor. All 3 breast cancer patients showed retention of {sup 18}F-fluoropaclitaxel and ultimately demonstrated a complete pathologic response (no invasive cancer in the breast or axillary nodes) to chemotherapy that included a taxane (either paclitaxel or docetaxel) at surgical resection. The tumor-to-background ratio increased with time to a maximum of 7.7 at 20 min. This study demonstrates the feasibility of using {sup 18}F-fluoropaclitaxel PET/CT tumor imaging and provides radiation dosimetry measurements in humans. Although further study is needed, it is hoped that the measured intratumoral {sup 18}F-fluoropaclitaxel distribution can serve as a surrogate for paclitaxel, and potentially other chemotherapeutic agent retention, in solid tumors.

18F-Fluoropaclitaxel (FPAC)↗
DOE OSTI2014

Microbial Activity and Precipitation at Solution-Solution Mixing Zones in Porous Media (Final Report)

The use of biological and chemical processes that degrade or immobilize contaminants in subsurface environments is a cornerstone of remediation technology. The enhancement of biological and chemical processes in situ, involves the transport, displacement, distribution and mixing of one or more reactive agents. Biological and chemical reactions all require diffusive transport of solutes to reaction sites at the molecular scale and accordingly, the success of processes at the meter-scale and larger is dictated by the success of phenomena that occur at the micron-scale. However, current understanding of scaling effects on the mixing and delivery of nutrients in biogeochemically dynamic porous media systems is limited, despite the limitations this imposes on the efficiency and effectiveness of the remediation challenges at hand. We therefore proposed to experimentally characterize and computationally describe the growth, evolution, and distribution of microbial activity and mineral formation as well as changes in transport processes in porous media that receive two or more reactive amendments. The model system chosen for this project was based on a method for immobilizing 90 Sr, which involves stimulating microbial urea hydrolysis with ensuing mineral precipitation (CaCO 3 ), and co-precipitation of Sr. Studies at different laboratory scales were used to visualize and quantitatively describe the spatial relationships between amendment transport and consumption that stimulate the production of biomass and mineral phases that subsequently modify the permeability and heterogeneity of porous media. Biomass growth, activity, and mass deposition in mixing zones was investigated using two-dimensional micro-model flow cells as well as flow cells that could be analyzed using synchrotron-based x-ray tomography. Larger-scale flow-cell experiments were conducted where the spatial distribution of media properties, flow, segregation of biological activity and impact on ancillary constituents (i.e., Sr) was determined. Model simulations accompanied the experimental efforts. The research contributed towards defining the key physical, chemical, and biological processes influencing the form and mobility of DOE priority contaminants (e.g., 60 Co, 90 Sr, U) in the subsurface. The work conducted and reported herein, will in the future (i) contribute to controlling the juxtaposition of microbial activity, contaminants and amendments, (ii) promote new strategies for delivering amendments, and (iii) allow new approaches for modifying permeability and flow in porous media. We feel that the work has already translated directly to improving the efficiency of amendment based remediation strategies. The results of the project have been published in a number of peer reviewed journal articles. The abstracts and citations to those articles, given in section 2.0 below, make up the bulk of this final report.

42 ENGINEERING↗
DOE OSTI2024

Implementation of the Glued Sphere Discrete Element Model for Non-Spherical Particles in MFiX Software

To enhance solver capabilities, simulation flexibility and model validation within the MFiX software, the U.S. Department of Energy (DOE) is funding efforts to develop and integrate the glued-sphere discrete element method into the latest version of MFiX as a dedicated computational module. The glued-sphere discrete element method is a numerical technique to depict the behavior of non-spherical particles in granular flows or particulate systems by representing them as a collection of component spheres. These spheres are bonded together to approximate the shape and mechanical/chemical properties of a more complex particle. The method effectively reuses the existing sphere-sphere collision algorithm, interphase momentum and heat transfer calculations utilized in the traditional discrete element method, extending these capabilities to non-spherical particles. Additionally, this method explicitly resolves intra-particle temperature and species distributions. The MFiX glued-sphere computational module includes tools for generating glued sphere configurations, a dedicated solver, and visualization capabilities in post-processing. More specifically within the computational module, collision detection and calculations were first performed on component spheres and then mapped onto non-spherical particles. The linear spring-dashpot model was utilized to simulate the sphere-sphere interactions.

Ke, Renjie↗
DOE OSTI2025

Linear Solver for Electromagnetic Simulation of General Distribution Feeders

High-fidelity electromagnetic transient (EMT) modeling is required for accurate simulation and analysis of power system dynamics in modern distribution feeders. However, the high-fidelity of EMT models often leads to significant computational challenges, particularly in terms of computational resources and simulation time. This paper investigates the development and application of a detailed EMT model for general distribution feeders, with a focus on improving computational efficiency. A direct linear solver is proposed for a bordered block diagonal (BBD) matrix structure commonly encountered in a EMT model of distribution feeders. The solver integrates the Schur complement method with the block tridiagonal matrix algorithm to enhance the computational performance. The proposed solver is validated using the primary feeder of the IEEE 342-node test system, demonstrating its accuracy and efficiency in EMT simulations. Furthermore, the solver’s performance is benchmarked against MATLAB’s built-in linear solvers, showing significant improvements in computation time while maintaining high fidelity and accuracy in simulation results.

Choi, Jongchan [ORNL] (ORCID:000000025952455X)↗
DOE OSTI2021

Decentralized Low-Rank State Estimation for Power Distribution Systems

This article considers the low-observability state estimation problem in power distribution networks and develops a decentralized state estimation algorithm leveraging the matrix completion methodology. Matrix completion has been shown to be an effective technique in state estimation that exploits the low dimensionality of the power system measurements to recover missing information. This technique can utilize an approximate (linear) load flow model, or it can be used with no physical models in a network where no information about the topology or line admittance is available. The direct application of matrix completion algorithms requires solving a semi-definite programming (SDP) problem, which becomes computationally challenging for large networks. We therefore develop a decentralized algorithm that capitalizes on the popular proximal alternating direction method of multipliers (proximal ADMM). The method allows us to distribute the computation among different areas of the network, leading to a scalable algorithm. By doing all computations at individual control areas and only communicating with neighboring areas, the algorithm eliminates the need for data to be sent to a central processing unit and thus increases efficiency and contributes to the goal of autonomous control of distribution networks. We illustrate the advantages of the proposed algorithm numerically using standard IEEE test cases.

41 EE - Solar Energy Technologies Office (EE-4S)↗
DOE OSTI2013

Manufacture of gradient micro-structures of magnesium alloys using two stage extrusion dies

This paper aims to manufacture magnesium alloy metals with gradient micro-structures using hot extrusion process. The extrusion die was designed to have a straight channel part combined with a conical part. Materials pushed through this specially-designed die generate a non-uniform velocity distribution at cross sections inside the die and result in different strain and strain rate distributions. Accordingly, a gradient microstructure product can be obtained. Using the finite element analysis, the forming temperature, effective strain, and effective strain rate distributions at the die exit were firstly discussed for various inclination angles in the conical die. Then, hot extrusion experiments with a two stage die were conducted to obtain magnesium alloy products with gradient micro-structures. The effects of the inclination angle on the grain size distribution at cross sections of the products were also discussed. Using a die of an inclination angle of 15°, gradient micro-structures of the grain size decreasing gradually from 17 μm at the center to 4 μm at the edge of product were achieved.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Consideration for the application of artificial neural networks to criticality code validation - 14270

The artificial neural network (ANN) is a consistent example for non-parametric parameter estimation and regression, with powerful universal approximation properties, that is distribution independent. ANN application to criticality code validation centers on aligning mathematical relationships between code validation parameters and neural network prediction for those parameter values. Two specific goals are: (1) an appropriate ANN architecture to determine validation parameter values, and (2) bootstrapping individual ANNs as a 'Prediction Interval' (PI) technique to estimate the Calculation Margin and Upper Subcritical Limit. The individual ANN workflow is segregated into Training, Validation (network), and Testing Sets, where the network predictive properties and performance are expressed in the latter set. The Bootstrap ANN replicates a large set of individual ANNs (taken from the original data distribution) where (1-α/2)100% PI estimates can be determined. Results from the ANN method are compared to a variety of validation examples: a 'toy example' from ANSI/ANS-8.24, Appendix D, a variety of methods noted in a Nuclear Criticality Safety Division 2009 paper, and several validation campaigns for highly enriched uranium and plutonium systems. The ANN method is robust and can be applied to single or multivariate inputs. A trained ANN can provide both an average estimate (e.g., single value) and an estimate over a range of independent variable values where the kUSL can be formed in a regression relationship suitable for extrapolation. The input dependence has been tested, ranging from 15 to over 400 sample size distributions. An expanded Shapiro-Wilk test (for sample sizes up to 5000) was implemented to confirm that the ANN method is transparent to distribution function. The Principle Component Analysis method was applied as an input preprocessing step embedded within the ANN in order to reduce dimensionality and correlation for multivariable inputs. This can potentially aid in determining accurate margins over applying arbitrary margins. The automatic optimization of the Bootstrap ANN PI is left for future consideration. (authors)

42 ENGINEERING↗
DOE OSTI2025

Viewing is understanding: Graphite microstructure effects on infiltrated molten salt distribution revealed by 3D neutron tomography

Molten salt infiltration in the pore network of nuclear graphite may cause unwanted changes to graphite's local structure and mechanical and thermal properties. A detailed and comprehensive understanding of molten salt intrusion (distribution across sample cross section and penetration depth) is needed to assess its effects. Here, in this work, we report on an improved methodology for the use of neutron imaging (computed tomography) to evaluate salt penetration and distribution of a wide range of graphite grades with diverse microstructures. Neutron tomography data were acquired on the same graphite sample before and after salt intrusion; the 3D reconstructed volumes were digitally co-registered and subtracted. The difference in neutron attenuation coefficient represents direct visualization of FLiNaK (LiF–NaF–KF) salt distribution in the salt-impregnated graphite samples. This improved methodology was applied to investigate the effect of exposure times (12 h and 336 h) and of graphite microstructure when exposed to FLiNaK at 750 °C and 3 bar (gauge) pressure, starting from flowing argon at near atmospheric pressure. The results show that medium-grained and fine-grained graphites evolve to equilibrium at significantly different rates: fast salt uptake in medium-grained graphites produces salt deposits throughout the volume of graphite specimens, whereas salt infiltration in fine-grained graphites is much slower and limited to exposed surfaces.

FLiNaK infiltration↗