Mesoscale simulations of particulate flows with parallel distributed Lagrange multiplier technique
Abstract not provided
SEARCH · Engineering Papers
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.
Abstract not provided
Abstract not provided.
Abstract not provided.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
In Fermilab's neutrino facilities such as the Neutrinos at the Main Injector (NuMI) and the upcoming Long Baseline Neutrino Facility (LBNF), a proton beam strikes high-power target, producing positively and negatively charged pions and kaons. There is a need for detailed simulations in order to capture all particle interactions and beam propagation from protons on target to short-lived mesons decaying into muons and neutrinos. The generation of individual beam simulations is a resource-intensive and time-consuming process. In this paper, we describe a method through which many simulation samples with high statistics can be generated to study the effects of beam scan across a target for given beam configurations.
This paper discusses the modeling, analysis, and testing in a real-time simulation environment of the Lanai power grid system for the integration and control of PhotoVoltaic (PV) distributed generation. The Lanai Island in Hawaii is part of the Hawaii Clean Energy Initiative (HCEI) to transition to 30% renewable green energy penetration by 2030. In Lanai the primary loads come from two Castle and Cook Resorts, in addition to residential needs. The total peak load profile is 12470 V, 5.5 MW. Currently there are several diesel generators that meet these loading requirements. As part of the HCEI, Lanai has initially installed 1.2 MW of PV generation. The goal of this study has been to evaluate the impact of the PV with respect to the conventional carbon-based diesel generation in real time simulation. For intermittent PV distributed generation, the overall stability and transient responses are investigated. A simple Lanai 'like' model has been developed in the Matlab/Simulink environment (see Fig. 1) and to accommodate real-time simulation of the hybrid power grid system the Opal-RT Technologies RT-Lab environment is used. The diesel generators have been modelled using the SimPowerSystems toolbox swing equations and a custom Simulink module has been developed for the High level PV generation. All of the loads have been characterized primarily as distribution lines with series resistive load banks with one VAR load bank. Three-phase faults are implemented for each bus. Both conventional and advanced control architectures will be used to evaluate the integration of the PV onto the current power grid system. The baseline numerical results include the stable performance of the power grid during varying cloud cover (PV generation ramping up/down) scenarios. The importance of assessing the real-time scenario is included.
In this paper, a normally-off AlGaN/GaN MIS-HEMT with improved threshold voltage (V{sub TH}) thermal stability is reported with investigations on its physical mechanism. The normally-off operation of the device is achieved from novel short argon plasma treatment (APT) prior to the fluorine plasma treatment (FPT) on Al{sub 2}O{sub 3} gate dielectrics. For the MIS-HEMT with FPT only, its V{sub TH} drops from 4.2 V at room temperature to 0.5 V at 200 °C. Alternatively, for the device with APT-then-FPT process, its V{sub TH} can retain at 2.5 V at 200 °C due to the increased amount of deep-level traps that do not emit electrons at 200 °C. This thermally stable V{sub TH} makes this device suitable for high power applications. The depth profile of the F atoms in Al{sub 2}O{sub 3}, measured by the secondary ion mass spectroscopy, reveals a significant increase in the F concentration when APT is conducted prior to FPT. The X-ray photoelectron spectroscopy (XPS) analysis on the plasma-treated Al{sub 2}O{sub 3} surfaces observes higher composition of Al-F bonds if APT was applied before FPT. The enhanced breaking of Al-O bonds due to Ar bombardment assisted in the increased incorporation of F radicals at the surface during the subsequent FPT process. The Schrödinger equation of Al{sub 2}O{sub x}F{sub y} cells, with the same Al-F compositions as obtained from XPS, was solved by Gaussian 09 molecular simulations to extract electron state distribution as a function of energy. The simulation results show creation of the deeper trap states in the Al{sub 2}O{sub 3} bandgap when APT is used before FPT. Finally, the trap distribution extracted from the simulations is verified by the gate-stress experimental characterization to confirm the physical mechanism described.
Pressure solution is an important process in the evolution of sedimentary rocks, which provide storage space for most of our petroleum resources. It directly influences the generation, migration, and storage of petroleum fluids in subsurface sedimentary rocks. Here, in this paper, we develop a pore-scale, mechanochemical model to demonstrate a possible positive feedback between the local porosity and pore surface stress, in which a higher local porosity causes a higher local pore surface stress, thus enhancing pressure solution and consequently further increasing the local porosity. Pore surface stress represents stress on a solid grain adjacent to a pore. Specifically, the pore-scale, mechanochemical model directly simulates the stress distribution over solid and pore surfaces using a finite element model. The dissolution of solids at the solid-pore interfaces under a far-from-equilibrium condition is simulated using a first-order kinetics model that accounts for the local stress distribution. The updated pore geometry, caused by pore surface dissolution, is then used in the stress simulation in the next numerical iteration. Two types of porous media, the Oriskany sandstone and an artificial porous medium with spherical pores, were tested in the mechanochemical simulation. The positive stress-porosity feedback during pressure solution was observed in both samples. In addition, the model quantitatively illustrated the distribution of local mineral dissolution rates on all pore surfaces, as well as its relation to the effective mineral dissolution rate of the entire sample. Based on the comparison between the two porous media, the local mineral dissolution was regulated by pore space distribution, geometry, and coalescence during pressure solution. This work is the first that uses direct, pore-scale numerical simulation to demonstrate the positive stress-porosity feedback during pressure solution, which has the potential to advance the understanding of the mechanical-chemical (MC) coupling in many geological processes that are relevant to subsurface energy systems, such as the recovery of petroleum hydrocarbons and geothermal energy.
Grid operators are now considering using distributed energy resources (DERs) to provide distribution voltage regulation rather than installing costly voltage regulation hardware. DER devices include multiple adjustable reactive power control functions, so grid operators have the difficult decision of selecting the best operating mode and settings for the DER. In this work, we develop a novel state estimation-based particle swarm optimization (PSO) for distribution voltage regulation using DER-reactive power setpoints and establish a methodology to validate and compare it against alternative DER control technologies (volt–VAR (VV), extremum seeking control (ESC)) in increasingly higher fidelity environments. Distribution system real-time simulations with virtualized and power hardware-in-the-loop (PHIL)-interfaced DER equipment were run to evaluate the implementations and select the best voltage regulation technique. Each method improved the distribution system voltage profile; VV did not reach the global optimum but the PSO and ESC methods optimized the reactive power contributions of multiple DER devices to approach the optimal solution.
Two sets of 48-hr hindcast experiments for Mei-yu season using non-hydrostatic global variable-resolution model (MPAS-Atmosphere) with regional refinement at 4km resolution (V4km) and uniform 60km resolution (U60km) are first-time conducted to attempt at simulating the Mei-yu rainfall in 2015 over East China. The analysis focuses on the seasonal average of diurnal variation of precipitation amount (PA), intensity (PI), and frequency (PF). Both simulations reasonably reproduce the spatial distribution of PA. V4km is more skillful in simulating the spatial distributions and magnitudes of PF and PI. The diurnal cycle of Mei-yu rainfall shows a major early morning peak and a minor afternoon peak, contributed by precipitation during two sub-periods with distinct synoptic circulations. With strong Mei-yu in the first sub-period, the diurnal variation of PA is controlled by nocturnal southwesterly jet. V4km overestimates the morning peak mainly due to its bias in simulating boundary layer inertial oscillation. Although differences in PA between the simulations are small, U60km overestimates PF and underestimates PI. With weak Mei-yu in the second sub-period, the diurnal variation of PA is controlled by both synoptic circulation and local convection. At both resolutions, deviations in the large-scale circulation modulated by a few typhoons lead to positive biases in the morning peak of PA. After removing the typhoon impacts, V4km captures the observed diurnal cycle of PA well, while U60km significantly underestimates PA and PI particularly in the afternoon. Future studies focusing on advancing modeling of southwesterly jet and typhoons may further improve convection permitting simulation of Mei-yu rainfall.
There is no consensus on how baryon feedback shapes the underlying matter distribution from either simulations or observations. We confront the uncertain landscape by jointly analyzing new measurements of the gas distribution around groups and clusters—DESI+ACT kinetic Sunyaev–Zel’dovich (kSZ) effect profiles and eROSITA X-ray gas masses—with mean halo masses characterized by galaxy–galaxy lensing. Across a wide range of halo masses ( M 500 = 10 13−14 M ⊙ ) and redshifts (0 < z < 1), we find evidence of more efficient gas expulsion beyond several R 500 than predicted by most state-of-the-art simulations. A like-with-like comparison reveals all kSZ and X-ray observations are inconsistent with the fiducial 1 Gpc 3 hydrodynamical FLAMINGO simulation, which was calibrated to reproduce pre-eROSITA X-ray gas fractions: eROSITA X-ray gas fractions are 2 × lower than the simulation, and the kSZ measurements are combined >8σ discrepant. The FLAMINGO simulation variant with the most gas expulsion, and therefore the most suppression of the matter power spectrum relative to a dark-matter-only simulation, provides a good description of how much gas is expelled and how far it extends; the enhanced gas depletion is achieved by more powerful but less frequent AGN outbursts. Joint kSZ, X-ray, and lensing measurements form a consistent picture of gas expulsion beyond several R 500 , implying a more suppressed matter power spectrum than predicted by most recent simulations. Complementary observables (e.g., thermal Sunyaev–Zel’dovich effect and fast radio bursts) and next-generation simulations are critical to understanding the physical mechanism behind this extreme gas expulsion and mapping its impact on the large-scale matter distribution.
Occupancy has significant impacts on building performance. However, in current building performance simulation programs, occupancy inputs are static and lack diversity, contributing to discrepancies between the simulated and actual building performance. This work presents an Occupancy Simulator that simulates the stochastic behavior of occupant presence and movement in buildings, capturing the spatial and temporal occupancy diversity. Each occupant and each space in the building are explicitly simulated as an agent with their profiles of stochastic behaviors. The occupancy behaviors are represented with three types of models: (1) the status transition events (e.g., first arrival in office) simulated with probability distribution model, (2) the random moving events (e.g., from one office to another) simulated with a homogeneous Markov chain model, and (3) the meeting events simulated with a new stochastic model. A hierarchical data model was developed for the Occupancy Simulator, which reduces the amount of data input by using the concepts of occupant types and space types. Finally, a case study of a small office building is presented to demonstrate the use of the Simulator to generate detailed annual sub-hourly occupant schedules for individual spaces and the whole building. The Simulator is a web application freely available to the public and capable of performing a detailed stochastic simulation of occupant presence and movement in buildings. Future work includes enhancements in the meeting event model, consideration of personal absent days, verification and validation of the simulated occupancy results, and expansion for use with residential buildings.
Understanding the nature of lanthanide transport in metallic fuels and its relation to fuel clad chemical interaction (FCCI) has been has been a topic that is long on interest and short on data. Recent experiments have shown a strong temperature-gradient based diffusion (Soret diffusion) and have thus provided vitally important information for modelers. However, details about the actual coefficients and many of the relevant processes that crucially participate in the cerium migration and precipitation process have not yet been disclosed. Using a combination of computing tools (primarily python and BISON) we are exploring methods that may extract useful information from visual data. In this specific example, we are using what is formally called 'blob detection' algorithms to analyzed the precipitate distribution of cerium in a U-Zr metal that has had a fixed temperature gradient enforced upon it for a specified period of time. These algorithms, once properly adapted to the problem at hand, may be useful at extracting information from radial and axial scans of metallic fuels which are much more abundant than direct measurements of lanthanides. We have developed a model that qualitatively works to move mass along a temperature gradient and precipitate heavily in the cold regions. Additionally, we have observed that the solubility of cerium in U-Zr held under this kind of temperature gradient is subject to some amount of U-Zr constituent redistribution. Our results predict movement of the Zr toward the hot end of the rod followed by a U-rich layer with the cooler 3/5 of the rod staying close the the initial U-Zr distribution. A more complete analysis will define groups of correlated parameters that correctly simulate the resulting distribution of precipitates. This restricted parameters space will give smaller bounds to modelers and will expedite the analysis of future works as new and more accurate data are made available. A follow-up experiment that started with 1/10 of the initial concentration of cerium may provide valuable additional information with regard to the effect of the rates of precipitation and dissolution as well as demonstrating the potential rapidity of migration in a lower burnup scenario. (authors)
Directional detection is an important way to detect dark matter. An input for these experiments is the dark matter velocity distribution. Recent hydrodynamical simulations have shown that the dark matter velocity distribution differs substantially from the Standard Halo Model. We study the impact of some of these updated velocity distributions in dark matter directional detection experiments. Here, we calculate the ratio of events required to confirm the forward-backward asymmetry and the existence of the ring of maximum recoil rate using different dark matter velocity distributions for 19 F and Xe targets. We show that with the use of updated dark matter velocity profiles, the forward-backward asymmetry and the ring of maximum recoil rate can be confirmed using a factor of ~ 2– 3 less events when compared to that using the Standard Halo Model.
Organic liquid scintillation detectors are valuable for the detection of special nuclear material since they are capable of detecting both neutrons and gamma rays. Scintillators can also provide energy information which is helpful in identification and characterization of the source. In order to design scintillation based measurement systems appropriate simulation tools are needed. MCNPX-PoliMi is capable of simulating scintillation detector response; however, simulations have traditionally been run in analog mode which leads to long computation times. In this paper, non-analog MCNPX-PoliMi mode which uses variance reduction techniques is applied and tested. The non-analog MCNPX-PoliMi simulation test cases use source biasing, geometry splitting and a combination of both variance reduction techniques to efficiently simulate pulse height distribution and then time-of-flight for a heavily shielded case with a {sup 252}Cf source. An improvement factor (I), is calculated for distributions in each of the three cases above to analyze the effectiveness of the non-analog MCNPX-PoliMi simulations in reducing computation time. It is found that of the three cases, the last case which uses a combination of source biasing and geometry splitting shows the most improvement in simulation run time for the same desired variance. For pulse height distributions speedup ranging from a factor 5 to 25 is observed, while for time-of-flights the speedup factors range from 3 to 10. (authors)
Cyclonic atmospheric vortices of varying intensity, collectively known as low-pressure systems (LPS), travel northwest across central India and produce more than half of the precipitation received by that fertile region and its ~600 million inhabitants. Yet, future changes in LPS activity are poorly understood, due in part to inadequate representation of these storms in current climate models. Using a high-resolution atmospheric general circulation model that realistically simulates the genesis distribution of LPS, here we show that Indian monsoon LPS activity declines about 45% by the late 21st century in simulations of a business-as-usual emission scenario. The distribution of LPS genesis shifts poleward as it weakens, with oceanic genesis decreasing by ~60% and continental genesis increasing by ~10%; over land the increase in storm counts is accompanied by a shift toward lower storm wind speeds. The weakening and poleward shift of the genesis distribution in a warmer climate are confirmed and attributed, via a statistical model, to the reduction and poleward shift of low-level absolute vorticity over the monsoon region, which in turn are robust features of most coupled model projections. The poleward shift in LPS activity results in an increased frequency of extreme precipitation events over northern India.
This paper presents a method for quantifying and enabling the resiliency of a power distribution system (PDS) using analytical hierarchical process and percolation theory. Using this metric, quantitative analysis can be done to analyze the impact of possible control decisions to pro-actively enable the resilient operation of distribution system with multiple microgrids and other resources. Developed resiliency metric can also be used in short term distribution system planning. The benefits of being able to quantify resiliency can help distribution system planning engineers and operators to justify control actions, compare different reconfiguration algorithms, develop proactive control actions to avert power system outage due to impending catastrophic weather situations or other adverse events. Validation of the proposed method is done using modified CERTS microgrids and a modified industrial distribution system. Furthermore, simulation results show topological and composite metric considering power system characteristics to quantify the resiliency of a distribution system with the proposed methodology, and improvements in resiliency using two-stage reconfiguration algorithm and multiple microgrids.