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,333 records · Page 74
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Analysis of reliable sub-ns spin-torque switching under transverse bias magnetic fields

The switching process of a magnetic spin-valve nanosystem subject to spin-polarized current pulses is considered. The dependence of the switching probability on the current pulse duration is investigated. The further application of a transverse field along the intermediate anisotropy axis of the particle is used to control the quasi-random relaxation of magnetization to the reversed magnetization state. The critical current amplitudes to realize the switching are determined by studying the phase portrait of the Landau-Lifshtz-Slonczewski dynamics. Macrospin numerical simulations are in good agreement with the theoretical prediction and demonstrate reliable switching even for very short (below 100 ps) current pulses.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Design of the new extraction grid for the NIO1 negative ion source

NIO1 is a compact source of negative ions jointly developed by RFX and INFN, to study the physics of production and acceleration of H{sup −} beams. Negative ions, up to 120 mA of current, are extracted from a radiofrequency driven plasma, by means of a gridded electrode (plasma grid, PG) featuring 9 apertures arranged in a 3x3 square lattice. The same aperture pattern is replicated in the following electrodes, allowing ion acceleration up to 60 keV. All electrodes are realized in copper, by electro-deposition technique, leaving empty slots in the metal to place magnets and to flow water for the grid cooling. The first set of electrodes was completed, installed in the source and tested. At the same time, an upgrade of the extraction system was carried out, in order to optimize the beam optics and to explore alternative electrostatic configurations. In particular, the accelerator will be modified by completely replacing the EG grid, exploiting the modularity of NIO1. The new electrode will feature other slots in between apertures, to place additional magnets. This allows testing different magnetic configurations, to optimize electron filtering and residual ion deflection. The present paper describes the theoretical activities driving the design of these new extractors, carried out with most updated numerical codes, and exploiting the synergy with the refined modeling of the 40 A ITER negative ion sources, under development at Consorzio RFX. Beam simulations are performed both with tracing codes (SLACCAD and OPERA) and with particle in cell codes (ACCPIC)

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Topology of classical molecular optimal control landscapes for multi-target objectives

This paper considers laser-driven optimal control of an ensemble of non-interacting molecules whose dynamics lie in classical phase space. The molecules evolve independently under control to distinct final states. We consider a control landscape defined in terms of multi-target (MT) molecular states and analyze the landscape as a functional of the control field. The topology of the MT control landscape is assessed through its gradient and Hessian with respect to the control. Under particular assumptions, the MT control landscape is found to be free of traps that could hinder reaching the objective. The Hessian associated with an optimal control field is shown to have finite rank, indicating an inherent degree of robustness to control noise. Both the absence of traps and rank of the Hessian are shown to be analogous to the situation of specifying multiple targets for an ensemble of quantum states. Numerical simulations are presented to illustrate the classical landscape principles and further characterize the system behavior as the control field is optimized.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Analysis of 440 GeV proton beam–matter interaction experiments at the High Radiation Materials test facility at CERN

In a previous paper [Schmidt et al., Phys. Plasmas 21, 080701 (2014)], we presented the first results on beam–matter interaction experiments that were carried out at the High Radiation Materials test facility at CERN. In these experiments, extended cylindrical targets of solid copper were irradiated with beam of 440 GeV protons delivered by the Super Proton Synchrotron (SPS). The beam comprised of a large number of high intensity proton bunches, each bunch having a length of 0.5 ns with a 50 ns gap between two neighboring bunches, while the length of this entire bunch train was about 7 μs. These experiments established the existence of the hydrodynamic tunneling phenomenon the first time. Detailed numerical simulations of these experiments were also carried out which were reported in detail in another paper [Tahir et al., Phys. Rev. E 90, 063112 (2014)]. Excellent agreement was found between the experimental measurements and the simulation results that validate our previous simulations done using the Large Hadron Collider (LHC) beam of 7 TeV protons [Tahir et al., Phys. Rev. Spec. Top.--Accel. Beams 15, 051003 (2012)]. According to these simulations, the range of the full LHC proton beam and the hadronic shower can be increased by more than an order of magnitude due to the hydrodynamic tunneling, compared to that of a single proton. This effect is of considerable importance for the design of machine protection system for hadron accelerators such as SPS, LHC, and Future Circular Collider. Recently, using metal cutting technology, the targets used in these experiments have been dissected into finer pieces for visual and microscopic inspection in order to establish the precise penetration depth of the protons and the corresponding hadronic shower. This, we believe will be helpful in studying the very important phenomenon of hydrodynamic tunneling in a more quantitative manner. The details of this experimental work together with a comparison with the numerical simulations are presented in this paper.

