SEARCH · Engineering Papers

Results for “Computerized 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 901 records · Page 50

All sources NASA NTRS DOE OSTI

DOE OSTI2016

Parametric resonance after hilltop inflation caused by an inhomogeneous inflaton field

We study preheating after hilltop inflation where the inflaton couples to another scalar field, e.g. a right-handed sneutrino, which provides a mechanism for generating the correct initial conditions for inflation and also a decay channel for the inflaton that allows for reheating and non-thermal leptogenesis. In the presence of such a coupling, we find that after the phases of tachyonic preheating and tachyonic oscillations, during which the inflaton field becomes inhomogeneous, there can be a subsequent preheating phase where the fluctuations of the other field get resonantly enhanced, from initial vacuum fluctuations up to amplitudes of the same order (and even larger) as the ones of the inflaton field. This resonant enhancement differs from the usual parametric resonance as the inflaton field is inhomogeneous at the time the enhancement takes place. We study this effect using lattice simulations as well as semi-analytically with a generalized Floquet analysis for inhomogeneous background fields.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2016

Efficient calculation of cosmological neutrino clustering in the non-linear regime

We study in detail how neutrino perturbations can be followed in linear theory by using only terms up to l =2 in the Boltzmann hierarchy. We provide a new approximation to the third moment and demonstrate that the neutrino power spectrum can be calculated to a precision of better than ∼ 5% for masses up to ∼ 1 eV and k ∼< 10 h /Mpc. The matter power spectrum can be calculated far more precisely and typically at least a factor of a few better than with existing approximations. We then proceed to study how the neutrino power spectrum can be reliably calculated even in the non-linear regime by using the non-linear gravitational potential, sourced by dark matter overdensities, as it is derived from semi-analytic methods based on N -body simulations in the Boltzmann evolution hierarchy. Our results agree extremely well with results derived from N -body simulations that include cold dark matter and neutrinos as independent particles with different properties.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2014

SU-E-T-586: Field Size Dependence of Output Factor for Uniform Scanning Proton Beams: A Comparison of TPS Calculation, Measurement and Monte Carlo Simulation

Purpose: Output dependence on field size for uniform scanning beams, and the accuracy of treatment planning system (TPS) calculation are not well studied. The purpose of this work is to investigate the dependence of output on field size for uniform scanning beams and compare it among TPS calculation, measurements and Monte Carlo simulations. Methods: Field size dependence was studied using various field sizes between 2.5 cm diameter to 10 cm diameter. The field size factor was studied for a number of proton range and modulation combinations based on output at the center of spread out Bragg peak normalized to a 10 cm diameter field. Three methods were used and compared in this study: 1) TPS calculation, 2) ionization chamber measurement, and 3) Monte Carlos simulation. The XiO TPS (Electa, St. Louis) was used to calculate the output factor using a pencil beam algorithm; a pinpoint ionization chamber was used for measurements; and the Fluka code was used for Monte Carlo simulations. Results: The field size factor varied with proton beam parameters, such as range, modulation, and calibration depth, and could decrease over 10% from a 10 cm to 3 cm diameter field for a large range proton beam. The XiO TPS predicted the field size factor relatively well at large field size, but could differ from measurements by 5% or more for small field and large range beams. Monte Carlo simulations predicted the field size factor within 1.5% of measurements. Conclusion: Output factor can vary largely with field size, and needs to be accounted for accurate proton beam delivery. This is especially important for small field beams such as in stereotactic proton therapy, where the field size dependence is large and TPS calculation is inaccurate. Measurements or Monte Carlo simulations are recommended for output determination for such cases.

07 ISOTOPE AND RADIATION SOURCES↗

DOE OSTI2014

TH-C-12A-10: Surface Dose Enhancement Using Novel Hybrid Electron and Photon Low-Z Therapy Beams: Monte Carlo Simulation

