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 1,135 records · Page 63
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Error Characteristics and Scale Dependence of Current Satellite Precipitation Estimates Products in Hydrological Modeling

Satellite precipitation estimates (SPEs) are promising alternatives to gauge observations for hydrological applications (e.g., streamflow simulation), especially in remote areas with sparse observation networks. However, the existing SPEs products are still biased due to imperfections in retrieval algorithms, data sources and post-processing, which makes the effective use of SPEs a challenge, especially at different spatial and temporal scales. In this study, we used a distributed hydrological model to evaluate the simulated discharge from eight quasi-global SPEs at different spatial scales and explored their potential scale effects of SPEs on a cascade of basins ranging from approximately 100 to 130,000 km 2 . The results indicate that, regardless of the difference in the accuracy of various SPEs, there is indeed a scale effect in their application in discharge simulation. Specifically, when the catchment area is larger than 20,000 km2, the overall performance of discharge simulation emerges an ascending trend with the increase of catchment area due to the river routing and spatial averaging. Whereas below 20,000 km 2 , the discharge simulation capability of the SPEs is more randomized and relies heavily on local precipitation accuracy. Our study also highlights the need to evaluate SPEs or other precipitation products (e.g., merge product or reanalysis data) not only at the limited station scale, but also at a finer scale depending on the practical application requirements. Here we have verified that the existing SPEs are scale-dependent in hydrological simulation, and they are not enough to be directly used in very fine scale distributed hydrological simulations (e.g., flash flood). More advanced retrieval algorithms, data sources and bias correction methods are needed to further improve the overall quality of SPEs.

DTVGM↗
DOE OSTI2016

Evaluating gyro-viscosity in the Kelvin-Helmholtz instability by kinetic simulations

In the present paper, the finite-Larmor-radius (gyro-viscous) term [K. V. Roberts and J. B. Taylor, Phys. Rev. Lett. 8, 197–198 (1962)] is evaluated by using a full kinetic Vlasov simulation result of the Kelvin-Helmholtz instability (KHI). The velocity field and the pressure tensor are calculated from the high-resolution data of the velocity distribution functions obtained by the Vlasov simulation, which are used to approximate the Finite-Larmor-Radius (FLR) term according to Roberts and Taylor [Phys. Rev. Lett. 8, 197–198 (1962)]. The direct comparison between the pressure tensor and the FLR term shows an agreement. It is also shown that the anisotropic pressure gradient enhanced the linear growth of the KHI when the inner product between the vorticity of the primary velocity shear layer and the magnetic field is negative, which is consistent with the previous FLR-magnetohydrodynamic simulation result. This result suggests that it is not sufficient for reproducing the kinetic simulation result by fluid simulations to include the FLR term (or the pressure tensor) only in the equation of motion for fluid.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2012

Aggregate Model for Heterogeneous Thermostatically Controlled Loads with Demand Response

Due to the potentially large number of Distributed Energy Resources (DERs) – demand response, distributed generation, distributed storage - that are expected to be deployed, it is impractical to use detailed models of these resources when integrated with the transmission system. Being able to accurately estimate the fast transients caused by demand response is especially important to analyze the stability of the system under different demand response strategies. On the other hand, a less complex model is more amenable to design feedback control strategies for the population of devices to provide ancillary services. The main contribution of this paper is to develop aggregated models for a heterogeneous population of Thermostatic Controlled Loads (TCLs) to accurately capture their collective behavior under demand response and other time varying effects of the system. The aggregated model efficiently includes statistical information of the population and accounts for a second order effect necessary to accurately capture the collective dynamic behavior. The developed aggregated models are validated against simulations of thousands of detailed building models using GridLAB-D (an open source distribution simulation software) under both steady state and severe dynamic conditions caused due to temperature set point changes.

Zhang, Wei↗
DOE OSTI2024

Mesoscale simulation of granular materials under weak shock compaction–pore size distribution effects

