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,171 records · Page 65
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

COSMOLOGICAL MAGNETOHYDRODYNAMIC SIMULATIONS OF GALAXY CLUSTER RADIO RELICS: INSIGHTS AND WARNINGS FOR OBSERVATIONS

Non-thermal radio emission from cosmic-ray electrons in the vicinity of merging galaxy clusters is an important tracer of cluster merger activity, and is the result of complex physical processes that involve magnetic fields, particle acceleration, gas dynamics, and radiation. In particular, objects known as radio relics are thought to be the result of shock-accelerated electrons that, when embedded in a magnetic field, emit synchrotron radiation in the radio wavelengths. In order to properly model this emission, we utilize the adaptive mesh refinement simulation of the magnetohydrodynamic evolution of a galaxy cluster from cosmological initial conditions. We locate shock fronts and apply models of cosmic-ray electron acceleration that are then input into radio emission models. We have determined the thermodynamic properties of this radio-emitting plasma and constructed synthetic radio observations to compare observed galaxy clusters. We find a significant dependence of the observed morphology and radio relic properties on the viewing angle of the cluster, raising concerns regarding the interpretation of observed radio features in clusters. We also find that a given shock should not be characterized by a single Mach number. We find that the bulk of the radio emission comes from gas with T > 5 Multiplication-Sign 10{sup 7} K, {rho} {approx} 10{sup -28}-10{sup -27} g cm{sup -3}, with magnetic field strengths of 0.1-1.0 {mu}G, and shock Mach numbers of M {approx} 3-6. We present an analysis of the radio spectral index which suggests that the spatial variation of the spectral index can mimic synchrotron aging. Finally, we examine the polarization fraction and position angle of the simulated radio features, and compare to observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

THE FORCE-FREE MAGNETOSPHERE OF A ROTATING BLACK HOLE

We revisit the Blandford-Znajek process and solve the fundamental equation that governs the structure of the steady-state force-free magnetosphere around a Kerr black hole. The solution depends on the distributions of the magnetic field angular velocity {omega} and the poloidal electric current I. These are not arbitrary. They are determined self-consistently by requiring that magnetic field lines cross smoothly the two singular surfaces of the problem: the inner ''light surface'' located inside the ergosphere and the outer ''light surface'' which is the generalization of the pulsar light cylinder. We find the solution for the simplest possible magnetic field configuration, the split monopole, through a numerical iterative relaxation method analogous to the one that yields the structure of the steady-state axisymmetric force-free pulsar magnetosphere. We obtain the rate of electromagnetic extraction of energy and confirm the results of Blandford and Znajek and of previous time-dependent simulations. Furthermore, we discuss the physical applicability of magnetic field configurations that do not cross both ''light surfaces''.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Numerical simulations of active region scale flux emergence: From spot formation to decay

We present numerical simulations of active region scale flux emergence covering a time span of up to 6 days. Flux emergence is driven by a bottom boundary condition that advects a semi-torus of magnetic field with 1.7 × 10{sup 22} Mx flux into the computational domain. The simulations show that, even in the absence of twist, the magnetic flux is able the rise through the upper 15.5 Mm of the convection zone and emerge into the photosphere to form spots. We find that spot formation is sensitive to the persistence of upflows at the bottom boundary footpoints, i.e., a continuing upflow would prevent spot formation. In addition, the presence of a torus-aligned flow (such flow into the retrograde direction is expected from angular momentum conservation during the rise of flux ropes through the convection zone) leads to a significant asymmetry between the pair of spots, with the spot corresponding to the leading spot on the Sun being more axisymmetric and coherent, but also forming with a delay relative to the following spot. The spot formation phase transitions directly into a decay phase. Subsurface flows fragment the magnetic field and lead to intrusions of almost field free plasma underneath the photosphere. When such intrusions reach photospheric layers, the spot fragments. The timescale for spot decay is comparable to the longest convective timescales present in the simulation domain. We find that the dispersal of flux from a simulated spot in the first two days of the decay phase is consistent with self-similar decay by turbulent diffusion.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

QUANTITATIVE EVALUATION OF THE HYPOTHESIS THAT BL LACERTAE OBJECTS ARE QSO REMNANTS

