Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulation of distributed computing systems”

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 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Multigrid treatment of implicit continuum diffusion

Implicit treatment of diffusive terms of various differential orders common in continuum mechanics modeling, such as computational fluid dynamics, is investigated with spectral and multigrid algorithms in non-periodic 2D domains. In doubly periodic time dependent problems these terms are handled efficiently by spectral methods, but in non-periodic systems solved with distributed memory parallel computing and 2D domain decomposition, this efficiency is lost for a large number of processors. Here, we built and present here a multigrid algorithm for these types of problems that outperforms a spectral solution employing the highly optimized FFTW library. This solver is suitable for high performance computing and may be able to efficiently treat implicit diffusion of arbitrary order by introducing auxiliary equations of lower order. We test these solvers for fourth and sixth order diffusion with harmonic test functions as well as turbulent 2D magnetohydrodynamic simulations. It is also shown that an anisotropic operator without mixed-derivative terms improves model accuracy and speed, and we examine the impact that the various diffusion operators have on the energy, the enstrophy, and the qualitative aspect of a simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Control over emergence and alignment of transient blisters in thermo-responsive gels using hierarchically patterned substrates

Surfaces with tunable microscale textures are vital in a large variety of technological applications, including heat transfer, antifouling and adhesion. To facilitate such broad-scale use, there is a need to create surfaces that undergo reconfigurable changes in topology and thus, enable switchable functionality. To date, there is a relative dearth of methods for engineering surfaces that can be actuated to change topography over a range of length scales, and hence, form tunable hierarchically structured layers. Combining modeling and experiments, we design a geometrically patterned, thermo-responsive poly (N-isopropylacrylamide) gel film that undergoes controllable hierarchical changes in topology with changes in temperature. At the bottom, the film is covalently bound to a solid, curved substrate; at the top, the film encompasses longitudinal rectangular ridges that are oriented perpendicular to the underlying cylindrical curves. At temperatures below lower critical solution temperature (LCST), the swollen gel exhibits 3D variations in polymer density and thickness defined by the gel's top and bottom topography. As the temperature rises above LCST, the interplay between the upper ridges and lower curves in the gel drives non-uniform, directional solvent transport, the nucleation and propagation of a phase-separated higher-density skin layer, and the resulting pressure buildup within the film. These different, interacting kinetic processes lead to an instability, which produces transient microscopic blisters in the film. Through simulations, we show how tuning the width of the ridges modifies the propagation of a skin layer and creates localized pressure build-up points, which enables control over the emergence, distribution, and alignment of the microscopic blisters. Futhermore, we provide a simple argument to predict the size of such microscopic features. Experiments confirm our predictions and further highlight how our computational model enables the rational design of topographical transitions in these tunable films. The development of actuatable, hierarchically structured films provides new routes for achieving switchable functionality in actuators, drug release systems and adhesives.

42 ENGINEERING↗
DOE OSTI2019

Estimated Fluid Force and Damping Characteristics of a Thin Film Damper Comparison Between Closed-Form Solutions and Numerical Analysis

Analytical solutions of thin film dampers are useful for determining critical speeds and stability of rotor systems. Most thin film dampers in use are of short axial length, and closed-form solutions to the Reynolds equations exist for estimating pressure, forces, and damping for these types of dampers. This article compares the fluid film forces and damping estimated by the short film bearing model form of the Reynolds equations to the calculated forces and damping of a transient computational fluid dynamic simulation. For this comparison, the fluid was assumed to be incompressible, laminar, and isoviscous. The fluid film forces and damping are calculated from integrating the pressure distribution over the surface of the damper due to small amplitude motions about a steady state static off-center circular orbit. In this case, no cavitation is assumed, and the journal has no angular velocity, so direct stiffness cannot be calculated from the closed-form solution. Radial clearance, journal length, and journal eccentricity have a significant effect on fluid force and damping within a thin film damper. Fluid density does not affect fluid force or damping substantially, while fluid viscosity does. Both the closed-form solutions and computational fluid dynamics simulation compare well with each other and reflect these trends.

Cook, Jason↗
DOE OSTI2020

