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 829 records · Page 46
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Characterization and Modeling of High-Intensity Evolution in the SNS Beam Test Facility

Modern high-power accelerators are charged with delivering reliable beam with low losses. Resolving the complex dynamics arising from space charge and nonlinear forces requires detailed models of the accelerator and particle-in-cell simulation. There has historically been discrepancy between simulated and measured beam distributions, particularly at the low-density halo level. The Beam Test Facility (BTF) at the Spallation Neutron Source is outfitted to study beam evolution in a high-power linear accelerator MEBT. This includes capability for high-dimensional measurements of the post-RFQ beam distribution, including interplane correlations that may be the key to accurate simulation. Beam is transported through a 4.6 m FODO channel (9.5 cells) to a second distribution measurement stage. Plans for validating simulations against BTF measurements of beam evolution in the FODO channel are discussed.

Ruisard, Kiersten↗
DOE OSTI2016

Garrotxa cosmological simulations of milky way-sized galaxies: general properties, hot-gas distribution, and missing baryons

We introduce a new set of simulations of Milky Way (MW)-sized galaxies using the AMR code ART + hydrodynamics in a Λ cold dark matter cosmogony. The simulation series is called GARROTXA and it follows the formation of a halo/galaxy from z = 60 to z = 0. The final virial mass of the system is ∼7.4 × 10{sup 11} M {sub ⊙}. Our results are as follows. (a) Contrary to many previous studies, the circular velocity curve shows no central peak and overall agrees with recent MW observations. (b) Other quantities, such as M-∗(6 × 10{sup 10} M {sub ⊙}) and R {sub d} (2.56 kpc), fall well inside the observational MW range. (c) We measure the disk-to-total ratio kinematically and find that D/T = 0.42. (d) The cold-gas fraction and star formation rate at z = 0, on the other hand, fall short of the values estimated for the MW. As a first scientific exploitation of the simulation series, we study the spatial distribution of hot X-ray luminous gas. We have found that most of this X-ray emitting gas is in a halo-like distribution accounting for an important fraction but not all of the missing baryons. An important amount of hot gas is also present in filaments. In all our models there is not a massive disk-like hot-gas distribution dominating the column density. Our analysis of hot-gas mock observations reveals that the homogeneity assumption leads to an overestimation of the total mass by factors of 3–5 or to an underestimation by factors of 0.7–0.1, depending on the used observational method. Finally, we confirm a clear correlation between the total hot-gas mass and the dark matter halo mass of galactic systems.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Rapid QSTS Simulations for High-Resolution Comprehensive Assessment of Distributed PV

The rapid increase in penetration of distributed energy resources on the electric power distribution system has created a need for more comprehensive interconnection modeling and impact analysis. Unlike conventional scenario-based studies, quasi-static time-series (QSTS) simulations can realistically model time-dependent voltage controllers and the diversity of potential impacts that can occur at different times of year. However, to accurately model a distribution system with all its controllable devices, a yearlong simulation at 1-second resolution is often required, which could take conventional computers a computational time of 10 to 120 hours when an actual unbalanced distribution feeder is modeled. This computational burden is a clear limitation to the adoption of QSTS simulations in interconnection studies and for determining optimal control solutions for utility operations. The solutions we developed include accurate and computationally efficient QSTS methods that could be implemented in existing open-source and commercial software used by utilities and the development of methods to create high-resolution proxy data sets. This project demonstrated multiple pathways for speeding up the QSTS computation using new and innovative methods for advanced time-series analysis, faster power flow solvers, parallel processing of power flow solutions and circuit reduction. The target performance level for this project was achieved with year-long high-resolution time series solutions run in less than 5 minutes within an acceptable error.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2014

Laboratory optical spectroscopy of the thiophenoxy radical and its profile simulation as a diffuse interstellar band based on rotational distribution by radiation and collisions

