Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Code”

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,585 records · Page 88
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

The FLUKA code: Overview and new developments

The FLUKA Monte Carlo Radiation Transport and Interaction code package is widely used to simulate the interaction of particles with matter in a variety of fields, including high energy physics, space radiation, medical applications, radiation protection and shielding assessments, accelerator studies, astrophysical studies and well logging. This paper gives a brief overview of the FLUKA program and describes recent developments, in particular, improvements in the modelling of particle interactions and transport are described in detail. In addition, an overview of selected applications is given.

Ballarini, Francesca↗
DOE OSTI2025

Status of GPU capabilities within the Shift Monte Carlo radiation transport code

Shift is a general-purpose Monte Carlo (MC) radiation transport code for fission, fusion, and national security applications. Shift has been adapted to efficiently run on GPUs in order to leverage leadership-class supercomputers. This work presents Shift’s current GPU capabilities. These include core radiation transport capabilities for eigenvalue and fixed-source simulations, and support for non-uniform domain decomposition, Doppler broadening, free-gas elastic scattering, general-purpose geometry, hybrid MC/deterministic transport, and depletion. Transport results demonstrate a 2–5× GPU-to-CPU speedup on a per-node basis for an eigenvalue problem on the Frontier supercomputer and a 28× speedup for a fixed-source problem on the Summit supercomputer.

Biondo, Elliott [ORNL] (ORCID:0000000290881360)↗
DOE OSTI2014

Spin wave based parallel logic operations for binary data coded with domain walls

We numerically investigate the feasibility of spin wave (SW) based parallel logic operations, where the phase of SW packet (SWP) is exploited as a state variable and the phase shift caused by the interaction with domain wall (DW) is utilized as a logic inversion functionality. A designed functional element consists of parallel ferromagnetic nanowires (6 nm-thick, 36 nm-width, 5120 nm-length, and 200 nm separation) with the perpendicular magnetization and sub-μm scale overlaid conductors. The logic outputs for binary data, coded with the existence (“1”) or absence (“0”) of the DW, are inductively read out from interferometric aspect of the superposed SWPs, one of them propagating through the stored data area. A practical exclusive-or operation, based on 2π periodicity in the phase logic, is demonstrated for the individual nanowire with an order of different output voltage V{sub out}, depending on the logic output for the stored data. The inductive output from the two nanowires exhibits well defined three different signal levels, corresponding to the information distance (Hamming distance) between 2-bit data stored in the multiple nanowires.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

Pressure distribution based optimization of phase-coded acoustical vortices

Based on the acoustic radiation of point source, the physical mechanism of phase-coded acoustical vortices is investigated with formulae derivations of acoustic pressure and vibration velocity. Various factors that affect the optimization of acoustical vortices are analyzed. Numerical simulations of the axial, radial, and circular pressure distributions are performed with different source numbers, frequencies, and axial distances. The results prove that the acoustic pressure of acoustical vortices is linearly proportional to the source number, and lower fluctuations of circular pressure distributions can be produced for more sources. With the increase of source frequency, the acoustic pressure of acoustical vortices increases accordingly with decreased vortex radius. Meanwhile, increased vortex radius with reduced acoustic pressure is also achieved for longer axial distance. With the 6-source experimental system, circular and radial pressure distributions at various frequencies and axial distances have been measured, which have good agreements with the results of numerical simulations. The favorable results of acoustic pressure distributions provide theoretical basis for further studies of acoustical vortices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

A model code for the radiative theta pinch

A model for the theta pinch is presented with three modelled phases of radial inward shock phase, reflected shock phase, and a final pinch phase. The governing equations for the phases are derived incorporating thermodynamics and radiation and radiation-coupled dynamics in the pinch phase. A code is written incorporating correction for the effects of transit delay of small disturbing speeds and the effects of plasma self-absorption on the radiation. Two model parameters are incorporated into the model, the coupling coefficient f between the primary loop current and the induced plasma current and the mass swept up factor f{sub m}. These values are taken from experiments carried out in the Chulalongkorn theta pinch.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Analysis of stimulated Raman backscatter and stimulated Brillouin backscatter in experiments performed on SG-III prototype facility with a spectral analysis code