PETSc DMNetwork: A Library for Scalable Network PDE-Based Multiphysics Simulations

We present DMNetwork, a high-level package included in the PETSc library for the simulation of multiphysics phenomena over large-scale networked systems. The library aims at applications that have networked structures such as those in electrical, gas, and water distribution systems. DMNetwork provides data and topology management, parallelization for multiphysics systems over a network, and hierarchical and composable solvers to exploit the problem structure. DMNetwork eases the simulation development cycle by providing the necessary infrastructure through simple abstractions to define and query the network components. This article presents the design of DMNetwork, illustrates its user interface, and demonstrates its ability to solve multiphysics systems, such as an electric circuit, a network of power grid and water subnetworks, and transient hydraulic systems over large networks with more than 2 billion variables on extreme-scale computers using up to 30,000 processors.

extreme-scale computing↗
DOE OSTI2021

Development of a Turbulent Liquid Spray Atomization Model for Diesel Engine Simulations (Final Technical Report)

This project addresses the systematic lack of predictive capabilities by spray models within engine CFD codes. We develop a new modeling approach to predict the breakup of diesel sprays based on recent literature showing that liquid turbulence plays a fundamental role in spray atomization. A new body of quantitative validation data is also developed as a critical element of the project, leveraging the joint capabilities of Georgia Tech’s high-pressure continuous-flow spray chamber and Argonne National Lab’s near-nozzle x-ray diagnostics at the Advanced Photon Source. This project contributes spatially-resolved measurements of drop size distribution within well-characterized diesel injectors, Spray A and D, from the Engine Combustion Network (ECN) to the engine combustion community for the first time. Utilizing this new body of measurements, we validate and demonstrate a new spray model for diesel sprays, termed the KH-Faeth model, that predicts global and local spray characteristic more accurately than the widely adopted and employed KH model. Predicted drop size distributions are seen to predict measured drops sizes both quantitatively and predictively, with accurate response in droplet size distributions over a wide range of ambient density, injection pressure, and injector nozzle size (Spray A and D) without model tuning. The KH-Faeth model can reduce error in the predicted centerline droplet size profile by up to 80% for ECN Spray D simulations when compared to use of the widely employed KH model.

33 ADVANCED PROPULSION SYSTEMS↗
NASA NTRS2003

North Alabama Lightning Mapping Array (LMA): VHF Source Retrieval Algorithm and Error Analyses

Two approaches are used to characterize how accurately the North Alabama Lightning Mapping Array (LMA) is able to locate lightning VHF sources in space and in time. The first method uses a Monte Carlo computer simulation to estimate source retrieval errors. The simulation applies a VHF source retrieval algorithm that was recently developed at the NASA Marshall Space Flight Center (MSFC) and that is similar, but not identical to, the standard New Mexico Tech retrieval algorithm. The second method uses a purely theoretical technique (i.e., chi-squared Curvature Matrix Theory) to estimate retrieval errors. Both methods assume that the LMA system has an overall rms timing error of 50 ns, but all other possible errors (e.g., multiple sources per retrieval attempt) are neglected. The detailed spatial distributions of retrieval errors are provided. Given that the two methods are completely independent of one another, it is shown that they provide remarkably similar results. However, for many source locations, the Curvature Matrix Theory produces larger altitude error estimates than the (more realistic) Monte Carlo simulation.

Koshak, W. J.↗
DOE OSTI2026

Hydrogen Dispersion Modeling for Development of Smart Distributed Monitoring

Studying hydrogen dispersion is crucial for ensuring the safe and effective deployment of hydrogen as an energy carrier. This study presents a comprehensive CFD modeling framework for simulating hydrogen dispersion at a real-world hydrogen production, storage, and utilization facility. Utilizing the Hydrogen Research Facility under the Advanced Research on Integrated Energy Systems (ARIES) at the National Renewable Energy Laboratory's (NREL) Flatirons campus, controlled hydrogen releases at 27 kg-H2/hr were simulated. The model incorporated site-specific atmospheric conditions, including hourly wind speeds and temperatures recorded between 8 AM and 8 PM from October to December 2023. To reduce computational demands, a statistical reduction technique was applied to condense the dataset to 100 representative scenarios, validated by statistical tests for wind speeds and power law coefficients. Simulations were conducted using the Reynolds-Averaged Navier-Stokes equations. Results demonstrated that wind speed substantially influences hydrogen dispersion, with low wind conditions forming concentrated clouds and higher wind speeds stretching the plume. Additionally, clustering analysis informed optimal sensor placement at various elevations with up to 10 sensor locations on each elevation. This framework offers a robust approach for understanding hydrogen behavior in ambient conditions and informing detection strategies.

