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 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Enhanced x-ray emissions from Au-Gd mixture targets ablated by a high-power nanosecond laser

As an important x-ray source, enhancement of x-ray emissions from laser-produced plasmas is imperative for various applications. High-Z Au-Gd mixture targets are proposed to enhance the laser to x-ray conversion efficiency compared to pure Au target. In the experiments, a 1 ns frequency-tripled (351 nm wavelength) laser light was used to obtain an intensity of 3×10{sup 14} W/cm{sup 2} on the targets. The x-ray spectra, total absolute x-ray emissions of all space, M-band fraction and backscattering from pure Au and Au-Gd mixture have been measured, respectively. It is shown that the absolute laser to x-ray conversion efficiency for the Au-Gd mixture containing 60% gold by atom is 47.7%, which has a 15% enhancement compared with that of the pure Au target. The experimental results are consistent with the radiation hydrodynamic simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Measurement of the energy, multiplicity and angular correlation of γ-rays from the thermal neutron capture reaction Gd(n, γ) using JPARC-ANNRI

We conducted an experiment using the JPARC-ANNRI spectrometer to measure the energy, multiplicity and correlation of γ-rays from the neutron capture of natural gadolinium. We incorporated the GEANT4 Monte Carlo (MC) simulation into the detector, and compared the data with the results of the MC simulation. We report our data analysis and compare our data with those obtained by the MC simulation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2014

Nonquasineutral electron vortices in nonuniform plasmas

Electron vortices are observed in the numerical simulation of current carrying plasmas on fast time scales where the ion motion can be ignored. In plasmas with nonuniform density n, vortices drift in the B × ∇n direction with a speed that is on the order of the Hall speed. This provides a mechanism for magnetic field penetration into a plasma. Here, we consider strong vortices with rotation speeds V{sub ϕ} close to the speed of light c where the vortex size δ is on the order of the magnetic Debye length λ{sub B}=|B|/4πen and the vortex is thus nonquasineutral. Drifting vortices are typically studied using the electron magnetohydrodynamic model (EMHD), which ignores the displacement current and assumes quasineutrality. However, these assumptions are not strictly valid for drifting vortices when δ ≈ λ{sub B}. In this paper, 2D electron vortices in nonuniform plasmas are studied for the first time using a fully electromagnetic, collisionless fluid code. Relatively large amplitude oscillations with periods that correspond to high frequency extraordinary modes are observed in the average drift speed. The drift speed W is calculated by averaging the electron velocity field over the vorticity. Interestingly, the time-averaged W from these simulations matches very well with W from the much simpler EMHD simulations even for strong vortices with order unity charge density separation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Illusion optics: Optically transforming the nature and the location of electromagnetic emissions

Complex electromagnetic structures can be designed by using the powerful concept of transformation electromagnetics. In this study, we define a spatial coordinate transformation that shows the possibility of designing a device capable of producing an illusion on an antenna radiation pattern. Indeed, by compressing the space containing a radiating element, we show that it is able to change the radiation pattern and to make the radiation location appear outside the latter space. Both continuous and discretized models with calculated electromagnetic parameter values are presented. A reduction of the electromagnetic material parameters is also proposed for a possible physical fabrication of the device with achievable values of permittivity and permeability that can be obtained from existing well-known metamaterials. Following that, the design of the proposed antenna using a layered metamaterial is presented. Full wave numerical simulations using Finite Element Method are performed to demonstrate the performances of such a device.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Simulation investigation of multipactor in metal components for space application with an improved secondary emission model

