Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulators”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

Electromagnetic instability of thin reconnection layers: Comparison of three-dimensional simulations with MRX observations

In this study, the influence of current-aligned instabilities on magnetic reconnection in weakly collisional regimes is investigated using experimental observations from Magnetic Reconnection Experiment (MRX) [M. Yamada et al., Phys. Plasmas 4, 1936 (1997)] and large-scale fully kinetic simulations. In the simulations as well as in the experiment, the dominant instability is localized near the center of the reconnection layer, produces large perturbations of the magnetic field, and is characterized by the wavenumber that is a geometric mean between electron and ion gyroradii k~(ρ e ρ i ) -1/2 . However, both the simulations and the experimental observations suggest the instability is not the dominant reconnection mechanism under parameters typical of MRX.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Parallel continuous simulated tempering and its applications in large-scale molecular simulations

In this paper, we introduce a parallel continuous simulated tempering (PCST) method for enhanced sampling in studying large complex systems. It mainly inherits the continuous simulated tempering (CST) method in our previous studies [C. Zhang and J. Ma, J. Chem. Phys. 130, 194112 (2009); C. Zhang and J. Ma, J. Chem. Phys. 132, 244101 (2010)], while adopts the spirit of parallel tempering (PT), or replica exchange method, by employing multiple copies with different temperature distributions. Differing from conventional PT methods, despite the large stride of total temperature range, the PCST method requires very few copies of simulations, typically 2–3 copies, yet it is still capable of maintaining a high rate of exchange between neighboring copies. Furthermore, in PCST method, the size of the system does not dramatically affect the number of copy needed because the exchange rate is independent of total potential energy, thus providing an enormous advantage over conventional PT methods in studying very large systems. The sampling efficiency of PCST was tested in two-dimensional Ising model, Lennard-Jones liquid and all-atom folding simulation of a small globular protein trp-cage in explicit solvent. The results demonstrate that the PCST method significantly improves sampling efficiency compared with other methods and it is particularly effective in simulating systems with long relaxation time or correlation time. We expect the PCST method to be a good alternative to parallel tempering methods in simulating large systems such as phase transition and dynamics of macromolecules in explicit solvent.

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

A geophysical shock and air blast simulator at the National Ignition Facility

The energy partitioning energy coupling experiments at the National Ignition Facility (NIF) have been designed to measure simultaneously the coupling of energy from a laser-driven target into both ground shock and air blast overpressure to nearby media. The source target for the experiment is positioned at a known height above the ground-surface simulant and is heated by four beams from the NIF. The resulting target energy density and specific energy are equal to those of a low-yield nuclear device. The ground-shock stress waves and atmospheric overpressure waveforms that result in our test system are hydrodynamically scaled analogs of full-scale seismic and air blast phenomena. This report summarizes the development of the platform, the simulations, and calculations that underpin the physics measurements that are being made, and finally the data that were measured. Agreement between the data and simulation of the order of a factor of two to three is seen for air blast quantities such as peak overpressure. Historical underground test data for seismic phenomena measured sensor displacements; we measure the stresses generated in our ground-surrogate medium. We find factors-of-a-few agreement between our measured peak stresses and predictions with modern geophysical computer codes.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2014

Kinetic simulations and reduced modeling of longitudinal sideband instabilities in non-linear electron plasma waves

Kinetic Vlasov simulations of one-dimensional finite amplitude Electron Plasma Waves are performed in a multi-wavelength long system. A systematic study of the most unstable linear sideband mode, in particular its growth rate γ and quasi- wavenumber δk, is carried out by scanning the amplitude and wavenumber of the initial wave. Simulation results are successfully compared against numerical and analytical solutions to the reduced model by Kruer et al. [Phys. Rev. Lett. 23, 838 (1969)] for the Trapped Particle Instability (TPI). A model recently suggested by Dodin et al. [Phys. Rev. Lett. 110, 215006 (2013)], which in addition to the TPI accounts for the so-called Negative Mass Instability because of a more detailed representation of the trapped particle dynamics, is also studied and compared with simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Study of the L-mode tokamak plasma “shortfall” with local and global nonlinear gyrokinetic δf particle-in-cell simulation

The δ f particle-in-cell code GEM is used to study the transport “shortfall” problem of gyrokinetic simulations. In local simulations, the GEM results confirm the previously reported simulation results of DIII-D [Holland et al., Phys. Plasmas 16, 052301 (2009)] and Alcator C-Mod [Howard et al., Nucl. Fusion 53, 123011 (2013)] tokamaks with the continuum code GYRO. Namely, for DIII-D the simulations closely predict the ion heat flux at the core, while substantially underpredict transport towards the edge; while for Alcator C-Mod, the simulations show agreement with the experimental values of ion heat flux, at least within the range of experimental error. Global simulations are carried out for DIII-D L-mode plasmas to study the effect of edge turbulence on the outer core ion heat transport. The edge turbulence enhances the outer core ion heat transport through turbulence spreading. However, this edge turbulence spreading effect is not enough to explain the transport underprediction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Verification of gyrokinetic particle simulation of current-driven instability in fusion plasmas. II. Resistive tearing mode