Purpose: The introduction of the TrueBeam linac platform provides access to an in-air target assembly making it possible to apply novel treatments using multiple target designs. One such novel treatment uses multiple low-Z targets to enhance surface dose replacing the use of synthetic tissue equivalent material (bolus). This treatment technique will decrease the common dosimetric and set up errors prevalent in using physical treatment accessories like bolus. The groundwork for a novel treatment beam used to enhance surface dose to within 80-100% of the dose at dmax by utilizing low-Z (Carbon) targets of various percent CSDA range thickness operated at 2.5–4 MeV used in conjunction with a clinical 6 MV beam is presented herein. Methods: A standard Monte Carlo model of a Varian Clinac accelerator was developed to manufacturers specifications. Simulations were performed using Be, C, AL, and C, as potential low-Z targets, placed in the secondary target position. The results determined C to be the target material of choice. Simulations of 15, 30 and 60% CSDA range C beams were propagated through slab phantoms. The resulting PDDs were weighted and combined with a standard 6 MV treatment beam. Versions of the experimental targets were installed into a 2100C Clinac and the models were validated. Results: Carbon was shown to be the low-Z material of choice for this project. Using combinations of 15, 30, 60% CSDA beams operated at 2.5 and 4 MeV in combination with a standard 6 MV treatment beam the surface dose was shown to be enhanced to within 80–100% the dose at dmax. Conclusion: The modeled low-Z beams were successfully validated using machined versions of the targets. Water phantom measurements and slab phantom simulations show excellent correlation. Patient simulations are now underway to compare the use of bolus with the proposed novel beams. NSERC.

07 ISOTOPE AND RADIATION SOURCES↗

DOE OSTI2015

SU-E-T-313: Dosimetric Deviation of Misaligned Beams for a 6 MV Photon Linear Accelerator Using Monte Carlo Simulations

Purpose: To quantify the dosimetric variations of misaligned beams for a linear accelerator by using Monte Carlo (MC) simulations. Method and Materials: Misaligned beams of a Varian 21EX Clinac were simulated to estimate the dosimetric effects. All the linac head components for a 6 MV photon beam were implemented in BEAMnrc/EGSnrc system. For incident electron beam parameters, 6 MeV with 0.1 cm full-width-half-max Gaussian beam was used. A phase space file was obtained below the jaw per each misalignment condition of the incident electron beam: (1) The incident electron beams were tilted by 0.5, 1.0 and 1.5 degrees on the x-axis from the central axis. (2) The center of the incident electron beam was off-axially moved toward +x-axis by 0.1, 0.2, and 0.3 cm away from the central axis. Lateral profiles for each misaligned beam condition were acquired at dmax = 1.5 cm and 10 cm depth in a rectangular water phantom. Beam flatness and symmetry were calculated by using the lateral profile data. Results: The lateral profiles were found to be skewed opposite to the angle of the incident beam for the tilted beams. For the displaced beams, similar skewed lateral profiles were obtained with small shifts of penumbra on the +x-axis. The variations of beam flatness were 3.89–11.18% and 4.12–42.57% for the tilted beam and the translated beam, respectively. The beam symmetry was separately found to be 2.95 −9.93% and 2.55–38.06% separately. It was found that the percent increase of the flatness and the symmetry values are approximated 2 to 3% per 0.5 degree tilt or per 1 mm displacement. Conclusion: This study quantified the dosimetric effects of misaligned beams using MC simulations. The results would be useful to understand the magnitude of the dosimetric deviations for the misaligned beams.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

Ribbon electron beam formation by a forevacuum plasma electron source

Results of the numerical analysis and experimental research on ribbon electron beam generation based on hollow cathode discharge at forevacuum gas pressure are presented. Geometry of the accelerating gap has modified. It lets us focus the ribbon electron beam and to transport it on a distance of several tens of centimeters in the absence of an axial magnetic field. The results of numerical simulations are confirmed by the experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2017

Study of the development of relativistic plasma bunches in a long mirror trap by optical and X-ray imaging and numerical simulations

The work presents experimental results demonstrating the feasibility of autoresonance acceleration of electrons in a long mirror trap with a reverse magnetic field. It is shown that gyromagnetic autoresonance results in the formation of a plasma bunch with average electron energy of several hundred keV, which is confined for a long time in the trap. The results of computer simulations of the regime of reverse gyromagnetic autoresonance agree well with the experimental data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Multiscale thermal hydraulic coupling methods for boiling water reactor simulation

During the last two years, the Virtual Environment for Reactor Applications (VERA) has been extended to simulate boiling water reactors (BWRs). The thermal hydraulic effects present in BWRs are far more complex than those in pressurized water reactors (PWRs). Therefore, the runtime is significantly increase and convergence behavior worse compared to PWRs. Most of the additional runtime is spent during the first few coupled iterations, when the power shape is still rapidly evolving, dramatically affecting the thermal hydraulics (TH). To alleviate the increased computational expense, a multiscale TH coupling approach was developed in VERA in which a highly efficient, simplified TH model is solved for several coupled iterations until the power shape is partially converged. The simplified solution is used to inform the assembly-wise flow distribution in the high-fidelity TH solver CTF, reducing the amount of work required to properly balance the pressure drop in each channel and reducing runtime for the couple calculation. Because the final TH calculations are performed with CTF, there is ultimately no impact on the accuracy of the converged solution. This paper presents the details of this multiscale TH coupling approach, along with results for a 4x4 array of GE-14 fuel bundles and whole-core coupled simulations of the Hatch core. These two cases show that the multiscale approach dramatically reduces the runtime of coupled BWR simulations. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2015