The gas-phase optical absorption spectrum of the thiophenoxy radical (C{sub 6}H{sub 5}S), a diffuse interstellar band (DIB) candidate molecule, was observed in the discharge of thiophenol using a cavity ringdown spectrometer. The ground-state rotational constants of the thiophenoxy radical were theoretically calculated, and the excited-state rotational constants were determined from the observed rotational profile. The rotational profile of a near prolate molecule having C {sub 2v} symmetry was simulated on the basis of a rotational distribution model by radiation and collisions. Although the simulated profile did not agree with the observed DIBs, the upper limit of the column density for the thiophenoxy radical in the diffuse clouds toward HD 204827 was evaluated to be 2 × 10{sup 13} cm{sup –2}. The profile simulation indicates that rotational distribution by radiation and collisions is important to reproduce a rotational profile for a DIB candidate and that the near prolate C {sub 2v} molecule is a possible candidate for DIB with a band width variation dependent on the line of sight.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Structural, dynamic, and vibrational properties during heat transfer in Si/Ge superlattices: A Car-Parrinello molecular dynamics study

The structural, dynamic, and vibrational properties during heat transfer process in Si/Ge superlattices are studied by analyzing the trajectories generated by the ab initio Car-Parrinello molecular dynamics simulation. The radial distribution functions and mean square displacements are calculated and further discussions are made to explain and probe the structural changes relating to the heat transfer phenomenon. Furthermore, the vibrational density of states of the two layers (Si/Ge) are computed and plotted to analyze the contributions of phonons with different frequencies to the heat conduction. Coherent heat conduction of the low frequency phonons is found and their contributions to facilitate heat transfer are confirmed. The Car-Parrinello molecular dynamics simulation outputs in the work show reasonable thermophysical results of the thermal energy transport process and shed light on the potential applications of treating the heat transfer in the superlattices of semiconductor materials from a quantum mechanical molecular dynamics simulation perspective.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Statistical analysis of spatial distribution of pores in metal additive manufacturing

We report mechanical properties of additively manufactured parts are sensitive to the presence of pores form during the manufacturing process. The impact of pores on the mechanical performance has been investigated extensively with respect to different parameters, such as pore volume fraction, shape, and size. However, statistical investigations focusing on the relationships between the spatial distribution of pores and process parameters; and consequently, the performance of the manufactured parts are scattered and limited. Also, these sparse investigations usually suffer from an ambiguous definition of terminologies. For instance, the required criteria to consider a point pattern as complete spatial randomness (CSR) are generally not clarified. Moreover, no numerical formalism is yet developed to show how much the observed spatial results are statistically significant. To address these shortcomings, the statistical definition of CSR in a point pattern and the procedure to quantitatively determine the deviation of a pattern from CSR were discussed. The explained statistical approach was used to investigate the effect of scanning speed parameter on the spatial distribution of gas pores in laser powder bed fusion manufactured stainless steel parts. Furthermore, and to highlight the impact of the spatial distribution of pores on mechanical properties, fatigue performances of parts with clustered and randomly distributed pores were simulated by finite element analysis. It is shown that by reducing the scanning speed, the spatial distribution of gas pores deviated more from CSR, and correspondingly fatigue performance deteriorated.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Computational modeling of graphite degradation in molten salt reactors: Role of infiltration

Molten salt reactors (MSRs) often employ graphite as a moderator and reflector. An important challenge for deploying graphite in these reactors is that, due to limited experimental data, our understanding of graphite’s structural integrity in molten salt environments remains incomplete. Here, this study addresses heat generation from fuel-bearing salt that has infiltrated open pores in the graphite, driven primarily by pressure differentials. This is one of multiple identified physical and chemical mechanisms through which molten salt could potentially degrade graphite. Thermally driven stresses are quantified using the Molten-Salt Reactor Experiment (MSRE) graphite moderator elements as a case study. Finite element simulations predict stress distributions at varying infiltration levels, indicating that thermal stresses increase with higher infiltration. Rare-event simulations using the parallel subset simulation framework identify the combinations and corresponding ranges of input parameters that lead to stresses above a specified threshold. In particular, combinations involving high infiltration amounts, high power density, and low thermal conductivity tend to induce the highest stresses. Under the inputs and assumptions considered in this work, the magnitudes of the thermally driven stresses are quite low, with a very low likelihood of causing failure due to exceeding the graphite’s tensile strength. Additionally, rare-event simulations were performed for two more scenarios: a scaled-up moderator geometry and a localized hotspot in the original geometry. Both cases resulted in increased susceptibility to failure, though not to a detrimental extent. Furthermore, the combined effects of irradiation and infiltration-induced thermal stresses were evaluated. The results showed that thermal stresses from infiltration were negligible compared to those caused by irradiation. The findings of such a study are inherently component-specific, but the methodology presented here could be used for similar assessments of salt-infiltration effects in other graphite components.

