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 1,297 records · Page 72
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

WE-D-BRF-01: FEATURED PRESENTATION - Investigating Particle Track Structures Using Fluorescent Nuclear Track Detectors and Monte Carlo Simulations

Purpose: To report on the efforts funded by the AAPM seed funding grant to develop the basis for fluorescent nuclear track detector (FNTD) based radiobiological experiments in combination with dedicated Monte Carlo simulations (MCS) on the nanometer scale. Methods: Two confocal microscopes were utilized in this study. Two FNTD samples were used to find the optimal microscope settings, one FNTD irradiated with 11.1 MeV/u Gold ions and one irradiated with 428.77 MeV/u Carbon ions. The first sample provided a brightly luminescent central track while the latter is used to test the capabilities to observe secondary electrons. MCS were performed using TOPAS beta9 version, layered on top of Geant4.9.6p02. Two sets of simulations were performed, one with the Geant4-DNA physics list and approximating the FNTDs by water, a second set using the Penelope physics list in a water-approximated FNTD and a aluminum-oxide FNTD. Results: Within the first half of the funding period, we have successfully established readout capabilities of FNTDs at our institute. Due to technical limitations, our microscope setup is significantly different from the approach implemented at the DKFZ, Germany. However, we can clearly reconstruct Carbon tracks in 3D with electron track resolution of 200 nm. A second microscope with superior readout capabilities will be tested in the second half of the funding period, we expect an improvement in signal to background ratio with the same the resolution.We have successfully simulated tracks in FNTDs. The more accurate Geant4-DNA track simulations can be used to reconstruct the track energy from the size and brightness of the observed tracks. Conclusion: We have achieved the goals set in the seed funding proposal: the setup of FNTD readout and simulation capabilities. We will work on improving the readout resolution to validate our MCS track structures down to the nanometer scales.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2017

Visualization of the microwave beam generated by a plasma relativistic microwave amplifier

A method based on the detection of emission of a dielectric screen with metal microinclusions in open air is applied to visualize the transverse structure of a high-power microwave beam. In contrast to other visualization techniques, the results obtained in this work provide qualitative information not only on the electric field strength, but also on the structure of electric field lines in the microwave beam cross section. The interpretation of the results obtained with this method is confirmed by numerical simulations of the structure of electric field lines in the microwave beam cross section by means of the CARAT code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Study of the anode plasma dynamics under the action of a high-power electron beam on epoxy resin

Results are presented from experimental studies of plasma dynamics in a diode gap under the action of a high-current relativistic electron beam on epoxy resin at energy densities in the range of 170–860 J/cm{sup 2}. The plasma expansion was studied by means of an optical streak camera. Three-dimensional numerical simulations in the one-temperature hydrodynamic approximation were also performed. The experimental data are compared with the results of numerical simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Effect of the J–J Interaction of Excited States of the Rare-Earth Ion Pr{sup 3+} on Magnetically Polarized Luminescence of Praseodymium-Yttrium Aluminum Garnet

The spectra of luminescence and magnetic circular polarization of luminescence of praseodymium–yttrium aluminum garnet Pr{sup 3+} : Y{sub 3}Al{sub 5}O{sub 12} (PrYAG) are studied in the visible spectral region at temperature T = 300 K. An analysis of spectral dependences of magnetooptical and optical spectra makes it possible to identify optical 4 f–4 f-transitions between Stark sublevels of multiplets {sup 3}P{sub 0}, {sup 3}P{sub 1}, {sup 3}Н{sub 5}, and {sup 3}Н{sub 6} in PrYAG. It was shown that an important role in the spectrum of the degree of magnetic circular polarization of luminescence of this paramagnetic garnet is played by the effect of quantum-mechanical J–J mixing of states of Stark singlets {sup 3}Н{sub 5} and {sup 3}Н{sub 6} of non-Kramer rare-earth ion Pr{sup 3+} in the “green” luminescence band related to forbidden 4 f → 4 f transition {sup 3}P{sub 0} → {sup 3}Н{sub 5} in the visible spectral region. To interpret the spectra of magnetic circular polarization of luminescence, the energy of experimentally determined Stark sublevels of multiplets under study, their irreducible representations and wave functions determined by numerical simulation of the energy spectrum of the rare-earth ion Pr{sup 3+} in the garnet structure are used.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2017

