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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
DOE OSTI2012

Test Cases for Wind Power Plant Dynamic Models on Real-Time Digital Simulator: Preprint

The objective of this paper is to present test cases for wind turbine generator and wind power plant models commonly used during commissioning of wind power plants to ensure grid integration compatibility. In this paper, different types of wind power plant models based on the Western Electricity Coordinating Council Wind Generator Modeling Group's standardization efforts are implemented on areal-time digital simulator, and different test cases are used to gauge their grid integration capability. The low-voltage ride through and reactive power support capability and limitations of wind turbine generators under different grid conditions are explored. Several types of transient events (e.g., symmetrical and unsymmetrical faults, frequency dips) are included in the test cases. Thedifferences in responses from different types of wind turbine are discussed in detail.

17 WIND ENERGY↗
DOE OSTI2015

Monte Carlo simulation methods in moment-based scale-bridging algorithms for thermal radiative-transfer problems

We present a moment-based acceleration algorithm applied to Monte Carlo simulation of thermal radiative-transfer problems. Our acceleration algorithm employs a continuum system of moments to accelerate convergence of stiff absorption–emission physics. The combination of energy-conserving tallies and the use of an asymptotic approximation in optically thick regions remedy the difficulties of local energy conservation and mitigation of statistical noise in such regions. We demonstrate the efficiency and accuracy of the developed method. We also compare directly to the standard linearization-based method of Fleck and Cummings [1]. A factor of 40 reduction in total computational time is achieved with the new algorithm for an equivalent (or more accurate) solution as compared with the Fleck–Cummings algorithm.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Numerical simulations of depressurization-induced gas production from gas hydrates using 3-D heterogeneous models of L-Pad, Prudhoe Bay Unit, North Slope Alaska

Gas production potential is estimated using a three-dimensional reservoir model based on gas hydrate deposits located in the Prudhoe Bay region of the Alaska's North Slope. The model incorporates two hydrate-bearing sand units using detailed reservoir geological and structural information obtained from past and recent drilling programs. Geostatistical porosity models, conditioned to log data from 78 wells drilled in the vicinity of the Prudhoe Bay “L-Pad,” were developed, providing 3D heterogeneity in porosity, porosity-dependent hydrate saturation, and intrinsic permeability. The simulations utilize both vertical and inclined wellbores to induce depressurization of the reservoir at a constant bottom-hole pressure. The results show the superior performance of the inclined well design. Average gas production rates during the first five years were ~6.0 × 10 4 ST m 3 /day (~2.1 MMSCF/day) and ~2.7 × 10 4 ST m 3 /day (~1.0 MMSCF/day) for incline and vertical wells, respectively. After 30 years 5.3 × 10 8 - 5.7 × 10 8 ST m 3 (18.6–20.2 BSCF) and 6.2 × 10 8 - 6.4 × 10 8 ST m 3 (22.0–22.7 BSCF) gas were produced using the vertical and inclined well configurations, respectively. The analysis reveals that 2D reservoir models with homogeneous representations for porosity and hydrate saturation significantly underestimate production potential. The heterogeneity implemented in this work provides complex porous network and preferential pathways for mobile phase flow to a producing well. Consequently, no secondary hydrate formation around a wellbore is predicted in contrast to models that utilize uniform porosity and hydrate saturation distributions.

03 NATURAL GAS↗
DOE OSTI2012

Hydration structures of U(III) and U(IV) ions from ab initio molecular dynamics simulations

We apply DFT+U-based ab initio molecular dynamics simulations to study the hydration structures of U(III) and U(IV) ions, pertinent to redox reactions associated with uranium salts in aqueous media. U(III) is predicted to be coordinated to 8 water molecules, while U(IV) has a hydration number between 7 and 8. At least one of the innershell water molecules of the hydrated U(IV) complex becomes spontaneously deprotonated. As a result, the U(IV)-O pair correlation function exhibits a satellite peak at 2.15 A associated with the shorter U(IV)-(OH{sup -}) bond. This feature is not accounted for in analysis of extended x-ray absorption fine structure and x-ray adsorption near edge structure measurements, which yield higher estimates of U(IV) hydration numbers. This suggests that it may be useful to include the effect of possible hydrolysis in future interpretation of experiments, especially when the experimental pH is close to the reported hydrolysis equilibrium constant value.

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