36 - MATERIALS SCIENCE↗
DOE OSTI2019

Cloud Phase and Relative Humidity Distributions over the Southern Ocean in Austral Summer Based on In Situ Observations and CAM5 Simulations

Cloud phase and relative humidity (RH) distributions at -67° to 0°C over the Southern Ocean during austral summer are compared between in situ airborne observations and global climate simulations. A scale-aware comparison is conducted using horizontally averaged observations from 0.1 to 50 km. Cloud phase frequencies, RH distributions, and liquid mass fraction are found to be less affected by horizontal resolutions than liquid and ice water content (LWC and IWC, respectively), liquid and ice number concentrations (Nc liq and Nc ice , respectively), and ice supersaturation (ISS) frequency. At -10° to 0°C, observations show 27%–34% and 17%–37% of liquid and mixed phases, while simulations show 60%–70% and 3%–4%, respectively. Simulations overestimate (underestimate) LWC and Nc liq in liquid (mixed) phase, overestimate Nc ice in mixed phase, underestimate IWC in ice and mixed phases, and underestimate (overestimate) liquid mass fraction below (above) -5°C, indicating that observational constraints are needed for different cloud phases. RH frequently occurs at liquid saturation in liquid and mixed phases for all datasets, yet the observed RH in ice phase can deviate from liquid saturation by up to 20%–40% at -20° to 0°C, indicating that the model assumption of liquid saturation for coexisting ice and liquid is inaccurate for low liquid mass fractions (<0.1). Simulations lack RH variability for partial cloud fractions (0.1–0.9) and underestimate (overestimate) ISS frequency for cloud fraction <0.1 (≥0.6), implying that improving RH subgrid-scale parameterizations may be a viable path to account for small-scale processes that affect RH and cloud phase heterogeneities. Lastly, two sets of simulations (nudged and free-running) show very similar results (except for ISS frequency) regardless of sample sizes, corroborating the statistical robustness of the model–observation comparisons.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2014

WE-A-17A-12: The Influence of Eye Plaque Design On Dose Distributions and Dose- Volume Histograms

Purpose: To investigate the effect of slot design of the model EP917 plaque on dose distributions and dose-volume histograms (DVHs). Methods: The dimensions and orientation of the slots in EP917 plaques were measured. In the MCNP5 radiation simulation geometry, dose distributions on orthogonal planes and DVHs for a tumor and sclera were generated for comparisons. 27 slot designs and 13 plaques were evaluated and compared with the published literature and the Plaque Simulator clinical treatment planning system. Results: The dosimetric effect of the gold backing composition and mass density was < 3%. Slot depth, width, and length changed the central axis (CAX) dose distributions by < 1% per 0.1 mm in design variation. Seed shifts in the slot towards the eye and shifts of the {sup 125} I-coated Ag rod within the capsule had the greatest impact on CAX dose distribution, increasing by 14%, 9%, 4%, and 2.5% at 1, 2, 5, and 10 mm, respectively, from the inner sclera. Along the CAX, dose from the full plaque geometry using the measured slot design was 3.4% ± 2.3% higher than the manufacturer-provided geometry. D{sub 10} for the simulated tumor, inner sclera, and outer sclera for the measured plaque was also higher, but 9%, 10%, and 20%, respectively. In comparison to the measured plaque design, a theoretical plaque having narrow and deep slots delivered 30%, 37%, and 62% lower D{sub 10} doses to the tumor, inner sclera, and outer sclera, respectively. CAX doses at −1, 0, 1, and 2 mm were also lower by a factor of 2.6, 1.4, 1.23, and 1.13, respectively. Conclusion: The study identified substantial sensitivity of the EP917 plaque dose distributions to slot design. However, it did not identify substantial dosimetric variations based on radionuclide choice ({sup 125}I, {sup 103}Pd, or {sup 131}Cs). COMS plaques provided lower scleral doses with similar tumor dose coverage.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2017

