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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

The preplasma effect on the properties of the shock wave driven by a fast electron beam

Strong shock wave generation by a mono-energetic fast electron beam in a plasma with an increasing density profile is studied theoretically. The proposed analytical model describes the shock wave characteristics for a homogeneous plasma preceded by a low density precursor. The shock pressure and the time of shock formation depend on the ratio of the electron stopping length to the preplasma areal density and on the initial energy of injected electrons. The conclusions of theoretical model are confirmed in numerical simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

On the reliability of microvariability tests in quasars

Microvariations probe the physics and internal structure of quasars. Unpredictability and small flux variations make this phenomenon elusive and difficult to detect. Variance-based probes such as the C and F tests, or a combination of both, are popular methods to compare the light curves of the quasar and a comparison star. Recently, detection claims in some studies have depended on the agreement of the results of the C and F tests, or of two instances of the F-test, for rejecting the non-variation null hypothesis. However, the C-test is a non-reliable statistical procedure, the F-test is not robust, and the combination of tests with concurrent results is anything but a straightforward methodology. A priori power analysis calculations and post hoc analysis of Monte Carlo simulations show excellent agreement for the analysis of variance test to detect microvariations as well as the limitations of the F-test. Additionally, the combined tests yield correlated probabilities that make the assessment of statistical significance unworkable. However, it is possible to include data from several field stars to enhance the power in a single F-test, increasing the reliability of the statistical analysis. This would be the preferred methodology when several comparison stars are available. An example using two stars and the enhanced F-test is presented. These results show the importance of using adequate methodologies and avoiding inappropriate procedures that can jeopardize microvariability detections. Power analysis and Monte Carlo simulations are useful tools for research planning, as they can demonstrate the robustness and reliability of different research approaches.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Conservative GRMHD simulations of moderately thin, tilted accretion disks

This paper presents our latest numerical simulations of accretion disks that are misaligned with respect to the rotation axis of a Kerr black hole. In this work, we use a new, fully conservative version of the Cosmos++ general relativistic magnetohydrodynamics (GRMHD) code, coupled with an ad hoc cooling function designed to control the thickness of the disk. Together these allow us to simulate the thinnest tilted accretion disks ever using a GRMHD code. In this way, we are able to probe the regime where the dimensionless stress and scale height of the disk become comparable. We present results for both prograde and retrograde cases. The simulated prograde tilted disk shows no sign of Bardeen-Petterson alignment even in the innermost parts of the disk. The simulated retrograde tilted disk, however, does show modest alignment. The implication of these results is that the parameter space associated with Bardeen-Petterson alignment for prograde disks may be rather small, only including very thin disks. Unlike our previous work, we find no evidence for standing shocks in our simulated tilted disks. We ascribe this to the black hole spin, tilt angle, and disk scale height all being small in these simulations. We also add to the growing body of literature pointing out that the turbulence driven by the magnetorotational instability in global simulations of accretion disks is not isotropic. Finally, we provide a comparison between our moderately thin, untilted reference simulation and other numerical simulations of thin disks in the literature.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

APPARENT CROSS-FIELD SUPERSLOW PROPAGATION OF MAGNETOHYDRODYNAMIC WAVES IN SOLAR PLASMAS

In this paper we show that the phase-mixing of continuum Alfvén waves and/or continuum slow waves in the magnetic structures of the solar atmosphere as, e.g., coronal arcades, can create the illusion of wave propagation across the magnetic field. This phenomenon could be erroneously interpreted as fast magnetosonic waves. The cross-field propagation due to the phase-mixing of continuum waves is apparent because there is no real propagation of energy across the magnetic surfaces. We investigate the continuous Alfvén and slow spectra in two-dimensional (2D) Cartesian equilibrium models with a purely poloidal magnetic field. We show that apparent superslow propagation across the magnetic surfaces in solar coronal structures is a consequence of the existence of continuum Alfvén waves and continuum slow waves that naturally live on those structures and phase-mix as time evolves. The apparent cross-field phase velocity is related to the spatial variation of the local Alfvén/slow frequency across the magnetic surfaces and is slower than the Alfvén/sound velocities for typical coronal conditions. Understanding the nature of the apparent cross-field propagation is important for the correct analysis of numerical simulations and the correct interpretation of observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Speeding up N-body simulations of modified gravity: Vainshtein screening models