Effects of the secondary electron emission (SEE) phenomenon of metal surface on the multipactor analysis of microwave components are investigated numerically and experimentally in this paper. Both the secondary electron yield (SEY) and the emitted energy spectrum measurements are performed on silver plated samples for accurate description of the SEE phenomenon. A phenomenological probabilistic model based on SEE physics is utilized and fitted accurately to the measured SEY and emitted energy spectrum of the conditioned surface material of microwave components. Specially, the phenomenological probabilistic model is extended to the low primary energy end lower than 20 eV mathematically, since no accurate measurement data can be obtained. Embedding the phenomenological probabilistic model into the Electromagnetic Particle-In-Cell (EM-PIC) method, the electronic resonant multipacting in microwave components can be tracked and hence the multipactor threshold can be predicted. The threshold prediction error of the transformer and the coaxial filter is 0.12 dB and 1.5 dB, respectively. Simulation results demonstrate that the discharge threshold is strongly dependent on the SEYs and its energy spectrum in the low energy end (lower than 50 eV). Multipacting simulation results agree quite well with experiments in practical components, while the phenomenological probabilistic model fit both the SEY and the emission energy spectrum better than the traditionally used model and distribution. The EM-PIC simulation method with the phenomenological probabilistic model for the surface collision simulation has been demonstrated for predicting the multipactor threshold in metal components for space application.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Threshold photoelectron spectroscopy of unstable N-containing compounds: Resolution of ΔK subbands in HNCO{sup +} and vibrational resolution in NCO{sup +}

The threshold photoelectron spectra (TPES) of two unstable nitrogen-containing species, HNCO and NCO, were recorded utilizing vacuum ultraviolet synchrotron radiation. Both are intermediates in combustion processes and play a role in the removal of nitrogen oxides from exhaust gases. The rovibronic structure of the first band in the TPES of HNCO{sup +} was analyzed within the framework of an orbital ionization model, and the resolved structure of the origin band was assigned to ΔK subbands. An ionization energy of 11.602 ± 0.005 eV was determined and the vibrational structure of the cationic ground state was analyzed by a Franck-Condon fit. Low lying electronically excited states of HNCO{sup +} were also observed. In a second series of experiments, the NCO radical was generated by flash pyrolysis from chlorine isocyanate. The ionization energy to the X{sup +} {sup 3}Σ{sup −} ground state was determined to be 11.76 ± 0.02 eV, while for the a{sup +1}Δ state, a value of 12.93 ± 0.02 eV was obtained. Vibrational structure was observed for both states, and bands were assigned by Franck-Condon simulations.

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

Silicon superconducting quantum interference device

We have studied a Superconducting Quantum Interference Device (SQUID) made from a single layer thin film of superconducting silicon. The superconducting layer is obtained by heavily doping a silicon wafer with boron atoms using the gas immersion laser doping technique. The SQUID is composed of two nano-bridges (Dayem bridges) in a loop and shows magnetic flux modulation at low temperature and low magnetic field. The overall behavior shows very good agreement with numerical simulations based on the Ginzburg-Landau equations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Hybrid stochastic simulation of reaction-diffusion systems with slow and fast dynamics

In this paper, we present a novel hybrid method to simulate discrete stochastic reaction-diffusion models arising in biochemical signaling pathways. We study moderately stiff systems, for which we can partition each reaction or diffusion channel into either a slow or fast subset, based on its propensity. Numerical approaches missing this distinction are often limited with respect to computational run time or approximation quality. We design an approximate scheme that remedies these pitfalls by using a new blending strategy of the well-established inhomogeneous stochastic simulation algorithm and the tau-leaping simulation method. The advantages of our hybrid simulation algorithm are demonstrated on three benchmarking systems, with special focus on approximation accuracy and efficiency.

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

Numerical simulation of transient moisture transfer into an electronic enclosure

Electronic systems are sometimes exposed to harsh environmental conditions of temperature and humidity. Moisture transfer into electronic enclosures and condensation can cause several problems such as corrosion and alteration in thermal stresses. It is therefore essential to study the local climate inside the enclosures to be able to protect the electronic systems. In this work, moisture transfer into a typical electronic enclosure is numerically studied using CFD. In order to reduce the CPU-time and make a way for subsequent factorial design analysis, a simplifying modification is applied in which the real 3D geometry is approximated by a 2D axial symmetry one. The results for 2D and 3D models were compared in order to calibrate the 2D representation. Furthermore, simulation results were compared with experimental data and good agreement was found.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Single laser pulse compression via strongly coupled stimulated Brillouin scattering in plasma

Laser amplification in plasma, including stimulated Raman scattering amplification and strongly coupled stimulated Brillouin scattering (sc-SBS) amplification, is very promising to generate ultrahigh-power and ultrashort laser pulses. But both are quite complex in experiments: at least three different laser pulses must be prepared; temporal delay and spatial overlap of these three pulses are difficult. We propose a single pulse compression scheme based on sc-SBS in plasma. Only one moderately long laser is applied, the front part of which ionizes the gas to produced plasma, and gets reflected by a plasma mirror at the end of the gas channel. The reflected front quickly depletes the remaining part of the laser by sc-SBS in the self-similar regime. The output laser is much stronger and shorter. This scheme is at first considered theoretically, then validated by using 1D PIC simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

