Engineering Papers⌕ Search

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.

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 OSTI2016

SU-F-T-322: A Comparison of Two Si Detectors for in Vivo Dosimetry

Purpose: To compare two types of semiconductor detectors for in vivo dosimetry by their dependence from various parameters in different conditions. Methods: QED yellow (Sun Nuclear) and EDP (Scanditronix) Si detectors were radiated by a Varian Clinac 2300 ix with 6 and 18 MV energies. 10 cm thickness water equivalent phantom consisted of 30×30 cm{sup 2} squared plates was used for experiments. Dose dependencies for different beam angles (0 – 180°), field size (3–40 cm), dose (50 – 300 MU), and dose rates (50 – 300 MU/min) were obtained and calibrated with Standard Farmer chamber (PTW). Results: Reproducibility, linearity, dose rate, angular dependence, and field size dependence were obtained for QED and EDP. They show no dose-rate dependence in available clinical dose rate range (100–600 MU/min). Both diodes have linear dependence with increasing the dose. Therefore even in case of high radiation therapy (including total body irradiation) it is not necessary to apply an additional correction during in vivo dosimetry. The diodes have different behavior for angular and field size dependencies. QED diode showed that dose value is stable for beam angles from 0 to 60°, for 60–180° correction factor has to be applied for each beam angle during in vivo measurements. For EDP diode dose value is sensitive to beam angle in whole range of angles. Conclusion: The study shows that QED diode is more suitable for in vivo dosimetry due to dose value independence from incident beam angle in the range 0–60°. There is no need in correction factors for increasing of dose and dose rate for both diodes. The next step will be to carry out measurements in non-standard conditions of total body irradiation. After this modeling of these experiments with Monte Carlo simulation for comparison calculated and obtained data is planned.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

TU-FG-BRB-07: GPU-Based Prompt Gamma Ray Imaging From Boron Neutron Capture Therapy

Purpose: The purpose of this research is to perform the fast reconstruction of a prompt gamma ray image using a graphics processing unit (GPU) computation from boron neutron capture therapy (BNCT) simulations. Methods: To evaluate the accuracy of the reconstructed image, a phantom including four boron uptake regions (BURs) was used in the simulation. After the Monte Carlo simulation of the BNCT, the modified ordered subset expectation maximization reconstruction algorithm using the GPU computation was used to reconstruct the images with fewer projections. The computation times for image reconstruction were compared between the GPU and the central processing unit (CPU). Also, the accuracy of the reconstructed image was evaluated by a receiver operating characteristic (ROC) curve analysis. Results: The image reconstruction time using the GPU was 196 times faster than the conventional reconstruction time using the CPU. For the four BURs, the area under curve values from the ROC curve were 0.6726 (A-region), 0.6890 (B-region), 0.7384 (C-region), and 0.8009 (D-region). Conclusion: The tomographic image using the prompt gamma ray event from the BNCT simulation was acquired using the GPU computation in order to perform a fast reconstruction during treatment. The authors verified the feasibility of the prompt gamma ray reconstruction using the GPU computation for BNCT simulations.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

TU-H-CAMPUS-IeP1-01: Bias and Computational Efficiency of Variance Reduction Methods for the Monte Carlo Simulation of Imaging Detectors

Purpose: Variance reduction techniques (VRTs) are employed in Monte Carlo simulations to obtain estimates with reduced statistical uncertainty for a given simulation time. In this work, we study the bias and efficiency of a VRT for estimating the response of imaging detectors. Methods: We implemented Directed Sampling (DS), preferentially directing a fraction of emitted optical photons directly towards the detector by altering the isotropic model. The weight of each optical photon is appropriately modified to maintain simulation estimates unbiased. We use a Monte Carlo tool called fastDETECT2 (part of the hybridMANTIS open-source package) for optical transport, modified for VRT. The weight of each photon is calculated as the ratio of original probability (no VRT) and the new probability for a particular direction. For our analysis of bias and efficiency, we use pulse height spectra, point response functions, and Swank factors. We obtain results for a variety of cases including analog (no VRT, isotropic distribution), and DS with 0.2 and 0.8 optical photons directed towards the sensor plane. We used 10,000, 25-keV primaries. Results: The Swank factor for all cases in our simplified model converged fast (within the first 100 primaries) to a stable value of 0.9. The root mean square error per pixel for DS VRT for the point response function between analog and VRT cases was approximately 5e-4. Conclusion: Our preliminary results suggest that DS VRT does not affect the estimate of the mean for the Swank factor. Our findings indicate that it may be possible to design VRTs for imaging detector simulations to increase computational efficiency without introducing bias.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2017

