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,549 records · Page 86

Investigations and advanced concepts on gyrotron interaction modeling and simulations

In gyrotron theory, the interaction between the electron beam and the high frequency electromagnetic field is commonly modeled using the slow variables approach. The slow variables are quantities that vary slowly in time in comparison to the electron cyclotron frequency. They represent the electron momentum and the high frequency field of the resonant TE modes in the gyrotron cavity. For their definition, some reference frequencies need to be introduced. These include the so-called averaging frequency, used to define the slow variable corresponding to the electron momentum, and the carrier frequencies, used to define the slow variables corresponding to the field envelopes of the modes. From the mathematical point of view, the choice of the reference frequencies is, to some extent, arbitrary. However, from the numerical point of view, there are arguments that point toward specific choices, in the sense that these choices are advantageous in terms of simulation speed and accuracy. In this paper, the typical monochromatic gyrotron operation is considered, and the numerical integration of the interaction equations is performed by the trajectory approach, since it is the fastest, and therefore it is the one that is most commonly used. The influence of the choice of the reference frequencies on the interaction simulations is studied using theoretical arguments, as well as numerical simulations. From these investigations, appropriate choices for the values of the reference frequencies are identified. In addition, novel, advanced concepts for the definitions of these frequencies are addressed, and their benefits are demonstrated numerically.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Design of high-gain, wideband antenna using microwave hyperbolic metasurface

In this work, we apply hyperbolic metasurfaces (HMSs) to design high-gain and wideband antennas. It is shown that HMSs formed by a single layer of split-ring resonators (SRRs) can be excited to generate highly directive beams. In particular, we suggest two types of the SRR-HMS: a capacitively loaded SRR (CLSRR)-HMS and a substrate-backed double SRR (DSRR)-HMS. Both configurations ensure that the periodicity of the structures is sufficiently small for satisfying the effective medium theory. For the antenna design, we propose a two-layer-stacked configuration for the 2.4 GHz frequency band based on the DSRR-HMS excited by a folded monopole. Measurement results confirm numerical simulations and demonstrate that an antenna gain of more than 5 dBi can be obtained for the frequency range of 2.1 - 2.6 GHz, with a maximum gain of 7.8 dBi at 2.4 GHz.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Hydraulic jumps in inhomogeneous strongly coupled toroidal dust flows

The inhomogeneous flow of strongly coupled dust particles in a toroidal particle trap with harmonic radial confinement is analyzed in the incompressible fluid limit. It is shown that the flow can spontaneously generate shock-like events, which are similar to the hydraulic jump in open channel flows. A definition of the Froude number for this model is given and the critical speed is recovered as the group velocity of surface waves. This hydraulic model is compared with molecular-dynamics simulations, which show that a sudden bifurcation of the flow lines and a localized temperature peak appear just at the point where the critical condition for the hydraulic jump is located.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Design and implementation of grid multi-scroll fractional-order chaotic attractors

This paper proposes a novel approach for generating multi-scroll chaotic attractors in multi-directions for fractional-order (FO) systems. The stair nonlinear function series and the saturated nonlinear function are combined to extend equilibrium points with index 2 in a new FO linear system. With the help of stability theory of FO systems, stability of its equilibrium points is analyzed, and the chaotic behaviors are validated through phase portraits, Lyapunov exponents, and Poincaré section. Choosing the order 0.96 as an example, a circuit for generating 2-D grid multiscroll chaotic attractors is designed, and 2-D 9 × 9 grid FO attractors are observed at most. Numerical simulations and circuit experimental results show that the method is feasible and the designed circuit is correct.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Rate theory of solvent exchange and kinetics of Li{sup +} − BF{sub 4}{sup −}/PF{sub 6}{sup −} ion pairs in acetonitrile