Synchrotron characterization of nanograined UO 2 grain growth

This activity is supported by the US Nuclear Energy Advanced Modeling and Simulation (NEAMS) Fuels Product Line (FPL) and aims at providing experimental data for the validation of the mesoscale simulation code MARMOT. MARMOT is a mesoscale multiphysics code that predicts the coevolution of microstructure and properties within reactor fuel during its lifetime in the reactor. It is an important component of the Moose-Bison-Marmot (MBM) code suite that has been developed by Idaho National Laboratory (INL) to enable next generation fuel performance modeling capability as part of the NEAMS Program FPL. In order to ensure the accuracy of the microstructure based materials models being developed within the MARMOT code, extensive validation efforts must be carried out. In this report, we summarize our preliminary synchrotron radiation experiments at APS to determine the grain size of nanograin UO 2 . The methodology and experimental setup developed in this experiment can directly apply to the proposed in-situ grain growth measurements. The investigation of the grain growth kinetics was conducted based on isothermal annealing and grain growth characterization as functions of duration and temperature. The kinetic parameters such as activation energy for grain growth for UO 2 with different stoichiometry are obtained and compared with molecular dynamics (MD) simulations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2016

THE CARNEGIE HUBBLE PROGRAM: THE DISTANCE AND STRUCTURE OF THE SMC AS REVEALED BY MID-INFRARED OBSERVATIONS OF CEPHEIDS

Using Spitzer observations of classical Cepheids we have measured the true average distance modulus of the Small Magellanic Cloud (SMC) to be 18.96±0.01{sub stat}±0.03{sub sys} mag (corresponding to 62±0.3 kpc), which is 0.48 ± 0.01 mag more distant than the LMC. This is in agreement with previous results from Cepheid observations, as well as with measurements from other indicators such as RR Lyrae stars and the tip of the red giant branch. Utilizing the properties of the mid-infrared Leavitt Law we measured precise distances to individual Cepheids in the SMC, and have confirmed that the galaxy is tilted and elongated such that its eastern side is up to 20 kpc closer than its western side. This is in agreement with the results from red clump stars and dynamical simulations of the Magellanic Clouds and Stream.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2016

REFLECTION OF PROPAGATING SLOW MAGNETO-ACOUSTIC WAVES IN HOT CORONAL LOOPS: MULTI-INSTRUMENT OBSERVATIONS AND NUMERICAL MODELING

Slow MHD waves are important tools for understanding coronal structures and dynamics. In this paper, we report a number of observations from the X-Ray Telescope (XRT) on board HINODE and Solar Dynamic Observatory /Atmospheric Imaging Assembly (AIA) of reflecting longitudinal waves in hot coronal loops. To our knowledge, this is the first report of this kind as seen from the XRT and simultaneously with the AIA. The wave appears after a micro-flare occurs at one of the footpoints. We estimate the density and temperature of the loop plasma by performing differential emission measure (DEM) analysis on the AIA image sequence. The estimated speed of propagation is comparable to or lower than the local sound speed, suggesting it to be a propagating slow wave. The intensity perturbation amplitude, in every case, falls very rapidly as the perturbation moves along the loop and eventually vanishes after one or more reflections. To check the consistency of such reflection signatures with the obtained loop parameters, we perform a 2.5D MHD simulation, which uses the parameters obtained from our observation as inputs, and perform forward modeling to synthesize AIA 94 Å images. Analyzing the synthesized images, we obtain the same properties of the observables as for the real observation. From the analysis we conclude that a footpoint heating can generate a slow wave which then reflects back and forth in the coronal loop before fading. Our analysis of the simulated data shows that the main agent for this damping is anisotropic thermal conduction.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2016

SPATIAL DISTRIBUTION OF PAIR PRODUCTION OVER THE PULSAR POLAR CAP

