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 487 records · Page 27

WE-AB-204-07: Spatiotemporal Distribution of the FDG PET Tracer in Solid Tumors: Contributions of Diffusion and Convection Mechanisms

Purpose: In this study, a mathematical model is utilized to simulate FDG distribution in tumor tissue. In contrast to conventional compartmental modeling, tracer distributions across space and time are directly linked together (i.e. moving beyond ordinary differential equations (ODEs) to utilizing partial differential equations (PDEs) coupling space and time). The diffusion and convection transport mechanisms are both incorporated to model tracer distribution. We aimed to investigate the contributions of these two mechanisms on FDG distribution for various tumor geometries obtained from PET/CT images. Methods: FDG transport was simulated via a spatiotemporal distribution model (SDM). The model is based on a 5K compartmental model. We model the fact that tracer concentration in the second compartment (extracellular space) is modulated via convection and diffusion. Data from n=45 patients with pancreatic tumors as imaged using clinical FDG PET/CT imaging were analyzed, and geometrical information from the tumors including size, shape, and aspect ratios were classified. Tumors with varying shapes and sizes were assessed in order to investigate the effects of convection and diffusion mechanisms on FDG transport. Numerical methods simulating interstitial flow and solute transport in tissue were utilized. Results: We have shown the convection mechanism to depend on the shape and size of tumors whereas diffusion mechanism is seen to exhibit low dependency on shape and size. Results show that concentration distribution of FDG is relatively similar for the considered tumors; and that the diffusion mechanism of FDG transport significantly dominates the convection mechanism. The Peclet number which shows the ratio of convection to diffusion rates was shown to be of the order of 10−{sup 3} for all considered tumors. Conclusion: We have demonstrated that even though convection leads to varying tracer distribution profiles depending on tumor shape and size, the domination of the diffusion phenomenon prevents these factors from modulating FDG distribution.

60 APPLIED LIFE SCIENCES↗

Evaluating Impacts of Battery Energy Storage System Functionalities on Distribution Systems Using Power Hardware in-the-Loop Simulation: Preprint

In this paper, we present results from a power hardware-in-the-loop (PHIL) simulation that was performed to test and demonstrate the impact of a battery energy storage system (BESS) on a distribution feeder. The PHIL platform includes a simulated distribution grid in a real-time digital simulator (RTDS), AC and DC power amplifiers, and a battery inverter as the device under test. Guidelines of specific procedures to set up a PHIL platform are provided on how to set up a PHIL platform to evaluate the impact of grid integrated systems. Accelerated time-series PHIL simulations of peak shaving and Volt/Var functionalities were performed using a 540 kVA BESS. These experimental results illustrate how PHIL simulations can be used to evaluate the impact of BESS functionalities on the distribution grid prior to installation in the field.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Simulation of the temperature distribution in crystals grown by Czochralski method

Production of perfect crystals, free of residual strain and dislocations and with prescribed dopant concentration, by the Czochralski method is possible only if the complex, interacting phenomena that affect crystal growth in a Cz-puller are fully understood and quantified. Natural and forced convection in the melt, thermocapillary effect and heat transfer in and around the crystal affect its growth rate, the shape of the crystal-melt interface and the temperature gradients in the crystal. The heat transfer problem in the crystal and between the crystal and all other surfaces present in the crystal pulling apparatus are discussed at length. A simulation and computer algorithm are used, based on the following assumptions: (1) only conduction occurs in the crystal (experimentally determined conductivity as a function of temperature is used), (2) melt temperature and the melt-crystal heat transfer coefficient are available (either as constant values or functions of radial position), (3) pseudo-steady state is achieved with respect to temperature gradients, (4) crystal radius is fixed, and (5) both direct and reflected radiation exchange occurs among all surfaces at various temperatures in the crystal puller enclosure.

Dudokovic, M. P.↗

Implementation of Parallel Dynamic Simulation on Shared-Memory vs. Distributed-Memory Environments

Power system dynamic simulation computes the system response to a sequence of large disturbance, such as sudden changes in generation or load, or a network short circuit followed by protective branch switching operation. It consists of a large set of differential and algebraic equations, which is computational intensive and challenging to solve using single-processor based dynamic simulation solution. High-performance computing (HPC) based parallel computing is a very promising technology to speed up the computation and facilitate the simulation process. This paper presents two different parallel implementations of power grid dynamic simulation using Open Multi-processing (OpenMP) on shared-memory platform, and Message Passing Interface (MPI) on distributed-memory clusters, respectively. The difference of the parallel simulation algorithms and architectures of the two HPC technologies are illustrated, and their performances for running parallel dynamic simulation are compared and demonstrated.

