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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

MODELING THE FREQUENCY DEPENDENCE OF RADIO BEAMS FOR CONE-DOMINANT PULSARS

Beam radii for cone-dominant pulsars follow a power-law relation with frequency, thetav = ({nu}/{nu}{sub 0}) {sup k} + thetav{sub 0}, which has not been well explained in previous works. We study this frequency dependence of beam radius (FDB) for cone-dominant pulsars by using the curvature radiation mechanism. Considering various density and energy distributions of particles in the pulsar open field-line region, we numerically simulate the emission intensity distribution across emission height and rotation phase, get integrated profiles at different frequencies, and obtain the FDB curves. For the density model of a conal-like distribution, the simulated profiles always shrink to one component at high frequencies. In the density model with two separated density patches, the profiles generally have two distinct components, and the power-law indices k are found to be in the range from -0.1 to -2.5, consistent with observational results. Energy distributions of streaming particles have significant influence on the frequency-dependence behavior. Radial energy decay of particles is desired to get proper thetav{sub 0} in models. We conclude that by using the curvature radiation mechanism, the observed FDB for the cone-dominant pulsars can only be explained by the emission model of particles in two density patches with a Gaussian energy distribution and a radial energy loss.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2011

Nanoscale silicon-on-insulator deformation induced by stressed liner structures

Rotation and strain fields were mapped across silicon-on-insulator (SOI) regions induced by overlying stressed Si{sub 3}N{sub 4} features using x-ray nanobeam diffraction. The distribution in SOI tilt exhibited an antisymmetric distribution with a maximum magnitude of 7.9 milliradians, representing one of the first direct measurements of the lattice tilt conducted in situ within buried layers using a spot size of less than 100 nm. The measured rotation distribution corresponds to simulated values generated by boundary element method modeling, indicating that the strain transfer into the underlying SOI primarily induces elastic deformation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

MO-E-18C-02: Hands-On Monte Carlo Project Assignment as a Method to Teach Radiation Physics

Purpose: To present a hands-on project on Monte Carlo methods (MC) recently added to the curriculum and to discuss the students' appreciation. Methods: Since 2012, a 1.5 hour lecture dedicated to MC fundamentals follows the detailed presentation of photon and electron interactions. Students also program all sampling steps (interaction length and type, scattering angle, energy deposit) of a MC photon transport code. A handout structured in a step-by-step fashion guides student in conducting consistency checks. For extra points, students can code a fully working MC simulation, that simulates a dose distribution for 50 keV photons. A kerma approximation to dose deposition is assumed. A survey was conducted to which 10 out of the 14 attending students responded. It compared MC knowledge prior to and after the project, questioned the usefulness of radiation physics teaching through MC and surveyed possible project improvements. Results: According to the survey, 76% of students had no or a basic knowledge of MC methods before the class and 65% estimate to have a good to very good understanding of MC methods after attending the class. 80% of students feel that the MC project helped them significantly to understand simulations of dose distributions. On average, students dedicated 12.5 hours to the project and appreciated the balance between hand-holding and questions/implications. Conclusion: A lecture on MC methods with a hands-on MC programming project requiring about 14 hours was added to the graduate study curriculum since 2012. MC methods produce “gold standard” dose distributions and slowly enter routine clinical work and a fundamental understanding of MC methods should be a requirement for future students. Overall, the lecture and project helped students relate crosssections to dose depositions and presented numerical sampling methods behind the simulation of these dose distributions. Research funding from governments of Canada and Quebec. PP acknowledges partial support by the CREATE Medical Physics Research Training Network grant of the Natural Sciences and Engineering Research Council (Grant number: 432290)

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2023

Structural water in amorphous carbonate minerals: ab initio molecular dynamics simulations of X-ray pair distribution experiments