In this paper, we describe our efforts to apply rate theories in studies of solvent exchange around Li{sup +} and the kinetics of ion pairings in lithium-ion batteries (LIBs). We report one of the first computer simulations of the exchange dynamics around solvated Li{sup +} in acetonitrile (ACN), which is a common solvent used in LIBs. We also provide details of the ion-pairing kinetics of Li{sup +}-[BF{sub 4}] and Li{sup +}-[PF{sub 6}] in ACN. Using our polarizable force-field models and employing classical rate theories of chemical reactions, we examine the ACN exchange process between the first and second solvation shells around Li{sup +}. We calculate exchange rates using transition state theory and weighted them with the transmission coefficients determined by the reactive flux, Impey, Madden, and McDonald approaches, and Grote-Hynes theory. We found the relaxation times changed from 180 ps to 4600 ps and from 30 ps to 280 ps for Li{sup +}-[BF{sub 4}] and Li{sup +}-[PF{sub 6}] ion pairs, respectively. These results confirm that the solvent response to the kinetics of ion pairing is significant. Our results also show that, in addition to affecting the free energy of solvation into ACN, the anion type also should significantly influence the kinetics of ion pairing. These results will increase our understanding of the thermodynamic and kinetic properties of LIB systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Monte Carlo simulation of a biological object with optical coherent tomography structural images using a voxel-based geometry of a medium

We describe a Monte Carlo algorithm for simulating an interference signal and constructing a structural image of a biological object using optical coherent tomography (OCT). The algorithm is specific in that the object geometry is described by using a three-dimensional array of voxels combined into parallelograms, and the border cross checking is carried out by using an improved Smith algorithm. The geometry of the simulated object is reconstructed based on the structure of real biological tissues obtained by OCT. The algorithm efficiency is estimated by comparing experimental and corresponding simulated images of a subcutaneous blood vessel and upper layer of human skin in vivo . (biophotonics)

60 APPLIED LIFE SCIENCES↗

The bispectrum in the Effective Field Theory of Large Scale Structure

We study the bispectrum in the Effective Field Theory of Large Scale Structure, consistently accounting for the effects of short-scale dynamics. We begin by proving that, as long as the theory is perturbative, it can be formulated to arbitrary order using only operators that are local in time. We then derive all the new operators required to cancel the UV-divergences and obtain a physically meaningful prediction for the one-loop bispectrum. In addition to new, subleading stochastic noises and the viscosity term needed for the one-loop power spectrum, we find three new effective operators. The three new parameters can be constrained by comparing with N-body simulations. The best fit is precisely what is suggested by the structure of UV-divergences, hence justifying a formula for the EFTofLSS bispectrum whose only fitting parameter is already fixed by the power spectrum. This result predicts the bispectrum of N-body simulations up to k{sub max}≈0.22 h Mpc{sup −1} at 0z=, an improvement by nearly a factor of two as compared to one-loop standard perturbation theory.

79 ASTRONOMY AND ASTROPHYSICS↗

MULTIDIMENSIONAL MODELING OF CORONAL RAIN DYNAMICS

We present the first multidimensional, magnetohydrodynamic simulations that capture the initial formation and long-term sustainment of the enigmatic coronal rain phenomenon. We demonstrate how thermal instability can induce a spectacular display of in situ forming blob-like condensations which then start their intimate ballet on top of initially linear force-free arcades. Our magnetic arcades host a chromospheric, transition region, and coronal plasma. Following coronal rain dynamics for over 80 minutes of physical time, we collect enough statistics to quantify blob widths, lengths, velocity distributions, and other characteristics which directly match modern observational knowledge. Our virtual coronal rain displays the deformation of blobs into V-shaped features, interactions of blobs due to mostly pressure-mediated levitations, and gives the first views of blobs that evaporate in situ or are siphoned over the apex of the background arcade. Our simulations pave the way for systematic surveys of coronal rain showers in true multidimensional settings to connect parameterized heating prescriptions with rain statistics, ultimately allowing us to quantify the coronal heating input.

79 ASTRONOMY AND ASTROPHYSICS↗