Verifying the reciprocity of interparticle interaction forces in strongly coupled systems

We propose a simple criterion for revealing the breaking of pair interaction symmetry in strongly coupled dissipative systems. The criterion is based on the analysis of correlations between the velocities of strongly interacting particles, which can be measured relatively easily in experiments with macroparticles in various media. We derive analytic relations that make it possible to calculate the derivatives of the interaction force between a pair of particles from the data on the correlations of their velocities and coordinates. The proposed criterion and relations are verified using the results of numerical simulation of the dynamics of dust particles in a plasma.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Parametric decay of an extraordinary electromagnetic wave in relativistic plasma

Parametric instability of an extraordinary electromagnetic wave in plasma preheated to a relativistic temperature is considered. A set of self-similar nonlinear differential equations taking into account the electron “thermal” mass is derived and investigated. Small perturbations of the parameters of the heated plasma are analyzed in the linear approximation by using the dispersion relation determining the phase velocities of the fast and slow extraordinary waves. In contrast to cold plasma, the evanescence zone in the frequency range above the electron upper hybrid frequency vanishes and the asymptotes of both branches converge. Theoretical analysis of the set of nonlinear equations shows that the growth rate of decay instability increases with increasing initial temperature of plasma electrons. This result is qualitatively confirmed by numerical simulations of plasma heating by a laser pulse injected from vacuum.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Synthesis of Electron-Optical Systems with Compression of Sheet Beam for Terahertz TWTs

Electron-optical systems are synthesized to form a converging sheet electron beam with a compression of 15 and 20, a cross section of 0.05 × 2 mm, and a current density of 100 A/cm{sup 2} in the presence of complete magnetic shielding of field-emission cathode. Deformation of low-perveance flux in the presence of magnetic field in the drift tunnel of slow-wave structure is analyzed with the aid of 3D computer simulation of the electron-optical system with the sheet electron beam.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Dakota Uncertainty Quantification Methods Applied to the CFD code Nek5000

This report presents the state of advancement of a Nuclear Energy Advanced Modeling and Simulation (NEAMS) project to characterize the uncertainty of the computational fluid dynamics (CFD) code Nek5000 using the Dakota package for flows encountered in the nuclear engineering industry. Nek5000 is a high-order spectral element CFD code developed at Argonne National Laboratory for high-resolution spectral-filtered large eddy simulations (LESs) and unsteady Reynolds-averaged Navier-Stokes (URANS) simulations.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2016

FORMATION OF GLOBULAR CLUSTERS IN ATOMIC-COOLING HALOS VIA RAPID GAS CONDENSATION AND FRAGMENTATION DURING THE EPOCH OF REIONIZATION

We investigate the formation of metal-poor globular clusters (GCs) at the center of two dark matter halos with M{sub halo}∼4×10{sup 7} M{sub ⊙} at z>10 using cosmological radiation-hydrodynamics simulations. We find that very compact (≲1 pc) and massive (∼6×10{sup 5} M{sub ⊙}) clusters form rapidly when pristine gas collapses isothermally with the aid of efficient Lyα emission during the transition from molecular-cooling halos to atomic-cooling halos. Because the local free-fall time of dense star-forming gas is very short (≪1 Myr), a large fraction of the collapsed gas is turned into stars before stellar feedback processes blow out the gas and shut down star formation. Although the early stage of star formation is limited to a small region of the central star-forming disk, we find that the disk quickly fragments due to metal enrichment from supernovae. Sub-clusters formed in the fragmented clouds eventually merge with the main cluster at the center. The simulated clusters closely resemble the local GCs in mass and size but show a metallicity spread that is much wider than found in the local GCs. We discuss a role of pre-enrichment by Pop III and II stars as a potential solution to the latter issue. Although not without shortcomings, it is encouraging that a naive blind (not tuned) cosmological simulation presents a possible channel for the formation of at least some massive GCs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Tracking the Magnetic Flux in and Around Sunspots

