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 865 records · Page 48
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Evaluating Domestic Hot Water Distribution System Options With Validated Analysis Models

A developing body of work is forming that collects data on domestic hot water consumption, water use behaviors, and energy efficiency of various distribution systems. A full distribution system developed in TRNSYS has been validated using field monitoring data and then exercised in a number of climates to understand climate impact on performance. This study builds upon previous analysis modelling work to evaluate differing distribution systems and the sensitivities of water heating energy and water use efficiency to variations of climate, load, distribution type, insulation and compact plumbing practices. Overall 124 different TRNSYS models were simulated. Of the configurations evaluated, distribution losses account for 13-29% of the total water heating energy use and water use efficiency ranges from 11-22%. The base case, an uninsulated trunk and branch system sees the most improvement in energy consumption by insulating and locating the water heater central to all fixtures. Demand recirculation systems are not projected to provide significant energy savings and in some cases increase energy consumption. Water use is most efficient with demand recirculation systems, followed by the insulated trunk and branch system with a central water heater. Compact plumbing practices and insulation have the most impact on energy consumption (2-6% for insulation and 3-4% per 10 gallons of enclosed volume reduced). The results of this work are useful in informing future development of water heating best practices guides as well as more accurate (and simulation time efficient) distribution models for annual whole house simulation programs.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2012

Minimization of Blast furnace Fuel Rate by Optimizing Burden and Gas Distribution

The goal of the research is to improve the competitive edge of steel mills by using the advanced CFD technology to optimize the gas and burden distributions inside a blast furnace for achieving the best gas utilization. A state-of-the-art 3-D CFD model has been developed for simulating the gas distribution inside a blast furnace at given burden conditions, burden distribution and blast parameters. The comprehensive 3-D CFD model has been validated by plant measurement data from an actual blast furnace. Validation of the sub-models is also achieved. The user friendly software package named Blast Furnace Shaft Simulator (BFSS) has been developed to simulate the blast furnace shaft process. The research has significant benefits to the steel industry with high productivity, low energy consumption, and improved environment.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

High-Frequency Signature-Based Fault Detection for Future MV Distribution Grids

Increasing penetration levels of inverter based distributed energy resources (DERs) impact the legacy distribution system protection. Inverter based DERs provide approximately 1.2-2 pu fault current. In systems with high penetration of inverter based DERs, it is difficult for over-current based protection schemes to differentiate between normal loading conditions and a fault. Directional, distance, and adaptive forms of protection schemes are also affected by low fault currents. This paper analyzes fault generated traveling wave (TW) based high-frequency signatures in the distribution system. In order to simulate such signatures, frequency-dependent distributed parameter line modeling approach is used in this research work to represent distribution lines and underground cables. Modified IEEE 13-bus medium voltage test system is modeled in electromagnetic transient simulation tool and multiple transient scenarios are simulated in this test system. The results are analyzed to understand the high-frequency signatures that can be used to detect and locate faults under high penetration of DERs.

41 EE - Solar Energy Technologies Office (EE-4S)↗
DOE OSTI2020

High-Frequency Signature-Based Fault Detection for Future MV Distribution Grids

Increasing penetration levels of inverter based distributed energy resources (DERs) impact the legacy distribution system protection. Inverter based DERs provide approximately 1.2-2 pu fault current. In systems with high penetration of inverter based DERs, it is difficult for over-current based protection schemes to differentiate between normal loading conditions and a fault. Directional, distance, and adaptive forms of protection schemes are also affected by low fault currents. This paper analyzes fault generated traveling wave (TW) based high-frequency signatures in the distribution system. In order to simulate such signatures, frequency-dependent distributed parameter line modeling approach is used in this research work to represent distribution lines and underground cables. Modified IEEE 13-bus medium voltage test system is modeled in electromagnetic transient simulation tool and multiple transient scenarios are simulated in this test system. The results are analyzed to understand the high-frequency signatures that can be used to detect and locate faults under high penetration of DERs.

41 EE - Solar Energy Technologies Office (EE-4S)↗
DOE OSTI2017

Significant inconsistency of vegetation carbon density in CMIP5 Earth system models against observational data