First-principles studies on vacancy-modified interstitial diffusion mechanism of oxygen in nickel, associated with large-scale atomic simulation techniques

This paper is concerned with the prediction of oxygen diffusivities in fcc nickel from first-principles calculations and large-scale atomic simulations. Considering only the interstitial octahedral to tetrahedral to octahedral minimum energy pathway for oxygen diffusion in fcc lattice, greatly underestimates the migration barrier and overestimates the diffusivities by several orders of magnitude. The results indicate that vacancies in the Ni-lattice significantly impact the migration barrier of oxygen in nickel. Incorporation of the effect of vacancies results in predicted diffusivities consistent with available experimental data. First-principles calculations show that at high temperatures the vacancy concentration is comparable to the oxygen solubility, and there is a strong binding energy and a redistribution of charge density between the oxygen atom and vacancy. Consequently, there is a strong attraction between the oxygen and vacancy in the Ni lattice, which impacts diffusion.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

Molecular simulation of thermodynamic and transport properties for the H{sub 2}O+NaCl system

Molecular dynamics and Monte Carlo simulations have been carried out to obtain thermodynamic and transport properties of the binary mixture H{sub 2}O+NaCl at temperatures from T = 298 to 473 K. In particular, vapor pressures, liquid densities, viscosities, and vapor-liquid interfacial tensions have been obtained as functions of pressure and salt concentration. Several previously proposed fixed-point-charge models that include either Lennard-Jones (LJ) 12-6 or exponential-6 (Exp6) functional forms to describe non-Coulombic interactions were studied. In particular, for water we used the SPC and SPC/E (LJ) models in their rigid forms, a semiflexible version of the SPC/E (LJ) model, and the Errington-Panagiotopoulos Exp6 model; for NaCl, we used the Smith-Dang and Joung-Cheatham (LJ) parameterizations as well as the Tosi-Fumi (Exp6) model. While none of the model combinations are able to reproduce simultaneously all target properties, vapor pressures are well represented using the SPC plus Joung-Cheathem model combination, and all LJ models do well for the liquid density, with the semiflexible SPC/E plus Joung-Cheatham combination being the most accurate. For viscosities, the combination of rigid SPC/E plus Smith-Dang is the best alternative. For interfacial tensions, the combination of the semiflexible SPC/E plus Smith-Dang or Joung-Cheatham gives the best results. Inclusion of water flexibility improves the mixture densities and interfacial tensions, at the cost of larger deviations for the vapor pressures and viscosities. The Exp6 water plus Tosi-Fumi salt model combination was found to perform poorly for most of the properties of interest, in particular being unable to describe the experimental trend for the vapor pressure as a function of salt concentration.

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

Development of an interpretive simulation tool for the proton radiography technique

Proton radiography is a useful diagnostic of high energy density (HED) plasmas under active theoretical and experimental development. In this paper, we describe a new simulation tool that interacts realistic laser-driven point-like proton sources with three dimensional electromagnetic fields of arbitrary strength and structure and synthesizes the associated high resolution proton radiograph. The present tool’s numerical approach captures all relevant physics effects, including effects related to the formation of caustics. Electromagnetic fields can be imported from particle-in-cell or hydrodynamic codes in a streamlined fashion, and a library of electromagnetic field “primitives” is also provided. This latter capability allows users to add a primitive, modify the field strength, rotate a primitive, and so on, while quickly generating a high resolution radiograph at each step. In this way, our tool enables the user to deconstruct features in a radiograph and interpret them in connection to specific underlying electromagnetic field elements. We show an example application of the tool in connection to experimental observations of the Weibel instability in counterstreaming plasmas, using ∼10{sup 8} particles generated from a realistic laser-driven point-like proton source, imaging fields which cover volumes of ∼10 mm{sup 3}. Insights derived from this application show that the tool can support understanding of HED plasmas.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2016