Experiments about the observations of stimulated Raman backscatter (SRS) and stimulated Brillouin backscatter (SBS) in Hohlraum were performed on Shenguang-III (SG-III) prototype facility for the first time in 2011. In this paper, relevant experimental results are analyzed for the first time with a one-dimension spectral analysis code, which is developed to study the coexistent process of SRS and SBS in Hohlraum plasma condition. Spectral features of the backscattered light are discussed with different plasma parameters. In the case of empty Hohlraum experiments, simulation results indicate that SBS, which grows fast at the energy deposition region near the Hohlraum wall, is the dominant instability process. The time resolved spectra of SRS and SBS are numerically obtained, which agree with the experimental observations. For the gas-filled Hohlraum experiments, simulation results show that SBS grows fastest in Au plasma and amplifies convectively in C{sub 5}H{sub 12} gas, whereas SRS mainly grows in the high density region of the C{sub 5}H{sub 12} gas. Gain spectra and the spectra of backscattered light are simulated along the ray path, which clearly show the location where the intensity of scattered light with a certain wavelength increases. This work is helpful to comprehend the observed spectral features of SRS and SBS. The experiments and relevant analysis provide references for the ignition target design in future.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Overcoming inherent magnetic instability, preventing spin canting and magnetic coding in an assembly of ferrimagnetic nanoparticles

The authors find that for mechanically milled Ni{sub 0.5}Zn{sub 0.5}Fe{sub 2}O{sub 4} (∼10 nm), the mechanical strain induced enhancement of anisotropy energy helps to retain stable magnetic order. The reduction of magnetization can be prevented by keeping the cation distribution of nanometric ferrites at its equilibrium ratio. Moreover, the sample can be used in coding, storing, and retrieving of binary bit (“0” and “1”) through magnetic field change.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Three dimensional nonlinear simulations of edge localized modes on the EAST tokamak using BOUT++ code

Experimental measurements of edge localized modes (ELMs) observed on the EAST experiment are compared to linear and nonlinear theoretical simulations of peeling-ballooning modes using the BOUT++ code. Simulations predict that the dominant toroidal mode number of the ELM instability becomes larger for lower current, which is consistent with the mode structure captured with visible light using an optical CCD camera. The poloidal mode number of the simulated pressure perturbation shows good agreement with the filamentary structure observed by the camera. The nonlinear simulation is also consistent with the experimentally measured energy loss during an ELM crash and with the radial speed of ELM effluxes measured using a gas puffing imaging diagnostic.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Development of a Monte Carlo code for the data analysis of the {sup 18}F(p,α){sup 15}O reaction at astrophysical energies

Novae are astrophysical events (violent explosion) occurring in close binary systems consisting of a white dwarf and a main-sequence star or a star in a more advanced stage of evolution. They are called 'narrow systems' because the two components interact with each other: there is a process of mass exchange with resulting in the transfer of matter from the companion star to the white dwarf, leading to the formation of this last of the so-called accretion disk, rich mainly of hydrogen. Over time, more and more material accumulates until the pressure and the temperature reached are sufficient to trigger nuclear fusion reactions, rapidly converting a large part of the hydrogen into heavier elements. The products of 'hot hydrogen burning' are then placed in the interstellar medium as a result of violent explosions. Studies on the element abundances observed in these events can provide important information about the stages of evolution stellar. During the outbursts of novae some radioactive isotopes are synthesized: in particular, the decay of short-lived nuclei such as {sup 13}N and {sup 18}F with subsequent emission of gamma radiation energy below 511 keV. The gamma rays from products electron-positron annihilation of positrons emitted in the decay of {sup 18}F are the most abundant and the first observable as soon as the atmosphere of the nova starts to become transparent to gamma radiation. Hence the importance of the study of nuclear reactions that lead both to the formation and to the destruction of {sup 18}F. Among these, the {sup 18}F(p,α){sup 15}O reaction is one of the main channels of destruction. This reaction was then studied at energies of astrophysical interest. The experiment done at Riken, Japan, has as its objective the study of the {sup 18}F(p,α){sup 15}O reaction, using a beam of {sup 18}F produced at CRIB, to derive important information about the phenomenon of novae. In this paper we present the experimental technique and the Monte Carlo code developed to be used in the data analysis process.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2015