Water is known to play a controlling role in directing mineralization pathways and stabilizing metastable amorphous intermediates in hydrous carbonate mineral MCO 3 ·nH 2 O systems, where M 2+ is a divalent metal cation. Despite this recognition, the nature of the controls on crystallization is poorly understood, largely owing to the difficulty in characterizing the dynamically disordered structures of amorphous intermediates at the atomic scale. Here, we present a series of atomistic models, derived from ab initio molecular dynamics simulation, across a range of experimentally relevant cations (M = Ca, Mg, Sr) and hydration levels (0 ≤ n ≤ 2). Theoretical simulations of the dependence of the X-ray pair distribution function on the hydration level n show good agreement with available experimental data and thus provide further evidence for a lack of significant nanoscale structure in amorphous carbonates. Upon dehydration, the metal coordination number does not change significantly, but the relative extent of water dissociation increases, indicating that a thermodynamic driving force exists for water dissociation to accompany dehydration. Mg strongly favors monodentate conformation of carbonate ligands and shows a marked preference to exchange monodentate carbonate O for water O upon hydration, whereas Ca and Sr exchange mono- and bidentate carbonate ligands with comparable frequency. Finally, water forms an extensive hydrogen bond network among both water and carbonate groups that exhibits frequent proton transfers for all three cations considered suggesting that proton mobility is likely predominantly due to water dissociation and proton transfer reactions rather than molecular water diffusion.

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

Evaluation of Cloud-Resolving Model Intercomparison Simulations Using TWP-ICE Observations: Precipitation and Cloud Structure

The Tropical Warm Pool – International Cloud Experiment (TWP-ICE) provided high quality model forcing and observational datasets through which detailed model and observational intercomparisons could be performed. In this first of a two part study, precipitation and cloud structures within nine cloud-resolving model simulations are compared with scanning radar reflectivity and satellite infrared brightness temperature observations during an active monsoon period from 19 to 25 January 2006. Most simulations slightly overestimate volumetric convective rainfall. Overestimation of simulated convective area by 50% or more in several simulations is somewhat offset by underestimation of mean convective rain rates. Stratiform volumetric rainfall is underestimated by 13% to 53% despite overestimation of stratiform area by up to 65% because stratiform rain rates in every simulation are much lower than observed. Although simulations match the peaked convective radar reflectivity distribution at low levels, they do not reproduce the peaked distributions observed above the melting level. Simulated radar reflectivity aloft in convective regions is too high in most simulations. In stratiform regions, there is a large spread in model results with none resembling observed distributions. Above the melting level, observed radar reflectivity decreases more gradually with height than simulated radar reflectivity. A few simulations produce unrealistically uniform and cold 10.8-μm infrared brightness temperatures, but several simulations produce distributions close to observed. Assumed ice particle size distributions appear to play a larger role than ice water contents in producing incorrect simulated radar reflectivity distributions aloft despite substantial differences in mean graupel and snow water contents across models.

Varble, Adam↗
DOE OSTI2024

Numerical study of implosion instability mitigation in magnetically driven solid liner dynamic screw pinches

Magnetized liner inertial fusion experiments on the Z accelerator suffer from magneto-Rayleigh–Taylor instabilities (MRTI) that compromise integrity of the imploding cylindrical liner, limiting achievable fusion fuel conditions and ultimately reducing magneto-inertial fusion target performance. Dynamic screw pinches (DSP) provide a method to reduce MRTI in-flight via application of magnetic field line tension to the imploding liner outer surface. In contrast with z-pinches that drive implosions with an azimuthal magnetic field, dynamic screw pinches enforce an additional axial drive magnetic field component, making the overall drive magnetic field helical. As the liner implodes, cumulative MRTI development is reduced by dynamically shifting the orientation of the fastest growing instability modes. Three-dimensional magnetohydrodynamic simulations show that the DSP mechanism effectively stabilizes initially solid cylindrical liner implosions driven by Z-scale current pulses, indicating that MRTI mitigation increases with the ratio of axial to azimuthal drive magnetic field components (i.e., the drive field ratio). We also performed a spectral analysis of the simulated imploding density distributions, extracting wavelength and pitch angle of the simulated MRTI structures to study their dynamics during the implosion. Simulations of liners initially perturbed with drive-field-aligned sinusoidal structures indicate that MRTI mitigation in DSP implosions decreases with perturbation wavelength, once again suggestive of magnetic field line tension effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Development of Estimation Method for Radioactivity Concentration Distribution of Radioactive Wastes in a Box-shaped Container - 17034

We are developing a new estimation method suitable for determining radioactivity concentration distribution of radioactive wastes stored in a square box-shaped container. This method uses the Maximum A Posteriori (MAP) method which is a way to solve the inverse problem, and it estimates the radioactivity concentration distribution by comparing simulated data with measured data. In the measurement, gamma rays emitted from the radioactive wastes are measured by several gamma ray spectrometers which are set at different positions in relation to each other. On the other hand, the gamma ray spectrum data of each spectrometer are simulated by a Monte Carlo simulation code. The method uses an iterative calculation to solve the concentration distribution with measured spectrum data and simulated spectrum data. We conducted a simulation test to check the algorithm of our estimation method. And we also conducted an experiment using activated objects. We found that the concentration of the wastes in a box-shaped container could be estimated within a 20 % error for the simulation, and within a 40 % error for the experiment and we confirmed that feasibility of our method was favorable. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2015