We have developed a procedure for tracking sunspots observed by the Helioseismic and Magnetic Imager on the Solar Dynamics Observatory and for making curvature-corrected space/time maps of the associated line-of-sight magnetic field and continuum intensity. We apply this procedure to 36 sunspots, each observed continuously for nine days around its central meridian passage time, and find that the proper motions separate into two distinct components depending on their speeds. Fast (∼3–5 km s{sup −1}) motions, comparable to Evershed flows, are produced by weak vertical fluctuations of the horizontal canopy field and recur on a timescale of 12–20 min. Slow (∼0.3–0.5 km s{sup −1}) motions diverge from a sunspot-centered ring whose location depends on the size of the sunspot, occurring in the mid-penumbra for large sunspots and at the outer edge of the penumbra for small sunspots. The slow ingoing features are contracting spokes of a quasi-vertical field of umbral polarity. These inflows disappear when the sunspot loses its penumbra, and may be related to inward-moving penumbral grain. The slow outgoing features may have either polarity depending on whether they originate from quasi-vertical fields of umbral polarity or from the outer edge of the canopy. When a sunspot decays, the penumbra and canopy disappear, and the moat becomes filled with slow outflows of umbral polarity. We apply our procedure to decaying sunspots, to long-lived sunspots, and to numerical simulations of a long-lived sunspot by Rempel.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Implementation and validation of a fully implicit accumulator model in RELAP-7

This paper presents the implementation and validation of an accumulator model in RELAP-7 under the framework of preconditioned Jacobian free Newton Krylov (JFNK) method, based on the similar model used in RELAP5. RELAP-7 is a new nuclear reactor system safety analysis code being developed at the Idaho National Laboratory (INL). RELAP-7 is a fully implicit system code. The JFNK and preconditioning methods used in RELAP-7 is briefly discussed. The slightly modified accumulator model is summarized for completeness. The implemented model was validated with LOFT L3-1 test and benchmarked with RELAP5 results. RELAP-7 and RELAP5 had almost identical results for the accumulator gas pressure and water level, although there were some minor difference in other parameters such as accumulator gas temperature and tank wall temperature. One advantage of the JFNK method is its easiness to maintain and modify models due to fully separation of numerical methods from physical models. It would be straightforward to extend the current RELAP-7 accumulator model to simulate the advanced accumulator design. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2019

Emplacement Room Mock-up Structure for Testing the Placement of Nuclear Waste Within a Deep Geologic Repository - 19271

The Nuclear Waste Management Organization (NWMO) has conceptualized a unique method of long-term storage for spent nuclear waste. The concept encompasses all aspects of the end-of-life process from the point the spent fuel is no longer in-service; including handling, shipping, containing, storing, monitoring, and finally decommissioning. The equipment developed for each process is unique and designed specifically for each operation. Since these methods have not been implemented before, it is required to prove each process and the related equipment through simulations and eventually practical tests to evaluate performance and to confirm effective and safe operation. The spent fuel will be stored in an underground repository several hundred meters below the surface using autonomous machinery. Considering the significant investment and labour involved in creating full-scale tests underground, it would be impractical at this project stage to attempt tests of the newly developed machinery and methods in other than a surface setting. Therefore, our scope was to provide a customizable mock-up structure to simulate the conditions and constraints an underground repository would present. (authors)

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

Development of a core coolant flow distribution calculation code for HTGRs - 18484

In order to ensure the thermal integrity of fuel in high temperature gas-cooled reactors (HTGRs), it is necessary that the maximum fuel temperature in the normal operation is to be lower than the thermal design target. In the core thermal-hydraulic design of block-type HTGRs, the maximum fuel temperature should be evaluated considering data such as thermal power, core geometry, power and neutron fluence distributions, and core coolant flow distribution. The core coolant flow distribution calculation code used in the design stage of High Temperature engineering Test Reactor (HTTR) presupposes to run on UNIX systems, and its operation and execution procedure are complicated and not user-friendly. Therefore, a new core coolant flow distribution calculation code with a user-friendly system such as simple and easy operations and execution procedures has been developed. In this paper, the outline of the new code is described and the simulation result of an out-of-pile test with one fuel column are shown as the first step of the code validation. The simulation results provide good agreement with the test one. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2018

Analytical Modeling of Triple-Metal Hetero-Dielectric DG SON TFET