THE MYSTERY OF THE σ-BUMP—A NEW SIGNATURE FOR MAJOR MERGERS IN EARLY-TYPE GALAXIES?

The stellar velocity dispersion as a function of the galactocentric radius of an early-type galaxy can generally be well approximated by a power law σ∝r {sup β}. However, some observed dispersion profiles show a deviation from this fit at intermediate radii, usually between one and three R {sub eff}, where the velocity dispersion remains constant with radius, showing a bump-like behavior, which we term the {sup σ-}bump. To understand the origin of this σ-bump, we study a set of simulated early-type galaxies formed in major mergers. We find the σ-bump in all of our simulated early-type galaxies, with the size and position of the bump slightly varying from galaxy to galaxy, suggesting that the bump is a characteristic of the major merger formation scenario. The feature can be seen both in the intrinsic and projected stellar velocity dispersions. In contrast to shells that form during the merger event but evolve with time and finally disappear, the σ-bump stays nearly constant with radius and is a permanent feature that is preserved until the end of the simulation. The σ-bump is not seen in the dark matter and gas components and we therefore conclude that it is a purely stellar feature of merger remnants.

79 ASTRONOMY AND ASTROPHYSICS↗

MAGNETOROTATIONAL CORE-COLLAPSE SUPERNOVAE IN THREE DIMENSIONS

We present results of new three-dimensional (3D) general-relativistic magnetohydrodynamic simulations of rapidly rotating strongly magnetized core collapse. These simulations are the first of their kind and include a microphysical finite-temperature equation of state and a leakage scheme that captures the overall energetics and lepton number exchange due to postbounce neutrino emission. Our results show that the 3D dynamics of magnetorotational core-collapse supernovae are fundamentally different from what was anticipated on the basis of previous simulations in axisymmetry (2D). A strong bipolar jet that develops in a simulation constrained to 2D is crippled by a spiral instability and fizzles in full 3D. While multiple (magneto-)hydrodynamic instabilities may be present, our analysis suggests that the jet is disrupted by an m = 1 kink instability of the ultra-strong toroidal field near the rotation axis. Instead of an axially symmetric jet, a completely new, previously unreported flow structure develops. Highly magnetized spiral plasma funnels expelled from the core push out the shock in polar regions, creating wide secularly expanding lobes. We observe no runaway explosion by the end of the full 3D simulation 185 ms after bounce. At this time, the lobes have reached maximum radii of ∼900 km.

79 ASTRONOMY AND ASTROPHYSICS↗

CORONAL HEATING BY THE INTERACTION BETWEEN EMERGING ACTIVE REGIONS AND THE QUIET SUN OBSERVED BY THE SOLAR DYNAMICS OBSERVATORY

The question of what heats the solar corona remains one of the most important puzzles in solar physics and astrophysics. Here we report Solar Dynamics Observatory Atmospheric Imaging Assembly observations of coronal heating by the interaction between emerging active regions (EARs) and the surrounding quiet Sun (QS). The EARs continuously interact with the surrounding QS, resulting in dark ribbons which appear at the boundary of the EARs and the QS. The dark ribbons visible in extreme-ultraviolet wavelengths propagate away from the EARs with speeds of a few km s{sup −1}. The regions swept by the dark ribbons are brightening afterward, with the mean temperature increasing by one quarter. The observational findings demonstrate that uninterrupted magnetic reconnection between EARs and the QS occurs. When the EARs develop, the reconnection continues. The dark ribbons may be the track of the interface between the reconnected magnetic fields and the undisturbed QS’s fields. The propagating speed of the dark ribbons reflects the reconnection rate and is consistent with our numerical simulation. A long-term coronal heating which occurs in turn from nearby the EARs to far away from the EARs is proposed.

79 ASTRONOMY AND ASTROPHYSICS↗

A multipurpose fast neutron beam capability at the MASURCA facility