Design and performance testing of a novel far field Gamma system to assay radioactive waste in 400 and 200 litre drums - 15339

A Far Field Gamma ray Measurement system has been developed for measuring and assaying both 400 and 200 litre waste drums with a wide range of activity at nuclear power stations. These include radioactive waste drums with surface dose-rates of less than 2 mSv/h (200 mrem/h) and, at the higher range of activity, 400 litre drums with surface dose-rates of up to 50 mSv/h (5 rem/h). The assay system incorporates an electro-mechanically cooled high purity Germanium profile detector and advanced digital spectroscopic analysis electronics. A novel optimised conical-trapezoidal collimator has been incorporated to reduce the effect of the radioactive background in addition to the conventional detector lead shield with tin and copper graded lining to reduce the effect of lead X-rays. Two features have been incorporated to reduce detector dead time and maintain detector resolution when measuring high activity drums. The first is the inclusion of two automated tungsten filters of different thickness as an integral part of the detector collimator. The filters are used in conjunction with an automated rail system on which the detector platform is mounted. The detector can be positioned at different distances from the drum surface in order to reduce detector dead time. The deployment of the filters and the adjustment of the drum - detector position are automatically controlled based on a user defined dose-rate limit table. The data is provided by two Geiger-Muller dose-rate detectors, which are used to measure the drum surface dose-rate. This data is subsequently used to adjust the detector position and in cases of higher dose-rate deploy one or other of the tungsten filters. The drum rotation platform incorporates a load cell to determine drum weight and hence the drum density. This information is used to calculate an attenuation correction based on the assumption of uniform drum density and uniform distribution of activity within the drum. Where the waste drums have a more complicated and known regular internal structure, such as a small radioactive region surrounded by an annulus of shielding material, the analysis algorithms incorporated into the far field geometry and attenuation correction spectroscopic data analysis code are able to make a more accurate determination of drum activity than would be possible assuming a completely uniform drum. The detector energy and efficiency calibration is achieved using a point source with multiple gamma ray energy peaks. The response of the high purity Germanium detector is also modelled using the Monte Carlo Neutron Photon code and the model is benchmarked using the measured point source calibration data. The performance of the instrument has been determined for a range of uniform density matrix drums containing a set of Eu-152 line or rod sources located in re-entrant tubes positioned on an equal volume basis. When rotated, the test drums with volume distributed line sources simulate waste drums of different density each with a uniformly distributed radioactive source. The measurements have been validated using Monte Carlo simulations. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2013

Electromagnetic Simulation of the Near-Field Distribution around a Wind Farm

An efficient approach to compute the near-field distribution around and within a wind farm under plane wave excitation is proposed. To make the problem computationally tractable, several simplifying assumptions are made based on the geometry problem. By comparing the approximations against full-wave simulations at 500 MHz, it is shown that the assumptions do not introduce significant errors into the resulting near-field distribution. The near fields around a 3 × 3 wind farm are computed using the developed methodology at 150 MHz, 500 MHz, and 3 GHz. Both the multipath interference patterns and the forward shadows are predicted by the proposed method.

Yang, Shang-Te↗
DOE OSTI2011

Fast-ions on Alcator C-Mod: Comparisons between Simulation and Experiment for Equilibrium and Evolving Distributions

Analysis of fast-ion distributions >100 T{sub e} is important for both ITER and a fusion reactor, as both will have a significant population of fusion-born alpha particles. Furthermore, energetic ions may display transport properties that differ from the bulk plasma. It is imperative to benchmark current simulation codes with measurements of highly energetic fast-ions on current devices. Experimental measurements of ICRF heated H-minority plasmas on Alcator C-Mod are performed with an upgraded Compact Neutral Particle Analyzer (CNPA). The CNPA consists of vertically-viewing photodiode detectors with viewing chords from r/a = {approx}0.18 to 0.6 and energies from 150 keV to 1.5 MeV. The CNPA measurements are compared to the output of a synthetic diagnostic that has been added to the Fokker-Planck solver CQL3D. CQL3D is coupled with the full-wave solver AORSA to evolve a self-consistent fast-ion distribution. Good agreement is obtained between simulation and experiment for fast ion distributions that are in steady-state (df/dt = 0). However, there are discrepancies between the experimental and simulated results for the time-dependent evolution of the fast-ion distribution. A possible explanation of the discrepancy is due to violation of certain assumptions of quasi-linear theory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Impact of bar resonances in the velocity–space distribution of the solar neighbourhood stars in a self-consistent N -body Galactic disc simulation