08 HYDROGEN↗
DOE OSTI2025

End-to-end microgrid protection using distributed data-driven methods

This paper introduces an end-to-end microgrid protection framework that offers real-time system monitoring, fault-related decision making, and circuit breaker control. This is achieved through the design of distributed data-driven techniques based on the support vector machine method, where each relay is responsible for distributed data collection, fault detection, fault localization, and fault isolation. Local communication is established among neighboring relays, fostering cooperative fault localization and isolation. This decentralized design not only reduces the computational and communication requirements but also enables the adaptability of each relay under varying operational dynamics. The proposed end-to-end protection framework was validated using MATLAB/Simulink simulations on a 100% renewable microgrid, achieving an accuracy of 93.1% with response time of 0.0523 s, in protecting against a range of fault scenarios that are characterized by various types, locations, impedances, load conditions, photovoltaic power levels, and microgrid operating modes.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Using kernel-based statistical distance to study the dynamics of charged particle beams in particle-based simulation codes

Measures of discrepancy between probability distributions (statistical distance) are widely used in the fields of artificial intelligence and machine learning. We describe how certain measures of statistical distance can be implemented as numerical diagnostics for simulations involving charged-particle beams. Related measures of statistical dependence are also described. The resulting diagnostics provide sensitive measures of dynamical processes important for beams in nonlinear or high-intensity systems, which are otherwise difficult to characterize. Here, the focus is on kernel-based methods such as maximum mean discrepancy, which have a well-developed mathematical foundation and reasonable computational complexity. Several benchmark problems and examples involving intense beams are discussed. While the focus is on charged-particle beams, these methods may also be applied to other many-body systems such as plasmas or gravitational systems.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Real-Time Hardware-in-The-Loop Implementation of Protection and Self-Healing of Microgrids

The design of a reliable protection scheme for microgrids often requires communication between protective devices and microgrid controllers. The authors have developed such a communication assisted scalable protection scheme with a self healing feature to protect a microgrid with 100% Inverter Based Resources (IBRs). The communication assisted scheme must be validated in real time with cyber physical co-simulation for successful demonstration. In this regard, the paper presents a co-simulation platform between a simulated power system model using RTDS and physical protective devices. The primary protection of the scheme is programmed in SEL 421-7 relay and backup protection is programmed in MATLAB on a generic computer acting as a microgrid controller. The IEC 61850 models are used to communicate between SEL-421-7 relay and RTDS, whereas TCP/IP communication connects the microgrid controller to RTDS. The paper’s focus is to demonstrate the co- simulation platform with communication links established using both protocols and validate the proposed scheme in real-time on the IEEE 123 node distribution feeder. The paper explains the configuration of the IEC 61850 and TCP/IP communications as the interface requires proper hardware and software setup. Furthermore, the real time performance indicates the Hardware In the Loop (HIL) framework as a competent testing environment for the developed protection scheme for microgrids.

42 ENGINEERING↗
DOE OSTI2019

Pulse Pileup Analysis for a Double-Sided Silicon Strip Detector Using Variable Pulse Shapes