We evaluate with numerical simulations the hypothesis that BL Lacertae objects (BLLs) are the remnants of quasi-stellar objects. This hypothesis is based on their highly peculiar redshift evolution. They have a comoving space density that increases with decreasing redshift, contrary to all other active galactic nuclei. We assume that relativistic jets are below detection in young radio-quiet quasars and increase in strength with cosmic time so that they eventually are detected as BLLs. Our numerical simulations fit very well the observed redshift distributions of BLLs. There are strong indications that only the high-synchrotron-peaked BLLs could be QSO remnants.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

MO-E-18C-04: Advanced Computer Simulation and Visualization Tools for Enhanced Understanding of Core Medical Physics Concepts

Purpose: Most medical physics programs emphasize proficiency in routine clinical calculations and QA. The formulaic aspect of these calculations and prescriptive nature of measurement protocols obviate the need to frequently apply basic physical principles, which, therefore, gradually decay away from memory. E.g. few students appreciate the role of electron transport in photon dose, making it difficult to understand key concepts such as dose buildup, electronic disequilibrium effects and Bragg-Gray theory. These conceptual deficiencies manifest when the physicist encounters a new system, requiring knowledge beyond routine activities. Methods: Two interactive computer simulation tools are developed to facilitate deeper learning of physical principles. One is a Monte Carlo code written with a strong educational aspect. The code can “label” regions and interactions to highlight specific aspects of the physics, e.g., certain regions can be designated as “starters” or “crossers,” and any interaction type can be turned on and off. Full 3D tracks with specific portions highlighted further enhance the visualization of radiation transport problems. The second code calculates and displays trajectories of a collection electrons under arbitrary space/time dependent Lorentz force using relativistic kinematics. Results: Using the Monte Carlo code, the student can interactively study photon and electron transport through visualization of dose components, particle tracks, and interaction types. The code can, for instance, be used to study kerma-dose relationship, explore electronic disequilibrium near interfaces, or visualize kernels by using interaction forcing. The electromagnetic simulator enables the student to explore accelerating mechanisms and particle optics in devices such as cyclotrons and linacs. Conclusion: The proposed tools are designed to enhance understanding of abstract concepts by highlighting various aspects of the physics. The simulations serve as virtual experiments that give deeper and long lasting understanding of core principles. The student can then make sound judgements in novel situations encountered beyond routine clinical activities.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

Fast multipurpose Monte Carlo simulation for proton therapy using multi- and many-core CPU architectures

Purpose: Accuracy in proton therapy treatment planning can be improved using Monte Carlo (MC) simulations. However the long computation time of such methods hinders their use in clinical routine. This work aims to develop a fast multipurpose Monte Carlo simulation tool for proton therapy using massively parallel central processing unit (CPU) architectures. Methods: A new Monte Carlo, called MCsquare (many-core Monte Carlo), has been designed and optimized for the last generation of Intel Xeon processors and Intel Xeon Phi coprocessors. These massively parallel architectures offer the flexibility and the computational power suitable to MC methods. The class-II condensed history algorithm of MCsquare provides a fast and yet accurate method of simulating heavy charged particles such as protons, deuterons, and alphas inside voxelized geometries. Hard ionizations, with energy losses above a user-specified threshold, are simulated individually while soft events are regrouped in a multiple scattering theory. Elastic and inelastic nuclear interactions are sampled from ICRU 63 differential cross sections, thereby allowing for the computation of prompt gamma emission profiles. MCsquare has been benchmarked with the GATE/GEANT4 Monte Carlo application for homogeneous and heterogeneous geometries. Results: Comparisons with GATE/GEANT4 for various geometries show deviations within 2%–1 mm. In spite of the limited memory bandwidth of the coprocessor simulation time is below 25 s for 10{sup 7} primary 200 MeV protons in average soft tissues using all Xeon Phi and CPU resources embedded in a single desktop unit. Conclusions: MCsquare exploits the flexibility of CPU architectures to provide a multipurpose MC simulation tool. Optimized code enables the use of accurate MC calculation within a reasonable computation time, adequate for clinical practice. MCsquare also simulates prompt gamma emission and can thus be used also for in vivo range verification.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2018