The Numerical Simulation of the Nanosecond Switching of a p-SOS Diode

Abrupt high-density reverse current interruption has been numerically simulated for switching from forward to reverse bias in a silicon p{sup +}P{sub 0}n{sup +} structure (p-SOS diode). It has been shown that the current interruption in this structure occurs as a result of the formation of two dynamic domains of a strong electric field in regions in which the free carrier concentration substantially exceeds the concentration of the doping impurity. The first domain is formed in the n{sup +} region at the n{sup +}P{sup 0} junction, while the second domain is formed in the P{sup 0} region at the interface with the p{sup +} layer. The second domain expands much faster, and this domain mainly determines the current interruption rate. Good agreement is achieved between the simulation results and the experimental data when the actual electric circuit determining the electron–hole plasma pumping in and out is accurately taken into account.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Putting multipoint kinetics to the test: preliminary analysis of the MUSE-4 experiment with APOLLO3

Modeling neutron kinetics is a challenging task and over the years, several approaches have been developed offering different tradeoffs between precision and cost. Probably the most popular one is point kinetics (PK), which, despite being very cheap, provides good results for many applications. One of the limitations of PK is that it is not able to model effects that happen before the prompt jump time scale. In practice, applications, where those time scales are interesting, are rare and those for which experimental data are available are even rarer. One of those cases is the MUSE-4 experiment for which the focus is on the very fast behavior and where we have no thermal feedback. As we have shown in a recent paper, a simple extension of PK, known as multipoint kinetics (MPK), can overcome this limitation. In this paper, which is intended as a preliminary study, we show how a simple MPK model can qualitatively reproduce the fast behavior that we cannot model with PK. Moreover, we will show how to exploit its low cost and compact form to model also slow effects as the precursors built up, behaviours that are not affordable to model with more expensive techniques, such as kinetic Monte Carlo simulations. For the periodic pulse transients with no feedbacks we will also provide an analytical solution in closed form. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Zero power experiment for molten salt reactor - code verification and validation request

Molten Salt Reactors (MSRs) are one of the six perspective Generation IV nuclear technologies that should fulfil sustainability, economics, safety, and proliferation resistance goals. In recent years these concepts have attracted increased interest in research as well as in industry. In the context of developing new, innovative technology, a zero-power experiment would be an ideal option for improving the understanding of MSR systems under low-risk and low-cost conditions. The given research was focused on cross-verification of the nuclear codes for MSR modelling and simulations as required for the development of the first design of a potential facility and a draft of an experimental programme. The analysis included the evaluation of lattice and Monte Carlo codes. The study showed that all considered nuclear tools were in good agreement with each other. Although this does not mean the results will match the real experimental data, which is essential for the regulatory approval of a novel technology. Currently, experimental data is almost absent for MSR technology. Hence, a zero power experiment is needed to assure we can rely on modelling and simulation studies for further development of MSR technology and deliver real-world experiments to support the licensing of future MSRs. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2015

Proteus-SN user manual

This user manual describes how to set up a neutron transport simulation with the PROTEUS-SN code. A companion methodology manual describes the theory and algorithms within PROTEUS-SN.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2015

Shape evolution in neutron-rich nuclei near the r-process path

To assess the resources needed to field a quantitative measurement to advance our knowledge on shape evolution, we have carried out a numerical simulation of the Coulomb excitation for 112 Ru at 460 MeV bombarding energy on a 208 Pb target of 1 mg/cm 2 using a semiclassical code, GOSIA with GRETINA/CHICO2.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2016

Development and Implementation of Mechanistic Terry Turbine Models in RELAP-7 to Simulate RCIC Normal Operation Conditions