dynamic simulation↗

Earthquake detection in a simulated lunar regolith using distributed acoustic sensing

Current models of inner lunar geology have largely been inferred from the seismic experiments and observations performed during the Apollo missions that comprised a relatively small number of seismic instruments. Refining constraints on fundamental lunar relationships such as crust-mantle and mantle-core boundaries in the future will require seismic arrays spanning larger epicentral distances. A promising technology for installing dense seismic arrays with minimal human effort is distributed acoustic sensing (DAS), an approach that allows a single length of fiber optic cable to act as hundreds or thousands of sensors when coupled with a DAS interrogator. While terrestrial uses of DAS technology for seismic monitoring rely on burying the cable to maximize fidelity of seismic signal transmission to the fiber, digging meters of trench to bury optical fiber on lunar or planetary surfaces is logistically infeasible. To evaluate DAS signal attenuation due to surface deployment of cable in lunar regolith, we completed earthquake detection analyses that evaluated the sensitivity of an optic-fiber DAS system to seismic signals at different burial depths. We deployed a single-mode fiber in a 10-m open-bottom wooden box filled with a lunar regolith simulant (LRS) with fiber buried at different depths within the LRS and recorded signals for four regional and local earthquakes. The results were used to identify and evaluate signal attenuation in surface-deployed fiber compared to buried fiber in the LRS. Burial depth responses to active-source signals were also evaluated similar to previous studies characterizing DAS sensitivity of surface-deployed fiber. Atmospheric noise was minimal as the cable was deployed in an indoor environment; however, where observed, atmospheric and anthropogenic noise was filtered out using the same bandpass filtering used to identify earthquake events. We found that signal attenuation of the surface-deployed fiber compared to buried fiber was relatively high in active-source experiments but was not consistently observed in earthquake signals. That burial depth is not highly correlated to attenuation of the observed earthquake signals indicates that in a noise-limited environment, placing DAS-interrogated fiber directly at the regolith surface may be a promising deployment strategy to consider for sensing remote seismic signals during lunar exploration.

58 GEOSCIENCES↗

Multi-band terahertz active device with complementary metamaterial

We describe a multi-band terahertz-active device using a composite structure made of complementary metamaterial and doped silicon that can be dynamically controlled. This special complementary metamaterial exhibits three resonances that produce three pass-bands. The pass-bands can be uniformly manipulated by exploiting the photoinduced characteristics of the doped silicon. Simulations were performed to analyze the magnetic field and surface current distributions. The simulation results agree well with experimental results obtained from terahertz time-domain spectroscopy. Using an 808-nm-wavelength laser beam, a modulation depth of up to 80% was obtained. In numerical simulations, we used a conductivity mode to characterize photoinduction. The development of multi-band terahertz-active devices has many potential applications, for example, in filters, modulators, switches, and sensors.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

LOG-NORMAL DISTRIBUTION OF COSMIC VOIDS IN SIMULATIONS AND MOCKS

Following up on previous studies, we complete here a full analysis of the void size distributions of the Cosmic Void Catalog based on three different simulation and mock catalogs: dark matter (DM), haloes, and galaxies. Based on this analysis, we attempt to answer two questions: Is a three-parameter log-normal distribution a good candidate to satisfy the void size distributions obtained from different types of environments? Is there a direct relation between the shape parameters of the void size distribution and the environmental effects? In an attempt to answer these questions, we find here that all void size distributions of these data samples satisfy the three-parameter log-normal distribution whether the environment is dominated by DM, haloes, or galaxies. In addition, the shape parameters of the three-parameter log-normal void size distribution seem highly affected by environment, particularly existing substructures. Therefore, we show two quantitative relations given by linear equations between the skewness and the maximum tree depth, and between the variance of the void size distribution and the maximum tree depth, directly from the simulated data. In addition to this, we find that the percentage of voids with nonzero central density in the data sets has a critical importance. If the number of voids with nonzero central density reaches ≥3.84% in a simulation/mock sample, then a second population is observed in the void size distributions. This second population emerges as a second peak in the log-normal void size distribution at larger radius.

79 ASTRONOMY AND ASTROPHYSICS↗

Benchmarking Data for the Proposed Signature of Used Fuel Casks