36 MATERIALS SCIENCE↗
DOE OSTI2015

CAMELOT: A machine learning approach for coarse-grained simulations of aggregation of block-copolymeric protein sequences

We report the development and deployment of a coarse-graining method that is well suited for computer simulations of aggregation and phase separation of protein sequences with block-copolymeric architectures. Our algorithm, named CAMELOT for Coarse-grained simulations Aided by MachinE Learning Optimization and Training, leverages information from converged all atom simulations that is used to determine a suitable resolution and parameterize the coarse-grained model. To parameterize a system-specific coarse-grained model, we use a combination of Boltzmann inversion, non-linear regression, and a Gaussian process Bayesian optimization approach. The accuracy of the coarse-grained model is demonstrated through direct comparisons to results from all atom simulations. We demonstrate the utility of our coarse-graining approach using the block-copolymeric sequence from the exon 1 encoded sequence of the huntingtin protein. This sequence comprises of 17 residues from the N-terminal end of huntingtin (N17) followed by a polyglutamine (polyQ) tract. Simulations based on the CAMELOT approach are used to show that the adsorption and unfolding of the wild type N17 and its sequence variants on the surface of polyQ tracts engender a patchy colloid like architecture that promotes the formation of linear aggregates. These results provide a plausible explanation for experimental observations, which show that N17 accelerates the formation of linear aggregates in block-copolymeric N17-polyQ sequences. The CAMELOT approach is versatile and is generalizable for simulating the aggregation and phase behavior of a range of block-copolymeric protein sequences.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

A new interpretation for performance improvement of high-efficiency vertical blue light-emitting diodes by InGaN/GaN superlattices

The effect of InGaN/GaN superlattices (SLs) on quantum efficiency and forward voltage of vertical blue InGaN/GaN multiple quantum well (MQW) light-emitting diodes (LED) grown on Si substrate has been experimentally and theoretically investigated. We have prepared two LED samples, in which the 30 and 45 periods of SLs are inserted between MQW active layers and n-GaN layer, respectively. Electroluminescence measurement shows that the LED with 45 periods of SLs has higher quantum efficiency but lower forward voltage. It is observed that V-shaped pits grow up in size with an increase in SLs period number by means of scan transmission electron microscope and secondary ion mass spectrometry. Further numerical simulations confirm that the performance improvement of LED by SLs is mainly ascribed to enhancing hole injection from the V-shaped pits.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Effect of flexibility on the growth of concentration fluctuations in a suspension of sedimenting fibers: Particle simulations

Three-dimensional numerical simulations are performed to study the stability of a sedimenting suspension of weakly flexible fibers. It is well known that a suspension of rigid rods sedimenting under gravity at low Reynolds number is unstable to concentration fluctuations owing to hydrodynamic interactions. Flexible fibers, however, reorient while settling and even weak flexibility can alter their collective dynamics. In our recent work [Manikantan et al., “The instability of a sedimenting suspension of weakly flexible fibres,” J. Fluid Mech. 756, 935–964 (2014)], we developed a mean-field theory to predict the linear stability of such a system. Here, we verify these predictions using accurate and efficient particle simulations based on a slender-body model. We also demonstrate the mechanisms by which flexibility-induced reorientation alters suspension microstructure, and through it, its stability. Specifically, we first show that the anisotropy of the base state in the case of a suspension of flexible fibers has a destabilizing effect compared to a suspension of rigid rods. Second, a conflicting effect of flexibility is also shown to suppress particle clustering and slow down the growth of the instability. The relative magnitude of filament flexibility and rotational Brownian motion dictates which effect dominates, and our simulations qualitatively follow theoretically predicted trends. The mechanism for either effects is tied to the flexibility-induced reorientation of particles, which we illustrate using velocity and orientation statistics from our simulations. Finally, we also show that, in the case of an initially homogeneous and isotropic suspension, flexibility always acts to suppress the growth of the instability.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Theoretical investigation of the dielectric-filled relativistic magnetron