Abstract Earth system models (ESMs) have been widely used for projecting global vegetation carbon dynamics, yet how well ESMs performed for simulating vegetation carbon density remains untested. We compiled observational data of vegetation carbon density from literature and existing data sets to evaluate nine ESMs at site, biome, latitude, and global scales. Three variables—root (including fine and coarse roots), total vegetation carbon density, and the root:total vegetation carbon ratios (R/T ratios), were chosen for ESM evaluation. ESM models performed well in simulating the spatial distribution of carbon densities in root ( r = 0.71) and total vegetation ( r = 0.62). However, ESM models had significant biases in simulating absolute carbon densities in root and total vegetation biomass across the majority of land ecosystems, especially in tropical and arctic ecosystems. Particularly, ESMs significantly overestimated carbon density in root (183%) and total vegetation biomass (167%) in climate zones of 10°S–10°N. Substantial discrepancies between modeled and observed R/T ratios were found: the R/T ratios from ESMs were relatively constant, approximately 0.2 across all ecosystems, along latitudinal gradients, and in tropic, temperate, and arctic climatic zones, which was significantly different from the observed large variations in the R/T ratios (0.1–0.8). There were substantial inconsistencies between ESM‐derived carbon density in root and total vegetation biomass and the R/T ratio at multiple scales, indicating urgent needs for model improvements on carbon allocation algorithms and more intensive field campaigns targeting carbon density in all key vegetation components.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Heating of aqueous suspensions of silicon nanoparticles by a 808-nm diode laser for application in local photohyperthermia

3D simulation of the distribution of thermal fields in aqueous suspensions of silicon nanoparticles upon irradiation by a high-power 808-nm laser diode operating in cw and quasi-cw repetitively pulsed regimes is performed. It is shown that irradiation by a repetitively pulsed laser with a peak power of 10 W and a pulse duration of 300 ms can form a region with a maximum temperature exceeding 42 °C, which is promising for application in cancer photohyperthermia for local heating of biotissues at a given depth. (paper)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2012

EVOLUTION OF POST-IMPACT COMPANION STARS IN SN Ia REMNANTS WITHIN THE SINGLE-DEGENERATE SCENARIO

The nature of the progenitor systems of Type Ia supernovae is still uncertain. One way to distinguish between the single-degenerate scenario and double-degenerate scenario is to search for the post-impact remnant star. To examine the characteristics of the post-impact remnant star, we have carried out three-dimensional hydrodynamic simulations of supernova impacts on main-sequence-like stars. We explore the evolution of the post-impact remnants using the stellar evolution code MESA. We find that the luminosity and radius of the remnant star dramatically increase just after the impact. After the explosion, post-impact companions continue to expand on a progenitor-dependent timescale of {approx}10{sup 2.5}-10{sup 3} years before contracting. It is found that the time evolution of the remnant star is dependent not only on the amount of energy absorbed but also on the depth of the energy deposition. We examine the viability of the candidate star Tycho G as the possible remnant companion in Tycho's supernova by comparing it to the evolved post-impact remnant stars in our simulations. The closest model in our simulations has a similar effective temperature, but the luminosity and radius are twice as large. By examining the angular momentum distribution in our simulations, we find that the surface rotational speed could drop to {approx}10 km s{sup -1} if the specific angular momentum is conserved during the post-impact evolution, implying that Tycho G cannot be completely ruled out because of its low surface rotation speed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Secrets of the Dark Universe: Simulating the Sky on the Blue Gene/Q, The Outer Rim Simulation

An astonishing 99.6% of our Universe is dark. Observations indicate that the Universe consists of 70% of a mysterious dark energy and 25% of a yet-unidentified dark matter component, and only 0.4% of the remaining ordinary matter is visible. Understanding the physics of this dark sector is the foremost challenge in cosmology today. Sophisticated simulations of the evolution of the Universe play a crucial task in this endeavor. This movie shows an intermediate stage in a large simulation of the distribution of matter in the Universe, the so-called cosmic web, accounting for the influence of dark energy. The simulation is evolving 1.1 trillion particles. The movie shows a snapshot of the Universe when it was 1.6 billion years old. Credits: Science: Hal Finkel, Salman Habib, Katrin Heitmann, Kalyan Kumaran, Vitali Morozov, Tom Peterka, Adrian Pope, Tim Williams, David Daniel, Patricia Fasel, Nicholas Frontiere and Zarija Lukic. Visualization: Mark Herald, Joseph A. Insley, Michael E. Papka and Venkatram Vishwanath

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2011