Using an analytic, axisymmetric approach that includes general relativity, coupled to a condition for pair production deduced from simulations, we derive general results about the spatial distribution of pair-producing field lines over the pulsar polar cap. In particular, we show that pair production on magnetic field lines operates over only a fraction of the polar cap for an aligned rotator for general magnetic field configurations, assuming the magnetic field varies spatially on a scale that is larger than the size of the polar cap. We compare our result to force-free simulations of a pulsar with a dipole surface field and find excellent agreement. Our work has implications for first-principles simulations of pulsar magnetospheres and for explaining observations of pulsed radio and high-energy emission.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2017

Toward a Greater Understanding of the Reduction of Drift Coefficients in the Presence of Turbulence

Drift effects play a significant role in the transport of charged particles in the heliosphere. A turbulent magnetic field is also known to reduce the effects of particle drifts. The exact nature of this reduction, however, is not clear. This study aims to provide some insight into this reduction and proposes a relatively simple, tractable means of modeling it that provides results in reasonable agreement with numerical simulations of the drift coefficient in a turbulent magnetic field.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2017

Buried AGNs in Advanced Mergers: Mid-infrared Color Selection as a Dual AGN Candidate Finder

A direct consequence of hierarchical galaxy formation is the existence of dual supermassive black holes, which may be preferentially triggered as active galactic nuclei (AGNs) during galaxy mergers. Despite decades of searching, however, dual AGNs are extremely rare, and most have been discovered serendipitously. Using the all-sky WISE survey, we identified a population of over 100 morphologically identified interacting galaxies or mergers that display red mid-infrared colors often associated in extragalactic sources with powerful AGNs. The vast majority of these advanced mergers are optically classified as star-forming galaxies, which suggests that they may represent an obscured population of AGNs that cannot be found through optical studies. In this work, we present Chandra/ACIS observations and near-infrared spectra with the Large Binocular Telescope of six advanced mergers with projected pair separations less than ∼10 kpc. The combined X-ray, near-infrared, and mid-infrared properties of these mergers provide confirmation that four out of the six mergers host at least one AGN, with four of the mergers possibly hosting dual AGNs with projected separations less than ∼10 kpc, despite showing no firm evidence for AGNs based on optical spectroscopic studies. Our results demonstrate that (1) optical studies miss a significant fraction of single and dual AGNs in advanced mergers, and (2) mid-infrared pre-selection is extremely effective in identifying dual AGN candidates in late-stage mergers. Our multi-wavelength observations suggest that the buried AGNs in these mergers are highly absorbed, with intrinsic column densities in excess of ∼N{sub H}>10{sup 24} cm{sup −2}, consistent with hydrodynamic simulations.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2019

A Study into the Benefits of Using Discrete Event Simulation in TRU and Mixed Low-Level Waste (MLLW) Disposal at Advanced Mixed Waste Treatment Plant (AMWTP) - 19536

As part of the 1995 Settlement agreement reached between the Idaho State Government and the Idaho National laboratory, it was agreed that the approximately 65,000 cubic metres of waste stored at the AMWTP site was to be shipped out of the state to the waste isolation pilot plant (WIPP) in New Mexico or a similar disposal facility by the end of December 2018. In June 2016, the operation of AMWTP was taken over by Fluor Idaho, a subsidiary of Fluor Inc. To improve their capability of processing all of the varied waste contained on the extensive site, Fluor Idaho requested the assistance of DBD Ltd. We were tasked with producing an all-encompassing operational research (OR) model/discrete event simulation (DES) model of the site, which is able to produce live end-date predictions based on current inventory and waste tracking against milestones. The DES model, is capable of predicting years worth of plant operation in a matter of minutes including the variation caused by such complexities as resource/staff availability, equipment breakdown and maintenance, prioritisation of different waste types, and varying shift patterns of both staff and facilities. The DES model can be run from its own bespoke dashboard, that can be customised to suit the audience it is being used to inform, producing clear key performance indicator (KPI) graphics and high-level overviews for senior management, down to drum by drum breakdowns of the predicted feed plans for each individual facility on the AMWTP site for a work planner to analyse and compare to the current operations. Simulation can produce these types of outputs for very low cost and low risk, as it avoids experimentation to find potential alternatives to the status-quo that could result in loss of production or time. These methods can also be used to produce cost benefit analyses for suggested processing options feeding into optioneering studies and cost reduction strategies. As a tool for plant optimisation and tracking, discrete event simulation can provide detailed results quickly and reliably. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2017

A dynamic model of once-through steam generator (OTSG) for prototypical advanced reactor - 374