A set of benchmarking measurements to test facets of the proposed extended storage signature was conducted on May 17, 2016. The measurements were designed to test the overall concept of how the proposed signature can be used to identify a used fuel cask based only on the distribution of neutron sources within the cask. To simulate the distribution, 4 Cf-252 sources were chosen and arranged on a 3x3 grid in 3 different patterns and raw neutron totals counts were taken at 6 locations around the grid. This is a very simplified test of the typical geometry studied previously in simulation with simulated used nuclear fuel.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Author Correction: US oil and gas system emissions from nearly one million aerial site measurements

Correction to: Naturehttps://doi.org/10.1038/s41586-024-07117-5 Published online 13 March 2024 In the version of the article initially published, several errors were present and have been corrected in the HTML and PDF versions of the article and Supplementary Information. The main results, conclusions, and our interpretations of the data remain unchanged. See the new Supplementary Information Section S15 for a more detailed description of the errors corrected and the resulting effects on the analysis. Data processing and methods corrections Overflight count correction: We previously used pre-computed source coverage data for some Carbon Mapper campaigns that was computed differently than was required for our analysis. We have re-computed Carbon Mapper source coverage based on flightline polygons and source coordinates. Transition point computation, well sites: The updated version now correctly compares the cumulative emissions distribution of simulated well site emissions with that of aerially detected sources (rather than plumes) when computing the transition point. Transition point computation, midstream: Additionally, the transition point calculation has been corrected to exclude aerially detected midstream emissions below the transition point, which was previously leading to double counting of these emissions. This error was not present for upstream (well site) emissions. Calculation errors Unit error: We corrected a specific unit conversion error affecting well site emissions in the Kairos Fort Worth dataset. Across all datasets, we also correct the conversion factor for converting from standard volume to mass for midstream emissions. Sorting error: We correct code that was applying incorrect sorting when computing correction factors to account for partial detection at well sites. Small typographical corrections were made in Fig. 1b and SI Section S4.1. Data processing and methods corrections Overflight count correction: We previously used pre-computed source coverage data for some Carbon Mapper campaigns that was computed differently than was required for our analysis. We have re-computed Carbon Mapper source coverage based on flightline polygons and source coordinates. Transition point computation, well sites: The updated version now correctly compares the cumulative emissions distribution of simulated well site emissions with that of aerially detected sources (rather than plumes) when computing the transition point. Transition point computation, midstream: Additionally, the transition point calculation has been corrected to exclude aerially detected midstream emissions below the transition point, which was previously leading to double counting of these emissions. This error was not present for upstream (well site) emissions. Calculation errors Unit error: We corrected a specific unit conversion error affecting well site emissions in the Kairos Fort Worth dataset. Across all datasets, we also correct the conversion factor for converting from standard volume to mass for midstream emissions. Sorting error: We correct code that was applying incorrect sorting when computing correction factors to account for partial detection at well sites. Small typographical corrections were made in Fig. 1b and SI Section S4.1. The following practices may help researchers conducting similar analyses avoid making similar errors: 1, Clear, accessible documentation explaining the interpretation of all columns in data input tables and all internal variables within the model, 2, Simple cross-check calculations computed before and after unit conversions.

Sherwin, Evan D↗

Integral functions of electron lateral distribution and their fluctuations in electron-photon cascades

Monte Carlo simulated lateral distribution functions for electrons of EPC developing in lead, at superhigh energies (.1-1 PeV) for depths t or = 60 c.u. delta t=1t. c.u. are presented. The higher moment characteristics, i.e., variation, asymmetry, excess, are presented along with analytical solutions for the same characteristics at fixed observation level calculated to theory approximations A and B by using numerical inversion of the Laplace transformation. The conclusion is made of a complex, usually non-Gaussian shape of the function of the particle number distribution within a circle of given radius at fixed depth.

Ivanenko, I. P.↗

Angular distribution of electrons directly accelerated by an intense tightly focused laser pulse

We report a study of spectral and angular distributions of electrons directly accelerated from an ultrathin nanofoil by a tightly focused, relativistically intense laser pulse. The approach applied is based on a realistic model describing the focusing of radiation by an off-axis parabolic mirror, the field distribution being simulated with the help of Stratton – Chu integrals. We have compared spectral and angular electron distributions for laser pulses having Gaussian transverse and rectangular intensity profiles on the mirror at the same laser pulse energy. It is shown that in the case of a pulse with a rectangular intensity profile, the energy of fast electrons is higher and the emission angles are smaller than those in the case of a pulse with a Gaussian profile. (paper)

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Counter-ion binding and mobility in the presence of hydrophobic polyions – combining molecular dynamics simulations and NMR