The fundamental mode frequency of a dielectric-filled relativistic magnetron is studied theoretically by the method of the equivalent circuit, and an exact fundamental mode frequency formula is derived. To prove the validity of the theoretical formula, simulation investigation is performed. The simulation results agree well with the theoretical formula, and the relative error does not exceed 3%. The comparative results verify the creditability of the theoretical formula.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Multi-scale modeling of the CD8 immune response

During the primary CD8 T-Cell immune response to an intracellular pathogen, CD8 T-Cells undergo exponential proliferation and continuous differentiation, acquiring cytotoxic capabilities to address the infection and memorize the corresponding antigen. After cleaning the organism, the only CD8 T-Cells left are antigen-specific memory cells whose role is to respond stronger and faster in case they are presented this very same antigen again. That is how vaccines work: a small quantity of a weakened pathogen is introduced in the organism to trigger the primary response, generating corresponding memory cells in the process, giving the organism a way to defend himself in case it encounters the same pathogen again. To investigate this process, we propose a non linear, multi-scale mathematical model of the CD8 T-Cells immune response due to vaccination using a maturity structured partial differential equation. At the intracellular scale, the level of expression of key proteins is modeled by a delay differential equation system, which gives the speeds of maturation for each cell. The population of cells is modeled by a maturity structured equation whose speeds are given by the intracellular model. We focus here on building the model, as well as its asymptotic study. Finally, we display numerical simulations showing the model can reproduce the biological dynamics of the cell population for both the primary response and the secondary responses.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Simulation of He{sup +} induced afterpulses in PMTs

We report on a Monte Carlo simulation study of afterpulses due to trace gases in EMI 9235QA photomultipliers that are widely used in many luminescence detection systems operating in single photon counting mode. The numerical simulation takes into account the nonuniform electric field distribution and processes including elastic scattering: e + He → e + He, excitation: e + He → e + He{sup ∗}, ionization: e + He → 2e + He{sup +}, elastic scattering: He{sup +} + He → He{sup +} + He, charge transfer: He{sup +} + He → He{sub f} + He{sup +} (He{sub f} indicates a fast neutral) and elastic scattering: He{sub f} + He → He{sub f} + He{sub (f)}. The simulated and the measured time of flight distributions agree well. In addition, the above simulated processes demonstrate mechanisms of the observed series of pulses brought about by a single helium atom ionization.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2016

Measurement of effect of electron cyclotron heating in a tandem mirror plasma using a semiconductor detector array and an electrostatic energy analyzer

Temporally and spatially resolved soft x-ray and end-loss-electron analyses of the electron cyclotron heated plasmas are carried out by using a semiconductor detector array and an electrostatic energy analyzer in the GAMMA 10 tandem mirror. The flux and the energy spectrum of the end loss electrons are measured by a multi-grid energy analyzer. Recently, the electron cyclotron heating power modulation experiments have been started in order to generate and control the high heat flux and to make the edge localized mode-like intermittent heat load pattern for the divertor simulation studies by the use of these detectors for electron properties.

43 PARTICLE ACCELERATORS↗
DOE OSTI2013

A MATHEMATICAL MODEL FOR AN HOURGLASS MAGNETIC FIELD

Starting with a mathematical boundary value problem for the magnetic vector potential in an axisymmetric cylindrical coordinate system, we derive a general solution for any arbitrary current distribution using the method of Green's functions. We use this to derive an analytic form for an hourglass magnetic field pattern created by electrical currents that are concentrated near (but not confined within) the equatorial plane of a cylindrical coordinate system. Our solution is not characterized by a cusp at the equatorial plane, as in previous solutions based on a current sheet. The pattern we derive provides a very good fit to hourglass magnetic field patterns emerging from three-dimensional numerical simulations of core formation, and can in principle be used for source-fitting of observed magnetic hourglass patterns.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

JOVIAN STRATOSPHERE AS A CHEMICAL TRANSPORT SYSTEM: BENCHMARK ANALYTICAL SOLUTIONS

