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

Comparison of neutron diffusion and Monte Carlo models for a fission wave

Many groups have used neutron diffusion simulations to study fission wave phenomena in natural or depleted uranium. However, few studies of fission wave phenomena have been published that use Monte Carlo simulations to confirm the results of diffusion models for this type of system. In the present work we show the results of a criticality and burnup simulation of a traveling wave reactor using MCNPX 2.7.0. The characteristics of the fission wave in this simulation are compared with those from a simple one-dimensional, one-group neutron diffusion model. The diffusion simulations produce a wave speed of 5.9 cm/yr versus 5.3 cm/yr for the Monte Carlo simulations. The axial flux profile in the Monte Carlo simulation is similar in shape to the diffusion results, but with different peak values, and the two profiles have an R2 value of 0.93. The {sup 238}U, {sup 239}Np and {sup 239}Pu burnup profiles from the diffusion simulation show good agreement with the Monte Carlo simulations, R values of 0.98, 0.93 and 0.97 respectively are observed. (authors)

97 MATHEMATICS AND COMPUTING↗
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 OSTI2013

Radio frequency cavity analysis, measurement, and calibration of absolute Dee voltage for K-500 superconducting cyclotron at VECC, Kolkata

Variable Energy Cyclotron Centre has commissioned a K-500 superconducting cyclotron for various types of nuclear physics experiments. The 3-phase radio-frequency system of superconducting cyclotron has been developed in the frequency range 9-27 MHz with amplitude and phase stability of 100 ppm and {+-}0.2{sup 0}, respectively. The analysis of the RF cavity has been carried out using 3D Computer Simulation Technology (CST) Microwave Studio code and various RF parameters and accelerating voltages ('Dee' voltage) are calculated from simulation. During the RF system commissioning, measurement of different RF parameters has been done and absolute Dee voltage has been calibrated using a CdTe X-ray detector along with its accessories and known X-ray source. The present paper discusses about the measured data and the simulation result.

43 PARTICLE ACCELERATORS↗
DOE OSTI2013

Modeling surface backgrounds from radon progeny plate-out

The next generation low-background detectors operating deep underground aim for unprecedented low levels of radioactive backgrounds. The surface deposition and subsequent implantation of radon progeny in detector materials will be a source of energetic background events. We investigate Monte Carlo and model-based simulations to understand the surface implantation profile of radon progeny. Depending on the material and region of interest of a rare event search, these partial energy depositions can be problematic. Motivated by the use of Ge crystals for the detection of neutrinoless double-beta decay, we wish to understand the detector response of surface backgrounds from radon progeny. We look at the simulation of surface decays using a validated implantation distribution based on nuclear recoils and a realistic surface texture. Results of the simulations and measured α spectra are presented.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2013

On collisional diffusion in a stochastic magnetic field

The effect of particle collisions on the transport in a stochastic magnetic field in tokamaks is investigated. The model of resonant magnetic perturbations generated by external coils at the plasma edge is used for the stochastic magnetic field. The particle collisions are simulated by a random walk process along the magnetic field lines and the jumps across the field lines at the collision instants. The dependencies of the local diffusion coefficients on the mean free path λ{sub mfp}, the diffusion coefficients of field lines D{sub FL}, and the collisional diffusion coefficients, χ{sub ⊥} are studied. Based on these numerical data and the heuristic arguments, the empirical formula, D{sub r}=χ{sub ⊥}+v{sub ||}D{sub FL}/(1+L{sub c}/λ{sub mfp}), for the local diffusion coefficient is proposed, where L{sub c} is the characteristic length of order of the connection length l{sub c}=πqR{sub 0}, q is the safety factor, R{sub 0} is the major radius. The formula quite well describes the results of numerical simulations. In the limiting cases, the formula describes the Rechester-Rosenbluth and Laval scalings.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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↗
DOE OSTI2014

The design and optimization of two low frequency energy harvesters employing 3C-SiC/AlN/Mo composite layers