This research established a systematic method to generate various pore-size distributions (PSDs) and studied the effect of PSDs on the shock compaction response of granular materials using two-dimensional mesoscale simulations under identical porosity. Simulations utilized various PSDs for three particle shapes (circle, ellipse, and square). Contacting particle configurations using three PSDs, characterized by spatially uniform distributed pores to heterogeneous distributed pores, and non-contacting particle configurations under a single case of PSD were tested. The PSD of generated particle sets was characterized using coordination number, mean diameter, and bimodality coefficient as statistical metrics. Mesoscale simulations showed that regardless of the conditions of pore distributions, shock compaction of granular materials consistently demonstrates a precursor, shock compaction front, and end. However, the shock compaction velocity of contacting particles was dependent on the PSDs despite the constant initial porosity. The compaction velocity was faster in particle configurations with relatively uniform pore distributions than in heterogeneous pore distributions, which our study demonstrated can be attributed to particle rearrangement during compaction. Circular-shaped particles had high sensitivity in shock compaction response to the various PSDs. Furthermore, a contacting particle configuration tended to propagate the shock compaction wave relatively faster than particles that were in a non-contact configuration. This study established the relative importance of considering PSD as a metric over the coordination number in studies of the shock compaction response of granular materials. Further, insights are provided on the evolving shock substructure to characterize the shock compaction response of granular materials.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Numerical and experimental studies on effects of moisture content on combustion characteristics of simulated municipal solid wastes in a fixed bed

Highlights: • The effects of moisture content on the burning process of MSW are investigated. • A two-dimensional mathematical model was built to simulate the combustion process. • Temperature distributions, process rates, gas species were measured and simulated. • The The conversion ratio of C/CO and N/NO in MSW are inverse to moisture content. - Abstract: In order to reveal the features of the combustion process in the porous bed of a waste incinerator, a two-dimensional unsteady state model and experimental study were employed to investigate the combustion process in a fixed bed of municipal solid waste (MSW) on the combustion process in a fixed bed reactor. Conservation equations of the waste bed were implemented to describe the incineration process. The gas phase turbulence was modeled using the k–ε turbulent model and the particle phase was modeled using the kinetic theory of granular flow. The rate of moisture evaporation, devolatilization rate, and char burnout was calculated according to the waste property characters. The simulation results were then compared with experimental data for different moisture content of MSW, which shows that the incineration process of waste in the fixed bed is reasonably simulated. The simulation results of solid temperature, gas species and process rate in the bed are accordant with experimental data. Due to the high moisture content of fuel, moisture evaporation consumes a vast amount of heat, and the evaporation takes up most of the combustion time (about 2/3 of the whole combustion process). The whole bed combustion process reduces greatly as MSW moisture content increases. The experimental and simulation results provide direction for design and optimization of the fixed bed of MSW.

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

Scale Mixtures of Frechet Distributions as Asymptotic Approximations of Extreme Precipitation

This paper is a further development of the results of [20] where, based on the negative binomial model for the duration of wet periods measured in days [16], an asymptotic approximation was proposed for the distribution of the maximum daily precipitation volume within a wet period. This approximation has the form of a scale mixture of the Fr´echet distribution with the gamma mixing distribution and coincides with the distribution of a positive power of a random variable having the Snedecor–Fisher distribution. Here we extend this result to the mth extremes, m ∈ ℕ, and sample quantiles. The proof of this result is based on the representation of the negative binomial distribution as a mixed geometric (and hence, mixed Poisson) distribution [17] and limit theorems for extreme order statistics in samples with random sizes having mixed Poisson distributions [10]. Some analytic properties of the obtained limit distribution are described. In particular, it is demonstrated that under certain conditions the limit distribution of the maximum precipitation is mixed exponential and hence, is infinitely divisible. It is shown that under the same conditions this limit distribution can be represented as a scale mixture of stable or Weibull or Pareto or folded normal laws. The corresponding product representations for the limit random variable can be used for its computer simulation. The results of fitting this distribution to real data are presented illustrating high adequacy of the proposed model. It is also shown that the limit distribution of sample quantiles is the well-known Student distribution. Several methods are proposed for the estimation of the parameters of the asymptotic distributions. The obtained mixture representations for the limit laws and the corresponding asymptotic approximations provide better insight into the nature of mixed probability (“Bayesian”) models.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Model Reduction for Fast Assessment of Grid Impact of High Penetration PV