We systematically investigated the solvable analytical benchmark cases in both one- and two-dimensional (1D and 2D) chemical-advective-diffusive systems. We use the stratosphere of Jupiter as an example but the results can be applied to other planetary atmospheres and exoplanetary atmospheres. In the 1D system, we show that CH{sub 4} and C{sub 2}H{sub 6} are mainly in diffusive equilibrium, and the C{sub 2}H{sub 2} profile can be approximated by modified Bessel functions. In the 2D system in the meridional plane, analytical solutions for two typical circulation patterns are derived. Simple tracer transport modeling demonstrates that the distribution of a short-lived species (such as C{sub 2}H{sub 2}) is dominated by the local chemical sources and sinks, while that of a long-lived species (such as C{sub 2}H{sub 6}) is significantly influenced by the circulation pattern. We find that an equator-to-pole circulation could qualitatively explain the Cassini observations, but a pure diffusive transport process could not. For slowly rotating planets like the close-in extrasolar planets, the interaction between the advection by the zonal wind and chemistry might cause a phase lag between the final tracer distribution and the original source distribution. The numerical simulation results from the 2D Caltech/JPL chemistry-transport model agree well with the analytical solutions for various cases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

THE SPACE MOTION OF LEO I: THE MASS OF THE MILKY WAY'S DARK MATTER HALO

We combine our Hubble Space Telescope measurement of the proper motion of the Leo I dwarf spheroidal galaxy (presented in a companion paper) with the highest resolution numerical simulations of Galaxy-size dark matter halos in existence to constrain the mass of the Milky Way's dark matter halo (M{sub vir,MW}). Despite Leo I's large Galactocentric space velocity (200 km s{sup -1}) and distance (261 kpc), we show that it is extremely unlikely to be unbound if Galactic satellites are associated with dark matter substructure, as 99.9% of subhalos in the simulations are bound to their host. The observed position and velocity of Leo I strongly disfavor a low-mass Milky Way: if we assume that Leo I is the least bound of the Milky Way's classical satellites, then we find that M{sub vir,MW} > 10{sup 12} M{sub Sun} at 95% confidence for a variety of Bayesian priors on M{sub vir,MW}. In lower mass halos, it is vanishingly rare to find subhalos at 261 kpc moving as fast as Leo I. Should an additional classical satellite be found to be less bound than Leo I, this lower limit on M{sub vir,MW} would increase by 30%. Imposing a mass-weighted {Lambda}CDM prior, we find a median Milky Way virial mass of M{sub vir,MW} = 1.6 Multiplication-Sign 10{sup 12} M{sub Sun }, with a 90% confidence interval of [1.0-2.4] Multiplication-Sign 10{sup 12} M{sub Sun }. We also confirm a strong correlation between subhalo infall time and orbital energy in the simulations and show that proper motions can aid significantly in interpreting the infall times and orbital histories of satellites.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

TOWARD A THEORY OF ASTROPHYSICAL PLASMA TURBULENCE AT SUBPROTON SCALES

We present an analytical study of subproton electromagnetic fluctuations in a collisionless plasma with a plasma beta of the order of unity. In the linear limit, a rigorous derivation from the kinetic equation is conducted focusing on the role and physical properties of kinetic-Alfvén and whistler waves. Then, nonlinear fluid-like equations for kinetic-Alfvén waves and whistler modes are derived, with special emphasis on the similarities and differences in the corresponding plasma dynamics. The kinetic-Alfvén modes exist in the lower-frequency region of phase space, ω << k v{sub Ti} , where they are described by the kinetic-Alfvén system. These modes exist both below and above the ion-cyclotron frequency. The whistler modes, which are qualitatively different from the kinetic-Alfvén modes, occupy a different region of phase space, k v{sub Ti} << ω << k{sub z}v{sub Te} , and they are described by the electron magnetohydrodynamics (MHD) system or the reduced electron MHD system if the propagation is oblique. Here, k{sub z} and k are the wavenumbers along and transverse to the background magnetic field, respectively, and v{sub Ti} and v{sub Te} are the ion and electron thermal velocities, respectively. The models of subproton plasma turbulence are discussed and the results of numerical simulations are presented. We also point out possible implications for solar-wind observations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