ELECTRON HEATING IN A RELATIVISTIC, WEIBEL-UNSTABLE PLASMA

The dynamics of two initially unmagnetized relativistic counter-streaming homogeneous ion–electron plasma beams are simulated in two dimensions (2D) using the particle-in-cell (PIC) method. It is shown that current filaments, which form due to the Weibel instability, develop a large-scale longitudinal electric field in the direction opposite to the current carried by the filaments as predicted by theory. This field, which is partially inductive and partially electrostatic, is identified as the main source of net electron acceleration, greatly exceeding that due to magnetic field decay at later stages. The transverse electric field, although larger than the longitudinal field, is shown to play a smaller role in heating electrons, contrary to previous claims. It is found that in one dimension, the electrons become strongly magnetized and are not accelerated beyond their initial kinetic energy. Rather, the heating of the electrons is enhanced by the bending and break up of the filaments, which releases electrons that would otherwise be trapped within a single filament and slow the development of the Weibel instability (i.e., the magnetic field growth) via induction as per Lenz’s law. In 2D simulations, electrons are heated to about one quarter of the initial kinetic energy of ions. The magnetic energy at maximum is about 4%, decaying to less than 1% by the end of the simulation. The ions are found to gradually decelerate until the end of the simulation, by which time they retain a residual anisotropy of less than 10%.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Simulating the galaxy cluster “El Gordo” and identifying the merger configuration

The observational features of the massive galaxy cluster “El Gordo” (ACT-CL J0102–4915), such as the X-ray emission, the Sunyaev–Zel’dovich (SZ) effect, and the surface mass density distribution, indicate that they are caused by an exceptional ongoing high-speed collision of two galaxy clusters, similar to the well-known Bullet Cluster. We perform a series of hydrodynamical simulations to investigate the merging scenario and identify the initial conditions for the collision in ACT-CL J0102–4915. By surveying the parameter space of the various physical quantities that describe the two colliding clusters, including their total mass (M), mass ratio (ξ), gas fractions (f{sub b}), initial relative velocity (V), and impact parameter (P), we find an off-axis merger with P∼800 h{sub 70}{sup −1} kpc, V∼2500 km s{sup −1}, M∼3×10{sup 15} M{sub ⊙}, and ξ=3.6 that can lead to most of the main observational features of ACT-CL J0102–4915. Those features include the morphology of the X-ray emission with a remarkable wake-like substructure trailing after the secondary cluster, the X-ray luminosity and the temperature distributions, and also the SZ temperature decrement. The initial relative velocity required for the merger is extremely high and rare compared to that inferred from currently available Λ cold dark matter (ΛCDM) cosmological simulations, which raises a potential challenge to the ΛCDM model, in addition to the case of the Bullet Cluster.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

DETECTION OF A MAGNETIZED DISK AROUND A VERY YOUNG PROTOSTAR

We present subarcsecond resolution polarimetric observations of the 878 μm thermal dust continuum emission obtained with the Submillimeter Array toward the IRAS 16293–2422 protostellar binary system. We report the detection of linearly polarized dust emission arising from the circumstellar disk associated with the IRAS 16293–2422 B protostar. The fractional polarization of ≅ 1.4% is only slightly lower than that expected from theoretical calculations in such disks. The magnetic field structure on the plane of the sky derived from the dust polarization suggests a complex magnetic field geometry in the disk, possibly associated with a rotating disk that is wrapping the field lines as expected from the simulations. The polarization around IRAS 16293–2422 A at subarcsecond angular resolution is only marginally detected.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

SU-F-T-156: Monte Carlo Simulation Using TOPAS for Synchrotron Based Proton Discrete Spot Scanning System