Counter-ion binding and mobility in aqueous solutions of partially hydrophobic ionene oligoions is studied here by a combination of all-atomic molecular dynamics (MD) simulations and NMR ({sup 19}F and {sup 81}Br nuclei) measurements. We present results for 12, 12–ionenes in the presence of different halide ions (F{sup −}, Cl{sup −}, Br{sup −} and I{sup −}), as well as their mixtures; the latter allowing us to probe counter-ion selectivity of these oligoions. We consolidate both structural and dynamic information, in particular simulated radial distribution functions and average residence times of counter-ions in the vicinity of ionenes and NMR data in the form of counter-ion chemical shift and self-diffusion coefficients. On one hand, previously reported enthalpy of dilution and mixing measurements show a reverse counter-ion sequence for 12, 12–ionenes with respect to their less hydrophobic 3, 3– and 6, 6– analogues. On the other hand, the current MD and NMR data, reflecting the counter-ion binding tendencies to the ionene chain, give evidence for the same ordering as that observed by MD for 3, 3–ionenes. This is not seen as a contradiction and can be rationalized on the basis of increasing chain hydrophobicity, which has different consequences for enthalpy and ion-binding. The latter is reflecting free energy changes and as such includes both enthalpic and entropic contributions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A Modular Framework for Modeling Hardware Elements in Distributed Engine Control Systems

Progress toward the implementation of distributed engine control in an aerospace application may be accelerated through the development of a hardware-in-the-loop (HIL) system for testing new control architectures and hardware outside of a physical test cell environment. One component required in an HIL simulation system is a high-fidelity model of the control platform: sensors, actuators, and the control law. The control system developed for the Commercial Modular Aero-Propulsion System Simulation 40k (C-MAPSS40k) provides a verifiable baseline for development of a model for simulating a distributed control architecture. This distributed controller model will contain enhanced hardware models, capturing the dynamics of the transducer and the effects of data processing, and a model of the controller network. A multilevel framework is presented that establishes three sets of interfaces in the control platform: communication with the engine (through sensors and actuators), communication between hardware and controller (over a network), and the physical connections within individual pieces of hardware. This introduces modularity at each level of the model, encouraging collaboration in the development and testing of various control schemes or hardware designs. At the hardware level, this modularity is leveraged through the creation of a Simulink(R) library containing blocks for constructing smart transducer models complying with the IEEE 1451 specification. These hardware models were incorporated in a distributed version of the baseline C-MAPSS40k controller and simulations were run to compare the performance of the two models. The overall tracking ability differed only due to quantization effects in the feedback measurements in the distributed controller. Additionally, it was also found that the added complexity of the smart transducer models did not prevent real-time operation of the distributed controller model, a requirement of an HIL system.

numerical simulation↗

A Modular Framework for Modeling Hardware Elements in Distributed Engine Control Systems

Progress toward the implementation of distributed engine control in an aerospace application may be accelerated through the development of a hardware-in-the-loop (HIL) system for testing new control architectures and hardware outside of a physical test cell environment. One component required in an HIL simulation system is a high-fidelity model of the control platform: sensors, actuators, and the control law. The control system developed for the Commercial Modular Aero-Propulsion System Simulation 40k (C-MAPSS40k) provides a verifiable baseline for development of a model for simulating a distributed control architecture. This distributed controller model will contain enhanced hardware models, capturing the dynamics of the transducer and the effects of data processing, and a model of the controller network. A multilevel framework is presented that establishes three sets of interfaces in the control platform: communication with the engine (through sensors and actuators), communication between hardware and controller (over a network), and the physical connections within individual pieces of hardware. This introduces modularity at each level of the model, encouraging collaboration in the development and testing of various control schemes or hardware designs. At the hardware level, this modularity is leveraged through the creation of a SimulinkR library containing blocks for constructing smart transducer models complying with the IEEE 1451 specification. These hardware models were incorporated in a distributed version of the baseline C-MAPSS40k controller and simulations were run to compare the performance of the two models. The overall tracking ability differed only due to quantization effects in the feedback measurements in the distributed controller. Additionally, it was also found that the added complexity of the smart transducer models did not prevent real-time operation of the distributed controller model, a requirement of an HIL system.

propulsion simulation↗