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 649 records · Page 36

Morphology of a columnar stack of coronene molecules encapsulated in a single-walled carbon nanotube

The morphology of stacked coronene molecules encapsulated in a single-walled carbon nanotube (SWCNT) is investigated using atomistic simulation. First, the minimum energy configuration of coronene molecules in a SWCNT is sought by means of conjugate gradient (CG) minimization. Secondly, encapsulation of coronene molecules into a SWCNT existing in a coronene atmosphere is simulated by means of molecular dynamics (MD). In both of the simulations, the diameter of the SWCNT ranges from 1.35 to 1.69 nm, and the final configurations of coronene molecules within a SWCNT are examined. In a thin SWCNT, coronene molecules tilt against the radial direction of the SWCNT and slide relative to each other, whereas in a thick SWCNT, they do not tilt but rotate relative to each other. In a SWCNT of the intermediate diameter, they tilt, slide, and rotate. For the SWCNT diameter less than or equal to 1.52 nm, the mean tilt angle of the stacked coronene molecules almost linearly decreases with increasing the diameter, whereas for the diameter above 1.52 nm, it is approximately 0{sup ∘}. To check the validity of the results, the MD simulations are performed changing the density of the coronene atmosphere and the length of the SWCNT; the results prove to be valid. Finally, the effects of temperature on the mean tilt angle and mean intermolecular distance of stacked coronene molecules are examined by a rather simplified simulation, which shows that both of them increase with increasing temperature.

36 MATERIALS SCIENCE↗

Combining molecular dynamics with mesoscopic Green’s function reaction dynamics simulations

In many reaction-diffusion processes, ranging from biochemical networks, catalysis, to complex self-assembly, the spatial distribution of the reactants and the stochastic character of their interactions are crucial for the macroscopic behavior. The recently developed mesoscopic Green’s Function Reaction Dynamics (GFRD) method enables efficient simulation at the particle level provided the microscopic dynamics can be integrated out. Yet, many processes exhibit non-trivial microscopic dynamics that can qualitatively change the macroscopic behavior, calling for an atomistic, microscopic description. We propose a novel approach that combines GFRD for simulating the system at the mesoscopic scale where particles are far apart, with a microscopic technique such as Langevin dynamics or Molecular Dynamics (MD), for simulating the system at the microscopic scale where reactants are in close proximity. This scheme defines the regions where the particles are close together and simulated with high microscopic resolution and those where they are far apart and simulated with lower mesoscopic resolution, adaptively on the fly. The new multi-scale scheme, called MD-GFRD, is generic and can be used to efficiently simulate reaction-diffusion systems at the particle level.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Smart darting diffusion Monte Carlo: Applications to lithium ion-Stockmayer clusters

In a recent investigation [K. Roberts et al., J. Chem. Phys. 136, 074104 (2012)], we have shown that, for a sufficiently complex potential, the Diffusion Monte Carlo (DMC) random walk can become quasiergodic, and we have introduced smart darting-like moves to improve the sampling. In this article, we systematically characterize the bias that smart darting moves introduce in the estimate of the ground state energy of a bosonic system. We then test a simple approach to eliminate completely such bias from the results. The approach is applied for the determination of the ground state of lithium ion-n–dipoles clusters in the n = 8–20 range. For these, the smart darting diffusion Monte Carlo simulations find the same ground state energy and mixed-distribution as the traditional approach for n < 14. In larger systems we find that while the ground state energies agree quantitatively with or without smart darting moves, the mixed-distributions can be significantly different. Some evidence is offered to conclude that introducing smart darting-like moves in traditional DMC simulations may produce a more reliable ground state mixed-distribution.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Output power fluctuations due to different weights of macro particles used in particle-in-cell simulations of Cerenkov devices