A simulation-based efficiency comparison of AC and DC power distribution networks in commercial buildings

Direct current (DC) power distribution has recently gained traction in buildings research due to the proliferation of on-site electricity generation and battery storage, and an increasing prevalence of internal DC loads. The research discussed in this paper uses Modelica-based simulation to compare the efficiency of DC building power distribution with an equivalent alternating current (AC) distribution. The buildings are all modeled with solar generation, battery storage, and loads that are representative of the most efficient building technology. A variety of paramet ric simulations determine how and when DC distribution proves advantageous. These simulations also validate previous studies that use simpler approaches and arithmetic efficiency models. This work shows that using DC distribution can be considerably more efficient: a medium sized office building using DC distribution has an expected baseline of 12% savings, but may also save up to 18%. In these results, the baseline simulation parameters are for a zero net energy (ZNE) building that can island as a microgrid. DC is most advantageous in buildings with large solar capacity, large battery capacity, and high voltage distribution.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Energy Resilient Recovery in Puerto Rico: Hurricane Maria Recovery Support

This report covers microgrids as a special application of distributed energy resources (DERs) in Puerto Rico. This is one of a series of reports describing the U.S. Department of Energy (DOE) multi-laboratory efforts undertaken. To ensure the sustainable, long-term recovery of Puerto Rico's electric power grid from hurricanes Maria and Irma and to build capacity to manage future potential natural disasters in the most secure and resilient way, DOE convened experts from multiple national laboratories to develop a comprehensive set of data, models, analytic tools, and studies, considering inputs from a wide variety of stakeholder groups, to support technically sound recommendations for Puerto Rico's energy investment decisions. In Phase 1 of the multi-laboratory effort to support Puerto Rico's recovery, the National Renewable Energy Laboratory (NREL) provided the utility company Puerto Rico Electric Power Authority (PREPA) recommendations for a new framework of interconnection standards to accelerate the integration of utility-scale, transmission-connected, renewable electrical generation and energy storage that ensure cross-technology compatibility and enable high deployment levels without compromising grid reliability, safety, or security (https://www.nrel.gov/docs/fy19osti/73848.pdf). In Phase 2, NREL published a report focused on the interconnection of distributed energy resources (DERs) to the electric distribution system in Puerto Rico (https://www.nrel.gov/docs/fy21osti/77127.pdf). This report familiarized the reader with Puerto Rico's distribution infrastructure and operational practices and procedures that are relevant to DER interconnection. The report also provided considerations for streamlining the interconnection process given the expected increase in deployments resulting from Puerto Rico's renewable portfolio standard goal of 100% renewables by 2050. Accordingly, the report identifies considerations and concerns associated with the increase in intermittent generation, strategies for DER interconnection best practices, and the potential use of the latest technological solutions identified in the latest revision of the Institute of Electrical and Electronics Engineers 1547-2018 interconnection standard. Additionally, the report identified ways to improve the physical resiliency of installed DERs. This task describes how distributed energy systems may participate in and contribute to improving the resiliency of Puerto Rico's energy infrastructure in the future. This report presents an analysis of the Puerto Rico grid to illustrate the possible different microgrid and minigrid scenarios. The capability to form sustained microgrids and minigrids has the potential to provide a high level of resiliency during disturbances on the grid. DERs, on the other hand, have the capability to aid in the microgrid operation by providing grid support as well as grid-forming functionality. Other tasks in Phase 3 include voltage regulation in distribution networks using DERs and transmission-and-distribution co-simulation.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

A statistical and simulation-informed model for estimating permeability from pore size distribution in saturated geomaterials