A poly-epoxy surface explored by Hartree-Fock ΔSCF simulations of C1s XPS spectra

Whereas poly-epoxy polymers represent a class of materials with a wide range of applications, the structural disorder makes them difficult to model. In the present work, we use good experimental model samples in the sense that they are pure, fully polymerized, flat and smooth, defect-free, and suitable for ultrahigh vacuum x-ray photoelectron spectroscopy, XPS, experiments. In parallel, we perform Hartree-Fock, HF, calculations of the binding energies, BEs, of the C1s electrons in a model molecule composed of the two constituents of the poly-epoxy sample. These C1s BEs were determined using the HF ΔSCF method, which is known to yield accurate values, especially for the shifts of the BEs, ΔBEs. We demonstrate the benefits of combining rigorous theory with careful XPS measurements in order to obtain correct assignments of the C1s XPS spectra of the polymer sample. Both the relative binding energies—by the ΔSCF method—and relative intensities—in the sudden approximation, SA, are calculated. It results in an excellent match with the experimental spectra. We are able to identify 9 different chemical environments under the C1s peak, where an exclusively experimental work would have found only 3 contributions. In addition, we observe that some contributions are localized at discrete binding energies, whereas others allow a much wider range because of the variation of their second neighbor bound polarization. Therefore, HF-ΔSCF simulations significantly increase the spectral resolution of XPS and thus offer a new avenue for the exploration of the surface of polymers.

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

System-size corrections for self-diffusion coefficients calculated from molecular dynamics simulations: The case of CO{sub 2}, n-alkanes, and poly(ethylene glycol) dimethyl ethers

Molecular dynamics simulations were carried out to study the self-diffusion coefficients of CO{sub 2}, methane, propane, n-hexane, n-hexadecane, and various poly(ethylene glycol) dimethyl ethers (glymes in short, CH{sub 3}O–(CH{sub 2}CH{sub 2}O){sub n}–CH{sub 3} with n = 1, 2, 3, and 4, labeled as G1, G2, G3, and G4, respectively) at different conditions. Various system sizes were examined. The widely used Yeh and Hummer [J. Phys. Chem. B 108, 15873 (2004)] correction for the prediction of diffusion coefficient at the thermodynamic limit was applied and shown to be accurate in all cases compared to extrapolated values at infinite system size. The magnitude of correction, in all cases examined, is significant, with the smallest systems examined giving for some cases a self-diffusion coefficient approximately 15% lower than the infinite system-size extrapolated value. The results suggest that finite size corrections to computed self-diffusivities must be used in order to obtain accurate results.

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

Atomistic simulations of tungsten surface evolution under low-energy neon implantation

Tungsten is a candidate material for the divertor of fusion reactors, where it will be subject to a high flux of particles coming from the fusion plasma as well as a significant heat load. Under helium plasma exposure in fusion-reactor-like conditions, a nanostructured morphology is known to form on the tungsten surface in certain temperature and incident energy ranges, although the formation mechanism is not fully established. A recent experimental study (Yajima et al 2013 Plasma Sci. Technol. 15 282-6) using neon or argon exposure did not produce similar nanostructure. Here, this article presents molecular dynamics simulations of neon implantation in tungsten aimed at investigating the surface evolution and elucidating the role of noble gas mass in fuzz formation. Finally, in contrast to helium, neon impacts can sputter both tungsten and previously implanted neon atoms. The shorter range of neon ions, along with sputtering, limit the formation of large bubbles and likely prevents nanostructure formation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