We introduce and demonstrate the power of a method to speed up current iterative techniques for N-body modified gravity simulations. Our method is based on the observation that the accuracy of the final result is not compromised if the calculation of the fifth force becomes less accurate, but substantially faster, in high-density regions where it is relatively weak due to screening. We focus on the nDGP model which employs Vainshtein screening, and test our method by running AMR simulations in which the fifth force on the finer levels of the mesh (high density) is not obtained iteratively, but instead interpolated from coarser levels. The calculation of the standard gravity component of the force still employs the full AMR structure. We show that the impact this has on the matter power spectrum is below 1% for k < 5h/Mpc at 0z = , and even smaller at higher redshift. The impact on halo properties is also small (∼< 3% for abundance, profiles, mass; and ∼< 0.05% for positions and velocities). The method can boost the performance of modified gravity simulations by more than a factor of 10. This allows them to run on timescales similar to GR simulations and to push them to resolution levels that were previously hard to achieve.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Analytic variance estimates of Swank and Fano factors

Purpose: Variance estimates for detector energy resolution metrics can be used as stopping criteria in Monte Carlo simulations for the purpose of ensuring a small uncertainty of those metrics and for the design of variance reduction techniques. Methods: The authors derive an estimate for the variance of two energy resolution metrics, the Swank factor and the Fano factor, in terms of statistical moments that can be accumulated without significant computational overhead. The authors examine the accuracy of these two estimators and demonstrate how the estimates of the coefficient of variation of the Swank and Fano factors behave with data from a Monte Carlo simulation of an indirect x-ray imaging detector. Results: The authors' analyses suggest that the accuracy of their variance estimators is appropriate for estimating the actual variances of the Swank and Fano factors for a variety of distributions of detector outputs. Conclusions: The variance estimators derived in this work provide a computationally convenient way to estimate the error or coefficient of variation of the Swank and Fano factors during Monte Carlo simulations of radiation imaging systems.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2016

SU-A-BRA-00: Education Council Symposium: Revitalizing Your Medical Physics Classroom: Some Examples and Thoughts from the Trenches

Vic Montemayor - No one has been more passionate about improving the quality and effectiveness of the teaching of Medical Physics than Bill Hendee. It was in August of 2008 that the first AAPM Workshop on Becoming a Better Teacher of Medical Physics was held, organized and run by Bill Hendee. This was followed up in July of 2010 with a summer school on the same topic, again organized by Bill. There has been continued interest in alternate approaches to teaching medical physics since those initial gatherings. The momentum established by these workshops is made clear each year in the annual Innovation in Medical Physics Education session, which highlights work being done in all forms of medical physics education, from one-on-one residencies or classroom presentations to large-scale program revisions and on-line resources for international audiences. This symposium, presented on behalf of the Education Council, highlights the work of three finalists from past Innovation in Education sessions. Each will be presenting their approaches to and innovations in teaching medical physics. It is hoped that audience members interested in trying something new in their teaching of medical physics will find some of these ideas and approaches readily applicable to their own classrooms. Rebecca Howell - The presentation will discuss ways to maximize classroom learning, i.e., increasing the amount of material covered while also enhancing students’ understanding of the broader implications of the course topics. Specifically, the presentation will focus on two teaching methodologies, project based learning and flip learning. These teaching methods will be illustrated using an example of graduate medical physics course where both are used in conjunction with traditional lectures. Additionally, the presentation will focus on our experience implementing these methods including challenges that were overcome. Jay Burmeister - My presentation will discuss the incorporation of active learning techniques into a traditional medical physics classroom course. I will describe these techniques and how they were implemented as well as student performance before and after implementation. Student feedback indicated that these course changes improved their ability to actively assimilate the course content, thus improving their understanding of the material. Shahid Naqvi - My talk will focus on ways to help students visualize crucial concepts that lie at the core of radiation physics. Although particle tracks generated by Monte Carlo simulations have served as an indispensable visualization tool, students often struggle to resolve the underlying physics from a simultaneous jumble of tracks. We can clarify the physics by “coding” the tracks, e.g., by coloring the tracks according to their “starting” or “crossing” regions. The regionally-coded tracks when overlaid with dose distributions help the students see the elusive connection between dose, kerma and electronic disequilibrium. Tracks coded according to local energy or energy-loss rate can illustrate the need for stopping power corrections in electron beams and explain the Bragg peak in a proton beam. Coding tracks according to parent interaction type and order can clarify the often misunderstood distinction between primary and scatter dose. The students can thus see the “whole” simultaneously with the “sum of the parts,” which enhances their physical insight and creates a sustainable foundation for further learning. After the presentations the speakers and moderator will be open to questions and discussion with the audience members. Learning Objectives: Be able to explain Project-Based Learning and how can it be incorporated into a Medical Physics classroom. Be able to explain Flipped Learning and how can it be incorporated into a Medical Physics classroom. Be able to explain active-learning strategies for the teaching of Medical Physics. Be able to explain how Monte Carlo simulations can be used to deepen a student’s understanding of radiation physics and dosimetry.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2016