As part of the efforts to understand the unexpected “self-regulating” mode of the RCIC (Reactor Core Isolation Cooling) systems in Fukushima accidents and extend BWR RCIC and PWR AFW (Auxiliary Feed Water) operational range and flexibility, mechanistic models for the Terry turbine, based on Sandia’s original work [1], have been developed and implemented in the RELAP-7 code to simulate the RCIC system. In 2016, our effort has been focused on normal working conditions of the RCIC system. More complex off-design conditions will be pursued in later years when more data are available. In the Sandia model, the turbine stator inlet velocity is provided according to a reduced-order model which was obtained from a large number of CFD (computational fluid dynamics) simulations. In this work, we propose an alternative method, using an under-expanded jet model to obtain the velocity and thermodynamic conditions for the turbine stator inlet. The models include both an adiabatic expansion process inside the nozzle and a free expansion process outside of the nozzle to ambient pressure. The combined models are able to predict the steam mass flow rate and supersonic velocity to the Terry turbine bucket entrance, which are the necessary input information for the Terry turbine rotor model. The analytical models for the nozzle were validated with experimental data and benchmarked with CFD simulations. The analytical models generally agree well with the experimental data and CFD simulations. The analytical models are suitable for implementation into a reactor system analysis code or severe accident code as part of mechanistic and dynamical models to understand the RCIC behaviors. The newly developed nozzle models and modified turbine rotor model according to the Sandia’s original work have been implemented into RELAP-7, along with the original Sandia Terry turbine model. A new pump model has also been developed and implemented to couple with the Terry turbine model. An input model was developed to test the Terry turbine RCIC system, which generates reasonable results. Both the INL RCIC model and the Sandia RCIC model produce results matching major rated parameters such as the rotational speed, pump torque, and the turbine shaft work for the normal operation condition. The Sandia model is more sensitive to the turbine outlet pressure than the INL model. The next step will be further refining the Terry turbine models by including two-phase flow cases so that off-design conditions can be simulated. The pump model could also be enhanced with the use of the homologous curves.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2016

INSIDE-OUT PLANET FORMATION. III. PLANET–DISK INTERACTION AT THE DEAD ZONE INNER BOUNDARY

The Kepler mission has discovered more than 4000 exoplanet candidates. Many of them are in systems with tightly packed inner planets. Inside-out planet formation (IOPF) has been proposed as a scenario to explain these systems. It involves sequential in situ planet formation at the local pressure maximum of a retreating dead zone inner boundary (DZIB). Pebbles accumulate at this pressure trap, which builds up a pebble ring and then a planet. The planet is expected to grow in mass until it opens a gap, which helps to both truncate pebble accretion and also induce DZIB retreat that sets the location of formation of the next planet. This simple scenario may be modified if the planet undergoes significant migration from its formation location. Thus, planet–disk interactions play a crucial role in the IOPF scenario. Here we present numerical simulations that first assess the degree of migration for planets of various masses that are forming at the DZIB of an active accretion disk, where the effective viscosity is undergoing a rapid increase in the radially inward direction. We find that torques exerted on the planet by the disk tend to trap the planet at a location very close to the initial pressure maximum where it formed. We then study gap opening by these planets to assess at what mass a significant gap is created. Finally, we present a simple model for DZIB retreat due to penetration of X-rays from the star to the disk midplane. Overall, these simulations help to quantify both the mass scale of first (“Vulcan”) planet formation and the orbital separation to the location of second planet formation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

TEMPERATURE SPECTRA OF INTERSTELLAR DUST GRAINS HEATED BY COSMIC RAYS. I. TRANSLUCENT CLOUDS