new design of Small Modular Reactor, the Prototypical Advanced Reactor (PAR), features two sodium primary systems each connected to dedicated steam generators (SGs) and sharing a balance of plant (BOP) system. PAR adopts once-through steam generators (OTSG) for its compact structure. The existence of superheated steam in a OTSG improves turbine design by eliminating the need of bulky moisture separator in SG. The two cores and SGs' design enable PAR to adjust power flexibly, as a result, introduce more challenges to the control. Therefore, a fast simulation model with high accuracy OTSG model is needed. In this paper, a dynamic model is built for OTSG of PAR based on the conservation of mass, momentum, and energy. The OTSG is simulated as lumped pipe and divided into subcooled, evaporating, and superheated regions by moving boundaries. SIMULINK is adopted to this OTSG simulation to be merged into existing PAR system models easily. The steady state results show that key parameters are rational. The 5% step increase of feedwater enthalpy transient demonstrates that the model is capable to respond to the perturbations. Since the basic laws and algebraic method used in this paper are suitable for other fluids, the model can be easily applied to any other OTSG. Future work will also focus on model optimization, merging and improved control. The final aim is to perform full scope realtime simulations for PAR under normal operational transients. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2016

Modeling the Resilience of Nuclear Power Plant for Unexpected Reactor Trip

Resilience can be defined as the intrinsic ability of a system to adjust to its functioning prior to, during, or following changes and disturbances, so that it can sustain required operations under both expected and unexpected conditions. Resilience engineering is a relatively new paradigm for safety management that focuses on how to cope with complexity under pressure or disturbance to achieve success. A resilient system combines inherent human abilities with engineered safety features and organizational features in a controlled and flexible manner both in expected and in unexpected situations. The concept of resilience has been suggested to complement conventional safety analyses. Conventional safety analysis methods such as PSA have several limitations: 1) they primarily focus on technical dimension, 2) the analysis are linear and sequential, 3) they are dominated by static models, 4) they do not take a systemic view into account, and 5) they focus primarily on why accidents happen and not how success is achieved. The concept of resilience provides a more integrated view for the safety analysis in safety-critical systems, e.g., nuclear power plants. The concept of resilience in the organization of NPPs has been highlighted since the Fukushima Daiichi Nuclear Power Plant accident. An IAEA report addresses that a resilient organization is one that quickly realizes deviation from normal operations and has the ability to make even the toughest and least popular decisions and to manage the margins in which it can maneuver. While traditional strategy to safety in NPPs is to identify what could go or has gone wrong, the concept of resilience is focused on what the organization does well and what it does to maintain its successful operation and to capitalize on those processes in the event of an unexpected situation. This study attempts to model the resilience of NPP dynamically, based on the operating experiences in Korean NPPs. First, this study revised the characterization of NPP. Then, a survey on the literature was conducted to identify the influences between the elements in the model of resilience in NPPs. These works provide a qualitative model for the resilience in NPPs. Then, this study reviewed 168 licensee events reports submitted to the Korean regulatory body since 2004 and quantified the relation between the attributes in the model by performing a Pearson correlation analysis. Finally, a few simulations using the System Dynamics were conducted to discuss which factors are most influencing the resilience. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2014

Uncertainty Quantification for Monte Carlo simulations - Paper 18

Uncertainty quantification (UQ) generally speaking refers to a wide variety of specific problems, also pertaining different scientific areas: then studies on Uncertainty Quantification involve many different methodologies, techniques and even terminologies: to be specific we are interested in studying how uncertainties on input physical data needed for Monte Carlo simulations transfer into uncertainties on the physical output of the simulation. We mention that other sources of uncertainty are be present in a simulation environment: they can be ascribed either to model inadequacy (for instance when one neglects atomic lattice at very low energies) or to the numerical procedure involved: for the case of MC simulation the latter corresponds to the statistical error affecting the simulation itself. So in general parameter uncertainties are mixed with other kind of uncertainties. The problem then remains seemingly quite complicated, even in this more specific domain: moreover we should better specify the meaning of the rather undetermined words physical data needed for Monte Carlo simulation. In fact they may refer both those data needed to specify the experimental setup that is going to be simulated - like, e.g., geometrical or composition data, but including also some external conditions applied to the system, like temperature, pressure, electromagnetic fields and so on - and physical data - cross sections or often physical models - required to describe particle transport. All these data are generally affected by uncertainties, which through the process of simulation transfer into uncertainty on its physical outcome. Then, from the MC user's own perspective, input data refers both to those that are under his direct experimental control and to many others, whose values and uncertainties will affect the result of the simulation. In this contribution we will summarize previous results, together with some new ones. (authors)

61 RADIATION PROTECTION AND DOSIMETRY↗