The output power fluctuations caused by weights of macro particles used in particle-in-cell (PIC) simulations of a backward wave oscillator and a travelling wave tube are statistically analyzed. It is found that the velocities of electrons passed a specific slow-wave structure form a specific electron velocity distribution. The electron velocity distribution obtained in PIC simulation with a relative small weight of macro particles is considered as an initial distribution. By analyzing this initial distribution with a statistical method, the estimations of the output power fluctuations caused by different weights of macro particles are obtained. The statistical method is verified by comparing the estimations with the simulation results. The fluctuations become stronger with increasing weight of macro particles, which can also be determined reversely from estimations of the output power fluctuations. With the weights of macro particles optimized by the statistical method, the output power fluctuations in PIC simulations are relatively small and acceptable.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Signal pulse emulation for scintillation detectors using Geant4 Monte Carlo with light tracking simulation

The anode pulse of a photomultiplier tube (PMT) coupled with a scintillator is used for pulse shape discrimination (PSD) analysis. We have developed a novel emulation technique for the PMT anode pulse based on optical photon transport and a PMT response function. The photon transport was calculated using Geant4 Monte Carlo code and the response function with a BC408 organic scintillator. The obtained percentage RMS value of the difference between the measured and simulated pulse with suitable scintillation properties using GSO:Ce (0.4, 1.0, 1.5 mol%), LaBr{sub 3}:Ce and BGO scintillators were 2.41%, 2.58%, 2.16%, 2.01%, and 3.32%, respectively. The proposed technique demonstrates high reproducibility of the measured pulse and can be applied to simulation studies of various radiation measurements.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

The charge imbalance in ultracold plasmas

Ultracold plasmas are regarded as quasineutral but not strictly neutral. The results of charge imbalance in the expansion of ultracold plasmas are reported. The calculations are performed by a full molecular-dynamics simulation. The details of the electron velocity distributions are calculated without the assumption of electron global thermal equilibrium and Boltzmann distribution. Spontaneous evolutions of the charge imbalance from the initial states with perfect neutrality are given in the simulations. The expansion of outer plasma slows down with the charge imbalance. The influences of plasma size and parameters on the charge imbalance are discussed. The radial profiles of electron temperature are given for the first time, and the self-similar expansion can still occur even if there is no global thermal equilibrium. The electron disorder induced heating is also found in the simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

HOW THE FIRST STARS SHAPED THE FAINTEST GAS-DOMINATED DWARF GALAXIES

Low-mass dwarf galaxies are very sensitive test-beds for theories of cosmic structure formation since their weak gravitational fields allow the effects of the relevant physical processes to clearly stand out. Up to now, no unified account has existed of the sometimes seemingly conflicting properties of the faintest isolated dwarfs in and around the Local Group, such as Leo T and the recently discovered Leo P and Pisces A systems. Using new numerical simulations, we show that this serious challenge to our understanding of galaxy formation can be effectively resolved by taking into account the regulating influence of the ultraviolet radiation of the first population of stars on a dwarf’s star formation rate while otherwise staying within the standard cosmological paradigm for structure formation. These simulations produce faint, gas-dominated, star-forming dwarf galaxies that lie on the baryonic Tully–Fisher relation and that successfully reproduce a broad range of chemical, kinematical, and structural observables of real late-type dwarf galaxies. Furthermore, we stress the importance of obtaining properties of simulated galaxies in a manner as close as possible to the typically employed observational techniques.

79 ASTRONOMY AND ASTROPHYSICS↗

SU-E-T-30: A Factor for Converting Dose to a Gold Nanoparticle Mixture to a Biologically-Relevant Dose