Heating of whole interstellar dust grains by cosmic-ray (CR) particles affects the gas–grain chemistry in molecular clouds by promoting molecule desorption, diffusion, and chemical reactions on grain surfaces. The frequency of such heating, f{sub T}, s{sup −1}, determines how often a certain temperature T{sub CR}, K, is reached for grains hit by CR particles. This study aims to provide astrochemists with a comprehensive and updated data set on CR-induced whole-grain heating. We present calculations of f{sub T} and T{sub CR} spectra for bare olivine grains with radius a of 0.05, 0.1, and 0.2 μ m and such grains covered with ice mantles of thickness 0.1 a and 0.3 a . Grain shape and structure effects are considered, as well as 30 CR elemental constituents with an updated energy spectrum corresponding to a translucent cloud with A{sub V} = 2 mag. Energy deposition by CRs in grain material was calculated with the srim program. We report full T{sub CR} spectra for all nine grain types and consider initial grain temperatures of 10 K and 20 K. We also provide frequencies for a range of minimum T{sub CR} values. The calculated data set can be simply and flexibly implemented in astrochemical models. The results show that, in the case of translucent clouds, the currently adopted rate for heating of whole grains to temperatures in excess of 70 K is underestimated by approximately two orders of magnitude in astrochemical numerical simulations. Additionally, grains are heated by CRs to modest temperatures (20–30 K) with intervals of a few years, which reduces the possibility of ice chemical explosions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Gravitational Waves from F-modes Excited by the Inspiral of Highly Eccentric Neutron Star Binaries

As gravitational wave instrumentation becomes more sensitive, it is interesting to speculate about subtle effects that could be analyzed using upcoming generations of detectors. One such effect that has great potential for revealing the properties of very dense matter is fluid oscillations of neutron stars. These have been found in numerical simulations of the hypermassive remnants of double neutron star mergers and of highly eccentric neutron star orbits. Here we focus on the latter and sketch out some ideas for the production, gravitational-wave detection, and analysis of neutron star oscillations. These events will be rare (perhaps up to several tens per year could be detected using third-generation detectors such as the Einstein Telescope or the Cosmic Explorer), but they would have unique diagnostic power for the analysis of cold, catalyzed, dense matter. Furthermore, these systems are unusual in that analysis of the tidally excited f-modes of the stars could yield simultaneous measurements of their masses, moments of inertia, and tidal Love numbers, using the frequency, damping time, and amplitude of the modes. They would thus present a nearly unique opportunity to test the I-Love-Q relation observationally. The analysis of such events will require significant further work in nuclear physics and general relativistic nonlinear mode coupling, and thus we discuss further directions that will need to be pursued. For example, we note that for nearly grazing encounters, numerical simulations show that the energy delivered to the f-modes may be up to two orders of magnitude greater than predicted in the linear theory.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

An Upper Limit on the Milky Way Mass from the Orbit of the Sagittarius Dwarf Satellite

As one of the most massive Milky Way (MW) satellites, the Sagittarius dwarf galaxy played an important role in shaping the Galactic disk and stellar halo morphologies. The disruption of Sagittarius over several close-in passages has populated the halo of our Galaxy with large-scale tidal streams and offers a unique diagnostic tool for measuring its gravitational potential. Here, we test different progenitor mass models for the MW and Sagittarius by modeling the full infall of the satellite. We constrain the mass of the Galaxy based on the kinematics of the satellite remnant and multiple tidal streams of Sagittarius. Our semianalytic modeling of the orbital dynamics agrees with full N-body simulations, and favors low values for the MW mass, ≲10{sup 12} M{sub ⊙}. This conclusion eases the tension between ΛCDM and the observed parameters of the MW satellites.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Research reactor plate fuel multi-physics hotspot model for the HFIR

Ongoing development of a multi-physics hotspot model for nuclear research reactor plate fuel is presented with a focus on the High Flux Isotope Reactor (HFIR) at Oak Ridge National Laboratory (ORNL). The background of plate fueled research reactors in general and the HFIR in particular are discussed along with the need to convert from high enriched uranium (HEU) to low enriched uranium (LEU) fuel, which drives the need to improve the analysis capability of fuel thermal-fluid performance. HFIR HEU fuel plate fabrication is described to explain the origin of fuel segregations and cladding non-bonds that currently contribute to the hotspot flux peaking factors. Previous efforts at modeling hotspots caused by these defects are reviewed and areas for improvement associated with high resolution multi-physics simulation are identified. Modes of heat transfer through non-bonds, including thermal contact conductance, are examined. Conjugate heat transfer over hotspots is also explored to high resolution up to and through the wall boundary. The challenges associated with multi-physics verification and validation (VV), particularly the requirements beyond VV 20, are presented. An approach to Validation within the limits of existing data for conduction heat transfer through the fuel plate, and convection into the cooling fluid is set forth. Preliminary simulations are presented assessing hot spot modeling approaches for contact conductance, and fluid turbulence modeling within COMSOL. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2017

