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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Finite-element analysis of transverse compressive and thermal loads on Nb 3 Sn wires with voids

High-field superconducting magnets play a very important role in many large-scale physics experiments, particularly particle colliders and fusion confinement devices such as Large Hadron Collider (LHC) and International Thermonuclear Experimental Reactor (ITER). The two most common superconductors used in these applications are NbTi and Nb 3 Sn. Nb 3 Sn wires are favored because of their significantly higher J c (critical current density) for higher field applications. The main disadvantage of Nb 3 Sn is that the superconducting performance of the wire is highly strain sensitive and it is very brittle. This strain sensitivity is strongly influenced by two factors: plasticity and cracked filaments. Cracks are induced by large stress concentrators that can be traced to the presence of voids in the wire. We develop detailed 2-D and 3-D finite-element models containing wire filaments and different possible distributions of voids in a bronze-route Nb 3 Sn wire. We apply compressive transverse loads for various cases of void distributions to simulate the stress and strain response of a Nb 3 Sn wire under the Lorentz force. Furthermore, this paper improves our understanding of the effect voids have on the Nb 3 Sn wire's mechanical properties, and in so, the connection between the distribution of voids and performance degradation such as the correlation between irreversible strain limit and the void-induced local stress concentrations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Simulation of local ion transport in lamellar block copolymer electrolytes based on electron micrographs

A method is presented in this paper to relate local morphology and ionic conductivity in a solid, lamellar block copolymer electrolyte for lithium batteries, by simulating conductivity through transmission electron micrographs. The electrolyte consists of polystyrene-block-poly(ethylene oxide) mixed with lithium bis(trifluoromethanesulfonyl) imide salt (SEO/LiTFSI), where the polystyrene phase is structural phase and the poly(ethylene oxide)/LiTFSI phase is ionically conductive. The electric potential distribution is simulated in binarized micrographs by solving the Laplace equation with constant potential boundary conditions. A morphology factor, f, is reported for each image by calculating the effective conductivity relative to a homogenous conductor. Images from two samples are examined, one annealed with large lamellar grains and one unannealed with small grains. The average value off is 0.45 ± 0.04 for the annealed sample, and 0.37 ± 0.03 for the unannealed sample, both close to the value predicted by effective medium theory, 1/2. Simulated conductivities are compared to published experimental conductivities. The value of f Unannealed /f Annealed is 0.82 for simulations and 6.2 for experiments. Simulation results correspond well to predictions by effective medium theory but do not explain the experimental measurements. Finally, observation of nanoscale morphology over length scales greater than the size of the micrographs (~1 μm) may be required to explain the experimental results.

36 MATERIALS SCIENCE↗
DOE OSTI2016

SU-F-R-36: Validating Quantitative Radiomic Texture Features for Oncologic PET: A Digital Phantom Study

Purpose: Radiomic texture features derived from the oncologic PET have recently been brought under intense investigation within the context of patient stratification and treatment outcome prediction in a variety of cancer types; however, their validity has not yet been examined. This work is aimed to validate radiomic PET texture metrics through the use of realistic simulations in the ground truth setting. Methods: Simulation of FDG-PET was conducted by applying the Zubal phantom as an attenuation map to the SimSET software package that employs Monte Carlo techniques to model the physical process of emission imaging. A total of 15 irregularly-shaped lesions featuring heterogeneous activity distribution were simulated. For each simulated lesion, 28 texture features in relation to the intensity histograms (GLIH), grey-level co-occurrence matrices (GLCOM), neighborhood difference matrices (GLNDM), and zone size matrices (GLZSM) were evaluated and compared with their respective values extracted from the ground truth activity map. Results: In reference to the values from the ground truth images, texture parameters appearing on the simulated data varied with a range of 0.73–3026.2% for GLIH-based, 0.02–100.1% for GLCOM-based, 1.11–173.8% for GLNDM-based, and 0.35–66.3% for GLZSM-based. For majority of the examined texture metrics (16/28), their values on the simulated data differed significantly from those from the ground truth images (P-value ranges from <0.0001 to 0.04). Features not exhibiting significant difference comprised of GLIH-based standard deviation, GLCO-based energy and entropy, GLND-based coarseness and contrast, and GLZS-based low gray-level zone emphasis, high gray-level zone emphasis, short zone low gray-level emphasis, long zone low gray-level emphasis, long zone high gray-level emphasis, and zone size nonuniformity. Conclusion: The extent to which PET imaging disturbs texture appearance is feature-dependent and could be substantial. It is thus advised that use of PET texture parameters for predictive and prognostic measurements in oncologic setting awaits further systematic and critical evaluation.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2013