In this paper, a 2-D analytical model of triple-metal hetero-dielectric DG TFET is presented by combining the concepts of triple material gate engineering and hetero-dielectric engineering. Three metals with different work functions are used as both front- and back gate electrodes to modulate the barrier at source/channel and channel/drain interface. In addition to this, front gate dielectric consists of high-K HfO{sub 2} at source end and low-K SiO{sub 2} at drain side, whereas back gate dielectric is replaced by air to further improve the ON current of the device. Surface potential and electric field of the proposed device are formulated solving 2-D Poisson’s equation and Young’s approximation. Based on this electric field expression, tunneling current is obtained by using Kane’s model. Several device parameters are varied to examine the behavior of the proposed device. The analytical model is validated with TCAD simulation results for proving the accuracy of our proposed model.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Influence of Thermal History on the Hot Ductility of a Continuously Cast Low Alloyed Cr-Mo Steel

The hot ductility of a low alloyed Cr-Mo steel has been investigated to evaluate the surface cracking sensitivity within the straightening or unbending regime during the continuous casting process. Tensile samples were subjected to various thermal treatments, including melting and solidification, and were tested at deforming temperatures ranging between 600 and 1100 °C using a strain rate of 10{sup −3} s{sup −1}. Hot ductility was evaluated based on reduction in area measurement and metallographic investigations. The investigated steel exhibits a drop in ductility at around 800 °C due to intergranular cracking. Microstructural examinations and supplementary thermokinetic computer simulations were carried out to describe the evolution of the microstructure during solidification and cooling.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Hopf bifurcation of a delayed chemostat model with general monotone response functions

In this paper, a chemostat model with general monotone response functions and two discrete time delays is proposed to describe the dynamical behavior about predator–prey system. First, by analyzing the characteristic equation associated with the model, we obtain the conditions of the existence and stability of extinction equilibrium and positive equilibrium. Choosing delays as bifurcation parameters, the existence of Hopf bifurcations is investigated in detail. Second, by virtue of the Poincaré normal form method and center manifold theorem, explicit formulas are derived to determine the direction of Hopf bifurcation and the stability of the bifurcating periodic solutions. Finally, some numerical simulations are carried out to illustrate the theoretical results and the biological significance.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Robust and non-fragile finite time H{sub ∞} synchronization control for complex networks with uncertain inner coupling

This paper addresses robust finite time H{sub ∞} synchronization problem for non-fragile complex networks with uncertain inner coupling. In order to achieve synchronization, a non-fragile controller subjected to randomly occurring perturbations will be considered. By constructing a suitable Lyapunov–Krasovskii functional and utilizing reciprocally convex approach, the criteria for the networks are established in terms of linear matrix inequalities (LMIs). Theoretical analysis and numerical simulations about robust finite time H{sub ∞} synchronization for uncertain inner coupling complex networks are presented.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Power spectrum and form factor in random diagonal matrices and integrable billiards

Triggered by a controversy surrounding a universal behavior of the power spectrum in quantum systems exhibiting regular classical dynamics, we focus on a model of random diagonal matrices (RDM), often associated with the Poisson spectral universality class, and examine how the power spectrum and the form factor get affected by two-sided truncations of RDM spectra. Having developed a nonperturbative description of both statistics, we perform their detailed asymptotic analysis to demonstrate explicitly how a traditional assumption (lying at the heart of the controversy) – that the power spectrum is merely determined by the spectral form factor – breaks down for truncated spectra. This observation has important consequences as we further argue that bounded quantum systems with integrable classical dynamics are described by heavily truncated rather than complete RDM spectra. High-precision numerical simulations of semicircular and irrational rectangular billiards lend independent support to these conclusions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Design of a backscatter 14-MeV neutron time-of-flight spectrometer for experiments at ITER

Neutron energy spectrometry diagnostics play an important role in present-day experiments related to fusion energy research. Measurements and thorough analysis of the neutron emission from the fusion plasma give information on a number of basic fusion performance quantities, on the condition of the neutron source and plasma behavior. Here we discuss the backscatter Time-of-Flight (bTOF) spectrometer concept as a possible instrument for performing high resolution measurements of 14 MeV neutrons. The instrument is based on two sets of scintillators, a first scatterer exposed to a collimated neutron beam and a second detector set placed in the backward direction. The scintillators of the first set are enriched in deuterium to achieve neutron backscattering. The energy resolution and efficiency of a bTOF instrument have been determined for various geometrical configurations. A preliminary design of optimal geometry for the two scintillator sets has been obtained by Monte Carlo simulations based on the MCNPX code.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