Purpose: Monte Carlo studies of gold nanoparticle (GNP) dose enhancement on macroscopic scales in radiotherapy have modeled GNPs in tissue as a homogeneous mixture of gold and tissue. Using an explicit model of GNPs randomly positioned in a small volume (1 µm{sup 3}) of tissue, this study aims to quantify the dose to the biologically relevant component of a goldtissue mixture, enabling a conversion from macroscopically-scored dose. Methods: Using the PENELOPE Monte Carlo code with the penEasy package, we modeled a 1 µm{sup 3} volume containing either a tissue-gold mixture or GNPs suspended in ICRU 4-component tissue at various gold concentrations (0, 5, 10, and 15 mg Au/g tissue) and GNP diameters (20, 30, 40, 50 nm). The volume was irradiated with monoenergetic photon and electron beams, ranging from 110 eV to 6 MeV. Interaction forcing was utilized to increase simulation efficiency. Energy deposition was scored in the tissue for each case and was converted to dose. For each scenario, we calculated a conversion factor, the ratio of dose-to-tissue to dose-to-mixture as a function of energy. Results: The conversion factor was plotted as a function of energy for both photons and electrons. For electrons, the conversion factor was relatively unaffected by any of the parameters, including energy, ranging between 0.98–1.02. For photons, the factor was very energy dependent, with a range of 0.49–1.02. The factor was lowest for 10–100 keV photons. The conversion factor generally decreased with increasing GNP concentration and increasing GNP size. Conclusion: With a large variation in the conversion factor with incident energy, dose deposition is dependent on the spectrum incident on a volume. By scoring the energy spectrum in a given volume, one can provide a scenario-specific conversion factor, allowing fast, detailed Monte Carlo simulations without the need for explicit GNP-definition.

60 APPLIED LIFE SCIENCES↗

Molecular dynamics studies of pathways of water movement in cyanobacterial photosystem II

Photosystem II (PSII) catalyzes the light-induced generation of oxygen from water. The oxygen-evolving complex is buried deep in the protein on the lumenal side of PSII, and water molecules need to pass through protein subunits to reach the active site—the manganese cluster. Previous studies on the elucidation of water channels in PSII were based on an analysis of the cavities in the static PSII structure determined by X-ray diffraction. In the present study, we perform molecular dynamics simulations of the water movement in the transport system of PSII.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Scalable Nonlinear Solvers for Fully Implicit Coupled Nuclear Fuel Modeling. Final Report

The focus of the project is on the development and customization of some highly scalable domain decomposition based preconditioning techniques for the numerical solution of nonlinear, coupled systems of partial differential equations (PDEs) arising from nuclear fuel simulations. These high-order PDEs represent multiple interacting physical fields (for example, heat conduction, oxygen transport, solid deformation), each is modeled by a certain type of Cahn-Hilliard and/or Allen-Cahn equations. Most existing approaches involve a careful splitting of the fields and the use of field-by-field iterations to obtain a solution of the coupled problem. Such approaches have many advantages such as ease of implementation since only single field solvers are needed, but also exhibit disadvantages. For example, certain nonlinear interactions between the fields may not be fully captured, and for unsteady problems, stable time integration schemes are difficult to design. In addition, when implemented on large scale parallel computers, the sequential nature of the field-by-field iterations substantially reduces the parallel efficiency. To overcome the disadvantages, fully coupled approaches have been investigated in order to obtain full physics simulations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Modeling investigation of the stability and irradiation-induced evolution of nanoscale precipitates in advanced structural materials