ABSTRACT The velocity–space distribution of the solar neighbourhood stars shows complex substructures. Most of the previous studies use static potentials to investigate their origins. Instead we use a self-consistent N-body model of the Milky Way, whose potential is asymmetric and evolves with time. In this paper, we quantitatively evaluate the similarities of the velocity–space distributions in the N-body model and that of the solar neighbourhood, using Kullback–Leibler divergence (KLD). The KLD analysis shows the time evolution and spatial variation of the velocity–space distribution. The KLD fluctuates with time, which indicates the velocity–space distribution at a fixed position is not always similar to that of the solar neighbourhood. Some positions show velocity–space distributions with small KLDs (high similarities) more frequently than others. One of them locates at $(R,\phi)=(8.2\,\,\rm{\mathrm{kpc}}, 30^\circ)$, where R and ϕ are the distance from the galactic centre and the angle with respect to the bar’s major axis, respectively. The detection frequency is higher in the inter-arm regions than in the arm regions. In the velocity maps with small KLDs, we identify the velocity–space substructures, which consist of particles trapped in bar resonances. The bar resonances have significant impact on the stellar velocity–space distribution even though the galactic potential is not static.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Probabilistic Cosmological Mass Mapping from Weak Lensing Shear

We infer gravitational lensing shear and convergence fields from galaxy ellipticity catalogs under a spatial process prior for the lensing potential. We demonstrate the performance of our algorithm with simulated Gaussian-distributed cosmological lensing shear maps and a reconstruction of the mass distribution of the merging galaxy cluster Abell 781 using galaxy ellipticities measured with the Deep Lens Survey. Given interim posterior samples of lensing shear or convergence fields on the sky, we describe an algorithm to infer cosmological parameters via lens field marginalization. In the most general formulation of our algorithm we make no assumptions about weak shear or Gaussian-distributed shape noise or shears. Because we require solutions and matrix determinants of a linear system of dimension that scales with the number of galaxies, we expect our algorithm to require parallel high-performance computing resources for application to ongoing wide field lensing surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Linking Deep and Shallow Convective Mass Fluxes via an Assumed Entrainment Distribution in CAM5‐CLUBB: Parameterization and Simulated Precipitation Variability

Abstract We modify the Zhang‐McFarlane (ZM) deep convection scheme in the Community Atmosphere Model version 5 to couple it with a unified parameterization for boundary‐layer turbulence and shallow convection, that is, Cloud Layers Unified by Binormals (CLUBB). By assuming a lognormal distribution of entrainment rate across the entire moist convective regimes, we link mass fluxes between shallow and deep convection, which are partitioned by the entrainment rate of the shallowest deep convective plume. Hence, a new deep convective closure is established which is coupled to the sub‐grid vertical motion variability in CLUBB. The convection feedback (or memory) effects are also considered to decrease the entrainment spectrum width and enhance the vertical velocity variability that further affect deep convection. Results show that the revised scheme improves the precipitation simulations in terms of the mean state and variability at various timescales, such as the alleviated double‐intertropical convergence zone and more realistic simulations of the seasonal variation of monsoon precipitation over East Asia, Madden‐Julian Oscillation, and precipitation diurnal phase propagations downstream of large terrains. The improvements are still seen in many aspects such as the mean‐state precipitation when turning off the convection feedback impacts in the revised scheme, emphasizing the benefits of using the modified mass‐flux closure. However, the convection feedbacks have considerable effects on the precipitation diurnal cycle simulations over regions with late‐afternoon precipitation peaks. Overall, the revised scheme provides a unified treatment for sub‐grid vertical motions across regimes of boundary‐layer turbulence, shallow convection, and deep convection, leading to better‐simulated precipitation at various timescales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

From Weakly Chaotic Dynamics to Deterministic Subdiffusion via Copula Modeling

