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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
DOE OSTI2012

DOES A ''STOCHASTIC'' BACKGROUND OF GRAVITATIONAL WAVES EXIST IN THE PULSAR TIMING BAND?

We investigate the effects of gravitational waves (GWs) from a simulated population of binary supermassive black holes (SMBHs) on pulsar timing array data sets. We construct a distribution describing the binary SMBH population from an existing semi-analytic galaxy formation model. Using realizations of the binary SMBH population generated from this distribution, we simulate pulsar timing data sets with GW-induced variations. We find that the statistics of these variations do not correspond to an isotropic, stochastic GW background. The ''Hellings and Downs'' correlations between simulated data sets for different pulsars are recovered on average, though the scatter of the correlation estimates is greater than expected for an isotropic, stochastic GW background. These results are attributable to the fact that just a few GW sources dominate the GW-induced variations in every Fourier frequency bin of a five-year data set. Current constraints on the amplitude of the GW signal from binary SMBHs will be biased. Individual binary systems are likely to be detectable in five-year pulsar timing array data sets where the noise is dominated by GW-induced variations. Searches for GWs in pulsar timing array data therefore need to account for the effects of individual sources of GWs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Baryon Acoustic Oscillations reconstruction with pixels

Gravitational non-linear evolution induces a shift in the position of the baryon acoustic oscillations (BAO) peak together with a damping and broadening of its shape that bias and degrades the accuracy with which the position of the peak can be determined. BAO reconstruction is a technique developed to undo part of the effect of non-linearities. We present and analyse a reconstruction method that consists of displacing pixels instead of galaxies and whose implementation is easier than the standard reconstruction method. We show that this method is equivalent to the standard reconstruction technique in the limit where the number of pixels becomes very large. This method is particularly useful in surveys where individual galaxies are not resolved, as in 21cm intensity mapping observations. We validate this method by reconstructing mock pixelated maps, that we build from the distribution of matter and halos in real- and redshift-space, from a large set of numerical simulations. We find that this method is able to decrease the uncertainty in the BAO peak position by 30-50% over the typical angular resolution scales of 21 cm intensity mapping experiments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2019

On benchmarking of simulations of particle transport in ITER

We present results of benchmarking of core particle transport simulations by a collection of codes widely used in transport modelling of tokamak plasmas. Our analysis includes formulation of transport equations, difference between electron and ion solvers, comparison of modules of the pellet and edge gas fuelling on the ITER baseline scenario. During the first phase of benchmarking we address the particle transport effects in the stationary phase. Firstly, simulations are performed with identical sources, sinks, transport coefficients, and boundary conditions prescribed in the flattop H-mode phase. The transformation of ion particle transport equations is introduced so to directly compare their results to electron transport solvers. Moreover, the pellet fuelling models are benchmarked in various conditions to evaluate the dependency of the pellet deposition on the pellet volume, injection side, pedestal and separatrix parameters. Thirdly, edge gas fuelling is benchmarked to assess sensitivities of source profile predictions to uncertainties in plasma conditions and detailed model assumptions. At the second phase, we address particle transport effects in the time-evolving plasma including the current ramp-up to the ramp-down phase. The ion and the electron solvers are benchmarked together. Differences between the simulation results of the solvers are investigated in terms of equilibrium, grid resolution, radial coordinate, radial grid distribution, and plasma volume evolution term. We found that the selection of the radial coordinate can yield prominent differences between the solvers mainly due to differences in the edge grid distribution. The simulations reveal that electron and ion solvers predict noticeably different density peaking for the same diffusion and pinch velocity while with the peaked profile of helium, expected in fusion reactors. The fuelling benchmarking demonstrates that gas puffing is not efficient for core fuelling in H-modes and density control should be done by the high field side pellet injection in contrast to present machines.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A Novel Approach to Fast Calculation of High-Order Q-Cumulants