Computer Simulation of Electron Thermalization in CsI and CsI(Tl)

A Monte Carlo (MC) model was developed and implemented to simulate the thermalization of electrons in inorganic scintillator materials. The model incorporates electron scattering with both longitudinal optical and acoustic phonons. In this paper, the MC model was applied to simulate electron thermalization in CsI, both pure and doped with a range of thallium concentrations. The inclusion of internal electric fields was shown to increase the fraction of recombined electron-hole pairs and to broaden the thermalization distance and thermalization time distributions. The MC simulations indicate that electron thermalization, following {gamma}-ray excitation, takes place within approximately 10 ps in CsI and that electrons can travel distances up to several hundreds of nanometers. Electron thermalization was studied for a range of incident {gamma}-ray energies using electron-hole pair spatial distributions generated by the MC code NWEGRIM (NorthWest Electron and Gamma Ray Interaction in Matter). These simulations revealed that the partition of thermalized electrons between different species (e.g., recombined with self-trapped holes or trapped at thallium sites) vary with the incident energy. Implications for the phenomenon of nonlinearity in scintillator light yield are discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2011

Quasilinear theory and simulation of Buneman instability

In a recently developed nonlinear theory of Buneman instability, a simplifying assumption of self-similarity was imposed for the electron distribution function, based upon which, a set of moment kinetic equations was derived and solved together with nonlinear wave kinetic equation [P. H. Yoon and T. Umeda, Phys. Plasmas 17, 112317 (2010)]. It was found that the theoretical result compared reasonably against one-dimensional electrostatic Vlasov simulation. In spite of this success, however, the simulated distribution deviated appreciably from the assumed self-similar form during the late stages of nonlinear evolution. In order to rectify this shortcoming, in this paper, the distribution function is computed on the basis of rigorous velocity space diffusion equation. A novel theoretical scheme is developed so that both the quasilinear particle diffusion equation and the adiabatic dispersion relation can be solved for an arbitrary particle distribution function. Comparison with Vlasov simulation over relatively early quasilinear phase of the instability shows a reasonable agreement, despite the fact that quasilinear theory lacks coherent nonlinear effects as well as mode-mode coupling effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Maximum entropy distributions of dark matter in ΛCDM cosmology

Context. Small-scale challenges to ΛCDM cosmology require a deeper understanding of dark matter physics. Aims. This paper aims to develop the maximum entropy distributions for dark matter particle velocity (denoted by X ), speed (denoted by Z ), and energy (denoted by E ) that are especially relevant on small scales where system approaches full virialization. Methods. For systems involving long-range interactions, a spectrum of halos of different sizes is required to form to maximize system entropy. While the velocity in halos can be Gaussian, the velocity distribution throughout the entire system, involving all halos of different sizes, is non-Gaussian. With the virial theorem for mechanical equilibrium, we applied the maximum entropy principle to the statistical equilibrium of entire system, such that the maximum entropy distribution of velocity (the X distribution) could be analytically derived. The halo mass function was not required in this formulation, but it did indeed result from the maximum entropy. Results. The predicted X distribution involves a shape parameter α and a velocity scale, v 0 . The shape parameter α reflects the nature of force ( α → 0 for long-range force or α → ∞ for short-range force). Therefore, the distribution approaches Laplacian with α → 0 and Gaussian with α → ∞. For an intermediate value of α , the distribution naturally exhibits a Gaussian core for v ≪ v 0 and exponential wings for v ≫ v 0 , as confirmed by N -body simulations. From this distribution, the mean particle energy of all dark matter particles with a given speed, v , follows a parabolic scaling for low speeds (∝ v 2 for v ≪ v 0 in halo core region, i.e., “Newtonian”) and a linear scaling for high speeds (∝ v for v ≫ v 0 in halo outskirt, i.e., exhibiting “non-Newtonian” behavior due to long-range gravity). We compared our results against N -body simulations and found a good agreement.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Influence of Intermittency on the Quasi-perpendicular Scaling in Three-dimensional Magnetohydrodynamic Turbulence