In this paper we investigate the possible future use of the CEA Cadarache MASURCA experimental fast reactor to generate a fairly high-intensity continuous beam of fast neutrons, having energies distributed in the 1 KeV to 5 MeV range. Such an extracted beam of fast neutrons, tailorable in intensity, size and energy, would be rather unique; it would be of interest to neutron-based research and could open a range of new applications at MASURCA. We report the results of numerical simulations which have been performed to evaluate the feasibility of such a beam port and to characterize it spectrally. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

TH-C-12A-08: New Compact 10 MV S-Band Linear Accelerator: 3D Finite-Element Design and Monte Carlo Dose Simulations

Purpose: To design a new compact S-band linac waveguide capable of producing a 10 MV x-ray beam, while maintaining the length (27.5 cm) of current 6 MV waveguides. This will allow higher x-ray energies to be used in our linac-MRI systems with the same footprint. Methods: Finite element software COMSOL Multiphysics was used to design an accelerator cavity matching one published in an experiment breakdown study, to ensure that our modeled cavities do not exceed the threshold electric fields published. This cavity was used as the basis for designing an accelerator waveguide, where each cavity of the full waveguide was tuned to resonate at 2.997 GHz by adjusting the cavity diameter. The RF field solution within the waveguide was calculated, and together with an electron-gun phase space generated using Opera3D/SCALA, were input into electron tracking software PARMELA to compute the electron phase space striking the x-ray target. This target phase space was then used in BEAM Monte Carlo simulations to generate percent depth doses curves for this new linac, which were then used to re-optimize the waveguide geometry. Results: The shunt impedance, Q-factor, and peak-to-mean electric field ratio were matched to those published for the breakdown study to within 0.1% error. After tuning the full waveguide, the peak surface fields are calculated to be 207 MV/m, 13% below the breakdown threshold, and a d-max depth of 2.42 cm, a D10/20 value of 1.59, compared to 2.45 cm and 1.59, respectively, for the simulated Varian 10 MV linac and brehmsstrahlung production efficiency 20% lower than a simulated Varian 10 MV linac. Conclusion: This work demonstrates the design of a functional 27.5 cm waveguide producing 10 MV photons with characteristics similar to a Varian 10 MV linac.

07 ISOTOPE AND RADIATION SOURCES↗

Sci—Fri PM: Dosimetry—02: A Nested Neutron Spectrometer to Measure Neutron Spectra in Radiotherapy

During high-energy radiotherapy treatments, neutrons are produced in the head of the linac through photonuclear interactions. This has been a concern for many years as photoneutrons contribute to the accepted, yet unwanted, out-of-field doses that pose an iatrogenic risk to patients and an occupational risk to personnel. Presently, in-room neutron measurements are difficult and time-consuming and have traditionally been carried out using Bonner spheres with activation foils and TLDs. In this work, a new detector, the Nested Neutron Spectrometer (NNS) is tested for use in radiotherapy bunkers. The NNS is designed for easy handling and is more practical than the traditional Bonner spheres. The NNS, operated in current mode, was used to measure the dose equivalent, average energy and energy spectrum at several positions in a radiotherapy bunker. The average energy and spectra were compared to Monte Carlo simulations while the dose equivalent was compared to bubble detector measurements. The average energies, as measured by the NNS and Monte Carlo simulations, differed by approximately 30% across the bunker. Measurements of the dose equivalent using the NNS and the bubble detectors agreed within 50% in the maze and less than 10% close to the linac head. Apart from some discrepancies at thermal energies, we also found reasonable agreement between NNS-measured and Monte Carlo-simulated spectra at a number of locations within our radiotherapy bunker. Our results demonstrate that the NNS is a suitable detector to be used in high dose-rate radiotherapy environments.

07 ISOTOPE AND RADIATION SOURCES↗

Ion-Optical System of an Ion Source with Energy Focusing in the Formed Beam