The method of Q-cumulants has been shown as a powerful tool to study the fine details of the azimuthal anisotropies in high-energy nucleus-nucleus collisions. A new method for the fast calculation of arbitrary order Q-cumulant $v_{n}\{2k\}$ values, based on the partition of a non-negative integer l $\le$ m for calculation of the 2m-particle azimuthal correlations is presented in this paper. Unlike the standard Q-cumulants method in which the calculation of high-order multi-particle calculations is impractical, the newly proposed method enables easy calculation. The validity of the method is proven via a toy model that uses the elliptic power distribution to simulate anisotropic emission of particles. The method enables the study of fine details of the $v_{2}$ distribution, such as higher-order central moments of the $v_{2}$ distribution, as well as the hydrodynamic behavior of the Quark-Gluon Plasma.

FOS: Physical sciences↗
DOE OSTI2017

Determining the strengths of HCP slip systems using harmonic analyses of lattice strain distributions

A robust methodology is presented to extract slip system strengths from lattice strain distributions for polycrystalline samples obtained from high-energy x-ray diffraction (HEXD) experiments with in situ loading. The methodology consists of matching the evolution of coefficients of a harmonic expansion of the distributions from simulation to the coefficients derived from measurements. Simulation results are generated via finite element simulations of virtual polycrystals that are subjected to the loading history applied in the HEXD experiments. Advantages of the methodology include: (1) its ability to utilize extensive data sets generated by HEXD experiments; (2) its ability to capture trends in distributions that may be noisy (both measured and simulated); and (3) its sensitivity to the ratios of the family strengths. The approach is used to evaluate the slip system strengths of Ti-6Al-4V using samples having relatively equiaxed grains. These strength estimates are compared to values in the literature.

42 ENGINEERING↗
DOE OSTI2011

3D modeling of effects of increased oxygenation and activity concentration in tumors treated with radionuclides and antiangiogenic drugs

Purpose: Formation of new blood vessels (angiogenesis) in response to hypoxia is a fundamental event in the process of tumor growth and metastatic dissemination. However, abnormalities in tumor neovasculature often induce increased interstitial pressure (IP) and further reduce oxygenation (pO{sub 2}) of tumor cells. In radiotherapy, well-oxygenated tumors favor treatment. Antiangiogenic drugs may lower IP in the tumor, improving perfusion, pO{sub 2} and drug uptake, by reducing the number of malfunctioning vessels in the tissue. This study aims to create a model for quantifying the effects of altered pO{sub 2}-distribution due to antiangiogenic treatment in combination with radionuclide therapy. Methods: Based on experimental data, describing the effects of antiangiogenic agents on oxygenation of GlioblastomaMultiforme (GBM), a single cell based 3D model, including 10{sup 10} tumor cells, was developed, showing how radionuclide therapy response improves as tumor oxygenation approaches normal tissue levels. The nuclides studied were {sup 90}Y, {sup 131}I, {sup 177}Lu, and {sup 211}At. The absorbed dose levels required for a tumor control probability (TCP) of 0.990 are compared for three different log-normal pO{sub 2}-distributions: {mu}{sub 1} = 2.483, {sigma}{sub 1} = 0.711; {mu}{sub 2} = 2.946, {sigma}{sub 2} = 0.689; {mu}{sub 3} = 3.689, and {sigma}{sub 3} = 0.330. The normal tissue absorbed doses will, in turn, depend on this. These distributions were chosen to represent the expected oxygen levels in an untreated hypoxic tumor, a hypoxic tumor treated with an anti-VEGF agent, and in normal, fully-oxygenated tissue, respectively. The former two are fitted to experimental data. The geometric oxygen distributions are simulated using two different patterns: one Monte Carlo based and one radially increasing, while keeping the log-normal volumetric distributions intact. Oxygen and activity are distributed, according to the same pattern. Results: As tumor pO{sub 2} approaches normal tissue levels, the therapeutic effect is improved so that the normal tissue absorbed doses can be decreased by more than 95%, while retaining TCP, in the most favorable scenario and by up to about 80% with oxygen levels previously achieved in vivo, when the least favourable oxygenation case is used as starting point. The major difference occurs in poorly oxygenated cells. This is also where the pO{sub 2}-dependence of the oxygen enhancement ratio is maximal. Conclusions: Improved tumor oxygenation together with increased radionuclide uptake show great potential for optimising treatment strategies, leaving room for successive treatments, or lowering absorbed dose to normal tissues, due to increased tumor response. Further studies of the concomitant use of antiangiogenic drugs and radionuclide therapy therefore appear merited.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
NASA NTRS1994