Electron energy balance and ionization in the channel of a stationary plasma thruster

The paper presents results of numerical simulations of the electron dynamics in the field of the azimuthal and longitudinal waves excited in the channel of a stationary plasma thruster (SPT). The simulations are based on the experimentally determined wave characteristics. The simulation results show that the azimuthal wave displayed as ionization instability enhances electron transport along the thruster channel. It is established that the electron transport rate in the azimuthal wave increases as compared to the rate of diffusion caused by electron scattering from neutral atoms in proportion to the ratio between the times of electron− neutral collisions responsible for ionization and elastic electron scattering, respectively. An expression governing the plasma conductivity is derived with allowance for electron interaction with the azimuthal wave. The Hall parameter, the electron component of the discharge current, and the electron heating power in the thruster channel are calculated for two model SPTs operating with krypton and xenon. The simulation results agree well with the results of experimental studies of these two SPTs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Numerical Study of the Voltage Waveform Effect on the Spatiotemporal Characteristics of a Dielectric Barrier Microdischarge in Argon

The effect of the shape of the feeding voltage on the spatiotemporal characteristics of an atmospheric- pressure barrier microdischarge in argon is demonstrated using numerical simulations based on an extended fluid model. Results of simulations performed for sinusoidal and square feeding voltages are analyzed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Experimental and Computational Study of a High-Speed Collision of an Aluminum Striker with a Solid and Mesh Obstacle

The results of computational and experimental studies of peculiarities of the process of impact between the striker and obstacle, which imitates shield protection of spacecrafts, have been considered and pictures of fragmentation of the striker and obstacle have been obtained. The numerical calculations were carried out using the method of smooth particle hydrodynamics (SPH). In the experiments, the striker was accelerated using a two-stage light-gas gun; the state of the obstacle and striker at the moment of their interaction was controlled using an X-ray and high-speed video camera. The recording equipment was switched on using induction cross sections. The comparative analysis of the results of numerical simulations and experiments have been considered.

42 ENGINEERING↗
DOE OSTI2017

Efficiency of Synchrotron Radiation from Rotation-powered Pulsars

Synchrotron radiation is widely considered to be the origin of the pulsed non-thermal emissions from rotation-powered pulsars in optical and X-ray bands. In this paper, we study the synchrotron radiation emitted by the created electron and positron pairs in the pulsar magnetosphere to constrain the energy conversion efficiency from the Poynting flux to the particle energy flux. We model two pair creation processes, two-photon collision, which efficiently works in young γ -ray pulsars (≲10{sup 6} year), and magnetic pair creation, which is the dominant process to supply pairs in old pulsars (≳10{sup 6} year). Using the analytical model, we derive the maximum synchrotron luminosity as a function of the energy conversion efficiency. From the comparison with observations, we find that the energy conversion efficiency to the accelerated particles should be an order of unity in the magnetosphere, even though we make a number of the optimistic assumptions to enlarge the synchrotron luminosity. In order to explain the luminosity of the non-thermal X-ray/optical emission from pulsars with low spin-down luminosity L {sub sd} ≲ 10{sup 34} erg s{sup −1}, non-dipole magnetic field components should be dominant at the emission region. For the γ -ray pulsars with L {sub sd} ≲ 10{sup 35} erg s{sup −1}, observed γ -ray to X-ray and optical flux ratios are much higher than the flux ratio between curvature and the synchrotron radiations. We discuss some possibilities such as the coexistence of multiple accelerators in the magnetosphere as suggested from the recent numerical simulation results. The obtained maximum luminosity would be useful to select observational targets in X-ray and optical bands.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Spatially Continuous Depletion Algorithm for Monte Carlo Simulations

To correctly predict reactor behavior during cycle operations, the evolution of nuclide number densities throughout the core must be accurately modeled. The time-varying spatial distribution of nuclide number densities is typically resolved by discretizing the Monte Carlo geometry into smaller cells over which number densities are assumed to be spatially invariant. The nuclide number densities in these smaller cells are integrated through time using reaction rate tallies on the same discretized geometry. However, detailed distributions of nuclide number densities in a full three dimensional simulation can require a prohibitive amount of tallies, and the spatial discretization of the base geometry makes coupling to external multi-physics tools difficult. In this paper a method for solving for spatially continuous number density distributions during depletion calculations will be described. The spatially continuous number densities can be used in the transport method proposed by Brown and Martin which allows for transporting neutrons through a material with continuously varying properties such as temperature and nuclide number densities. Coupled with the ability of Functional Expansion Tallies (FETs) to represent tallied quantities as continuous functions, it is possible to both solve for and make use of spatially continuous nuclide number densities. The need for this capability was alluded to by Brown et. al., but no solution has yet been proposed. With a continuous depletion method, recent work which utilized FETs and continuous material tracking to incorporate multi-physics feedback in Monte Carlo simulations can be extended to simulations that include depletion analysis. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2016