Solar wind fluctuations reveal the ubiquity of intermittency, which is believed to affect the spectral signatures of turbulence. In this work, based on simulation of driven compressible MHD turbulence, we apply the wavelet technique to the magnetic field and velocity to identify intermittency, and we analyze the influence of the intermittency on the quasi-perpendicular scaling in the inertial range. The numerical results show that the original magnetic and velocity fluctuations are anisotropic, and have a power anisotropy with a spectral index approaching the Iroshnikov–Kraichnan −3/2 scaling in the direction quasi-perpendicular to the local mean magnetic field. As in observations of the solar wind fluctuations, as the scale decreases in the simulation, the calculated probability distribution functions (pdfs) of the wavelet coefficients become extended on both tails of the non-Gaussian distribution, with a rapid increase in flatness. After intermittency has been removed from the driven turbulence, at each scale, the pdfs approach a Gaussian distribution, with the flatness being ∼3. Meanwhile, the quasi-perpendicular scaling for both fluctuations becomes steeper and close to a Kolmogorov −5/3 scaling, which may be a result of the stronger intermittency in the quasi-perpendicular direction and at the smaller scales. These results suggest that there is intermittency superposed on the “background” turbulence that seems to have the Kolmogorov scaling, whereby the overall slope is getting flatter with the involvement of intermittency.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Petascale Application of a Coupled CPU-GPU Algorithm for Simulation and Analysis of Multiphase Flow Solutions in Porous Medium Systems, In: 2014 IEEE 28th International Parallel and Distributed Processing Symposium

Large-scale simulation can provide a wide range of information needed to develop and validate theoretical models for multiphase flow in porous medium systems. In this paper, we consider a coupled solution in which a multiphase flow simulator is coupled to an analysis approach used to extract the interfacial geometries as the flow evolves. This has been implemented using MPI to target heterogeneous nodes equipped with GPUs. The GPUs evolve the multiphase flow solution using the lattice Boltzmann method while the CPUs compute up scaled measures of the morphology and topology of the phase distributions and their rate of evolution. Our approach is demonstrated to scale to 4,096 GPUs and 65,536 CPU cores to achieve a maximum performance of 244,754 million-lattice-node updates per second (MLUPS) in double precision execution on Titan. In turn, this approach increases the size of systems that can be considered by an order of magnitude compared with previous work and enables detailed in situ tracking of averaged flow quantities at temporal resolutions that were previously impossible. Furthermore, it virtually eliminates the need for post-processing and intensive I/O and mitigates the potential loss of data associated with node failures.

Computer Science↗
DOE OSTI2014

Nucleation and stabilization of carbon-rich structures in interstellar media

We study the conditions under which carbon clusters of different sizes form and stabilize. We describe the approach to equilibrium by simulating tenuous carbon gas dynamics to long times. First, we use reactive molecular dynamics simulations to describe the nucleation of long chains, large clusters, and complex cage structures in carbon- and hydrogen-rich interstellar gas phases. We study how temperature, particle density, the presence of hydrogen, and carbon inflow affect the nucleation of molecular moieties with different characteristics, in accordance with astrophysical conditions. We extend the simulations to densities that are orders of magnitude lower than current laboratory densities, to temperatures that are relevant to circumstellar environments of planetary nebulae, and microsecond formation times. We correlate cluster size distributions from the simulations with thermodynamic equilibrium at low temperatures and gas densities, where entropy plays a significant role.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Predictive Inverse Model for Advective Heat Transfer in a Short–Circuited Fracture: Dimensional Analysis, Machine Learning, and Field Demonstration