Materials used in extremely hostile environment such as nuclear reactors are subject to a high flux of neutron irradiation, and thus vast concentrations of vacancy and interstitial point defects are produced because of collisions of energetic neutrons with host lattice atoms. The fate of these defects depends on various reaction mechanisms which occur immediately following the displacement cascade evolution and during the longer-time kinetically dominated evolution such as annihilation, recombination, clustering or trapping at sinks of vacancies, interstitials and their clusters. The long-range diffusional transport and evolution of point defects and self-defect clusters drive a microstructural and microchemical evolution that are known to produce degradation of mechanical properties including the creep rate, yield strength, ductility, or fracture toughness, and correspondingly affect material serviceability and lifetimes in nuclear applications. Therefore, a detailed understanding of microstructural evolution in materials at different time and length scales is of significant importance. The primary objective of this work is to utilize a hierarchical computational modeling approach i) to evaluate the potential for nanoscale precipitates to enhance point defect recombination rates and thereby the self-healing ability of advanced structural materials, and ii) to evaluate the stability and irradiation-induced evolution of such nanoscale precipitates resulting from enhanced point defect transport to and annihilation at precipitate interfaces. This project will utilize, and as necessary develop, computational materials modeling techniques within a hierarchical computational modeling approach, principally including molecular dynamics, kinetic Monte Carlo and spatially-dependent cluster dynamics modeling, to identify and understand the most important physical processes relevant to promoting the “selfhealing” or radiation resistance in advanced materials containing nanoscale precipitates. In particular, the interfacial structure of embedded nanoscale precipitates will be evaluated by electronic- and atomic-scale modeling methods, and the efficiency of the validated interfaces for trapping point defects will next be evaluated by atomic-scale modeling (e.g., determining the sink strength of the precipitates), addressing key questions related to the optimal interface characteristics to attract point defects and enhance their recombination. Kinetic models will also be developed to simulate microstructural evolution of the nanoscale features and irradiation produced defect clusters, and compared with observed microstructural changes.

36 MATERIALS SCIENCE↗

( Chandra ) view of magnetically confined wind in HD 191612: theory versus observations

High-resolution spectra of the magnetic star HD 191612 were acquired using the Chandra X-ray Observatory at both maximum and minimum emission phases. We confirm the flux and hardness variations previously reported with XMM-Newton, demonstrating the high repeatability of the behavior of HD 191612 over a decade. The line profiles appear typical for magnetic massive stars: no significant line shift, relatively narrow lines for high-Z elements, and formation radius at about 2 R&lowast;. Line ratios confirm the softening of the X-ray spectrum at the minimum emission phase. Shift or width variations appear of limited amplitude at most (slightly lower velocity and slightly increased broadening at minimum emission phase, but within 1–2σ of values at maximum). In addition, a fully self-consistent 3D magnetohydrodynamic simulation of the confined wind in HD 191612 was performed. The simulation results were directly fitted to the data, leading to a remarkable agreement overall between them.

79 ASTRONOMY AND ASTROPHYSICS↗

Modeling Turbulence in Galactic Centers

Turbulence is a prevalent phenomenon in the interstellar medium, and in particular, the environment at the centers of galaxies. For example, detailed observations of the Milky Way’s Central Molecular Zone (CMZ) revealed that it has a complex and turbulent structure. Turbulence on galactic scales is often modeled using star formation and feedback. However, these effects do not appear to be sufficient for explaining the high-velocity dispersion observed in the CMZ, indicating that additional gas-stirring processes are likely to be operating. Here we introduce a proof-of-concept method to drive turbulence in gas that orbits under the influence of a galactic potential. Instead of relying on a particular physical mechanism, we have adopted a Fourier forcing module and have applied it using a smoothed particle hydrodynamics code. To test our method, we performed simulations of a simplistic model of the CMZ. Our turbulence injection method is capable of balancing the self-gravity of the gas, which allows us to run the simulations for long timescales and thereby follow the evolution of the CMZ. Our results show that turbulence induces a flocculent spiral pattern in our model, analogous to that found in galactic-scale simulations. Furthermore, we find that our turbulence injection method induces inward migration of gas, a result consistent with previous numerical simulations. We submit that this injection method is a promising new tool to simulate turbulence in galactic centers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Neutronics and computational fluid dynamics analyses for a natural circulation loop of molten fuel salt

A conceptual reactor system of molten fuel salt loop designed for Mo-99 production was analyzed by coupling TRITON neutronics and Fluent computational fluid dynamics (CFD) simulations. The loop system operated in 1 MW thermal power features natural circulation and potential extraction of non-soluble fission products (FP) in fuel salt. The behavior of natural circulation was studied and the power and temperature distributions in equilibrium were obtained. A long-term fuel depletion analysis of the system considering fuel circulation and reprocessing was then performed using an approximate stepwise sequence of simulations. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DNS Study of Supercritical CO{sub 2} Flow Using Spectral Element Code Nek5000