Observations of stratospheric hydrogen fluoride by halogen occultation experiment (HALOE)

The Halogen Occultation Experiment (HALOE) Hydrogen Fluoride (HF) channel on the Upper Atmospheric Research Satellite (UARS) is providing the first global measurements of stratospheric HF, the dominant flourine reservoir in the atmosphere. This paper describes the latitudinal and seasonal variations of HALOE-observed HF in terms of vertical profiles, altitude/latitude cross sections, and column abundances. The HF global distribution shows a 'tracerlike' structure and its column amount increases with latitude, in agreement with previous aircraft measurements of the HF column amount. A comparison between the HALOE HF column above 20 km and the ATMOS 1985 measurements is used to estimate the annual rate of increase of stratospheric HF. Exponential rates of 4.9-6.6%/yr and linear growth rates of 6-8.6%/yr in 1985 and 4.3-5.5%/yr in 1992-1993 are found. HALOE HF measurements during the 1993 Antarctic spring are briefly described. This species behaves like a conserved tracer and its distribution shows an area of enhanced mixing ratios correlated with the polar vortex that has a clear latitude boundary. Finally, simulated HF distributions by the National Center for Atmospheric Research (NCAR) two-dimensional model are used to compare with HALOE observations of HF. Reasonable agreements in the global structure and the absolute amount of HF are found. The differences between the model and the observed results indicate the need for improving treatment of atmospheric dynamics and fluorine-related chemical parameters in the model simulations.

Luo, M.↗
DOE OSTI2011

Modeling pulsed-laser melting of embedded semiconductor nanoparticles

Pulsed-laser melting (PLM) is commonly used to achieve a fast quench rate in both thin films and nanoparticles. A model for the size evolution during PLM of nanoparticles confined in a transparent matrix, such as those created by ion-beam synthesis, is presented. A self-consistent mean-field rate equations approach that has been used successfully to model ion beam synthesis of germanium nanoparticles in silica is extended to include the PLM process. The PLM model includes classical optical absorption, multiscale heat transport by both analytical and finite difference methods, and melting kinetics for confined nanoparticles. The treatment of nucleation and coarsening behavior developed for the ion beam synthesis model is modified to allow for a non-uniform temperature gradient and for interacting liquid and solid particles with different properties. The model allows prediction of the particle size distribution after PLM under various laser fluences, starting from any particle size distribution including as-implanted or annealed simulated samples. A route for narrowing the size distribution of embedded nanoparticles is suggested, with simulated distribution widths as low as 15% of the average size.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Inner Belt Electron Decay Timescales: A Comparison of Van Allen Probes and DREAM3D Losses Following the June 2015 Storm

