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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
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 OSTI2011

Online Robot Dead Reckoning Localization Using Maximum Relative Entropy Optimization With Model Constraints

The principle of Maximum relative Entropy optimization was analyzed for dead reckoning localization of a rigid body when observation data of two attached accelerometers was collected. Model constraints were derived from the relationships between the sensors. The experiment's results confirmed that accelerometers each axis' noise can be successfully filtered utilizing dependency between channels and the dependency between time series data. Dependency between channels was used for a priori calculation, and a posteriori distribution was derived utilizing dependency between time series data. There was revisited data of autocalibration experiment by removing the initial assumption that instantaneous rotation axis of a rigid body was known. Performance results confirmed that such an approach could be used for online dead reckoning localization.

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↗
NASA NTRS2011

Background and Pickup Ion Velocity Distribution Dynamics in Titan's Plasma Environment: 3D Hybrid Simulation and Comparison with CAPS T9 Observations

In this report we discuss the ion velocity distribution dynamics from the 3D hybrid simulation. In our model the background, pickup, and ionospheric ions are considered as a particles, whereas the electrons are described as a fluid. Inhomogeneous photoionization, electron-impact ionization and charge exchange are included in our model. We also take into account the collisions between the ions and neutrals. The current simulation shows that mass loading by pickup ions H(+); H2(+), CH4(+) and N2(+) is stronger than in the previous simulations when O+ ions are introduced into the background plasma. In our hybrid simulations we use Chamberlain profiles for the atmospheric components. We also include a simple ionosphere model with average mass M = 28 amu ions that were generated inside the ionosphere. The moon is considered as a weakly conducting body. Special attention will be paid to comparing the simulated pickup ion velocity distribution with CAPS T9 observations. Our simulation shows an asymmetry of the ion density distribution and the magnetic field, including the formation of the Alfve n wing-like structures. The simulation also shows that the ring-like velocity distribution for pickup ions relaxes to a Maxwellian core and a shell-like halo.

Lipatov, A. S.↗
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↗
NASA NTRS2012

A Comparison of Filter-based Approaches for Model-based Prognostics

Model-based prognostics approaches use domain knowledge about a system and its failure modes through the use of physics-based models. Model-based prognosis is generally divided into two sequential problems: a joint state-parameter estimation problem, in which, using the model, the health of a system or component is determined based on the observations; and a prediction problem, in which, using the model, the stateparameter distribution is simulated forward in time to compute end of life and remaining useful life. The first problem is typically solved through the use of a state observer, or filter. The choice of filter depends on the assumptions that may be made about the system, and on the desired algorithm performance. In this paper, we review three separate filters for the solution to the first problem: the Daum filter, an exact nonlinear filter; the unscented Kalman filter, which approximates nonlinearities through the use of a deterministic sampling method known as the unscented transform; and the particle filter, which approximates the state distribution using a finite set of discrete, weighted samples, called particles. Using a centrifugal pump as a case study, we conduct a number of simulation-based experiments investigating the performance of the different algorithms as applied to prognostics.

prognostics↗
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 OSTI2011

Cosmic string induced CMB maps

We compute maps of CMB temperature fluctuations seeded by cosmic strings using high resolution simulations of cosmic strings in a Friedmann-Robertson-Walker universe. We create full-sky, 18 deg. and 3 deg. CMB maps, including the relevant string contribution at each resolution from before recombination to today. We extract the angular power spectrum from these maps, demonstrating the importance of recombination effects. We briefly discuss the probability density function of the pixel temperatures, their skewness, and kurtosis.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
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↗
DOE OSTI2020

Spectral solar irradiance on inclined surfaces: A fast Monte Carlo approach

Estimating spectral plane-of-array (POA) solar irradiance on inclined surfaces is an important step in the design and performance evaluation of both photovoltaic and concentrated solar plants. This work introduces a fast, line-by-line spectral, Monte Carlo (MC) radiative transfer model approach to simulate anisotropic distributions of shortwave radiation through the atmosphere as photon bundles impinge on inclined surfaces. This fast Monte Carlo approach reproduces the angular distribution of solar irradiance without the undesirable effects of spatial discretization and thus computes detailed POA irradiance values on surfaces at any orientation and also when surfaces are subjected to the anisotropic ground and atmospheric scattering. Polarization effects are also easily incorporated into this approach that can be considered as direct numerical simulation of the physics involved. Here, we compare our Monte Carlo radiative transfer model with the most widely used empirical transposition model, Perez4, under various conditions. The results show that the Perez4 model reproduces the more detailed Monte Carlo simulations with less than 10% deviation under clear skies for all relevant surface tilt and azimuth angles. When optically thin clouds are present, observed deviations are larger, especially when the receiving surface is strongly tilted. Deviations are also observed for large azimuth angle differences between the receiving surface and the solar position. When optically thick clouds are present, the two models agree within 15% deviation for nearly all surface orientation and tilt angles. The overall deviations are smaller when compared with cases for optically thin clouds. The Perez4 model performs very well (6.0% deviation) in comparison with the detailed MC simulations for all cases, thus validating its widespread use for practical solar applications. When detailed atmospheric profiles and cloud optical properties are available, the proposed fast Monte Carlo radiative model reproduces accurate spectral and angular POA irradiance levels for various atmospheric and cloud cover conditions, surface orientations, and different surface and ground properties.

42 ENGINEERING↗