Accurate permeability estimation is essential across subsurface engineering applications but remains challenging due to the complex pore structures of natural geomaterials. Traditional empirical methods and simplified theoretical models often inadequately capture the role of pore size distribution and connectivity. Here, this study develops a statistical and simulation-informed permeability model that collapses pore-scale complexity into a compact scaling of the form k = αϕμ d 2 , where ϕ is porosity, μ d is mean pore size, and α is a weakly varying coefficient. By combining pore network simulations with statistical analysis of unimodal and bimodal pore size distributions, we identify three key findings: (i) permeability is much more sensitive to mean pore size than to porosity; (ii) across extensive datasets, the ratio σ d /μ d (standard deviation to mean) clusters around a characteristic value ∼0.4, allowing the effects of the full pore size distribution to be represented by μ d and a narrowly varying α ≈ 0.05; and (iii) for bimodal systems, there exists a critical fraction of small pores ∼0.78 above which flow becomes small-pore dominated, enabling the definition of an effective flow-controlling pore population and facilitating simplified permeability estimation for such systems. The resulting model, which requires only porosity and a representative mean pore size as inputs, is validated against comprehensive experimental datasets (>1700 samples) spanning diverse soils and rocks and achieves good predictive accuracy. Overall, this work provides a physically grounded yet practically simple permeability estimator suitable for subsurface engineering, environmental protection, and resource management applications.

Permeability↗
DOE OSTI2016

Evidence for a liquid-liquid critical point in supercooled water within the E3B3 model and a possible interpretation of the kink in the homogeneous nucleation line

Supercooled water exhibits many thermodynamic anomalies, and several scenarios have been proposed to interpret them, among which the liquid-liquid critical point (LLCP) hypothesis is the most commonly discussed. We investigated Widom lines and the LLCP of deeply supercooled water, by using molecular dynamics simulation with a newly reparameterized water model that explicitly includes three-body interactions. Seven isobars are studied from ambient pressure to 2.5 kbar, and Widom lines are identified by calculating maxima in the coefficient of thermal expansion and the isothermal compressibility (both with respect to temperature). From these data we estimate that the LLCP of the new water model is at 180 K and 2.1 kbar. The oxygen radial distribution function is calculated along the 2 kbar isobar. It shows a steep change in the height of its second peak between 180 and 185 K, which indicates a transition between the high-density liquid and low-density liquid phases and which is consistent with the ascribed location of the critical point. The good agreement of the height of the second peak of the radial distribution function between simulation and experiment at 1 bar, as a function of temperature, supports the validity of the model. The location of the LLCP within the model is close to the kink in the experimental homogeneous nucleation line. We use existing experimental data to argue that the experimental LLCP is at 168 K and 1.95 kbar and speculate how this LLCP and its Widom line might be responsible for the kink in the homogeneous nucleation line.

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

Cluster infall for mass calibration in the stage-IV era

The outskirts of galaxy clusters present a promising avenue for constraining cluster masses in a way that is robust to the impact of baryonic physics. We assess the accuracy to which the cluster infall regions can be used for cluster mass calibration. Building on previous work, we parametrize the velocity distribution 𝑃⁡(𝑣r,𝑣tan|𝑟,𝑀) of dark matter halos on scales 𝑟 ≥ 5⁢ℎ −1 Mpc as the product of the marginalized distribution 𝑃⁡(𝑣 r |𝑟,𝑀) and the conditional distribution 𝑃⁡(𝑣 tan |𝑣 r ,𝑟,𝑀), calibrating the radial and mass dependence of these distributions in numerical simulations. We then project our model along the line of sight to obtain accurate predictions for the distributions of line-of-sight velocities at a given projected radius and cluster mass 𝑃⁡(𝑣 LOS |𝑅,𝑀), which we can observe with spectroscopic survey data. Furthermore, with our model, we forecast that spectra from the Dark Energy Spectroscopic Instrument can constrain cluster masses with subpercent-level precision, comparable to that of stage-IV weak lensing surveys.

79 ASTRONOMY AND ASTROPHYSICS↗