A key challenge in properly analyzing the impact of solar photovoltaics (PV) during planning, interconnection and operation in the distribution grid is computational complexity. Such studies require the use of Quasi-Steady-Time-Series (QSTS) power flow to capture the complex interaction of PV with voltage regulating equipment and protection devices and a large number of scenario studies are needed to capture location specific constraints and uncertainties. In addition, a typical distribution feeder can have thousands of nodes and there can be thousands of such feeders in a given utility system. To solve this problem, a new and systematic method for distribution network reduction is proposed. The objective is to reduce the size of distribution systems by more than 80% while preserving the accuracy of power flow results. To this end, an extended Ward-equivalent method has been proposed to create an equivalent circuit. Moreover, load and PV aggregation based on the backward sweep method has been developed to further simplify the feeder. The proposed method is tested with a modified IEEE 123-node distribution test system. The simulation results demonstrate the high efficiency and accuracy of the proposed approach for distribution network reduction.

reduced order systems, load flow, load modeling, v↗
DOE OSTI2014

Bringing simulation and observation together to better understand the intergalactic medium

The methods by which one characterizes the distribution of matter in cosmological simulations is intrinsically different from how one performs the same task observationally. In this paper, we make substantial steps toward comparing simulations and observations of the intergalactic medium (IGM) in a more sensible way. We present a pipeline that generates and fits synthetic QSO absorption spectra using sight lines cast through a cosmological simulation, and simultaneously identifies structure by directly analyzing the variations in H I and O VI number density. We compare synthetic absorption spectra with a less observationally motivated, but more straightforward density threshold-based method for finding absorbers. Our efforts focus on H I and O VI to better characterize the warm/hot IGM, a subset of the IGM that is challenging to conclusively identify observationally. We find that the two methods trace roughly the same quantities of H I and O VI above observable column density limits, but the synthetic spectra typically identify more substructure in absorbers. We use both methods to characterize H I and O VI absorber properties. We find that both integrated and differential column density distributions from both methods generally agree with observations. The distribution of Doppler parameters between the two methods are similar for Lyα and compare reasonably with observational results, but while the two methods agree with each other with O VI systems, they both are systematically different from observations. We find a strong correlation between O VI baryon fraction and O VI column density. We also discuss a possible bimodality in the temperature distribution of the gas traced by O VI.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

The Effect of Atmospheric Cooling on Vertical Velocity Dispersion and Density Distribution of Brown Dwarfs

We present a Monte Carlo simulation designed to predict the vertical velocity dispersion of brown dwarfs in the Milky Way. We show that since these stars are constantly cooling, the velocity dispersion has a noticeable trend with the spectral type. With realistic assumptions for the initial mass function, star formation history, and the cooling models, we show that the velocity dispersion is roughly consistent with what is observed for M dwarfs, decreases to cooler spectral types, and increases again for the coolest types in our study (∼T9). We predict a minimum in the velocity dispersions for L/T transition objects, however, the detailed properties of the minimum predominately depend on the star formation history. Since this trend is due to brown dwarf cooling, we expect that the velocity dispersion as a function of spectral type should deviate from the constancy around the hydrogen-burning limit. We convert from velocity dispersion to vertical scale height using standard disk models and present similar trends in disk thickness as a function of spectral type. We suggest that future, wide-field photometric and/or spectroscopic missions may collect sizable samples of distant (∼1 kpc) dwarfs that span the hydrogen-burning limit. As such, we speculate that such observations may provide a unique way of constraining the average spectral type of hydrogen burning.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Reduced Toxicity With Intensity Modulated Radiation Therapy (IMRT) for Desmoplastic Small Round Cell Tumor (DSRCT): An Update on the Whole Abdominopelvic Radiation Therapy (WAP-RT) Experience