Purpose: This study provides an overview of the design and commissioning of the Monte Carlo (MC) model of the spot-scanning proton therapy nozzle and its implementation for the patient plan simulation. Methods: The Hitachi PROBEAT V scanning nozzle was simulated based on vendor specifications using the TOPAS extension of Geant4 code. FLUKA MC simulation was also utilized to provide supporting data for the main simulation. Validation of the MC model was performed using vendor provided data and measurements collected during acceptance/commissioning of the proton therapy machine. Actual patient plans using CT based treatment geometry were simulated and compared to the dose distributions produced by the treatment planning system (Varian Eclipse 13.6), and patient quality assurance measurements. In-house MATLAB scripts are used for converting DICOM data into TOPAS input files. Results: Comparison analysis of integrated depth doses (IDDs), therapeutic ranges (R90), and spot shape/sizes at different distances from the isocenter, indicate good agreement between MC and measurements. R90 agreement is within 0.15 mm across all energy tunes. IDDs and spot shapes/sizes differences are within statistical error of simulation (less than 1.5%). The MC simulated data, validated with physical measurements, were used for the commissioning of the treatment planning system. Patient geometry simulations were conducted based on the Eclipse produced DICOM plans. Conclusion: The treatment nozzle and standard option beam model were implemented in the TOPAS framework to simulate a highly conformal discrete spot-scanning proton beam system.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

Study of a macrodefect in a silicon carbid single crystal by means of X-ray phase contrast

The morphology of a macrodefect in a single-crystal silicon carbide wafer has been investigated by the computer simulation of an experimental X-ray phase-contrast image. A micropipe, i.e., a long cavity with a small (elliptical in the general case) cross section, in a single crystal has been considered as a macrodefect. A far-field image of micropipe has been measured with the aid of synchrotron radiation without a monochromator. The parameters of micropipe elliptical cross section are determined based on one projection in two directions: parallel and perpendicular to the X-ray beam propagation direction, when scanning along the pipe axis. The results demonstrate the efficiency of the phase contrast method supplemented with computer simulation for studying such macrodefects when the defect position in the sample volume is unknown beforehand.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Studying Nanosecond Discharge in Argon at Atmospheric Pressure with Preionization

The influence of initial conditions on specific features of the formation and development of a cathode-directed ionization wave between two flat electrodes in argon at atmospheric pressure was studied at nanosecond time resolution using a high-speed photoelectron monitor and numerical simulations based on the two-dimensional axisymmetric diffusion–drift model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

CEFR simulation using diffusion code system RAST-F

In this study, the RAST-F code system, which is based on the two-step approach that couples a multi-group cross-section generation Monte-Carlo (MC) code and a multi-group nodal diffusion solver, was used for the neutronic simulation of the CEFR start-up experiments. The numerical solution of the RAST-F system was verified against the full core MC solution MCS at all control rods fully inserted and withdrawn states. The RAST-F solution of the selected experimental simulations was compared against the measurement data. A good agreement between RAST-F and MCS solutions was observed with less than 120 pcm discrepancies and 1.2% root-mean-square error in terms of k{sub eff} and power distribution, respectively. Meanwhile, the RAST-F result agreed well with the experimental data within two-sigma of experimental uncertainty. The good agreement of these results indicates that RAST-F can be used to neutronic steady-state simulations for small core-size SFR, which was challenged to deterministic code systems. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2015

Supplying materials needed for grain growth characterizations of nano-grained UO 2

This activity is supported by the US Nuclear Energy Advanced Modeling and Simulation (NEAMS) Fuels Product Line (FPL) and aims at providing experimental data for the validation of the mesoscale simulation code MARMOT. MARMOT is a mesoscale multiphysics code that predicts the coevolution of microstructure and properties within reactor fuel during its lifetime in the reactor. It is an important component of the Moose-Bison-Marmot (MBM) code suite that has been developed by Idaho National Laboratory (INL) to enable next generation fuel performance modeling capability as part of the NEAMS Program FPL. In order to ensure the accuracy of the microstructure based materials models being developed within the MARMOT code, extensive validation efforts must be carried out. In this report, we summarize our preliminary synchrotron radiation experiments at APS to determine the grain size of nanograin UO 2 . The methodology and experimental setup developed in this experiment can directly apply to the proposed in-situ grain growth measurements. The investigation of the grain growth kinetics was conducted based on isothermal annealing and grain growth characterization as functions of duration and temperature. The kinetic parameters such as activation energy for grain growth for UO 2 with different stoichiometry are obtained and compared with molecular dynamics (MD) simulations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