On account of pulse pileup, photon counting detectors (PCDs) suffer from count loss and energy distortion when operating in high count rate environments. Here, we researched the pulse pileup of a double-sided silicon strip detector (DSSSD) to evaluate its potential application in a mammography system. We analyzed the pulse pileup using pulses of varied shapes, where the shape of the pulse depends on the location of photon interaction within the detector. To obtain the shaped pulses, first, transient currents for photons interacting at different locations were simulated using a Technology Computer-Aided Design (TCAD) software. Next, the currents were shaped by a CR- RC 2 shaping circuit, calculated using Simulink. After obtaining these pulses, both the different orders of pileup and the energy spectrum were calculated by taking into account the following two factors: 1) spatial distribution of photon interactions within the detector and 2) time interval distribution between successive photons under a given photon flux. We found that for a DSSSD with a thickness of 300 μm, a pitch of 25 μm, and a strip length of 1 cm, under a bias voltage of 50 V, the variable pulse shape model predicts that the fraction free of pileup can be >90% under a photon flux of 3.75 Mcps/mm 2 (million counts per second per square mm). The DSSSD is a promising candidate for digital mammography applications.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Development of a Multi-Robot System for Autonomous Inspection of Nuclear Waste Tank Pits

This paper introduces the overall design plan, development timeline, and preliminary progress of the Autonomous Pit Exploration System project. This project aims to develop an advanced multi-robot system for the efficient inspection of nuclear waste-storage tank pits. The project is structured into three phases: Phase 1 involves data collection and interface definition in collaboration with Hanford Site experts and university partners, focusing on tank riser geometry and hardware solutions. Phase 2 includes the selection of sensors and robot components, detailed mechanical design, and prototyping. Phase 3 integrates all components into a cohesive system managed by a master control package which also incorporates digital twin and surrogate models, and culminates in comprehensive testing and validation at a simulated tank pit at the Idaho National Laboratory. Additionally, the system’s communication design ensures coordinated operation through shared data, power, and control signals. For transportation and deployment, an electric vehicle (EV) is chosen to support the system for a full 10 h shift with better regulatory compliance for field deployment. A telescopic arm design is selected for its simple configuration and superior reach capability and controllability. Preliminary testing utilizes an educational robot to demonstrate the feasibility of splitting computational tasks between edge and cloud computers. Successful simultaneous localization and mapping (SLAM) tasks validate our distributed computing approach. More design considerations are also discussed, including radiation hardness assurance, SLAM performance, software transferability, and digital twinning strategies.

Nuclear waste management↗
NASA NTRS1986

Correlation of QSO absorption lines in universes dominated by cold dark matter

Theoretical predictions for the redshift correlations between QSO absorption-line systems are investigated in the context of 'cold dark matter' cosmological models. Particles in 'particle-mesh' N-body simulations are interpreted as absorbing clouds at epochs corresponding to mean redshifts, z, of 0.0, 1.25, and 3.0. The velocity correlation function for absorbing clouds is found by passing lines-of-sight through the systems and computing velocity differences for those particles which lie close to the lines. It depends strongly on z and Omega but only weakly, if at all, on the number density, diameter or mass of the clouds. Two interpretations are possible: (1) the heavy element absorption systems are associated with galaxies which are an unbiased sample of the mass distribution in an Omega(0) = 0.2 universe or (2) the Lyman-alpha absorbers are an unbiased sample of the mass in an Omega(0) = 1 universe and the heavy-element absorption systems, like galaxies, are more strongly clustered than the mass.

Salmon, J.↗
NASA NTRS1987

Performance of a three-dimensional Navier-Stokes code on CYBER 205 for high-speed juncture flows

A vectorized 3D Navier-Stokes code has been implemented on CYBER 205 for solving the supersonic laminar flow over a swept fin/flat plate junction. The code extends MacCormack's predictor-corrector finite volume scheme to a generalized coordinate system in a locally one dimensional time split fashion. A systematic parametric study is conducted to examine the effect of fin sweep on the computed flow field. Calculated results for the pressure distribution on the flat plate and fin leading edge are compared with the experimental measurements of a right angle blunt fin/flat plate junction. The decrease in the extent of the separated flow region and peak pressure on the fin leading edge, and weakening of the two reversed supersonic zones with increase in fin sweep have been clearly observed in the numerical simulation.

Lakshmanan, B.↗
NASA NTRS2024

Development of a Data Fusion Methodology for Lineload Aerodynamic Databases for a Launch Vehicle during Liftoff and Transition

