Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “distributed simulation”

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 127 records · Page 7

PythonFOAM: In-situ data analyses with OpenFOAM and Python

Here, we outline the development of a general-purpose Python-based data analysis tool for OpenFOAM. Our implementation relies on the construction of OpenFOAM applications that have bindings to data analysis libraries in Python. Double precision data in OpenFOAM is cast to a NumPy array using the NumPy C-API and Python modules may then be used for arbitrary data analysis and manipulation on flow-field information. We highlight how the proposed wrapper may be used for an in-situ online singular value decomposition (SVD) implemented in Python and accessed from the OpenFOAM solver PimpleFOAM. Here, 'in-situ' refers to a programming paradigm that allows for a concurrent computation of the data analysis on the same computational resources utilized for the partial differential equation solver. In addition, to demonstrate parallel deployments, we deploy a distributed SVD, which collects snapshot data across the ranks of a distributed simulation to compute the global left singular vectors. Crucially, both OpenFOAM and Python share the same message passing interface (MPI) communicator for this deployment which allows Python objects and functions to exchange NumPy arrays across ranks. Subsequently, we provide scaling assessments of this distributed SVD on multiple nodes of Intel Broadwell and KNL architectures for canonical test cases such as the large eddy simulations of a backward facing step and a channel flow at friction Reynolds number of 395. Finally, we demonstrate the deployment of a deep neural network for compressing the flow-field information using an autoencoder to demonstrate an ability to use state-of-the-art machine learning tools in the Python ecosystem.

97 MATHEMATICS AND COMPUTING↗

SU-E-J-175: Comparison of the Treatment Reproducibility of Tumors Affected by Breathing Motion

Purpose: The aim of the dose distribution simulations was to form a global idea of intensity-modulated radiation therapy (IMRT) realization, by its comparison to three-dimensional conformal radiation therapy (3DCRT) delivery for tumors affected by respiratory motion. Methods: In the group of 10patients both 3DCRT and IMRT plans were prepared.For each field the motion kernel was generated with the largest movement amplitude of 4;6 and 8mm.Additionally,the sets of reference measurements were made in no motion conditions(0 mm).The evaluation of plan delivery,using a diode array placed on moving platform,was based on the Gamma Index analysis with distance to agreement of 3mm and dose difference of 3%. Results: IMRT plans tended to spare doses delivered to lungs compared to 3DCRT.Nonetheless,analyzed volumes showed no significant difference between the static and dynamic techniques,except for the volumes of both lungs receiving 10 and 15Gy.After adding the components associated with the respiratory movement,all IMRT lung parameters evaluated for the ipsilateral,contralateral and both lungs together,revealed considerable differences between the 0vs.6, 0vs.8 and 4vs.8-mm amplitudes.Similar results were obtained for the 3DCRT lung measurements,but without significance between the 0vs.6-mm amplitude.Taking into account the CTV score parameter in 3DCRT and IMRT plans,there was no statistically significant difference between the motion patterns with the smallest amplitudes.The differences were found for the 8-mm amplitude when it was compared both with static conditions and 4-mm amplitude (for 3DCRT) and between 0vs.6, 0vs.8 and 4vs.8-mm amplitudes (for IMRT).All accepted and measured 3DCRT and IMRT doses to spinal cord,esophagus and heart were always below the QUANTEC limits. Conclusion: The application of IMRT technique in lung radiotherapy affords possibilities for reducing the lung doses.For maximal amplitudes of breathing trajectory below 4mm,the disagreement between CTV planned and delivered dose distribution can be neglected.Whereas respiratory motion with an amplitude above 5mm leads to significant changes in IMRT delivery.

60 APPLIED LIFE SCIENCES↗

Technical analysis of prospective photovoltaic systems in Utah

This report explores the technical feasibility of prospective utility-scale photovoltaic system (PV) deployments in Utah. Sandia National Laboratories worked with Rocky Mountain Power (RMP), a division of PacifiCorp operating in Utah, to evaluate prospective 2-megawatt (MW) PV plants in different locations with respect to energy production and possible impact on the RMP system and customers. The study focused on 2-MW{sub AC} nameplate PV systems of different PV technologies and different tracking configurations. Technical feasibility was evaluated at three different potential locations in the RMP distribution system. An advanced distribution simulation tool was used to conduct detailed time-series analysis on each feeder and provide results on the impacts on voltage, demand, voltage regulation equipment operations, and flicker. Annual energy performance was estimated.

42 ENGINEERING↗