Abstract NASA's Van Allen Probes observed significant, long‐lived fluxes of inner belt electrons up to ∼1 MeV after geomagnetic storms in March and June 2015. Reanalysis of Magnetic Electron Ion Spectrometer (MagEIS) data with improved background correction showed a clearer picture of the relativistic electron population that persisted through 2016 and into 2017 above the Fennell et al. (2015, https://doi.org/10.1002/2014gl062874 ) limit. The intensity and duration of these enhancements allow estimation of decay timescales for comparison with simulated decay rates and theoretical lifetimes. We compare decay timescales from these data and DREAM3D simulations based on them using geomagnetic activity‐dependent pitch angle diffusion coefficients derived from plasmapause‐indexed wave data (Malaspina et al., 2016, https://doi.org/10.1002/2016gl069982 , 2018, https://doi.org/10.1029/2018gl078564 ) and phase space densities derived from MagEIS observations. Simulated decay rates match observed decay rates more closely than the theoretical lifetime due to significantly nonequilibrium pitch angle distributions in simulation and data. We conclude that nonequilibrium effects, rather than a missing diffusion or loss process, account for observed short decay rates.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

SU-E-T-496: Influence of the Intrafraction Patient Repositioning On Spot Scanning Proton Therapy for Prostate Cancer

Purpose To verify the effect of intrafraction patient repositioning based on real time image guidance on the proton dose distribution for prostate cancer. Methods Dose distributions were simulated using CT images of prostate cancer patients. Treatment plans were prepared using proton beam from left side of the patient with 3 and 6 mm lateral margins. During the beam delivery, prostate position was assumed to be shifted by 5 mm in anterior, posterior, superior and inferior directions. The dose distributions were calculated for two scenarios; (i) the patient position was readjusted by 5 mm just after the prostate shift occurs (ii) the patient position was not adjusted. Such a real time readjustment can be realized with real time image guidance. The differences of the values D{sub 5} and D{sub 95} in CTV (D{sub 5–95}) were compared between different sizes of the lateral margins and different scenarios. Results D{sub 5–95} values in the case of scenario (i) are smaller than scenario (ii). The differences between scenario (i) and (ii) were not depend on the directions of the prostate shift. In the case of scenario (ii), margin expansion makes D{sub 5–95} values small. In the case of scenario (i), margin size is uncorrelated to D{sub 5–95} values. Conclusion These results indicate that the intrafraction patient repositioning is effective to minimize dose deterioration when the prostate moves during the irradiation. In the case without intrafraction patient repositioning, margin expansion is effective to reduce dose deterioration. In the case with intrafraction patient repositioning, 3 mm margin is enough to minimize dose deterioration.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

SU-D-19A-03: Monte Carlo Investigation of the Mobetron to Perform Modulated Electron Beam Therapy

Purpose: Modulated electron radiotherapy (MERT) has been proposed as a mean of delivering conformal dose to shallow tumors while sparing distal structures and surrounding tissues. In intraoperative radiotherapy (IORT) utilizing Mobetron, an applicator is placed as closely as possible to the suspected cancerous tissues to be treated. In this study we investigate the characteristics of Mobetron electron beams collimated by an in-house prospective electron multileaf collimator (eMLC) and its feasibility for MERT. Methods: IntraOp Mobetron™ dedicated to perform radiotherapy during surgery was used in the study. It provides several energies (6, 9 and 12 MeV). Dosimetry measurements were performed to obtain percentage depth dose curves (PDD) and profiles for a 10-cm diameter applicator using the PTW MP3/XS 3D-scanning system and the semiflex ion chamber. MCBEAM/MCSIM Monte Carlo codes were used for the treatment head simulation and phantom dose calculation. The design of an electron beam collimation by an eMLC attached to the Mobetron head was also investigated using Monte Carlo simulations. Isodose distributions resulting from eMLC collimated beams were compared to that collimated using cutouts. The design for our Mobetron eMLC is based on our previous experiences with eMLCs designed for clinical linear accelerators. For Mobetron the eMLC is attached to the end of a spacer-mounted rectangular applicator at 50 cm SSD. Steel will be used as the leaf material because other materials would be toxic and will not be suitable for intraoperative applications. Results: Good agreement (within 2%) was achieved between measured and calculated PDD curves and profiles for all available energies. Dose distributiosn provided by the eMLC showed reasonable agreement (∼3%/1mm) with those obtained by conventional cutouts. Conclusion: Monte Carlo simulations are capable of modeling Mobetron electron beams with a reliable accuracy. An eMLC attached to the Mobteron treatment head will allow better treatment options with those machines.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Observed and Bin Model Simulated Evolution of Drop Size Distributions in High-Based Cumulus Congestus Over the United Arab Emirates