Study of Lasing in Liquid-Crystal Systems with Microgratings

Waveguide lasing in a layer of nematic liquid crystal (NLC) doped with a dye has been experimentally observed. A chromium micrograting with a period of 1.76 μm has been deposited on the surface one of the cell substrates to provide a distributed feedback (DFB) and partially extract laser radiation to a quartz substrate. The DFB provides lasing in the ninth diffraction order for waveguide TM modes. Laser radiation is observed at the output of the substrate end face at an angle of 67.0° ± 1.5° with respect to the normal to the waveguide plane. An increase in voltage across the micrograting electrodes leads to red-shifted multimode lasing. A numerical simulation of a structure imitating the experimental cell has shown good agreement between the calculation results and experimental data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

THREE-DIMENSIONAL SIMULATIONS OF LONG DURATION GAMMA-RAY BURST JETS: TIMESCALES FROM VARIABLE ENGINES

Gamma-ray burst (GRB) light curves are characterized by marked variability, each showing unique properties. The origin of this variability, at least for a fraction of long GRBs, may be the result of an unsteady central engine. It is thus important to study the effects that an episodic central engine has on the jet propagation and, eventually, on the prompt emission within the collapsar scenario. Thus, in this study we follow the interaction of pulsed outflows with their progenitor stars with hydrodynamic numerical simulations in both two and three dimensions. We show that the propagation of unsteady jets is affected by the interaction with the progenitor material well after the break-out time, especially for jets with long quiescent times comparable to or larger than a second. We also show that this interaction can lead to an asymmetric behavior in which pulse durations and quiescent periods are systematically different. After the pulsed jets drill through the progenitor and the interstellar medium, we find that, on average, the quiescent epochs last longer than the pulses (even in simulations with symmetrical active and quiescent engine times). This could explain the asymmetry detected in the light curves of long quiescent time GRBs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Nonuniversal star formation efficiency in turbulent ISM

We present a study of a star formation prescription in which star formation efficiency (SFE) depends on local gas density and turbulent velocity dispersion, as suggested by direct simulations of SF in turbulent giant molecular clouds (GMCs). We test the model using a simulation of an isolated Milky-Way-sized galaxy with a self-consistent treatment of turbulence on unresolved scales. We show that this prescription predicts a wide variation of local SFE per free-fall time, ϵ{sub ff} ∼ 0.1%–10%, and gas depletion time, t{sub dep} ∼ 0.1–10 Gyr. In addition, it predicts an effective density threshold for star formation due to suppression of ϵ{sub ff} in warm diffuse gas stabilized by thermal pressure. We show that the model predicts star formation rates (SFRs) in agreement with observations from the scales of individual star-forming regions to the kiloparsec scales. This agreement is nontrivial, as the model was not tuned in any way and the predicted SFRs on all scales are determined by the distribution of the GMC-scale densities and turbulent velocities σ in the cold gas within the galaxy, which is shaped by galactic dynamics. The broad agreement of the star formation prescription calibrated in the GMC-scale simulations with observations both gives credence to such simulations and promises to put star formation modeling in galaxy formation simulations on a much firmer theoretical footing.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Common Envelope Light Curves. I. Grid-code Module Calibration

The common envelope (CE) binary interaction occurs when a star transfers mass onto a companion that cannot fully accrete it. The interaction can lead to a merger of the two objects or to a close binary. The CE interaction is the gateway of all evolved compact binaries, all stellar mergers, and likely many of the stellar transients witnessed to date. CE simulations are needed to understand this interaction and to interpret stars and binaries thought to be the byproduct of this stage. At this time, simulations are unable to reproduce the few observational data available and several ideas have been put forward to address their shortcomings. The need for more definitive simulation validation is pressing and is already being fulfilled by observations from time-domain surveys. In this article, we present an initial method and its implementation for post-processing grid-based CE simulations to produce the light curve so as to compare simulations with upcoming observations. Here we implemented a zeroth order method to calculate the light emitted from CE hydrodynamic simulations carried out with the 3D hydrodynamic code Enzo used in unigrid mode. The code implements an approach for the computation of luminosity in both optically thick and optically thin regimes and is tested using the first 135 days of the CE simulation of Passy et al., where a 0.8 M {sub ⊙} red giant branch star interacts with a 0.6 M {sub ⊙} companion. This code is used to highlight two large obstacles that need to be overcome before realistic light curves can be calculated. We explain the nature of these problems and the attempted solutions and approximations in full detail to enable the next step to be identified and implemented. We also discuss our simulation in relation to recent data of transients identified as CE interactions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Optimization design of helium circulator diffuser in HTR-PM - 18488