Simulation and control of a passively cooled small modular reactor Samet - 142

This work regulates a small modular reactor (SMR) with traditional control methodologies and analyzes its response to typical disturbances. An analytical model for the passively cooled NuScale SMR is developed to represent the essential dynamics. The MATLAB/Simulink based model consists of reactor core, steam generator, pressurizer, hot leg riser and downcomer, and turbine representations. The point kinetics equations with a single combined neutron precursor group and the lumped parameter representations with single-phase natural circulation account for the neutronics and thermohydraulics in the reactor core. A lumped parameter, moving-boundary approach forms the basis of the nonlinear model adopted for the once-through helical-coil heat exchangers in which boundaries between regions of different fluid states (i.e., subcooled, boiling, and superheated) vary over time. A proportional-integral-derivative controller employing a sliding-average temperature strategy commonly used for pressurized water reactors is implemented. The responses of the SMR to two different disturbances are analyzed and discussed in detail. For the first scenario, a 5% step increase in the load is applied to the system, which results in a change in the steam valve opening. For comparison, two different simulations under this same disturbance are run: one including and one excluding the control systems. The second scenario tests the effectiveness of the control system to increase (or decrease) the reactor thermal power to a certain level within a desired time period. The simulation results show that important system variables are kept at the desired values by the proposed control strategy. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2017

Multi-Modular Coordination Control of HTR-PM600 Plant - 190

HTR-PM600 is a high temperature gas-cooled reactor plant based on the multi-modular scheme, i.e. six modular high temperature gas-cooled reactor (MHTGR) driving a common secondary loop system. Six MHTGR-based nuclear steam supplying system (NSSS) modules are tightly coupled by the secondary fluid flow network (FFN) which contains the secondary side of steam generators, feedwater system, main steam system, turbine generator system. The secondary FFN control, which focuses on the feedwater flowrates and main steam pressure, plays a central role in multi-modular coordination of HTR-PM600 plant. In this paper, a differential algebraic model of this secondary FFN is first proposed. Then, the corresponding flowrate-pressure control strategy is proposed, which takes the simple form of proportional-integral feedback laws. Numerical simulation results show both the feasibility and satisfactory transient performance of this controller. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2017

High Performance-Multidisciplinary Simulation for Regional Scale Earthquake Hazard and Risk Assessments

Moderate-to-large magnitude earthquakes represent a significant risk around the world and a major issue across the Department of Energy (DOE) mission interests, ranging from DOE's one-of-a-kind mission critical facilities to national energy systems (nuclear power plants, electric/gas distribution systems, and renewable energy production facilities). As a result, for the DOE there is a strong motivation in promoting and developing a risk-informed and performance-based approach for both the assessment of existing facilities and the design of new structures. DOE National Labs have historically made multidisciplinary contributions to the solution of challenging seismic hazard and risk problems and have been at the forefront for the development of design standards for major facilities. The current approach for estimating the ground motion hazard for critical facilities is based upon a Probabilistic Seismic Hazard Analysis (PSHA), whereby measured ground motion data are collected and homogenized into a probabilistic framework to estimate the hazard for future events. However, recent studies have highlighted the site-specific nature of earthquakes and shown the limitations of a purely empirical approach. Moreover, the design of infrastructure systems has been moving toward a performance-based approach, in which specific performance goals are specified at different levels of excitation and for different limit states. Such an approach requires the availability of robust and reliable computational tools for the nonlinear dynamic analysis of complex structural systems. Ongoing research at Lawrence Berkeley National Laboratory (LBNL) and Lawrence Livermore National Laboratory (LLNL) aims at developing an advanced simulation environment at the exascale for a coupled assessment of earthquake hazard (ground motions) and earthquake risk (structural demands)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