Evaluation of the use of color-set geometry during lattice physics constants generation for boiling water reactor simulation

Current methods for BWR nuclear design and analysis consist of using lattice physics neutron transport methods to generate the two-group homogenized cross-sections that are then used in a nodal diffusion theory code. The lattice transport solutions are performed for a single assembly with reflective boundary conditions, which is a practical approximation. A method is developed to account for assembly exposure distributions (environment) in the core within the lattice transport calculations with the use of color-sets (2x2) geometry. The loading pattern is examined and an appropriate number of characteristic color-set cells are selected for analysis. Treatment of the co-resident exposed fuel within this method is also presented. The calculation process was followed for a recent BWR cycle design with comparisons being performed on both a lattice and core-wide basis to evaluate the proposed method. The lattice based comparisons show noticeable differences in the pin power distribution predictions, which require further investigation to see how this translates into core performance calculations. The core-wide comparisons show minor differences and are generally in a good agreement, which is expected with this small perturbation. A slight improvement was noticed in the reduction of the power distribution uncertainty. However, given the additional amount of work and computer run time increase, further evaluation, especially of core pin power predictions, is needed to consider this method for production level design and safety analysis calculations. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2011

Modeling Impacts of Climate Change Mitigation Technologies on Power Grids

This paper describes a modeling approach that simulate the impacts of different climate change mitigation technologies on power grids for power system planning purposes. Because the historical data is less credible when new technologies are being deployed to the system, it is then critical to model them to address their impacts. This paper illustrated how to integrate modeling results obtained from different modeling tools to give a reasonable forecast of the future. Building simulation tools, distribution power grid modeling tools, and power system planning tools are used to model and aggregate the impacts from the end-use to the system level. Electricity generation, production cost, emission, and transmission congestions are used to quantify the influence of different mitigation technologies. Modeling results have shown that the cross-discipline modeling approach provided the modeler with the necessary time resolution and input details to address the variables that influence the modeling results. Different modeling issues are also addressed in the paper.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2014

New approach for absolute fluence distribution calculations in Monte Carlo simulations of light propagation in turbid media

A novel way to attain three dimensional fluence rate maps from Monte-Carlo simulations of photon propagation is presented in this work. The propagation of light in a turbid medium is described by the radiative transfer equation and formulated in terms of radiance. For many applications, particularly in biomedical optics, the fluence rate is a more useful quantity and directly derived from the radiance by integrating over all directions. Contrary to the usual way which calculates the fluence rate from absorbed photon power, the fluence rate in this work is directly calculated from the photon packet trajectory. The voxel based algorithm works in arbitrary geometries and material distributions. It is shown that the new algorithm is more efficient and also works in materials with a low or even zero absorption coefficient. The capabilities of the new algorithm are demonstrated on a curved layered structure, where a non-scattering, non-absorbing layer is sandwiched between two highly scattering layers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2012

Implications for the electron distribution from the stationary hydrodynamic model of a one-dimensional plasma expansion into vacuum

It is shown that the hydrodynamic model of a one-dimensional collisionless plasma expansion is contained in the kinetic description as a special case. This belongs to a specific choice for the electron distribution function. Moreover, the consequences of the use of the hydrodynamic approach regarding the temporal evolution of the electron phase space density are investigated. It turns out that only the case of a hydrodynamic description with the adiabatic constant {kappa}=3 is physically self-consistent. Numerical simulations confirm this argumentation. The analysis for the case {kappa}=3 is extended to the kinetics of a relativistic electron gas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Measurements of the populations of metastable and resonance levels in the plasma of an RF capacitive discharge in argon

The behavior of the populations of two metastable and two lower resonance levels of argon atoms in the plasma of an RF capacitive discharge was studied. The populations were measured by two methods: the method of emission self-absorption and the method based on measurements of the intensity ratios of spectral lines. It is shown that the populations of resonance levels increase with increasing power deposited in the discharge, whereas the populations of metastable levels is independent of the RF power. The distribution of the populations over energy levels is not equilibrium under these conditions. The population kinetics of argon atomic levels in the discharge plasma is simulated numerically. The distribution function of plasma electrons recovered from the measured populations of atomic levels and numerical simulations is found to be non-Maxwellian.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