Magnetic enhancement and coding in mechanosynthesized Ni{sub 0.3}Zn{sub 0.7}Fe{sub 2}O{sub 4} nanoparticles

The nanosized Ni{sub 0.3}Zn{sub 0.7}Fe{sub 2}O{sub 4} (∼ 15 nm) has been synthesized by high energy ball milling of the bulk powder sample. The sample has been characterized by powder x-ray diffraction, field emission and transmission electron microscopic and dc magnetic measurement techniques. The dc magnetic measurement on the sample indicates that the sample exhibit enhancement of magnetization compared to its counterparts synthesized by chemical methods. Moreover, the system stores the memory of either decrease or increase of magnetic field enabling a magnetic coding of “0” and “1” which can be profitably used in magnetic storage and sensing devices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Spontaneous emergence of autocatalytic information-coding polymers

Self-replicating systems based on information-coding polymers are of crucial importance in biology. They also recently emerged as a paradigm in material design on nano- and micro-scales. We present a general theoretical and numerical analysis of the problem of spontaneous emergence of autocatalysis for heteropolymers capable of template-assisted ligation driven by cyclic changes in the environment. Our central result is the existence of the first order transition between the regime dominated by free monomers and that with a self-sustaining population of sufficiently long chains. We provide a simple, mathematically tractable model supported by numerical simulations, which predicts the distribution of chain lengths and the onset of autocatalysis in terms of the overall monomer concentration and two fundamental rate constants. Another key result of our study is the emergence of the kinetically limited optimal overlap length between a template and each of its two substrates. The template-assisted ligation allows for heritable transmission of the information encoded in chain sequences thus opening up the possibility of long-term memory and evolvability in such systems.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2015

Modeling of EAST ICRF antenna performance using the full-wave code TORIC

Access to advanced operating regimes in the EAST tokamak will require a combination of electron-cyclotron resonance heating (ECRH), neutral beam injection (NBI) and ion cyclotron range frequency heating (ICRF), with the addition of lower-hybrid current drive (LHCD) for current profile control. Prior experiments at the EAST tokamak facility have shown relatively weak response of the plasma temperature to application of ICRF heating, with typical coupled power about 2 MW out of 12 MW source. The launched spectrum, at n{sub φ} = 34 for 0-π -0-π phasing and 27 MHz, is largely inaccessible at line-averaged densities of approximately 2 × 10{sup 19} m{sup −3}. However, with variable antenna phasing and frequency, this system has considerable latitude to explore different heating schemes. To develop an ICRF actuator control model, we have used the full-wave code TORIC to explore the physics of ICRF wave propagation in EAST. The results presented from this study use a spectrum analysis using a superposition of n{sub φ} spanning −50 to +50. The low density regime typical of EAST plasmas results in a perpendicular wavelength comparable to the minor radius which results in global cavity resonance effects and eigenmode formation when the single-pass absorption is low. This behavior indicates that improved performance can be attained by lowering the peak of the k{sub ||} spectrum by using π/3 phasing of the 4-strap antenna. Based on prior studies conducted at Alcator C-Mod, this phasing is also expected to have the advantage of nearly divergence-free box currents, which should result in reduced levels of impurity production. Significant enhancements of the loading resistance may be achieved by using low k{sub ||} phasing and a combination of magnetic field and frequency to vary the location of the resonance and mode conversion regions. TORIC calculations indicate that the significant power may be channeled to the electrons and deuterium majority. We expect that implementation of these recommendations in EAST will yield substantial improvements in the net absorbed power that will greatly assist in the attempt to access advanced tokamak operating regimes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