High Performance Interactive System Dynamics Visualization

This brochure describes a system dynamics simulation (SD) framework that supports an end-to-end analysis workflow that is optimized for deployment on ESIF facilities(Peregrine and the Insight Center). It includes (I) parallel and distributed simulation of SD models, (ii) real-time 3D visualization of running simulations, and (iii) comprehensive database-oriented persistence of simulation metadata, inputs, and outputs.

97 MATHEMATICS AND COMPUTING↗

High Performance Interactive System Dynamics Visualization

This presentation describes a system dynamics simulation (SD) framework that supports an end-to-end analysis workflow that is optimized for deployment on ESIF facilities(Peregrine and the Insight Center). It includes (I) parallel and distributed simulation of SD models, (ii) real-time 3D visualization of running simulations, and (iii) comprehensive database-oriented persistence of simulation metadata, inputs, and outputs.

97 MATHEMATICS AND COMPUTING↗

A Sparse Distributed Gigascale Resolution Material Point Method

In this paper, we present a four-layer distributed simulation system and its adaptation to the Material Point Method (MPM). The system is built upon a performance portable C++ programming model targeting major High-Performance-Computing (HPC) platforms. A key ingredient of our system is a hierarchical block-tile-cell sparse grid data structure that is distributable to an arbitrary number of Message Passing Interface (MPI) ranks. We additionally propose strategies for efficient dynamic load balance optimization to maximize the efficiency of MPI tasks. Our simulation pipeline can easily switch among backend programming models, including OpenMP and CUDA, and can be effortlessly dispatched onto supercomputers and the cloud. Finally, we construct benchmark experiments and ablation studies on supercomputers and consumer workstations in a local network to evaluate the scalability and load balancing criteria. We demonstrate massively parallel, highly scalable, and gigascale resolution MPM simulations of up to 1.01 billion particles for less than 323.25 seconds per frame with 8 OpenSSH-connected workstations.

97 MATHEMATICS AND COMPUTING↗

A New Capability of E4D For 3D Parallel Joint Inversion of DC Resistivity And Traveltime Data on Unstructured Mesh

A major challenge in interpreting geophysical data is how to derive consistent three-dimensional (3D) earth models of different physical properties from spatially and temporally limited measurements. Joint inversion with cross-gradient constraints is an approach to find such models by imposing structural similarities between different physical parameters. We have developed a parallel distributed-memory joint inversion code for direct-current (DC) resistivity and traveltime data using the cross-gradient constraint on unstructured mesh. The code utilizes existing E4D framework for parallel forward simulation, distributed storage and computation of the Jacobian matrix of forward operator, and parallel execution of matrix-vector multiplication during inversion. Besides, the joint inversion is solved by nonlinear conjugate gradient algorithm parallelized for DC resistivity and traveltime data. The joint inversion capability of E4D was tested using synthetic data from cross-borehole DC resistivity and traveltime data. The results indicate that the shape and size of the anomalies from the joint inversion are more reliable than those from separate inversions.

58 GEOSCIENCES↗

RETIRED A STARS: THE EFFECT OF STELLAR EVOLUTION ON THE MASS ESTIMATES OF SUBGIANTS

Doppler surveys have shown that the occurrence rate of Jupiter-mass planets appears to increase as a function of stellar mass. However, this result depends on the ability to accurately measure the masses of evolved stars. Recently, Lloyd called into question the masses of subgiant stars targeted by Doppler surveys. Lloyd argues that very few observable subgiants have masses greater than 1.5 M {sub Sun }, and that most of them have masses in the range 1.0-1.2 M {sub Sun }. To investigate this claim, we use Galactic stellar population models to generate an all-sky distribution of stars. We incorporate the effects that make massive subgiants less numerous, such as the initial mass function and differences in stellar evolution timescales. We find that these effects lead to negligibly small systematic errors in stellar mass estimates, in contrast to the Almost-Equal-To 50% errors predicted by Lloyd. Additionally, our simulated target sample does in fact include a significant fraction of stars with masses greater than 1.5 M {sub Sun }, primarily because the inclusion of an apparent magnitude limit results in a Malmquist-like bias toward more massive stars, in contrast to the volume-limited simulations of Lloyd. The magnitude limit shifts the mean of our simulated distribution toward higher masses and results in a relatively smaller number of evolved stars with masses in the range 1.0-1.2 M {sub Sun }. We conclude that, within the context of our present-day understanding of stellar structure and evolution, many of the subgiants observed in Doppler surveys are indeed as massive as main-sequence A stars.