Purpose: Desmoplastic small round cell tumor (DSRCT) is a rare malignancy typically involving the peritoneum in young men. Whole abdominopelvic radiation therapy (WAP-RT) using conventional 2-dimensional (2D) radiation therapy (RT) is used to address local recurrence but has been limited by toxicity. Our objectives were to assess the benefit of intensity modulated radiation therapy (IMRT) on toxicity and to update the largest series on radiation for DSRCT. Methods and Materials: The records of 31 patients with DSRCT treated with WAP-RT (22 with 2D-RT and 9 with IMRT) between 1992 and 2011 were retrospectively reviewed. All received multi-agent chemotherapy and maximal surgical debulking followed by 30 Gy of WAP-RT. A further focal boost of 12 to 24 Gy was used in 12 cases. Boost RT and autologous stem cell transplantation were nearly exclusive to patients treated with 2D-RT. Toxicities were assessed with the Common Terminology Criteria for Adverse Events. Dosimetric analysis compared IMRT and simulated 2D-RT dose distributions. Results: Of 31 patients, 30 completed WAP-RT, with a median follow-up after RT of 19 months. Acute toxicity was reduced with IMRT versus 2D-RT: P=.04 for gastrointestinal toxicity of grade 2 or higher (33% vs 77%); P=.02 for grade 4 hematologic toxicity (33% vs 86%); P=.01 for rates of granulocyte colony-stimulating factor; and P=.04 for rates of platelet transfusion. Post treatment red blood cell and platelet transfusion rates were also reduced (P=.01). IMRT improved target homogeneity ([D05-D95]/D05 of 21% vs 46%) and resulted in a 21% mean bone dose reduction. Small bowel obstruction was the most common late toxicity (23% overall). Updated 3-year overall survival and progression-free survival rates were 50% and 24%, respectively. Overall survival was associated with distant metastasis at diagnosis on multivariate analysis. Most failures remained intraperitoneal (88%). Conclusions: IMRT for consolidative WAP-RT in DSRCT improves hematologic toxicity in particular. Although the long-term efficacy of current treatment options remains disappointing, the improved therapeutic index of IMRT may aid in generalizing its use and allowing the addition of novel approaches such as intraperitoneal immunotherapy.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2023

An Ensemble of Neural Networks for Moist Physics Processes, Its Generalizability and Stable Integration