Circulator is very important in HTR-PM. It needs to supply enough pressure rise to circulate cooled helium in the primary loop of HTR-PM to bring heat out. Its pressure rise requirement is the highest among the worldwide circulators. Diffuser is one main component to produce pressure rise for circulator. And the most particular for HTR-PM is that the helium must flow out from two round tubes after flowing through the diffuser, which is very different from other circulators. This special export structure increases flow loss and decreases diffuser efficiency. So it is necessary to optimize diffuser to satisfy high pressure rise requirement of HTR-PM circulator with the special export structure. Four geometrical variables of diffuser have been studied and optimized. Based on the flow and performance analysis, the optimum structure was chosen. Then the detailed simulation was done by combining it with impeller to obtain performance characteristics of the whole circulator. Finally, the simulation results were compared with experimental results. It was obtained that simulation pressure rise was 201 kPa and experimental result was 210 kPa, which verified that optimization design of diffuser was effective. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2016

A Neutron Streaming Problem to Test Rattlesnake Methods for TREAT

This work demonstrates the progress made at Idaho National Laboratory in the neutronic modeling of neutron streaming in the TREAT reactor core. The streaming in the slotted elements is important because of its effect on the axial distribution and energy deposition in the experiment vehicle. Two methods implemented in Rattlesnake: 1) tensor diffusion coefficients and, 2) super-homogenization, yield much improved power distributions using a diffusion solution with homogenized fuel elements. The RMS difference in the power distribution decreases below 0.661% compared to a reference Monte Carlo simulation. The maximum differences are within 1.5% of this reference. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2016

Experimental and CFD Analysis of Condensation with Non-Condensable Gas in Multi-Tubes

After the Fukushima accident, the importance of safety system that relies on passive cooling for decay heat removal has been reemphasized. The passive containment cooling system (PCCS) is the system for cooling down a containment and lowering a containment pressure even under the station blackout. In investigating a PCCS, steam condensation in the presence of non-condensable gas is the dominant heat and mass transfer phenomenon. The non-condensable gases degrade the heat removal performance through the PCCS heat exchangers (HXs) because they are accumulated near the cooling surface and deteriorate the approach of steam by diffusion to the surface. Therefore, in order to understand steam condensation phenomenon in the PCCS, we conducted multi-tube condensation experiment with air. And we also analyzed the results by CFD simulation for the CFD-based approach to condensation phenomena prediction. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2017

Using the MCNP6.2 Correlated Fission Multiplicity Models, CGMF and FREYA

The culmination of a multi-year, collaborative effort by researchers at Los Alamos National Laboratory (LANL), Lawrence Berkeley National Laboratory (LBNL), Lawrence Livermore National Laboratory (LLNL) and the University of Michigan, led to the delivery of two new correlated fission multiplicity models available within the latest release of MCNP{sup R} version 6.2. The needs within the nuclear nonproliferation community, particularly in simulating and predicting the neutron and gamma-ray signatures of special nuclear materials (SNM), provided the motivation for this major collaborative effort. From LANL, the CGMF code has been implemented as an upgrade to the original Cascading Gamma-ray Multiplicity (CGM) code. Additionally, the LLNL Fission Library has been upgraded from the original implementation released in prior versions of MCNP6/X to include the LBNL/LLNL Fission Reaction Event Yield Algorithm (FREYA) code. Both new codes are capable of simulating fission events with correlated neutron and gamma-ray emissions. Due to the release of these new models within MCNP6.2, the purpose of this paper is to describe the differences between these new models and existing capabilities, demonstrate their use in the context of simulation inputs and outputs, and to discuss some of the ongoing work directed toward the next generation of these features.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