CTF Parallel Performance Improvements

CTF is a modernized version of COBRA-TF (Coolant Boiling in Rod Arrays-Two Fluid), which is a thermal hydraulic (T/H) sub-channel code based on a two-phase, two-fluid model that formulates the conservation equations of mass, energy, and momentum for three fields of vapor, continuous liquid, and entrained liquid droplets. Under the Consortium for Advanced Simulation of Light Water Reactors (CASL) program sponsored by the U.S. Department of Energy, CTF is being further developed and improved for modeling and simulation of light water reactors. The code improvements for predicting T/H responses of pressurized water reactor (PWR) cores are focused on its applications to specific challenge problems deemed significant to the industry, such as departure from nucleate boiling (DNB) and crud-induced power shift. The code is also being used for T/H feedback in normal operating condition depletions in PWR and boiling water reactor designs in the reactor core simulator, Virtual Environment for Reactor Applications Core Simulator (VERA-CS), being developed by CASL. CTF modifications under the CASL program include software optimization, new closure models, and parallelization for modeling full reactor core T/H responses. The purpose of this work was to improve the parallel performance of CTF. CTF was originally parallelized via a domain-decomposition strategy using the Message Passing Interface (MPI), in which the basis of a solution domain was a physical assembly in the model. In other words, a user who wanted to perform a parallel solution with CTF had to provide the code one processor per assembly in the model (e.g., 56 cores for a 56 assembly quarter-core model, or 193 cores for a 193 assembly full-core model). This design was rigid, demanding a very specific number of processors for a model; furthermore, when CTF was run in VERA-CS-which typically uses thousands of processors- many processors sat idle during the CTF solve because CTF was unable to use any more processors than there are physical assemblies. This was not necessarily a problem if the CTF execution speed was very fast, but the average time spent in a CTF solve for prototypical PWR operating conditions was significant (i.e., over 20 % of the simulation wall time was being spent in CTF for a standard core depletion simulation). Larger, more computationally intensive simulations performed in support of the CASL DNB challenge problem demonstrated a need for improvements in this original parallelization approach. The specific goal of this work was to allow CTF to use more processors than there are physical assemblies in the model, which would not only speed up CTF and improve the overall efficiency of VERA-CS but also give users more-flexible options for executing their simulations. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2017

Tersoff Benchmarking of Be-C-H Interatomic Potential

To demonstrate prolonged fusion power in ITERs (International Thermonuclear Experimental Reactors), the selection of appropriate plasma facing materials (PFMs) is a vital issue that requires the use of interatomic potentials to simulate the plasma interactions with the PFMs. So it is no surprise that the use of computer codes such as LAMMPS to perform studies in order to determine proper PFMs is necessary based on the complexity of such studies. In recent Tokamak advancements, beryllium was chosen as a first wall material due to its low Z (lower energy loss due to plasma contamination), its low tritium inventory and its oxygen gettering quality. With its high melting point and ability to withstand high heat fluxes, beryllium makes a fine choice for a PFM. In this work, Be-C-H (beryllium-carbon-hydrogen) bonds were investigated in a LAMMPS Tersoff simulation, where the interatomic potentials for Be-C-H were converted to LAMMPS Tersoff format, and were validated by calculating the cohesive energy (E{sub coh}), lattice constant (a), and atomic volume (V{sub 0}) for comparison to Density Functional Theory (DFT), Analytical Bond-Order Potentials (ABOP) and experiment potentials. This benchmark serves as the validation of performing Be-C-H bombardment simulations using LAMMPS.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Uncertainty Quantification and Sensitivity Analysis for Coupled VERA-CS and FRAPCON