MAGNETIC FIELD COMPONENTS ANALYSIS OF THE SCUPOL 850 μm POLARIZATION DATA CATALOG

We present an extensive analysis of the 850 μm polarization maps of the SCUBA Polarimeter Legacy (SCUPOL) Catalogue produced by Matthews et al., focusing exclusively on the molecular clouds and star-forming regions. For the sufficiently sampled regions, we characterize the depolarization properties and the turbulent-to-mean magnetic field ratio of each region. Similar sets of parameters are calculated from two-dimensional synthetic maps of dust-emission polarization produced with three-dimensional magnetohydrodynamics (MHD) numerical simulations scaled to the S106, OMC-2/3, W49, and DR21 molecular cloud polarization maps. For these specific regions, the turbulent MHD regimes retrieved from the simulations, as described by the turbulent Alfvén and Sonic Mach numbers, are consistent within a factor one to two with the values of the same turbulent regimes estimated from the analysis of Zeeman measurements data provided by Crutcher. Constraints on the values of the inclination angle α of the mean magnetic field with respect to the line of sight are also given. The values obtained from the comparison of the simulations with the SCUPOL data are consistent with the estimates made by using two observational methods provided by other authors. Our main conclusion is that simple, ideal, isothermal, and non-self-gravitating MHD simulations are sufficient in order to describe the large-scale observed physical properties of the envelopes of this set of regions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Small-scale gradients of charged particles in the heliospheric magnetic field

Using numerical simulations of charged-particles propagating in the heliospheric magnetic field, we study small-scale gradients, or 'dropouts,' in the intensity of solar energetic particles seen at 1 AU. We use two turbulence models, the foot-point random motion model and the two-component model, to generate fluctuating magnetic fields similar to spacecraft observations at 1 AU. The turbulence models include a Kolmogorov-like magnetic field power spectrum containing a broad range of spatial scales from those that lead to large-scale field-line random walk to small scales leading to resonant pitch-angle scattering of energetic particles. We release energetic protons (20 keV-10 MeV) from a spatially compact and instantaneous source. The trajectories of energetic charged particles in turbulent magnetic fields are numerically integrated. Spacecraft observations are mimicked by collecting particles in small windows when they pass the windows at a distance of 1 AU. We show that small-scale gradients in the intensity of energetic particles and velocity dispersions observed by spacecraft can be reproduced using the foot-point random motion model. However, no dropouts are seen in simulations using the two-component magnetic turbulence model. We also show that particle scattering in the solar wind magnetic field needs to be infrequent for intensity dropouts to form.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

KINEMATICS OF SUPERBUBBLES AND SUPERSHELLS IN THE IRREGULAR GALAXY, NGC 1569

We present observations in the optical lines of Hα and [S II] (λλ6717, 6731 Å) and in X-rays of the irregular galaxy, NGC 1569. The observations in Hα and [S II] were made with the UNAM scanning Fabry-Perot interferometer (PUMA) and the X-ray data were obtained from the Chandra data archive. We detected several superbubbles, filaments, and supershells in NGC 1569 for which we determined size as well as their kinematic properties. We present a catalog of expansion velocities of 12 superbubbles, listing their positions, diameters, and physical parameters. Likewise, we present a catalog of 15 filaments and 4 supershells. In order to identify possible X-ray emission from the superbubbles in this galaxy, we analyzed the X-ray emission of NGC 1569 in two energy bands: 0.2-2.0 keV (soft X-rays) and 2.0-8.0 keV (hard X-rays). Based on X-ray images, we detected X-ray emission that could possibly be related to some of the superbubbles. The spectrum of the X-ray superbubbles can be described by an optically thin thermal plasma model. In order to identify the possible coexistence of galactic super winds and superbubbles we have performed adiabatic three-dimensional N-body/smoothed particle hydrodynamics simulations to follow the evolution of the most important stellar clusters in this galaxy, SSC A and SSC B, using the GADGET-2 code. Those simulations demonstrate that depending on the specific initial conditions, the formation of superbubbles or a galactic superwind can result in NGC 1569.

79 ASTRONOMY AND ASTROPHYSICS↗