Abstract With the recent advances in data science, machine learning has been increasingly applied to convection and cloud parameterizations in global climate models (GCMs). This study extends the work of Han et al. (2020, https://doi.org/10.1029/2020MS002076 ) and uses an ensemble of 32‐layer deep convolutional residual neural networks, referred to as ResCu‐en, to emulate convection and cloud processes simulated by a superparameterized GCM, SPCAM. ResCu‐en predicts GCM grid‐scale temperature and moisture tendencies, and cloud liquid and ice water contents from moist physics processes. The surface rainfall is derived from the column‐integrated moisture tendency. The prediction uncertainty inherent in deep learning algorithms in emulating the moist physics is reduced by ensemble averaging. Results in 1‐year independent offline validation show that ResCu‐en has high prediction accuracy for all output variables, both in the current climate and in a warmer climate with +4K sea surface temperature. The analysis of different neural net configurations shows that the success to generalize in a warmer climate is attributed to convective memory and the 1‐dimensional convolution layers incorporated into ResCu‐en. We further implement a member of ResCu‐en into CAM5 with real world geography and run the neural‐network‐enabled CAM5 (NCAM) for 5 years without encountering any numerical integration instability. The simulation generally captures the global distribution of the mean precipitation, with a better simulation of precipitation intensity and diurnal cycle. However, there are large biases in temperature and moisture in high latitudes. These results highlight the importance of convective memory and demonstrate the potential for machine learning to enhance climate modeling.

Meteorology & Atmospheric Sciences↗
DOE OSTI2025

Dynamics-based halo model for large scale structure

Accurate modeling of the one-to-two halo transition has long been difficult to achieve. Here, we demonstrate that physically motivated halo definitions that respect the bimodal phase-space distribution of dark matter particles near halos resolves this difficulty. Specifically, the two phase-space components are overlapping and correspond to (1) particles orbiting the halo and (2) particles infalling into the halo for the first time. Motivated by this decomposition, García et al. [Mon. Not. R. Astron. Soc. 521, 2464 (2023)] advocated for defining halos as the collection of particles orbiting their self-generated potential. This definition identifies the traditional one-halo term of the halo-mass correlation function with the distribution of orbiting particles around a halo, while the two-halo term governs the distribution of infalling particles. We use dark matter simulations to demonstrate that the distribution of orbiting particles is finite and can be characterized by a single physical scale 𝑟 h , which we refer to as the halo radius. The two-halo term is described using a simple yet accurate empirical model based on the Zel’dovich correlation function. We further demonstrate that the halo radius imprints itself on the distribution of infalling particles at small scales. Our final model for the halo-mass correlation function is accurate at the ≈ 2% level for 𝑟∈ [0.1, 50] ℎ −1 Mpc. The Fourier transform of our best-fit model describes the halo-mass power spectrum with comparable accuracy for 𝑘 ∈ [0.06, 6.0] ℎ Mpc −1 .

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

OPTIMIZATION OF COOLING DISTRIBUTION OF THE EIC SHC COOLER ERL

The Electron-Ion Collider (EIC) Hadron Storage Ring (HSR) will use strong hadron cooling to maintain the beam brightness and high luminosity during long collision ex?periments. An Energy Recovery Linac is used to deliver the high-current high-brightness electron beam for cooling. For the best cooling effect, the electron beam requires low emittance, small energy spread, and uniform longitudinal distribution. In this work, we simulate and optimize the longitudinal laser-beam distribution shaping at the photo?cathode, modeling space charge forces accurately. Machine parameters such as RF cavity phases are optimized in con?junction with the beam distribution using a genetic optimizer. We demonstrate improvement of the cooling distribution in key parameters.

Wang, N.↗
DOE OSTI2022

Evaluating Basin-Scale Forest Adaptation Scenarios: Wildfire, Streamflow, Biomass, and Economic Recovery Synergies and Trade-Offs

Active forest management is applied in many parts of the western United States to reduce wildfire severity, mitigate vulnerability to drought and bark beetle mortality, and more recently, to increase snow retention and late-season streamflow. A rapidly warming climate accelerates the need for these restorative treatments, but the treatment priority among forest patches varies considerably. We simulated four treatment scenarios across the 3,450 km 2 Wenatchee River basin in eastern Washington, United States. We used a decision support tool (DST) to assess trade-offs and synergies within and among treatments on wildfire risk and smoke emissions, water yield and snow retention, biomass production, and economic return. Treatment scenarios emphasized prescribed burning ( BurnOnly ), biomass production ( MaxBiomass ), gap-based thinning to optimize water yield ( IdealWater ), and a principle-based restoration scenario ( RA1 ). Fire hazard, smoke emissions, and biomass production metrics were evaluated across scenarios using the Forest Vegetation Simulator, and water yields were modeled using the Distributed Hydrology Soil Vegetation Model. Simulations were summarized to both patch- (10 1 –10 2 ha) and subwatershed- (10 3 –10 4 ha) scales, and treatment effects were evaluated against an untreated baseline landscape. We used logic models to rank effect sizes by scenario across metrics along a continuum between −1 (no or weak effect) to +1 (large effect). All treatments produced benefits across one or more ecosystem services and led to synergistic benefits to water yield and wildfire hazard reduction. Tradeoffs among resource benefits were clear in wilderness where reliance on prescribed burning without mechanical treatment increased costs and eliminated the potential for biomass recovery. The BurnOnly scenario improved fire risk metrics and streamflow, but effect sizes were lower compared to other treatments. IdealWater showed the strongest benefits overall, demonstrating the ability to capture multiple resource benefits through spatially explicit thinning. Our study provides a framework for integrating strategic and tactical models that evaluate tradeoffs and synergies gained through varied management approaches. We demonstrate the utility of decision support modeling to enhance management synergies across large landscapes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

DIFFUSE AURORA ON GANYMEDE DRIVEN BY ELECTROSTATIC WAVES

The role of electrostatic electron cyclotron harmonic (ECH) waves in producing diffuse auroral emission O i 1356 Å on Ganymede is investigated. Electron precipitation flux entering the atmosphere of Ganymede due to pitch-angle diffusion by ECH waves into the atmospheric loss-cone is calculated. The analytical yield spectrum approach for electron energy degradation in gases is used for calculating diffuse auroral intensities. It is found that calculated O i 1356 Å intensity resulting from the precipitation of magnetospheric electrons observed near Ganymede is insufficient to account for the observed diffuse auroral intensity. This is in agreement with estimates made in earlier works. Heating and acceleration of ambient electrons by ECH wave turbulence near the magnetic equator on the field line connecting Ganymede and Jupiter are considered. Two electron distribution functions are used to simulate the heating effect by ECH waves. Use of a Maxwellian distribution with temperature 100 eV can produce about 50–70 Rayleigh O i 1356 Å intensities, and the kappa distribution with characteristic energy 50 eV also gives rise to intensities with similar magnitude. Numerical experiments are performed to study the effect of ECH wave spectral intensity profile, ECH wave amplitude, and temperature/characteristic energy of electron distribution functions on the calculated diffuse auroral intensities. The proposed missions, joint NASA/ESA Jupiter Icy Moon Explorer and the present JUNO mission to Jupiter, would provide new data to constrain the ECH wave and other physical parameters near Ganymede. These should help confirm the findings of the present study.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2011

The impact of size distribution assumptions in a bulk one-moment microphysics scheme on simulated surface precipitation and storm dynamics during a low-topped supercell case in Belgium

In this research the impact of modifying the size distribution assumptions of the precipitating hydrometeors in a bulk one-moment microphysics scheme on simulated surface precipitation and storm dynamics has been explored for long-lived low-topped supercells in Belgium. It was shown that weighting the largest precipitating ice species of the microphysics scheme to small graupel results in an increase of surface precipitation because of counteracting effects. On the one hand, the precipitation formation process slowed down, resulting in lower precipitation efficiency. On the other hand, latent heat release associated with freezing favored more intense storms. In contrast to previous studies finding decreased surface precipitation when graupel was present in the microphysics parameterization, storms were rather shallow in the authors simulations. This left little time for graupel sublimation. The impact of size distribution assumptions of snow was found to be small, but more realistic size distribution assumptions of rain led to the strongest effect on surface precipitation. Cold pools shrunk because of weaker rain evaporation at the cold pool boundaries, leading to a decreased surface rain area.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Acceleration of Oxygen Ions in Dipolarization Events: 1. CPS Distributions

Abstract Using an MHD simulation of near tail reconnection associated with a flow burst and the collapse (dipolarization) of the inner tail in combination with test particle tracing we study the acceleration and flux increases of energetic oxygen ions (O + ). The characteristic orbits, distributions, and acceleration mechanisms are governed by the dimensionless parameter σ = ω ci t n , where ω ci is the ion gyro frequency and t n a characteristic Alfvén time of the MHD simulation. For σ < 1, central plasma sheet (CPS) populations after the passage of the dipolarization front are found to resemble half‐shells in velocity space oriented toward dusk. They originate from within the CPS and are energized typically by a single encounter of the region of enhanced cross‐tail electric field associated with the flow burst. For larger σ values ( σ > 1) the O + distributions resemble more closely those of protons, consisting of two counter‐streaming field‐aligned beams and an, albeit more tenuous and irregular, ring population perpendicular to the magnetic field. The existence of the beams, however, depends on suitable earthward moving source populations in the plasma sheet boundary layer or the adjacent lobes. The acceleration to higher energies is found to indicate a charge dependence, consistent with a dominance of more highly charged ions at energies of a few hundred keV. As in earlier simulations, the simulated fluxes show large anisotropies and nongyrotropic effects, phase bunching, and spatially and temporally localized beams.

Birn, J.↗