Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SIMULATOR”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Hybrid simulations of magnetic reconnection with kinetic ions and fluid electron pressure anisotropy

We present the first hybrid simulations with kinetic ions and recently developed equations of state for the electron fluid appropriate for reconnection with a guide field. The equations of state account for the main anisotropy of the electron pressure tensor.Magnetic reconnection is studied in two systems, an initially force-free current sheet and a Harris sheet. The hybrid model with the equations of state is compared to two other models, hybrid simulations with isothermal electrons and fully kinetic simulations. Including the anisotropicequations of state in the hybrid model provides a better match to the fully kinetic model. In agreement with fully kinetic results, the main feature captured is the formation of an electron current sheet that extends several ion inertial lengths. This electron current sheet modifies the Hall magnetic field structure near the X-line, and it is not observed in the standard hybrid model with isotropic electrons. The saturated reconnection rate in this regime nevertheless remains similar in all three models. Here, implications for global modeling are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Gyrokinetic simulation of driftwave instability in field-reversed configuration

Following the recent remarkable progress in magnetohydrodynamic (MHD) stability control in the C-2U advanced beam driven field-reversed configuration (FRC), turbulent transport has become one of the foremost obstacles on the path towards an FRC-based fusion reactor. Significant effort has been made to expand kinetic simulation capabilities in FRC magnetic geometry. The recently upgraded Gyrokinetic Toroidal Code (GTC) now accommodates realistic magnetic geometry from the C-2U experiment at Tri Alpha Energy, Inc. and is optimized to efficiently handle the FRC's magnetic field line orientation. Initial electrostatic GTC simulations find that ion-scale instabilities are linearly stable in the FRC core for realistic pressure gradient drives. Estimated instability thresholds from linear GTC simulations are qualitatively consistent with critical gradients determined from experimental Doppler backscattering fluctuation data, which also find ion scale modes to be depressed in the FRC core. Beyond GTC, A New Code (ANC) has been developed to accurately resolve the magnetic field separatrix and address the interaction between the core and scrape-off layer regions, which ultimately determines global plasma confinement in the FRC. The current status of ANC and future development targets are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Nonlinear electromagnetic formulation for particle-in-cell simulation of lower hybrid waves in toroidal geometry

An electromagnetic particle simulation model has been formulated and verified for nonlinear processes of lower hybrid (LH) waves in fusion plasmas. Electron dynamics are described by the drift kinetic equation using either kinetic momentum or canonical momentum. Ion dynamics are treated as the fluid system or by the Vlasov equation. Compressible magnetic perturbation is retained to simulate both the fast and slow LH waves. Numerical properties are greatly improved by using the electron continuity equation to enforce the consistency between electrostatic potential and vector potential, and by using the importance sampling scheme. The simulation model has been implemented in the gyrokinetic toroidal code (GTC), and verified for the dispersion relation and nonlinear particle trapping of the electromagnetic LH waves.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Speeding up simulations by slowing down particles: Speed-limited particle-in-cell simulation

Based on the particle-in-cell (PIC) plasma simulation method, the speed-limited PIC (SLPIC) method delivers faster kinetic plasma simulation in cases where the particle distributions evolve slowly compared with the maximum stable PIC timestep. SLPIC thus offers more feasible, fully kinetic simulation in regimes that historically have required fluid approaches, such as magnetohydrodynamic (MHD), two-fluid, or Boltzmann electron treatments. In particular, SLPIC allows an explicit time advance with steps much larger than the inverse plasma frequency, avoiding the instability explicit PIC faces with large timesteps. SLPIC avoids this instability by slowing down fast particles (e.g., electrons) in a way that is rigorously underpinned by an approximate Vlasov equation; unlike MHD, two-fluid, and Boltzmann electron approaches, SLPIC does not fundamentally neglect any first-principles plasma physics, although the choices of the grid cell size, timestep, and number of macroparticles per cell naturally limit the physical phenomena that can be accurately represented. SLPIC can be implemented with minor modifications of a standard PIC code and does not require an implicit time advance. It enables large timesteps in first-principles kinetic plasma simulation of appropriately slow phenomena, and it can handle many of the same complications as PIC, such as boundary conditions and collisions. Finally, in an argon plasma sheath test problem, a SLPIC simulation achieved a speed-up of a factor of 160 over the corresponding PIC simulation, without loss of accuracy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Assessment and Commissioning of Electrical Substation Grid Testbed with a Real-Time Simulator and Protective Relays/Power Meters in the Loop