RAM pressure stripping of the Large Magellanic Cloud’s disk as a probe of the Milky Way’s circumgalactic medium

Recent observations have constrained the orbit and structure of the Large Magellanic Cloud (LMC), implying a well-constrained pericentric passage about the Milky Way (MW) ∼50 Myr ago. In this scenario, the LMC’s gaseous disk has recently experienced stripping, suggesting the extent of its H{sub I} disk directly probes the medium in which it is moving. From the observed stellar and H{sub I} distributions of the system we find evidence of a truncated gas profile along the windward “leading edge” of the LMC disk, despite a far more extended stellar component. We explore the implications of this ram pressure stripping signature, using both analytic prescriptions and full three-dimensional hydrodynamic simulations of the LMC. Our simulations can match the radial extent of the LMC’s leading (windward) edge only in scenarios where the MW CGM density at pericentric passage is n{sub p}(R=48.2±5 kpc)=1.1{sub −.45}{sup +.44}×10{sup −4} cm{sup −3}. The implied pericentric density proves to be insensitive to both the broader CGM and temperature profiles, thus providing a model-independent constraint on the local gas density. This result imposes an important constraint on the density profile of the MW’s CGM at 50 kpc, and thus also on the total baryon content of the MW. From our work, assuming a β-profile valid to ∼r {sub vir}, we infer a total diffuse CGM mass M(300 kpc)=2.7±1.4×10{sup 10}M{sub ⊙} or approximately 15% of a 10{sup 12}M{sub ⊙} MW’s baryonic mass budget.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Optimal Sizing and Placement of Battery Energy Storage in Distribution System Based on Solar Size for Voltage Regulation

This paper proposes a new strategy to achieve voltage regulation in distributed power systems in the presence of solar energy sources and battery storage systems. The goal is to find the minimum size of battery storage and its corresponding location in the network based on the size and place of the integrated solar generation. The proposed method formulates the problem by employing the network impedance matrix to obtain an analytical solution instead of using a recursive algorithm such as power flow. The required modifications for modeling the slack and PV buses (generator buses) are utilized to increase the accuracy of the approach. The use of reactive power control to regulate the voltage regulation is not always an optimal solution as in distribution systems R/X is large. In this paper the minimum size and the best place of battery storage is achieved by optimizing the amount of both active and reactive power exchanged by battery storage and its gridtie inverter (GTI) based on the network topology and R/X ratios in the distribution system. Simulation results for the IEEE 14-bus system verify the effectiveness of the proposed approach.

14 SOLAR ENERGY↗
DOE OSTI2023

DUET - Distribution of Understory using Elliptical Transport: A mechanistic model of leaf litter and herbaceous spatial distribution based on tree canopy structure

Heterogeneity in surface fuels produced by overstory trees and understory vegetation is a major driver of fire behavior and ecosystem dynamics. Previous attempts at predicting tree leaf and needle litter accumulation over time have been constrained in scope to probabilistic models that consider a limited number of key factors influencing tree litter dispersal patterns and decomposition processes. Here we present a mechanistic model for estimating variation in surface fuels called the Distribution of Understory using Elliptical Transport (DUET). DUET uses a pre-generated voxelated canopy array and distributes the leaf litter from the trees in an elliptical shape based on tree species characteristics, wind data, and location-specific features. DUET then calculates grass growth between trees and around the litter patterns and sets moisture levels based on relative humidity. In this study, we use DUET to simulate spatial distributions of leaf litter and illustrate model sensitivity to several parameters, providing inputs for testing fuel spatial heterogeneity impacts on fire behavior using process-based models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Findings from Design and Operation of Connected Neighborhoods

Distributed energy resources and load management are an important emerging part of smart grids. Concurrent management of homeowner preferences and utility objectives requires an intelligent negotiation strategy supported by physical infrastructure that learns, optimizes, and controls the system. While there is a wide body of research related to modelling and simulation of distributed resources, relatively little is known about their practical feasibility. One of the first attempts to address this knowledge gap was through a 62-house connected neighborhood in Alabama. This study makes one more step in advancing this knowledge through a 46-townhome demonstration neighborhood located in Atlanta, GA. It reports hardware design, system architecture, and the results from the summer experimental work. The findings are discussed in the context of earlier experience, and analysis is expanded to areas which were not the focus of earlier research.

Tsybina, Eve↗
DOE OSTI2016

An improved multilevel Monte Carlo method for estimating probability distribution functions in stochastic oil reservoir simulations