Copula modeling consists in finding a probabilistic distribution, called copula, whereby its coupling with the marginal distributions of a set of random variables produces their joint distribution. The present work aims to use this technique to connect the statistical distributions of weakly chaotic dynamics and deterministic subdiffusion. More precisely, we decompose the jumps distribution of Geisel–Thomae map into a bivariate one and determine the marginal and copula distributions respectively by infinite ergodic theory and statistical inference techniques. We verify therefore that the characteristic tail distribution of subdiffusion is an extreme value copula coupling Mittag–Leffler distributions. We also present a method to calculate the exact copula and joint distributions in the case where weakly chaotic dynamics and deterministic subdiffusion statistical distributions are already known. Numerical simulations and consistency with the dynamical aspects of the map support our results.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2012

Electron scattering in helium for Monte Carlo simulations

An analytical approximation for differential cross-section of electron scattering on helium atoms is introduced. It is intended for Monte Carlo simulations, which, instead of angular distributions based on experimental data (or on first-principle calculations), usually rely on approximations that are accurate yet numerically efficient. The approximation is based on the screened-Coulomb differential cross-section with energy-dependent screening. For helium, a two-pole approximation of the screening parameter is found to be highly accurate over a wide range of energies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Impact of GaN cap on charges in Al₂O₃/(GaN/)AlGaN/GaN metal-oxide-semiconductor heterostructures analyzed by means of capacitance measurements and simulations

Oxide/semiconductor interface trap density (D{sub it}) and net charge of Al₂O₃/(GaN)/AlGaN/GaN metal-oxide-semiconductor high-electron mobility transistor (MOS-HEMT) structures with and without GaN cap were comparatively analyzed using comprehensive capacitance measurements and simulations. D{sub it} distribution was determined in full band gap of the barrier using combination of three complementary capacitance techniques. A remarkably higher D{sub it} (∼5–8 × 10¹²eV⁻¹ cm⁻²) was found at trap energies ranging from E C -0.5 to 1 eV for structure with GaN cap compared to that (D{sub it} ∼ 2–3 × 10¹²eV⁻¹ cm⁻²) where the GaN cap was selectively etched away. D{sub it} distributions were then used for simulation of capacitance-voltage characteristics. A good agreement between experimental and simulated capacitance-voltage characteristics affected by interface traps suggests (i) that very high D{sub it} (>10¹³eV⁻¹ cm⁻²) close to the barrier conduction band edge hampers accumulation of free electron in the barrier layer and (ii) the higher D{sub it} centered about E C -0.6 eV can solely account for the increased C-V hysteresis observed for MOS-HEMT structure with GaN cap. Analysis of the threshold voltage dependence on Al₂O₃ thickness for both MOS-HEMT structures suggests that (i) positive charge, which compensates the surface polarization, is not necessarily formed during the growth of III-N heterostructure, and (ii) its density is similar to the total surface polarization charge of the GaN/AlGaN barrier, rather than surface polarization of the top GaN layer only. Some constraints for the positive surface compensating charge are discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Measurement of muon-induced high-energy neutrons from rock in an underground Gd-doped water detector

In this paper, we present a measurement of the rate of correlated neutron captures in the WATCHBOY detector, deployed at a depth of approximately 390 meters water equivalent (m.w.e.) in the Kimballton Underground Research Facility (KURF). WATCHBOY consists of a cylindrical 2 ton water target doped with 0.1% gadolinium, surrounded by a 40 ton undoped water hermetic shield. We present a comparison of our results with the expected rate of correlated neutron captures arising from high-energy neutrons incident on the outside of the WATCHBOY shield, predicted by a hybrid FLUKA/GEANT4-based simulation. The incident neutron energy distribution used in the simulation was measured by a fast neutron spectrometer, the 1.8-ton Multiplicity and Recoil Spectrometer (MARS) detector, at the same depth. We nd that the measured detection rate of two correlated neutrons is consistent with that predicted by simulation. The result lends additional confidence in the detection technique used by MARS, and therefore in the MARS spectra as measured at three different depths. Confirmation of the fast neutron flux and spectrum is important as it helps validate the scaling models used to predict the fast neutron fluxes at different overburdens.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2019

Bragg Deflectors of Wave Fluxes for High-Power Relativistic Masers

Quasi-optical Bragg structures that represent sections of planar waveguides with tilted (with respect to the beam propagation direction) corrugation are studied. It is shown that such structures may serve as efficient deflectors for high-power wave fluxes and, hence, can be used for separation of microwave radiation and electron beam in relativistic masers. The corrugation configuration is optimized to increase the efficiency of transformation of wave fluxes to the transverse direction and improve uniformity of the spatial distribution of scattered radiation. The simulated results are verified using cold electrodynamic tests.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