This paper presents the design and simulation of twocantilever-based energy harvesters that employs cubic silicon carbide on silicon (3C-SiC-on-Si) wafer as the base material and bottom electrode. Aluminum Nitride (AlN) is employed as the piezoelectric/middle layer due to its excellent material properties and high stability in varying temperature and harsh environment. Molybdenum (Mo) serves as the top layer/electrode. The thickness of the structural layers are optimized through MATLAB and also analyzed via Finite Element Analysis using Intellisuite. Two designs are proposed at low resonant frequency, one with conventional cantilever beam, the other being a T-shaped cantilever beam. Both structures are simulated and their performances are compared.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

Post Fukushima tsunami simulations for Malaysian coasts

The recent recurrences of mega tsunamis in the Asian region have rekindled concern regarding potential tsunamis that could inflict severe damage to affected coastal facilities and communities. The 11 March 2011 Fukushima tsunami that crippled nuclear power plants in Northern Japan has further raised the level of caution. The recent discovery of petroleum reserves in the coastal water surrounding Malaysia further ignites the concern regarding tsunami hazards to petroleum facilities located along affected coasts. Working in a group, federal government agencies seek to understand the dynamics of tsunami and their impacts under the coordination of the Malaysian National Centre for Tsunami Research, Malaysian Meteorological Department. Knowledge regarding the generation, propagation and runup of tsunami would provide the scientific basis to address safety issues. An in-house tsunami simulation models known as TUNA has been developed by the authors to assess tsunami hazards along affected beaches so that mitigation measures could be put in place. Capacity building on tsunami simulation plays a critical role in the development of tsunami resilience. This paper aims to first provide a simple introduction to tsunami simulation towards the achievement of tsunami simulation capacity building. The paper will also present several scenarios of tsunami dangers along affected Malaysia coastal regions via TUNA simulations to highlight tsunami threats. The choice of tsunami generation parameters reflects the concern following the Fukushima tsunami.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2014

Constant pressure and temperature discrete-time Langevin molecular dynamics

We present a new and improved method for simultaneous control of temperature and pressure in molecular dynamics simulations with periodic boundary conditions. The thermostat-barostat equations are built on our previously developed stochastic thermostat, which has been shown to provide correct statistical configurational sampling for any time step that yields stable trajectories. Here, we extend the method and develop a set of discrete-time equations of motion for both particle dynamics and system volume in order to seek pressure control that is insensitive to the choice of the numerical time step. The resulting method is simple, practical, and efficient. The method is demonstrated through direct numerical simulations of two characteristic model systems—a one-dimensional particle chain for which exact statistical results can be obtained and used as benchmarks, and a three-dimensional system of Lennard-Jones interacting particles simulated in both solid and liquid phases. The results, which are compared against the method of Kolb and Dünweg [J. Chem. Phys. 111, 4453 (1999)], show that the new method behaves according to the objective, namely that acquired statistical averages and fluctuations of configurational measures are accurate and robust against the chosen time step applied to the simulation.

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

Maximizing semi-active vibration isolation utilizing a magnetorheological damper with an inner bypass configuration

A single-degree-of-freedom (SDOF) semi-active vibration control system based on a magnetorheological (MR) damper with an inner bypass is investigated in this paper. The MR damper employing a pair of concentric tubes, between which the key structure, i.e., the inner bypass, is formed and MR fluids are energized, is designed to provide large dynamic range (i.e., ratio of field-on damping force to field-off damping force) and damping force range. The damping force performance of the MR damper is modeled using phenomenological model and verified by the experimental tests. In order to assess its feasibility and capability in vibration control systems, the mathematical model of a SDOF semi-active vibration control system based on the MR damper and skyhook control strategy is established. Using an MTS 244 hydraulic vibration exciter system and a dSPACE DS1103 real-time simulation system, experimental study for the SDOF semi-active vibration control system is also conducted. Simulation results are compared to experimental measurements.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