Ways to reduce ion beam intensity losses in a mass spectrometric ion source, which are caused by chromatic aberration of its immersion ion-optical system, are considered. These losses are rather significant in forming a beam from ions with a large energy spread. The reduction of aberration losses is especially important when ion sources are used jointly with mass analyzers with energy focusing. It is shown that these losses can be reduced significantly by using a new type of ion-optical system of the ion source, which includes achromatic elements. A special method for calculating such elements is given. Computer simulations have shown high efficiency of such elements in the ion-optical path of ion sources of mass spectrometers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Application of computational fluid dynamics for the simulation of cryogenic molecular sieve bed absorber of hydrogen isotopes recovery system for Indian LLCB-TBM

One of the most challenging tasks in the design of the fuel cycle system lies in the effective design of Tritium Extraction System (TES) which involves proper extraction and purification of tritium in the fuel cycle of the fusion reactor. Indian Lead Lithium cooled Ceramic Breeder Test Blanket Module (LLCB-TBM) would extract hydrogen isotopes through Cryogenic Molecular Sieve Bed (CMSB) adsorber system. A prototype Hydrogen Isotopes Recovery System (HIRS) is being developed to validate the concepts for tritium extraction by adsorption mass transfer mechanism. In this study, a design model has been developed and analyzed to simulate the adsorption mass transfer kinetics in a fixed bed adsorption column. The simulation leads primarily to effective design of HIRS, which is a state-of-the-art technology. The paper describes the process simulation approach and the results of Computational Fluid Dynamics (CFD) analysis. The effects of different operating conditions are studied to investigate their influence on the hydrogen isotopes adsorption capacity. The results of the present simulation study would be used to understand the best optimized transport phenomenon before realizing the TES as a system for LLCB-TBM. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Neutronics Modeling of the pulsed plasma rocket reactor using rattlesnake

In the pulsed plasma rocket (PPR) reactor, neutron induced fission processes are utilized to implement a series of pulsed micro-explosions of very high power and eject plasma as a propellant. More specifically, projectiles (bullets) composed of moderated uranium are sent through the chamber of an unmoderated uranium barrel. By inducing rapid fission within the bullets, a plasma can be generated with appropriate delivery of neutrons. This study aims to assess the neutronics performance of the PPR reactor during normal operations using the MOOSE-based Rattlesnake code through the evaluation of the impact of the movement of the fuel bullet and the rotation of control drums on the criticality of the system. The Monte Carlo (MC) code Serpent 2 was employed to generate material-based cross sections for use in Rattlesnake and the reference neutronics solution. Cubit was used to generate the mesh for the Rattlesnake model. Parametric studies were conducted to evaluate the best approaches for cross section and mesh generation to ensure accurate results from Rattlesnake. As part of the verification process, the eigenvalue results of the system at various fuel bullet positions were obtained using Rattlesnake and compared with the reference solutions. The acceptable differences show that the Rattlesnake model with appropriate cross section and mesh generation procedure is a sufficiently accurate approximation of the continuous energy (CE) MC model for micro-sized reactors like the PPR reactor. Next, the impact of the rotation of the control drums on the system criticality was evaluated using the verified model. It was shown that the high-fidelity simulation using the deterministic code Rattlesnake can produce sufficiently accurate results for the evaluation of the reactor's neutronics performance in different phases of the normal operation with acceptable computational cost.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Status of the phenomena representation, 3D modeling, and cloud-based software architecture development

Early in 2013, researchers at the Idaho National Laboratory outlined a technical framework to support the implementation of state-of-the-art probabilistic risk assessment to predict the safety performance of advanced small modular reactors. From that vision of the advanced framework for risk analysis, specific tasks have been underway in order to implement the framework. This report discusses the current development of a several tasks related to the framework implementation, including a discussion of a 3D physics engine that represents the motion of objects (including collision and debris modeling), cloud-based analysis tools such as a Bayesian-inference engine, and scenario simulations. These tasks were performed during 2015 as part of the technical work associated with the Advanced Reactor Technologies Program.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