79 ASTRONOMY AND ASTROPHYSICS↗

Development of a Laser-Electron-Interaction Region for Compact High Repetition-Rate Mono-Energetic Laser Based Gamma Sources

Inverse-Compton scattering (ICS) is a well-established technique capable of producing monochromatic radiation over a large range of the electromagnetic spectrum. Most recently, ICS has attracted attention, owing to its compactness and possible coupling with novel acceleration methods to yield compact X- and 𝛄-ray sources. Under this award, we carried out a systematic optimization of the ICS and demonstrate via numerical simulation that combining ICS with a conventional superconducting linac could produce Watt-level average-power 𝛄 rays with high spectral brightness. The study is performed with idealized laser and electron-beam distributions and confirmed with the realistic distributions simulated from start-to-end numerical modeling of the electron accelerator. The studies led to the design and engineering design of an ICS interaction point for the FAST facility along with some beam measurement. Due to the cancellation of the program in Year 3 long with delay in the commissioning of the FAST facility, we were not able to demonstrate the generation of 𝛄 rays, we, however, procured some of the hardware, installed some of the required infrastructures and plan to eventually generate 𝛄 ray pending future interest on the topics and availability of funding.

43 PARTICLE ACCELERATORS↗

Chapter 16 - Simulating Discrete-Event Models in the Classic Worldviews: A New Approach to Simultaneous Events and Parallel Execution

This chapter describes a new approach to the activity scanning worldview of discrete-event simulation. With this new approach it becomes possible to realize activity scanning simulations with a parallel computer while resolving simultaneous events in a repeatable manner. The motivation for this approach is the two-phase simulation procedure intrinsic to cellular automata, coupled difference equations, and the discrete-event system specification. This new approach is introduced via an intuitively appealing construction of a new simulation procedure for activity scanning models, and this introduction is followed by an embedding of the new approach into the discrete-event system specification. In doing so, we illustrate a practical approach to distributed simulation and demonstrate how the proposed approach can be applied to the event scheduling and process-oriented worldviews.

Nutaro, James J.↗

Computational Fluid Dynamics (CFD) Simulations of Taylor Bubbles in Vertical and Inclined Pipes with Upward and Downward Liquid Flow

Summary Two-phase flow is a common occurrence in pipes of oil and gas developments. Current predictive tools are based on the mechanistic two-fluid model, which requires the use of closure relations to predict integral flow parameters such as liquid holdup (or void fraction) and pressure gradient. However, these closure relations carry the highest uncertainties in the model. In particular, significant discrepancies have been found between experimental data and closure relations for the Taylor bubble velocity in slug flow, which has been determined to strongly affect the mechanistic model predictions (Lizarraga-García 2016). In this work, we study the behavior of Taylor bubbles in vertical and inclined pipes with upward and downward flow using a validated 3D computational fluid dynamics (CFD) approach with level set method implemented in a commercial code. A total of 56 cases are simulated, covering a wide range of fluid properties, pipe diameters, and inclination angles: Eo ∈ [10, 700]; Mo ∈ [1×10–6, 5×103]; ReSL ∈ [–40, 10]; θ ∈ [5°, 90°]. For bubbles in vertical upward flows, the simulated distribution parameter, C0, is successfully compared with an existing model. However, the C0 values of downward and inclined slug flows where the bubble becomes asymmetric are shown to be significantly different from their respective vertical upward flow values, and no current model exists for the fluids simulated here. The main contributions of this work are (1) the relatively large 3D numerical database generated for this type of flow, (2) the study of the asymmetric nature of inclined and some vertical downward slug flows, and (3) the analysis of its impact on the distribution parameter, C0.

Engineering↗

A modeling study of coastal inundation induced by storm surge, sea-level rise, and subsidence in the Gulf of Mexico