In this work, global gyrokinetic particle simulation of resistive tearing modes has been developed and verified in the gyrokinetic toroidal code (GTC). GTC linear simulations in the fluid limit of the kink-tearing and resistive tearing modes in the cylindrical geometry agree well with the resistive magnetohydrodynamic eigenvalue and initial value codes. Ion kinetic effects are discovered to reduce the radial width of the tearing modes. GTC simulations of the resistive tearing modes in the toroidal geometry find that the toroidicity reduces the growth rates.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Comparisons between tokamak fueling of gas puffing and supersonic molecular beam injection in 2D simulations

Plasma fueling with high efficiency and deep injection is very important to enable fusion power performance requirements. It is a powerful and efficient way to study neutral transport dynamics and find methods of improving the fueling performance by doing large scale simulations. Two basic fueling methods, gas puffing (GP) and supersonic molecular beam injection (SMBI), are simulated and compared in realistic divertor geometry of the HL-2A tokamak with a newly developed module, named trans-neut, within the framework of BOUT++ boundary plasma turbulence code [Z. H. Wang et al., Nucl. Fusion 54, 043019 (2014)]. The physical model includes plasma density, heat and momentum transport equations along with neutral density, and momentum transport equations. Transport dynamics and profile evolutions of both plasma and neutrals are simulated and compared between GP and SMBI in both poloidal and radial directions, which are quite different from one and the other. It finds that the neutrals can penetrate about four centimeters inside the last closed (magnetic) flux surface during SMBI, while they are all deposited outside of the LCF during GP. It is the radial convection and larger inflowing flux which lead to the deeper penetration depth of SMBI and higher fueling efficiency compared to GP.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Self-consistent simulation of CdTe solar cells with active defects

We demonstrate a self-consistent numerical scheme for simulating an electronic device which contains active defects. As a specific case, we consider copper defects in cadmium telluride solar cells. The presence of copper has been shown experimentally to play a crucial role in predicting device performance. The primary source of this copper is migration away from the back contact during annealing, which likely occurs predominantly along grain boundaries. We introduce a mathematical scheme for simulating this effect in 2D and explain the numerical implementation of the system. Finally, we will give numerical results comparing our results to known 1D simulations to demonstrate the accuracy of the solver and then show results unique to the 2D case.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Simulation of magnetic hysteresis loops and magnetic Barkhausen noise of α-iron containing nonmagnetic particles

Here, the magnetic hysteresis loops and Barkhausen noise of a single α-iron with nonmagnetic particles are simulated to investigate into the magnetic hardening due to Cu-rich precipitates in irradiated reactor pressure vessel (RPV) steels. Phase field method basing Landau-Lifshitz-Gilbert (LLG) equation is used for this simulation. The results show that the presence of the nonmagnetic particle could result in magnetic hardening by making the nucleation of reversed domains difficult. The coercive field is found to increase, while the intensity of Barkhausen noise voltage is decreased when the nonmagnetic particle is introduced. Simulations demonstrate the impact of nucleation field of reversed domains on the magnetization reversal behavior and the magnetic properties.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Three- and two-dimensional simulations of counter-propagating shear experiments at high energy densities at the National Ignition Facility

Three- and two-dimensional numerical studies have been carried out to simulate recent counter-propagating shear flow experiments on the National Ignition Facility. A multi-physics three-dimensional, time-dependent radiation hydrodynamics simulation code is used. Using a Reynolds Averaging Navier-Stokes model, we show that the evolution of the mixing layer width obtained from the simulations agrees well with that measured from the experiments. A sensitivity study is conducted to illustrate a 3D geometrical effect that could confuse the measurement at late times, if the energy drives from the two ends of the shock tube are asymmetric. Implications for future experiments are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Vlasov simulations of electron-ion collision effects on damping of electron plasma waves

Collisional effects can play an essential role in the dynamics of plasma waves by setting a minimum damping rate and by interfering with wave-particle resonances. Kinetic simulations of the effects of electron-ion pitch angle scattering on Electron Plasma Waves (EPWs) are presented here. In particular, the effects of such collisions on the frequency and damping of small-amplitude EPWs for a range of collision rates and wave phase velocities are computed and compared with theory. Both the Vlasov simulations and linear kinetic theory find the direct contribution of electron-ion collisions to wave damping significantly reduced from that obtained through linearized fluid theory. To our knowledge, this simple result has not been published before. Simulations have been carried out using a grid-based (Vlasov) approach, based on a high-order conservative finite difference method for discretizing the Fokker-Planck equation describing the evolution of the electron distribution function. Details of the implementation of the collision operator within this framework are presented. Such a grid-based approach, which is not subject to numerical noise, is of particular interest for the accurate measurements of the wave damping rates.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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↗