Electrical utility substations are wired with intelligent electronic devices (IEDs), such as protective relays, power meters, and communication switches. Substation engineers commission these IEDs to assess the appropriate measurements for monitoring, control, power system protection, and communication applications. Like real electrical utility substations, complex electrical substation grid testbeds (ESGTs) need to be assessed for measuring current and voltage signals in monitoring, power system protection, control (synchro check), and communication applications that are limited by small measurement percentage errors. In the process of setting an ESGT with real-time simulators and IEDs in the loop, protective relays, power meters, and communication devices must be commissioned before running experiments. In this study, an ESGT with IEDs and distributed ledger technology was developed. The ESGT with a real-time simulator and IEDs in the loop was satisfactorily assessed and commissioned. The commissioning and problem-solving tasks of the testbed are described to define a method with flowcharts to assess possible trouble-shooting in ESGTs. This method was based on comparing the simulations versus IED measurements for the phase current and voltage magnitudes, three-phase phasor diagrams, breaker states, protective relay times with selectivity coordination at electrical faults, communication data points, and time-stamp sources.

Piesciorovsky, Emilio↗
DOE OSTI2016

The development of explosions in axisymmetric ab initio core-collapse supernova simulations of 12–25 M{sub ⊙} STARS

We present four ab initio axisymmetric core-collapse supernova simulations initiated from 12, 15, 20, and 25 M{sub ⊙} zero-age main sequence progenitors. All of the simulations yield explosions and have been evolved for at least 1.2 s after core bounce and 1 s after material first becomes unbound. These simulations were computed with our Chimera code employing RbR spectral neutrino transport, special and general relativistic transport effects, and state-of-the-art neutrino interactions. Continuing the evolution beyond 1 s after core bounce allows the explosions to develop more fully and the processes involved in powering the explosions to become more clearly evident. We compute explosion energy estimates, including the negative gravitational binding energy of the stellar envelope outside the expanding shock, of 0.34, 0.88, 0.38, and 0.70 Bethe (B ≡ 10{sup 51} erg) and increasing at 0.03, 0.15, 0.19, and 0.52 B s{sup −1}, respectively, for the 12, 15, 20, and 25 M{sub ⊙} models at the endpoint of this report. We examine the growth of the explosion energy in our models through detailed analyses of the energy sources and flows. We discuss how the explosion energies may be subject to stochastic variations as exemplfied by the effect of the explosion geometry of the 20 M{sub ⊙} model in reducing its explosion energy. We compute the proto-neutron star masses and kick velocities. We compare our results for the explosion energies and ejected {sup 56}Ni masses against some observational standards despite the large error bars in both models and observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Pair-instability Supernova Simulations: Progenitor Evolution, Explosion, and Light Curves

In recent years, the viability of the pair-instability supernova (PISN) scenario for explaining superluminous supernovae has all but disappeared except for a few slowly-evolving examples. However, PISNe are not predicted to be superluminous throughout the bulk of their mass range. In fact, it is more likely that the first PISN we see (if we have not seen one already) will not be superluminous. Here, we present hydrodynamic simulations of PISNe for four stellar models with unique envelope properties spanning the PISN mass range. In addition, we compute synthetic light curves (LCs) for comparison with current and future observations. We also investigate, in the context of our most massive model, the prospect of mixing in the supernova ejecta, alleviating discrepancies between current PISN models and the remaining superluminous candidate events. To this end, we present the first published 3D hydrodynamic simulations of PISNe. After achieving convergence between 1D, 2D, and 3D simulations, we examine mixing in the supernova ejecta and its affect on the bolometric LC. We observe slight deviations from spherical symmetry, which increase with the number of dimensions. We find no significant effects on the bolometric LC; however, we conclude that mixing between the silicon and oxygen rich layers caused by the Rayleigh–Taylor instability may affect spectra.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2003

HYDRA : High-speed simulation architecture for precision spacecraft formation simulation