Identifying fluid flow maldistribution in planar geometries is a well–established problem in subsurface science/engineering. Of particular importance to the thermal performance of enhanced (or “engineered”) geothermal systems is identifying the existence of nonuniform (i.e., heterogeneous) permeability and subsequently predicting advective heat transfer. Here, machine learning via a genetic algorithm (GA) identifies the spatial distribution of an unknown permeability field in a two–dimensional Hele–Shaw geometry (i.e., parallel plates). The inverse problem is solved by minimizing the L2 norm between simulated residence time distribution (RTD) and measurements of an inert tracer breakthrough curve (BTC) (C–Dot nanoparticle). Principal component analysis (PCA) of spatially correlated permeability fields enabled reduction of the parameter space by more than a factor of 10 and restricted the inverse search to reservoir–scale permeability variations. Thermal experiments and tracer tests conducted at the mesoscale Altona Field Laboratory (AFL) demonstrate that the method accurately predicts the effects of extreme flow channeling on heat transfer in a single bedding–plane rock fracture. However, this is only true when the permeability distributions provide adequate matches to both tracer RTD and frictional pressure loss. Without good agreement to frictional pressure loss, it is still possible to match a simulated RTD to measurements, but subsequent predictions of heat transfer are grossly inaccurate. Here, the results of this study suggest that it is possible to anticipate the thermal effects of flow maldistribution, but only if both simulated RTDs and frictional pressure loss between fluid inlets and outlets are in good agreement with measurements.

42 ENGINEERING↗
DOE OSTI2010

Wide field array calibration dependence on the stability of measured dose distributions

Purpose: The aim of this work was to simulate the effect of dose distribution changes on detector array calibrations and to explore compensatory methods that are used during calibration measurements. Methods: The array calibration technique that was investigated is known as wide field (WF) calibration. Using this method, a linear array [y-axis (65 detectors) of the IC PROFILER (Sun Nuclear Corporation, Melbourne, FL)] is calibrated with three measurements ({alpha}, {theta}, and {lambda}); each measurement uses the same radiation field, which is larger than the array. For measurement configuration {theta}, the array is rotated by 180 deg. from its position in {alpha}; for {lambda}, the array is shifted by one detector from its position in {theta}. The relative detector sensitivities are then determined through ratios of detector readings at the same field locations (using {theta} and {lambda}). This method results in error propagation that is proportional to the number of detectors in the array. During the procedure, the calibration protocol operates under three postulates, which state that (a) the beam shape does not change between measurements; (b) the relative sensitivities of the detectors do not change; and (c) the scatter to the array does not change as the array is moved. The WF calibration's sensitivity to a postulate (a) violation was quantified by applying a sine shaped perturbation (of up to 0.1%) to {alpha}, {theta}, or {lambda}, and then determining the change relative to a baseline calibration. Postulate (a) violations were minimized by using a continuous beam and mechanized array movement during {theta} and {lambda}. A continuously on beam demonstrated more stable beam symmetry as compared to cycling the beam on and off between measurements. Additional side-scatter was also used to satisfy postulate (c). Results: Simulated symmetry perturbations of 0.1% to {theta} or {lambda} resulted in calibration errors of up to 2%; {alpha} was relatively immune to perturbation (<0.1% error). Wide field calibration error on a linear accelerator with similar symmetry variations was {+-}1.6%. Using a continuous beam during {theta} and {lambda} with additional side-scatter reduced the calibration error from {+-}1.6% to {+-}0.48%. Conclusions: This work increased the reproducibility of WF calibrations by limiting the effect of measurement perturbations primarily due to linear accelerator symmetry variations. The same technique would work for any array using WF calibration.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

Midplane neutral density profiles in the National Spherical Torus Experiment

Atomic and molecular density data in the outer midplane of NSTX [Ono et al., Nucl. Fusion 40, 557 (2000)] are inferred from tangential camera data via a forward modeling procedure using the DEGAS 2 Monte Carlo neutral transport code. The observed Balmer-β light emission data from 17 shots during the 2010 NSTX campaign display no obvious trends with discharge parameters such as the divertor Balmer-α emission level or edge deuterium ion density. Simulations of 12 time slices in 7 of these discharges produce molecular densities near the vacuum vessel wall of 2–8 × 10{sup 17 }m{sup −3} and atomic densities ranging from 1 to 7 × 10{sup 16 }m{sup −3}; neither has a clear correlation with other parameters. Validation of the technique, begun in an earlier publication, is continued with an assessment of the sensitivity of the simulated camera image and neutral densities to uncertainties in the data input to the model. The simulated camera image is sensitive to the plasma profiles and virtually nothing else. The neutral densities at the vessel wall depend most strongly on the spatial distribution of the source; simulations with a localized neutral source yield densities within a factor of two of the baseline, uniform source, case. The uncertainties in the neutral densities associated with other model inputs and assumptions are ≤50%.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