The need for databases for the distributed loading on launch vehicles during the early portion of flight necessitates the use of expensive computational flows in regimes where wake effects dominate. While also being expensive, this is a regime that computational tools tend to historically have problems simulating accurately. To help tackle this problem, a method of data fusion to combine computational results to wind tunnel derived force and moment data is developed. Using this method, significant reduction in computational costs and increases in confidence of the final product is possible and has been used to generate several databases for the Space Launch System (SLS) at NASA. While the full details of database generation are not part of this work, the crucial method at its core is developed here. Two SLS geometries are used throughout the work to demonstrate the techniques. These are two of the larger geometries and represent both planned crewed missions to the Moon as well as potential cargo missions to deep space. The method uses principal component analysis (PCA) to generate a reduced ordered model (ROM) to help fill in the full parameter space. Other similar techniques are explored, but were not found to have a significant result on the predictions of the ROM. Because the full number of components are kept to generate the model, this lack of difference is expected. This method is then extended to ensure that predicted surfaces match trusted force and moment data derived from wind tunnel testing. This extension is done by setting up a constrained optimization problem in order to minimize the deviation from the surface resolved computational data while still integrating to the desired values. When generating the constrained optimization problem, a weighting factor to balance these competing needs is introduced. The work compares previously introduced weighting terms from similar work to the proposed terms and shows that the previously used terms do not have as desirable behavior in this flow regime. This method is then expanded by developing a technique to incorporate uncertainty quantification into the developed data fusion methodology. This expansion takes a two pronged approach. One examines transferring the uncertainties in the force and moment database and characterizes how those adjustments change the predicted lineloads. The second looks at model form error and looks how rebuilding the model using slightly different data changes the predictions. These two terms are then combined in order to create an uncertainty model that takes both effects into account. The limitations of the proposed methods is then discussed as well as possible techniques to address these shortcomings.

Launch Vehicles↗
DOE OSTI2020

Thermal Overloading Risk Mitigation With a Semi-Analytical Probabilistic Model on Branch Current

A semi-analytical formulation is presented in this paper for the probability computation of branch current in multiphase systems. The developed formula is derived based on the linear power flow model in rectangular coordinates. The system uncertainty injections can be renewable energy resources or loads and are modeled using a Gaussian mixture model (GMM). The developed formula can be used to compute the line current violation probability as well as integrate into optimal power flow problem as chance-constraint relaxation. Here, the proposed formula is first compared with the Matlab embedded numerical integration function to show its performance. Besides, the semi-analytical formula is validated and compared with the Monte Carlo simulation method using the IEEE 123-bus system, EPRI Ckt5, and Ckt7 systems.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Flexible silicon photonic architecture for accelerating distributed deep learning

The increasing size and complexity of deep learning (DL) models have led to the wide adoption of distributed training methods in datacenters (DCs) and high-performance computing (HPC) systems. However, communication among distributed computing units (CUs) has emerged as a major bottleneck in the training process. In this study, we propose Flex-SiPAC, a flexible silicon photonic accelerated compute cluster designed to accelerate multi-tenant distributed DL training workloads. Flex-SiPAC takes a co-design approach that combines a silicon photonic hardware platform with a tailored collective algorithm, optimized to leverage the unique physical properties of the architecture. The hardware platform integrates a novel wavelength-reconfigurable transceiver design and a micro-resonator-based wavelength-reconfigurable switch, enabling the system to achieve flexible bandwidth steering in the wavelength domain. The collective algorithm is designed to support reconfigurable topologies, enabling efficient all-reduce communications that are commonly used in DL training. The feasibility of the Flex-SiPAC architecture is demonstrated through two testbed experiments. First, an optical testbed experiment demonstrates the flexible routing of wavelengths by shuffling an array of input wavelengths using a custom-designed spatial-wavelength selective switch. Second, a four-GPU testbed running two DL workloads shows a 23% improvement in job completion time compared to a similarly sized leaf-spine topology. We further evaluate Flex-SiPAC using large-scale simulations, which show that Flex-SiPAC is able to reduce the communication time by 26% to 29% compared to state-of-the-art compute clusters under representative collective operations.

Wu, Zhenguo (ORCID:0000000322847985)↗