The current nuclear power plant fleet in the U.S. is facing numerous challenges economically ranging from low natural gas price to the rapidly growing renewable energy sources such as wind and solar power. In order for the nuclear industry to keep the viability of its operating fleet, new innovations have to be introduced to lower the costs to operate the fleet. One area that could contribute to the cost reduction is to lower the fuel cycle cost and improve the plant safety through initiatives such as higher burnup fuel and accident tolerant fuel. At the same time, the nuclear industry has to address future challenges on the horizon such as load following and other flexible operating strategies as well as those associated with new regulations such as 10 CFR 50.69 and 10 CFR50.46c etc. All these require coupled analysis of core neutronics, thermal hydraulics and fuel performance, which necessitates the development of a multi-physics toolkit. Since different models, assumptions, and fidelities have gone into various computer codes, how to consistently propagate uncertainties in a multi-physics environment warrants careful evaluation. The U.S. DOE's Light Water Reactor Sustainability (LWRS) Program initiated a project to develop analytical capabilities to support the industry to address these challenges. The general idea behind the initiative was the development of an integrated multiphysics evaluation model. This integrated evaluation model is called LOTUS which stands for LOCA Toolkit for the U.S. Light Water Reactors. The focus of LOTUS is to establish the automation of the interface between the five disciplines. LOTUS framework ensures that uncertainties can be easily propagated to solve the problems forward and backward (inverse) with large numbers of simulations. The objective of this work is to perform uncertainty quantification (UQ) and sensitivity analysis (SA) on a coupling of CASL's Virtual Environment for Reactor Applications-Core Simulator (VERA-CS) and NRC's fuel performance code FRAPCON. The work presented herein is a subset of the overall LOTUS framework. The figures of merit (FOM) selected for analysis are minimum departure from nucleate boiling ratio (MDNBR), maximum fuel centerline temperature (MFCT), and the gap conductance at peak power (GCPP). All FOM are obtained at the end of the of a depletion case for a half cycle lasting 263 days. The minimum is taken over all spatial location for each pin at the end of the half cycle. This value is calculated in VERACS alone and has been analyzed in previous work by Brown and Zhang. The MFCT is available in both VERA-CS and FRAPCON, but is calculated via different methods (see Model Inconsistencies section for more details). For this work, only the FRAPCON values are reported due to the code's maturity from extensive validation in fuel performance. Previous work includes SA on a standalone nuclear fuel performance code VTT-modified FRAPCON- 3.4 by Ikonen. The GCPP is the gap conductance at the location corresponding to maximum power in the fuel rod. This value is available in FRAPCON alone and highlights the capability of LOTUS to obtain unique fuel performance metrics for each pin in a core simulation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2018

Numerical Modeling of Fiber-Reinforced Metal Matrix Composite Processing by the Liquid Route: Literature Contribution

One of the technologies used to produce metal matrix composites (MMCs) is liquid route processing. One solution is to inject a liquid metal under pressure or at constant rate through a fibrous preform. This foundry technique overcomes the problem of the wettability of ceramic fibers by liquid metal. The liquid route can also be used to produce semiproducts by coating a filament with a molten metal. These processes involve physical phenomena combined with mass and heat transfer and phase change. The phase change phenomena related to solidification and also to the melting of the metal during the process notably result in modifications to the permeability of porous media, in gaps in impregnation, in the appearance of defects (porosities), and in segregation in the final product. In this article, we provide a state-of-the-art review of numerical models and simulation developed to study these physical phenomena involved in MMC processing by the liquid route.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Promising Hard Carbon Coatings on Cu Substrates: Corrosion and Tribological Performance with Theoretical Aspect

Protecting the surface of metals and alloys against corrosion and wear is of abundant importance owing to their widespread applications. In the present work, we report the improved anticorrosion and tribo-mechanical performance of copper (Cu) by a hard carbon (HC) coating synthesized in different pyrolysis temperature. Structural and surface characterization with roughness measurements was systematically investigated using various techniques. Effect of pyrolysis temperature on the corrosion behavior of coated Cu substrates in 0.6 M NaCl solution was evaluated via electrochemical impedance spectroscopy, potentiodynamic polarization. Pin-on-disk wear test of coated Cu substrate showed the influence of the pyrolysis temperature on the wear resistance performance of the HC coatings. According to the obtained results, it could be concluded that the HC coatings synthesized at 1100 °C revealed an enhanced comprehensive performance, revealing their possible utilization as a protective coating for Cu substrates in chloride environment. Monte Carlo simulations have been utilized to elucidate the interaction between the Cu surface and HC coatings.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Strain Rate Dependency of Simulated Welding Residual Stresses

In the present work, the influence of taking into consideration the strain rate-dependent nature of steel S355 during a calculation of the welding residual stresses through a finite element simulation is investigated. The Perzyna material model is calibrated using experimental values found in the literature and is introduced to a validated weld simulation model, where the strain rate dependency had not been taken into consideration before this study. The calculated profiles of the welding residual stresses, for strain rate-dependent and independent behavior, are then compared with experimentally measured profiles. The results of this first-step analysis show that taking into consideration strain rate dependency during a welding simulation of the investigated S355 has non-negligible influence on the calculated welding residual stresses.

36 MATERIALS SCIENCE↗