In this paper, we develop an improved multilevel Monte Carlo (MLMC) method for estimating cumulative distribution functions (CDFs) of a quantity of interest, coming from numerical approximation of large-scale stochastic subsurface simulations. Compared with Monte Carlo (MC) methods, that require a significantly large number of high-fidelity model executions to achieve a prescribed accuracy when computing statistical expectations, MLMC methods were originally proposed to significantly reduce the computational cost with the use of multifidelity approximations. The improved performance of the MLMC methods depends strongly on the decay of the variance of the integrand as the level increases. However, the main challenge in estimating CDFs is that the integrand is a discontinuous indicator function whose variance decays slowly. To address this difficult task, we approximate the integrand using a smoothing function that accelerates the decay of the variance. In addition, we design a novel a posteriori optimization strategy to calibrate the smoothing function, so as to balance the computational gain and the approximation error. The combined proposed techniques are integrated into a very general and practical algorithm that can be applied to a wide range of subsurface problems for high-dimensional uncertainty quantification, such as a fine-grid oil reservoir model considered in this effort. The numerical results reveal that with the use of the calibrated smoothing function, the improved MLMC technique significantly reduces the computational complexity compared to the standard MC approach. Finally, we discuss several factors that affect the performance of the MLMC method and provide guidance for effective and efficient usage in practice.

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

Effects of graded distribution of scattering centers on ballistic transport

The transmission coefficient of a two dimensional scattering region connected to ideal leads was calculated for the case of electrons interacting with an inhomogeneous distribution of repulsive or attractive scattering centers. The scattering centers with Gaussian profiles were positioned at regular intervals perpendicular to the transport direction, but were spaced according to a power law along this direction. The transmission function was obtained using a scattering formalism based on the R-matrix method. The simulations revealed that although, overall, the transmission coefficient decreases and becomes almost monotonously dependent on energy as the inhomogeneity of both attractive and repulsive scattering centers increases, the redistribution of transmission between open channels depends on the type of scattering centers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

ReaxFF molecular dynamics simulation of intermolecular structure formation in acetic acid-water mixtures at elevated temperatures and pressures

The intermolecular structure formation in liquid and supercritical acetic acid-water mixtures was explored using ReaxFF-based molecular dynamics simulations. The microscopic structures of acetic acid-water mixtures with different acetic acid mole fractions (1.0 ≥ x HAc ≥ 0.2) at ambient and critical conditions were examined. The potential energy surface associated with the dissociation of acetic acid molecules was calculated using a metadynamics procedure to optimize the dissociation energy of ReaxFF potential. At ambient conditions, depending on the acetic acid concentration, either acetic acid clusters or water clusters are dominant in the liquid mixture. When acetic acid is dominant (0.4 ≤ x HAc ), cyclic dimers and chain structures between acetic acid molecules are present in the mixture. Both structures disappear at increased water content of the mixture. It was found by simulations that the acetic acid molecules released from these dimer and chain structures tend to stay in a dipole-dipole interaction. These structural changes are in agreement with the experimental findings. When switched to critical conditions, the long-range interactions (e.g., second or fourth neighbor) disappear and the water-water and acetic acid-acetic acid structural formations become disordered. The simulated radial distribution function for water-water interactions is in agreement with experimental and computational studies. The first neighbor interactions between acetic acid and water molecules are preserved at relatively lower temperatures of the critical region. As higher temperatures are reached in the critical region, these interactions were observed to weaken. These simulations suggest that ReaxFF molecular dynamics simulations are an appropriate tool for studying supercritical water/organic acid mixtures.

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

Adaptive stratified importance sampling: hybridization of extrapolation and importance sampling Monte Carlo methods for estimation of wind turbine extreme loads

Wind turbine extreme load estimation is especially difficult because turbulent inflow drives nonlinear turbine physics and control strategies; thus there can be huge differences in turbine response to essentially equivalent environmental conditions. The two main current approaches, extrapolation and Monte Carlo sampling, are both unsatisfying: extrapolation-based methods are dangerous because by definition they make predictions outside the range of available data, but Monte Carlo methods converge too slowly to routinely reach the desired 50-year return period estimates. Thus a search for a better method is warranted. We introduce an adaptive stratified importance sampling approach that allows for treating the choice of environmental conditions at which to run simulations as a stochastic optimization problem that minimizes the variance of unbiased estimates of extreme loads. Furthermore, the framework, built on the traditional bin-based approach used in extrapolation methods, provides a close connection between sampling and extrapolation, and thus allows the solution of the stochastic optimization (i.e., the optimal distribution of simulations in different wind speed bins) to guide and recalibrate the extrapolation. Results show that indeed this is a promising approach, as the variance of both the Monte Carlo and extrapolation estimates are reduced quickly by the adaptive procedure. We conclude, however, that due to the extreme response variability in turbine loads to the same environmental conditions, our method and any similar method quickly reaches its fundamental limits, and that therefore our efforts going forward are best spent elucidating the underlying causes of the response variability.