The northern coasts of the Gulf of Mexico are highly vulnerable to the direct threats of climate change, such as hurricane-induced storm surge, and such risks can be potentially exacerbated by land subsidence and global sea level rise. This paper presents an application of a coastal storm surge model to study the coastal inundation process induced by tide and storm surge, and its response to the effects of land subsidence and sea level rise in the northern Gulf coast. An unstructured-grid Finite Volume Coastal Ocean Model was used to simulate tides and hurricane-induced storm surges in the Gulf of Mexico. Simulated distributions of co-amplitude and co-phase of semi-diurnal and diurnal tides are in good agreement with previous modeling studies. The storm surges induced by four historical hurricanes (Rita, Katrina, Ivan and Dolly) were simulated and compared to observed water levels at National Oceanic and Atmospheric Administration tide stations. Effects of coastal subsidence and future global sea level rise on coastal inundation in the Louisiana coast were evaluated using a parameter “change of inundation depth” through sensitivity simulations that were based on a projected future subsidence scenario and 1-m global sea level rise by the end of the century. Model results suggested that hurricane-induced storm surge height and coastal inundation could be exacerbated by future global sea level rise and subsidence, and that responses of storm surge and coastal inundation to the effects of sea level rise and subsidence are highly nonlinear and vary on temporal and spatial scales.

Yang, Zhaoqing↗

CyDER – an FMI-based co-simulation platform for distributed energy resources

The increased integration of distributed energy resources (DERs) is bringing a number of challenges to the power grid. These include reverse power flows in distribution systems and potentially transmission systems and grid stability. So far, specialized tools have been developed to capture some of the impact of DERs at the distribution level. However, distribution system operators lack visibility into the overall system conditions. Furthermore, the impact of increasing DERs is not limited to the distribution level but also influences the transmission grid. To support the planning and operation of the grid, we developed a co-simulation platform called CyDER (A Cyber Physical Co-simulation Platform for Distributed Energy Resources in Smart Grids) that integrates various domain-specific simulation tools. CyDER is based on the functional mock-up interface standard. This paper gives an overview of CyDER and demonstrates its use based on two applications.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Regional Analysis of Building Distributed Energy Costs and CO2 Abatement: A U.S. - China Comparison

The following paper conducts a regional analysis of the U.S. and Chinese buildings? potential for adopting Distributed Energy Resources (DER). The expected economics of DER in 2020-2025 is modeled for a commercial and a multi-family residential building in different climate zones. The optimal building energy economic performance is calculated using the Distributed Energy Resources Customer Adoption Model (DER CAM) which minimizes building energy costs for a typical reference year of operation. Several DER such as combined heat and power (CHP) units, photovoltaics, and battery storage are considered. The results indicate DER have economic and environmental competitiveness potential, especially for commercial buildings in hot and cold climates of both countries. In the U.S., the average expected energy cost savings in commercial buildings from DER CAM?s suggested investments is 17percent, while in Chinese buildings is 12percent. The electricity tariffs structure and prices along with the cost of natural gas, represent important factors in determining adoption of DER, more so than climate. High energy pricing spark spreads lead to increased economic attractiveness of DER. The average emissions reduction in commercial buildings is 19percent in the U.S. as a result of significant investments in PV, whereas in China, it is 20percent and driven by investments in CHP. Keywords: Building Modeling and Simulation, Distributed Energy Resources (DER), Energy Efficiency, Combined Heat and Power (CHP), CO2 emissions 1. Introduction The transition from a centralized and fossil-based energy paradigm towards the decentralization of energy supply and distribution has been a major subject of research over the past two decades. Various concerns have brought the traditional model into question; namely its environmental footprint, its structural inflexibility and inefficiency, and more recently, its inability to maintain acceptable reliability of supply. Under such a troubled setting, distributed energy resources (DER) comprising of small, modular, electrical renewable or fossil-based electricity generation units placed at or near the point of energy consumption, has gained much attention as a viable alternative or addition to the current energy system. In 2010, China consumed about 30percent of its primary energy in the buildings sector, leading the country to pay great attention to DER development and its applications in buildings. During the 11th Five Year Plan (FYP), China has implemented 371 renewable energy building demonstration projects, and 210 photovoltaics (PV) building integration projects. At the end of the 12th FYP, China is targeting renewable energy to provide 10percent of total building energy, and to save 30 metric tons of CO2 equivalents (mtce) of energy with building integrated renewables. China is also planning to implement one thousand natural gas-based distributed cogeneration demonstration projects with energy utilization rates over 70percent in the 12th FYP. All these policy targets require significant DER systems development for building applications. China?s fast urbanization makes building energy efficiency a crucial economic issue; however, only limited studies have been done that examine how to design and select suitable building energy technologies in its different regions. In the U.S., buildings consumed 40percent of the total primary energy in 2010 [1] and it is estimated that about 14 billion m2 of floor space of the existing building stock will be remodeled over the next 30 years. Most building?s renovation work has been on building envelope, lighting and HVAC systems. Although interest has emerged, less attention is being paid to DER for buildings. This context has created opportunities for research, development and progressive deployment of DER, due to its potential to combine the production of power and heat (CHP) near the point of consumption and delivering multiple benefits to customers, such as cost