Radiolytic Gas Bubbles' Model and its Application to AHR's Transient Neutronics Analysis

Aqueous Homogeneous Reactor (AHR) receives much attention recently because its advantages in medical isotopes production. Fissile material and moderator is mixed in the fuel solution, leads to big negative feedback and generation of radiolytic gas bubbles; the fluidity of fuel and movement of bubbles bring out big challenge of the neutronics analysis methods. This paper aims to establish a reliable and effective gas bubbles' model, firstly a coupled Neutronics-CFD method is introduced to recognize gas bubbles' behavior characteristics, then it is simplified into mathematical expression and applied to Transient Calculation Code of Aqueous Homogeneous Reactor (TCCAHR), finally a typical AHR's transient process is simulated and compared with the experiment.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2018

Effect of Low Temperature on a 4 W/60 K Pulse-Tube Cryocooler for Cooling HgCdTe Detector

Temperature is an extremely important parameter for the material of the space-borne infrared detector. To cool an HgCdTe-infrared detector, a Stirling-type pulse-tube cryocooler (PTC) has been developed based on a great deal of numerical simulations, which are performed to investigate the thermodynamic behaviors of the PTC. The effects of different low temperatures are presented to analyze different energy flows, losses, phase shifts, and impedance matching of the PTC at a temperature range of 40–120 K, where woven wire screens are used. Finally, a high-efficiency coaxial PTC has been designed, built, and tested, operating around 60 K after a number of theoretical and experimental studies. The PTC can offer a no-load refrigeration temperature of 40 K with an input electric power of 150 W, and a cooling power of 4 W at 60 K is obtained with Carnot efficiency of 12%. In addition, a comparative study of simulation and experiment has been carried out, and some studies on reject temperatures have been presented for a thorough understanding of the PTC system.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Capturing Hidden Geochemical Anomalies in Scarce Data by Fractal Analysis and Stochastic Modeling

Fractal/multifractal modeling is a widely used geomathematical approach to capturing different populations in geochemical mapping. The rationale of this methodology is based on empirical frequency density functions attained from global or local distributions. This approach is quite popular because of its simplicity and versatility; it accounts for the frequency and spatial distribution of geochemical data considering self-similarity across a range of scales. Using this technique for detection of geochemical anomalies in scarce data, however, is problematic and can lead to systematic bias in the characterization of the underlying populations. In this paper, an innovative technique is presented that provides good results without a priori assumptions. A simulation approach is adopted for fractal analysis by generating different possible distribution scenarios for the variable under study to reveal the underlying populations that are frequently hidden due to lack of data. The proposed technique is called the global simulated size–number method, and it is validated in a case study with two synthetic datasets and another case study with real dataset from the Ushtagan gold deposit in northeast Kazakhstan.

58 GEOSCIENCES↗
DOE OSTI2017

Parameter Calibration of GTN Damage Model and Formability Analysis of 22MnB5 in Hot Forming Process

Hot forming of high strength steel at elevated temperatures is an attractive technology to achieve the lightweight of vehicle body. The mechanical behavior of boron steel 22MnB5 strongly depends on the variation of temperature which makes the process design more difficult. In this paper, the Gurson–Tvergaard–Needleman (GTN) model is used to study the formability of 22MnB5 sheet at different temperatures. Firstly, the rheological behavior of 22MnB5 is analyzed through a series of hot tensile tests at a temperature range of 600-800 °C. Then, a detailed process to calibrate the damage parameters is given based on the response surface methodology and genetic algorithm method. The GTN model together with the damage parameters calibrated is then implemented to simulate the deformation and damage evolution of 22MnB5 in the process of high-temperature Nakazima test. The capability of the GTN model as a suitable tool to evaluate the sheet formability is confirmed by comparing experimental and calculated results. Finally, as a practical application, the forming limit diagram of 22MnB5 at 700 °C is constructed using the Nakazima simulation and Marciniak–Kuczynski (M–K) model, respectively. And the simulation integrated GTN model shows a higher reliability by comparing the predicted results of these two approaches with the experimental ones.

36 MATERIALS SCIENCE↗