The study of supercritical CO{sub 2} (s-CO{sub 2}) heat transfer has become popular in recent decades. The driving factor is the growing interest of using s-CO{sub 2} in future power conversion cycles. The supercritical CO{sub 2} Brayton cycle has gained the most focus for its high efficiency, compact design, economics, etc, . In this cycle, heat exchangers are used to heat up or cool down CO{sub 2}. Thus, the heat transfer characteristics of s-CO{sub 2} need to be well understood for heat exchanger design. In this research, we use the open-source spectral-element code Nek5000, developed by Argonne National Laboratory, to perform a direct numerical simulation (DNS) for s-CO{sub 2} flow through a straight heated tube. The pioneering paper by Bae et al. is used to validate the results from the Nek5000 simulation. (authors)

74 ATOMIC AND MOLECULAR PHYSICS↗

An Investigation of Chaotic Diffusion in a Family of Hamiltonian Mappings Whose Angles Diverge in the Limit of Vanishingly Action

The chaotic diffusion for a family of Hamiltonian mappings whose angles diverge in the limit of vanishingly action is investigated by using the solution of the diffusion equation. The system is described by a two-dimensional mapping for the variables action, I, and angle, θ and controlled by two control parameters: (i) ϵ, controlling the nonlinearity of the system, particularly a transition from integrable for ϵ=0 to non-integrable for ϵ≠0 and; (ii) γ denoting the power of the action in the equation defining the angle. For ϵ≠0 the phase space is mixed and chaos is present in the system leading to a finite diffusion in the action characterized by the solution of the diffusion equation. The analytical solution is then compared to the numerical simulations showing a remarkable agreement between the two procedures.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A Simple Kinetic Model for the Phase Transition of the van der Waals Fluid

A simple kinetic model, which is presumably minimum, for the phase transition of the van der Waals fluid is presented. In the model, intermolecular collisions for a dense gas has not been treated faithfully. Instead, the expected interactions as the non-ideal gas effect are confined in a self-consistent force term. Collision term plays just a role of thermal bath. Accordingly, it conserves neither momentum nor energy, even globally. It is demonstrated that (i) by a natural separation of the mean-field self-consistent potential, the potential for the non-ideal gas effect is determined from the equation of state for the van der Waals fluid, with the aid of the balance equation of momentum, (ii) a functional which monotonically decreases in time is identified by the H theorem and is found to have a close relation to the Helmholtz free energy in thermodynamics, and (iii) the Cahn–Hilliard type equation is obtained in the continuum limit of the present kinetic model. Numerical simulations based on the Cahn–Hilliard type equation are also performed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Accounting for Melt Flow Pattern and Solid Fraction Evolution in DC Casting of Al-Cu Alloy Using v{sup 2}−f Turbulence Model

Experimental techniques remain reliable approaches for investigating the effects of parameters of the direct chill (DC) casting process on macrosegregation, grain structure formation, and surface quality. However, experimental analysis methods of these effects are expensive and can often prove difficult to implement. Understanding of the complex physics of the solidification process remains limited, and this makes it indispensable to use accurate numerical modeling tools as a complement to experiments. Numerical modeling without consideration of turbulent flow is unsuitable for application to solidification occurring during DC casting, since significant turbulent flow exists in the upper section of the cast, where the liquid pool exists. The present work employed a low-Re v{sup 2}−f turbulence model and a dual-zone solidification flow model to predict the melt flow pattern, turbulence level, sump depth, transition zone thickness, and mushy zone thickness in the DC casting of Al-Cu alloy. The obtained results provided evidence of significant damping of the flow in the slurry zone. Furthermore, results of different simulation cases revealed a narrower transition zone and shallower sump.

36 MATERIALS SCIENCE↗