New developments of the CARTE thermochemical code: A two-phase equation of state for nanocarbons

We developed a new equation of state (EOS) for nanocarbons in the thermodynamic range of high explosives detonation products (up to 50 GPa and 4000 K). This EOS was fitted to an extensive database of thermodynamic properties computed by molecular dynamics simulations of nanodiamonds and nano-onions with the LCBOPII potential. We reproduced the detonation properties of a variety of high explosives with the CARTE thermochemical code, including carbon-poor and carbon-rich explosives, with excellent accuracy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

On a stochastic approach to a code performance estimation

The main goal of an efficient profiling of software is to minimize the runtime overhead under certain constraints and requirements. The traces built by a profiler during the work, affect the performance of the system itself. One of important aspect of an overhead arises from the randomness of variability in the context in which the application is embedded, e.g., due to possible cache misses, etc. Such uncertainty needs to be taken into account in the design phase. In order to overcome these difficulties we propose to investigate this issue through the analysis of the probability distribution of the difference between profiler’s times for the same code. The approximating model is based on the finite normal mixtures within the framework of the method of moving separation of mixtures. We demonstrate some results for the MATLAB profiler using plotting of 3D surfaces by the function surf. The idea can be used for an estimating of a program efficiency.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Alfvén ionization in an MHD-gas interactions code

A numerical model of partially ionized plasmas is developed in order to capture their evolving ionization fractions as a result of Alfvén ionization (AI). The mechanism of, and the parameter regime necessary for, AI is discussed and an expression for the AI rate based on fluid parameters, from a gas-MHD model, is derived. This AI term is added to an existing MHD-gas interactions' code, and the result is a linear, 2D, two-fluid model that includes momentum transfer between charged and neutral species as well as an ionization rate that depends on the velocity fields of both fluids. The dynamics of waves propagating through such a partially ionized plasma are investigated, and it is found that AI has a significant influence on the fluid dynamics as well as both the local and global ionization fraction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Verification of long wavelength electromagnetic modes with a gyrokinetic-fluid hybrid model in the XGC code

As an alternative option to kinetic electrons, the gyrokinetic total-f particle-in-cell (PIC) code XGC1 has been extended to the MHD/fluid type electromagnetic regime by combining gyrokinetic PIC ions with massless drift-fluid electrons. Here, two representative long wavelength modes, shear Alfven waves and resistive tearing modes, are verified in cylindrical and toroidal magnetic field geometries.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Extension of the SIESTA MHD equilibrium code to free-plasma-boundary problems

Here, SIESTA is a recently developed MHD equilibrium code designed to perform fast and accurate calculations of ideal MHD equilibria for three-dimensional magnetic configurations. Since SIESTA does not assume closed magnetic surfaces, the solution can exhibit magnetic islands and stochastic regions. In its original implementation SIESTA addressed only fixed-boundary problems. That is, the shape of the plasma edge, assumed to be a magnetic surface, was kept fixed as the solution iteratively converges to equilibrium. This condition somewhat restricts the possible applications of SIESTA. In this paper we discuss an extension that will enable SIESTA to address free-plasma-boundary problems, opening up the possibility of investigating problems in which the plasma boundary is perturbed either externally or internally. As an illustration, SIESTA is applied to a configuration of the W7-X stellarator.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Validation of extended magnetohydrodynamic simulations of the HIT-SI3 experiment using the NIMROD code

The HIT-SI3 experiment uses a set of inductively driven helicity injectors to apply a non-axisymmetric current drive on the edge of the plasma, driving an axisymmetric spheromak equilibrium in a central confinement volume. These helicity injectors drive a non-axisymmetric perturbation that oscillates in time, with relative temporal phasing of the injectors modifying the mode structure of the applied perturbation. A set of three experimental discharges with different perturbation spectra are modelled here using the NIMROD extended magnetohydrodynamics code, and comparisons are made to both magnetic and fluid measurements. These models successfully capture the bulk dynamics of both the perturbation and the equilibrium, though disagreements related to the pressure gradients experimentally measured exist.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