This study examines microphysical processes in developing high-based cumulus congestus over the United Arab Emirates using aircraft observations and a large-eddy-simulation model with bin microphysics. A notable feature of this case is the lack of mm-sized drops despite having a liquid cloud layer >3.5 km deep, contrasting with copious large drops observed in maritime tropical cumulus congestus having a similar vertical extent. Modeled drop size distributions are similar to observations at various temperatures between 9.5 and –12°C, including the lack of mm-sized drops. Cloud dilution leads to low-to-moderate liquid water contents (~0.5–1.5 g m –3 ) in most of the cloud core, several times smaller than adiabatic values. Dilution is enhanced in the inflowing branch of the toroidal circulations associated with individual cloud thermals, which are favored regions for secondary droplet activation (activation above cloud base). Secondary activation in general contributes substantially to the droplet population. Turning it off leads to a sharp decrease in droplet concentration and increase in mean size aloft, but does little to increase rain drop production. Warm rain generation (or lack thereof) in this case is therefore determined more by the sub-cloud aerosol and the cloud base droplet size distribution (DSD) than DSD evolution aloft from secondary droplet activation. Decreasing the aerosol concentration by a factor of 10 greatly increases production of large drops via collision-coalescence. Thus, despite its high base (low temperatures) and substantial dilution, the simulated cloud is thermodynamically and dynamically capable of rapidly producing copious mm-sized drops from collision-coalescence under pristine aerosol conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Three dimensional particle-in-cell simulations of electron beams created via reflection of intense laser light from a water target

We present 3D Particle-in-Cell (PIC) modeling of an ultra-intense laser experiment by the Extreme Light group at the Air Force Research Laboratory using the Large Scale Plasma (LSP) PIC code. This is the first time PIC simulations have been performed in 3D for this experiment which involves an ultra-intense, short-pulse (30 fs) laser interacting with a water jet target at normal incidence. The laser-energy-to-ejected-electron-energy conversion efficiency observed in 2D(3v) simulations were comparable to the conversion efficiencies seen in the 3D simulations, but the angular distribution of ejected electrons in the 2D(3v) simulations displayed interesting differences with the 3D simulations' angular distribution; the observed differences between the 2D(3v) and 3D simulations were more noticeable for the simulations with higher intensity laser pulses. An analytic plane-wave model is discussed which provides some explanation for the angular distribution and energies of ejected electrons in the 2D(3v) simulations. We also performed a 3D simulation with circularly polarized light and found a significantly higher conversion efficiency and peak electron energy, which is promising for future experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2012

DISTRIBUTION OF MAGNETIC DISCONTINUITIES IN THE SOLAR WIND AND IN MAGNETOHYDRODYNAMIC TURBULENCE

The statistical properties of magnetic discontinuities in the solar wind are investigated by measuring fluctuations in the magnetic field direction, given by the rotation {Delta}{theta} that the magnetic field vector undergoes during time interval {Delta}t. We show that the probability density function, P({Delta}{theta}), can be described by a simple model in which the magnetic field vector purely rotates with a relative increment {Delta}B/B that is lognormally distributed. We find that the probability density function of increments, P({Delta}B/B), has a remarkable scaling property: the normalized variable x = ({Delta}B/B) {center_dot} ({Delta}t/{Delta}t{sub 0}){sup -{alpha}} has a universal lognormal distribution for all time intervals {Delta}t. We then compare measurements from the solar wind with those from direct numerical simulations of magnetohydrodynamic (MHD) turbulence. We find good agreement for P({Delta}{theta}) obtained in the two cases when the magnetic guide field to fluctuations ratio B{sub 0}/b{sub rms} is chosen accordingly. However, the scale invariance of P({Delta}B/B) is broken in the MHD simulations with relatively limited inertial interval, which causes P({Delta}{theta}) to scale with measurement interval differently than in the solar wind.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