25 ENERGY STORAGE↗

GRIDAPPSD/GOSS-GridAPPS-D

GridAPPS-D provides a method for developers to run their new applications on a real-time simulator with extensive modeling and tool support. GridAPPS-D is built around standard data models like the CIM. It readily interfaces to existing software products, which may also 1) use components of GridAPPS-D and 2) supplement or replace the built-in distribution simulator, facilitating the deployment of new ADMS applications to existing software products.

Allwardt, Craig↗

Systematic study of target localization for bioluminescence tomography guided radiation therapy

Purpose: To overcome the limitation of CT/cone-beam CT (CBCT) in guiding radiation for soft tissue targets, the authors developed a spectrally resolved bioluminescence tomography (BLT) system for the small animal radiation research platform. The authors systematically assessed the performance of the BLT system in terms of target localization and the ability to resolve two neighboring sources in simulations, tissue-mimicking phantom, and in vivo environments. Methods: Multispectral measurements acquired in a single projection were used for the BLT reconstruction. The incomplete variables truncated conjugate gradient algorithm with an iterative permissible region shrinking strategy was employed as the optimization scheme to reconstruct source distributions. Simulation studies were conducted for single spherical sources with sizes from 0.5 to 3 mm radius at depth of 3–12 mm. The same configuration was also applied for the double source simulation with source separations varying from 3 to 9 mm. Experiments were performed in a standalone BLT/CBCT system. Two self-illuminated sources with 3 and 4.7 mm separations placed inside a tissue-mimicking phantom were chosen as the test cases. Live mice implanted with single-source at 6 and 9 mm depth, two sources at 3 and 5 mm separation at depth of 5 mm, or three sources in the abdomen were also used to illustrate the localization capability of the BLT system for multiple targets in vivo. Results: For simulation study, approximate 1 mm accuracy can be achieved at localizing center of mass (CoM) for single-source and grouped CoM for double source cases. For the case of 1.5 mm radius source, a common tumor size used in preclinical study, their simulation shows that for all the source separations considered, except for the 3 mm separation at 9 and 12 mm depth, the two neighboring sources can be resolved at depths from 3 to 12 mm. Phantom experiments illustrated that 2D bioluminescence imaging failed to distinguish two sources, but BLT can provide 3D source localization with approximately 1 mm accuracy. The in vivo results are encouraging that 1 and 1.7 mm accuracy can be attained for the single-source case at 6 and 9 mm depth, respectively. For the 2 sources in vivo study, both sources can be distinguished at 3 and 5 mm separations, and approximately 1 mm localization accuracy can also be achieved. Conclusions: This study demonstrated that their multispectral BLT/CBCT system could be potentially applied to localize and resolve multiple sources at wide range of source sizes, depths, and separations. The average accuracy of localizing CoM for single-source and grouped CoM for double sources is approximately 1 mm except deep-seated target. The information provided in this study can be instructive to devise treatment margins for BLT-guided irradiation. These results also suggest that the 3D BLT system could guide radiation for the situation with multiple targets, such as metastatic tumor models.

60 APPLIED LIFE SCIENCES↗

Toward the Development of a Cold Regions Regional-Scale Hydrologic Model, Final Project Report

This project improves meso-scale hydrologic modeling in the boreal forest by: (1) demonstrating the importance of capturing the heterogeneity of the landscape using small scale datasets for parameterization for both small and large basins; (2) demonstrating that in drier parts of the landscape and as the boreal forest dries with climate change, modeling approaches must consider the sensitivity of simulations to soil hydraulic parameters - such as residual water content - that are usually held constant. Thus, variability / flexibility in residual water content must be considered for accurate simulation of hydrologic processes in the boreal forest; (3) demonstrating that assessing climate change impacts on boreal forest hydrology through multiple model integration must account for direct effects of climate change (temperature and precipitation), and indirect effects from climate impacts on landscape characteristics (permafrost and vegetation distribution). Simulations demonstrated that climate change will increase runoff, but will increase ET to a greater extent and result in a drying of the landscape; and (4) vegetation plays a significant role in boreal hydrologic processes in permafrost free areas that have deciduous trees. This landscape type results in a decoupling of ET and precipitation, a tight coupling of ET and temperature, low runoff, and overall soil drying.

54 ENVIRONMENTAL SCIENCES↗