e Hierarchical Distributed Reconfigurable Architecture- is a scalable simulation architecture that provides flexibility and ease-of-use which take advantage of modern computation and communication hardware. It also provides the ability to implement distributed - or workstation - based simulations and high-fidelity real-time simulation from a common core. Originally designed to serve as a research platform for examining fundamental challenges in formation flying simulation for future space missions, it is also finding use in other missions and applications, all of which can take advantage of the underlying Object-Oriented structure to easily produce distributed simulations. Hydra automates the process of connecting disparate simulation components (Hydra Clients) through a client server architecture that uses high-level descriptions of data associated with each client to find and forge desirable connections (Hydra Services) at run time. Services communicate through the use of Connectors, which abstract messaging to provide single-interface access to any desired communication protocol, such as from shared-memory message passing to TCP/IP to ACE and COBRA. Hydra shares many features with the HLA, although providing more flexibility in connectivity services and behavior overriding.

formation flying↗
DOE OSTI2015

Probing the mechanism of CO 2 capture in diamine-appended metal–organic frameworks using measured and simulated X-ray spectroscopy

Diamine-appended metal–organic frameworks display great promise for carbon capture applications, due to unusual step-shaped adsorption behavior that was recently attributed to a cooperative mechanism in which the adsorbed CO 2 molecules insert into the metal–nitrogen bonds to form ordered ammonium carbamate chains. We present a detailed study of this mechanism by in situ X-ray absorption spectroscopy and density functional theory calculations. Distinct spectral changes at the N and O K-edges are apparent upon CO 2 adsorption in both mmen-Mg 2 (dobpdc) and mmen-Mn 2 (dobpdc), and these are evaluated based upon computed spectra from three potential adsorption structures. The computations reveal that the observed spectral changes arise from specific electronic states that are signatures of a quasi-trigonal planar carbamate species that is hydrogen bonded to an ammonium cation. This eliminates two of the three structures studied, and confirms the insertion mechanism. Here, we note the particular sensitivity of X-ray absorption spectra to the insertion step of this mechanism, underpinning the strength of the technique for examining subtle chemical changes upon gas adsorption.

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

Monte Carlo Simulation of Electron Thermalization in Scintillator Materials: Implications for Scintillator Nonproportionality

The lack of reliable quantitative estimates of the length and time scales associated with hot electron thermalization after a gamma-ray induced energy cascade obscures the interplay of various microscopic processes controlling scintillator performance and hampers the search for improved detector materials. We apply a detailed microscopic kinetic Monte Carlo model of the creation and subsequent thermalization of hot electrons produced by gamma irradiation of six important scintillating crystals to determine the spatial extent of the cloud of excitations produced by gamma rays and the time required for the cloud to thermalize with the host lattice. The main ingredients of the model are ensembles of microscopic track structures produced upon gamma excitation (including the energy distribution of the excited carriers), numerical estimates of electron-phonon scattering rates, and a calculated particle dispersion to relate the speed and energy of excited carriers. All these ingredients are based on first-principles density functional theory calculations of the electronic and phonon band structures of the materials. Details of the Monte Carlo model are presented along with results for thermalization time and distance distributions. These results are discussed in light of previous work. It is found that among the studied materials, calculated thermalization distances are positively correlated with measured nonproportionality. In the important class of halide scintillators, the particle dispersion is found to be more influential than the largest phonon energy in determining the thermalization distance.

Phonons, radiation detection, scintillators, thema↗
DOE OSTI2014

Combined molecular dynamics-spin dynamics simulations of bcc iron

Using a classical model that treats translational and spin degrees of freedom on an equal footing, we study phonon-magnon interactions in BCC iron with combined molecular and spin dynamics methods. The atomic interactions are modeled via an empirical many-body potential while spin dependent interactions are established through a Hamiltonian of the Heisenberg form with a distance dependent magnetic exchange interaction obtained from first principles electronic structure calculations. The temporal evolution of translational and spin degrees of freedom was determined by numerically solving the coupled equations of motion, using an algorithm based on the second order Suzuki-Trotter decomposition of the exponential operators. By calculating Fourier transforms of space- and time-displaced correlation functions, we demonstrate that the the presence of lattice vibrations leads to noticeable softening and damping of spin wave modes. As a result of the interplay between lattice and spin subsystems, we also observe additional longitudinal spin wave excitations, with frequencies which coincide with that of the longitudinal lattice vibrations.

Perera, Meewanage Dilina↗