The Effect of the Velocity Distribution on Kilonova Emission

The electromagnetic emission from the nonrelativistic ejecta launched in neutron star mergers (either dynamically or through a disk wind) has the potential to probe both the total mass and composition of this ejecta. These observations are crucial in understanding the role of these mergers in the production of r-process elements in the Universe. However, many properties of the ejecta can alter the light curves and we must both identify which properties play a role in shaping this emission and understand the effects these properties have on the emission before we can use observations to place strong constraints on the amount of r-process elements produced in the merger. This paper focuses on understanding the effect of the velocity distribution (amount of mass moving at different velocities) for lanthanide-rich ejecta on the light curves and spectra. The simulations use distributions guided by recent calculations of disk outflows and compare the velocity-distribution effects to those of ejecta mass, velocity, and composition. Our comparisons show that uncertainties in the velocity distribution can lead to a factor of 2–4 uncertainties in the inferred ejecta mass based on peak infrared luminosities. We also show that early-time UV or optical observations may be able to constrain the velocity distribution, reducing the uncertainty in the ejecta mass.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

ON THE CONTRIBUTION OF FLUORESCENCE TO LYα HALOS AROUND STAR-FORMING GALAXIES

We quantify the contribution of Lyα fluorescence to observed spatially extended Lyα halos around Lyα emitters at redshift z = 3.1. The key physical quantities that describe the fluorescent signal include (i) the distribution of cold gas in the circumgalactic medium (CGM); we explore simple analytic models and fitting functions to recent hydrodynamical simulations; and (ii) local variations in the ionizing background due to ionizing sources that cluster around the central galaxy. We account for clustering by boosting the observationally inferred volumetric production rate of ionizing photons, ϵ{sub LyC}, by a factor of 1+ξ{sub LyC}(r), in which ξ{sub LyC}(r) quantifies the clustering of ionizing sources around the central galaxy. We compute ξ{sub LyC}(r) by assigning an “effective” bias parameter to the ionizing sources. This novel approach allows us to quantify our ignorance of the population of ionizing sources in a simple parametrized form. We find a maximum enhancement in the local ionizing background in the range 50–200 at r ∼ 10 physical kpc. For spatially uncorrelated ionizing sources and fluorescing clouds we find that fluorescence can contribute up to ∼ 50%–60% of the observed spatially extended Lyα emission. We briefly discuss how future observations can shed light on the nature of Lyα halos around star-forming galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

The Effect of the Velocity Distribution on Kilonova Emission

Abstract The electromagnetic emission from the nonrelativistic ejecta launched in neutron star mergers (either dynamically or through a disk wind) has the potential to probe both the total mass and composition of this ejecta. These observations are crucial in understanding the role of these mergers in the production of r -process elements in the Universe. However, many properties of the ejecta can alter the light curves and we must both identify which properties play a role in shaping this emission and understand the effects these properties have on the emission before we can use observations to place strong constraints on the amount of r -process elements produced in the merger. This paper focuses on understanding the effect of the velocity distribution (amount of mass moving at different velocities) for lanthanide-rich ejecta on the light curves and spectra. The simulations use distributions guided by recent calculations of disk outflows and compare the velocity-distribution effects to those of ejecta mass, velocity, and composition. Our comparisons show that uncertainties in the velocity distribution can lead to a factor of 2–4 uncertainties in the inferred ejecta mass based on peak infrared luminosities. We also show that early-time UV or optical observations may be able to constrain the velocity distribution, reducing the uncertainty in the ejecta mass.

79 ASTRONOMY AND ASTROPHYSICS↗