17 WIND ENERGY↗
DOE OSTI2019

Remanent Magnetic Flux Distribution in Superconducting-Ferromagnetic Layered Heterostructures

Stacks of superconducting tape, made of second-generation coated conductors, present an interesting alternative to bulk superconductors. The layered geometry allows incorporation of other functional, namely ferromagnetic, materials. The ferromagnetic material changes the field configuration within the superconductor and can be used to shape or divert the flux density from its interior. This can be used to optimise the remanent flux density after magnetisation. FEM simulations are employed to test the performance of a variety of superconductor/ferromagnet structures and explain the role of the ferromagnetic material in the composite structure. Cases where trapped field is much higher and lower than remanent magnetisation of the ferromagnetic material are considered, showing that performance can be enhanced or maintained by using ferromagnetic material while reducing the usage of costly superconducting tape.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Transport and Dynamics in Toroidal Fusion Systems

The study entitled, "Transport and Dynamics in Toroidal Fusion Systems," (TDTFS) applied analytical theory and numerical computation to investigate topics of importance to confining plasma, the fourth state of matter, with magnetic fields. A central focus of the work is how non-thermal components of the ion particle distribution affect the "sawtooth" collective oscillation in the core of the tokamak magnetic configuration. Previous experimental and analytical research had shown and described how the oscillation frequency decreases and amplitude increases, leading to "monster" or "giant" sawteeth, when the non-thermal component is increased by injecting particle beams or by exciting ions with imposed electromagnetic waves. The TDTFS study applied numerical computation to self-consistently simulate the interaction between macroscopic collective plasma dynamics and the non-thermal particles. The modeling used the NIMROD code [Sovinec, Glasser, Gianakon, et al., J. Comput. Phys. 195, 355 (2004)] with the energetic component represented by simulation particles [Kim, Parker, Sovinec, and the NIMROD Team, Comput. Phys. Commun. 164, 448 (2004)]. The computations found decreasing growth rates for the instability that drives the oscillations, but they were ultimately limited from achieving experimentally relevant parameters due to computational practicalities. Nonetheless, this effort provided valuable lessons for integrated simulation of macroscopic plasma dynamics. It also motivated an investigation of the applicability of fluid-based modeling to the ion temperature gradient instability, leading to the journal publication [Schnack, Cheng, Barnes, and Parker, Phys. Plasmas 20, 062106 (2013)]. Apart from the tokamak-specific topics, the TDTFS study also addressed topics in the basic physics of magnetized plasma and in the dynamics of the reversed-field pinch (RFP) configuration. The basic physics work contributed to a study of two-fluid effects on interchange dynamics, where "two-fluid" refers to modeling independent dynamics of electron and ion species without full kinetic effects. In collaboration with scientist Ping Zhu, who received separate support, it was found that the rule-of-thumb criteria on stabilizing interchange has caveats that depend on the plasma density and temperature profiles. This work was published in [Zhu, Schnack, Ebrahimi, et al., Phys. Rev. Lett. 101, 085005 (2008)]. An investigation of general nonlinear relaxation with fluid models was partially supported by the TDTFS study and led to the publication [Khalzov, Ebrahimi, Schnack, and Mirnov, Phys. Plasmas 19, 012111 (2012)]. Work specific to the RFP included an investigation of interchange at large plasma pressure and support for applications [for example, Scheffel, Schnack, and Mirza, Nucl. Fusion 53, 113007 (2013)] of the DEBS code [Schnack, Barnes, Mikic, Harned, and Caramana, J. Comput. Phys. 70, 330 (1987)]. Finally, the principal investigator over most of the award period, Dalton Schnack, supervised a numerical study of modeling magnetic island suppression [Jenkins, Kruger, Hegna, Schnack, and Sovinec, Phys. Plasmas 17, 12502 (2010)].

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