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,405 records · Page 78
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Properties of a soft-core model of methanol: An integral equation theory and computer simulation study

Thermodynamic and structural properties of a coarse-grained model of methanol are examined by Monte Carlo simulations and reference interaction site model (RISM) integral equation theory. Methanol particles are described as dimers formed from an apolar Lennard-Jones sphere, mimicking the methyl group, and a sphere with a core-softened potential as the hydroxyl group. Different closure approximations of the RISM theory are compared and discussed. The liquid structure of methanol is investigated by calculating site-site radial distribution functions and static structure factors for a wide range of temperatures and densities. Results obtained show a good agreement between RISM and Monte Carlo simulations. The phase behavior of methanol is investigated by employing different thermodynamic routes for the calculation of the RISM free energy, drawing gas-liquid coexistence curves that match the simulation data. Preliminary indications for a putative second critical point between two different liquid phases of methanol are also discussed.

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

Energy resolved actinometry for simultaneous measurement of atomic oxygen densities and local mean electron energies in radio-frequency driven plasmas

A diagnostic method for the simultaneous determination of atomic oxygen densities and mean electron energies is demonstrated for an atmospheric pressure radio-frequency plasma jet. The proposed method is based on phase resolved optical emission measurements of the direct and dissociative electron-impact excitation dynamics of three distinct emission lines, namely, Ar 750.4 nm, O 777.4 nm, and O 844.6 nm. The energy dependence of these lines serves as basis for analysis by taking into account two line ratios. In this frame, the method is highly adaptable with regard to pressure and gas composition. Results are benchmarked against independent numerical simulations and two-photon absorption laser-induced fluorescence experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Looking for Auger signatures in III-nitride light emitters: A full-band Monte Carlo perspective

Recent experiments of electron emission spectroscopy (EES) on III-nitride light-emitting diodes (LEDs) have shown a correlation between droop onset and hot electron emission at the cesiated surface of the LED p-cap. The observed hot electrons have been interpreted as a direct signature of Auger recombination in the LED active region, as highly energetic Auger-excited electrons would be collected in long-lived satellite valleys of the conduction band so that they would not decay on their journey to the surface across the highly doped p-contact layer. We discuss this interpretation by using a full-band Monte Carlo model based on first-principles electronic structure and lattice dynamics calculations. The results of our analysis suggest that Auger-excited electrons cannot be unambiguously detected in the LED structures used in the EES experiments. Additional experimental and simulative work are necessary to unravel the complex physics of GaN cesiated surfaces.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Shock wave acceleration of protons in inhomogeneous plasma interacting with ultrashort intense laser pulses

The acceleration of protons, triggered by solitary waves in expanded solid targets is investigated using particle-in-cell simulations. The near-critical density plasma is irradiated by ultrashort high power laser pulses, which generate the solitary wave. The transformation of this soliton into a shock wave during propagation in plasma with exponentially decreasing density profile is described analytically, which allows to obtain a scaling law for the proton energy. The high quality proton bunch with small energy spread is produced by reflection from the shock-front. According to the 2D simulations, the mechanism is stable only if the laser pulse duration is shorter than the characteristic development time of the parasitic Weibel instability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

The microscopic Z-pinch process of current-carrying rarefied deuterium plasma shell

For insight into the microscopic mechanism of Z-pinch dynamic processes, a code of two-dimensional particle-in-cell (PIC) simulation has been developed in cylindrical coordinates. In principle, the Z-pinch of current-carrying rarefied deuterium plasma shell has been simulated by means of this code. Many results related to the microscopic processes of the Z-pinch are obtained. They include the spatio-temporal distributions of electromagnetic field, current density, forces experienced by the ions and electrons, positions and energy distributions of particles, and trailing mass and current. In radial direction, the electric and magnetic forces exerted on the electrons are comparable in magnitude, while the forces exerted on the ions are mainly the electric forces. So in the Z-pinch process, the electrons are first accelerated in Z direction and get higher velocities; then, they are driven inwards to the axis at the same time by the radial magnetic forces (i.e., Lorentz forces) of them. That causes the separations between the electrons and ions because the ion mass is much larger than the electron's, and in turn a strong electrostatic field is produced. The produced electrostatic field attracts the ions to move towards the electrons. When the electrons are driven along the radial direction to arrive at the axis, they shortly move inversely due to the static repellency among them and their tiny mass, while the ions continue to move inertially inwards, and later get into stagnation, and finally scatter outwards. Near the stagnation, the energies of the deuterium ions mostly range from 0.3 to 6 keV, while the electron energies are mostly from 5 to 35 keV. The radial components, which can contribute to the pinched plasma temperature, of the most probable energies of electron and ion at the stagnation are comparable to the Bennett equilibrium temperature (about 1 keV), and also to the highest temperatures of electron and ion obtained in one dimensional radiation magnetohydrodynamic simulation of the plasma shell Z-pinch. The trailing mass is about 20% of the total mass of the shell, and the maximum trailing current is about 7% of the driven current under our trailing definition. Our PIC simulation also demonstrates that the plasma shell first experiences a snow-plow like implosion process, which is relatively stable.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Muon neutrino charged current inclusive charged pion (CCπ{sup ±}) production in MINERνA

The production of charged pions by neutrinos interacting on nuclei is of great interest in nuclear physics and neutrino oscillation experiments. The MINERνA experiment is working towards releasing the world’s first high statistics neutrino pion production measurements in a few-GeV neutrino beam. We describe MINERνA’s CCπ{sup ±} analysis event selection in both the neutrino and antineutrino beams, noting reconstruction resolutions and kinematic limits. We also show area-normalized data-simulation comparisons of the reconstructed muon and charged pion kinetic energy distributions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

The effect of particle size on the morphology and thermodynamics of diblock copolymer/tethered-particle membranes

A combination of self-consistent field theory and density functional theory was used to examine the effect of particle size on the stable, 3-dimensional equilibrium morphologies formed by diblock copolymers with a tethered nanoparticle attached either between the two blocks or at the end of one of the blocks. Particle size was varied between one and four tenths of the radius of gyration of the diblock polymer chain for neutral particles as well as those either favoring or disfavoring segments of the copolymer blocks. Phase diagrams were constructed and analyzed in terms of thermodynamic diagrams to understand the physics associated with the molecular-level self-assembly processes. Typical morphologies were observed, such as lamellar, spheroidal, cylindrical, gyroidal, and perforated lamellar, with the primary concentration region of the tethered particles being influenced heavily by particle size and tethering location, strength of the particle-segment energetic interactions, chain length, and copolymer radius of gyration. The effect of the simulation box size on the observed morphology and system thermodynamics was also investigated, indicating possible effects of confinement upon the system self-assembly processes.

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

Wavelength dependence of femtosecond laser-induced damage threshold of optical materials

An experimental and numerical study of the laser-induced damage of the surface of optical material in the femtosecond regime is presented. The objective of this work is to investigate the different processes involved as a function of the ratio of photon to bandgap energies and compare the results to models based on nonlinear ionization processes. Experimentally, the laser-induced damage threshold of optical materials has been studied in a range of wavelengths from 1030 nm (1.2 eV) to 310 nm (4 eV) with pulse durations of 100 fs with the use of an optical parametric amplifier system. Semi-conductors and dielectrics materials, in bulk or thin film forms, in a range of bandgap from 1 to 10 eV have been tested in order to investigate the scaling of the femtosecond laser damage threshold with the bandgap and photon energy. A model based on the Keldysh photo-ionization theory and the description of impact ionization by a multiple-rate-equation system is used to explain the dependence of laser-breakdown with the photon energy. The calculated damage fluence threshold is found to be consistent with experimental results. From these results, the relative importance of the ionization processes can be derived depending on material properties and irradiation conditions. Moreover, the observed damage morphologies can be described within the framework of the model by taking into account the dynamics of energy deposition with one dimensional propagation simulations in the excited material and thermodynamical considerations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Feasibility study of archaeological structures scanning by muon tomography

One of the main concerns in archaeology is to find of a method to study precisely archaeological structures in the least invasive way possible to avoid damage. The requirement of preserving the structures integrity prevents, in the case of pyramids or tumuli, the study of any internal structure (halls or tombs) which are not reachable by existing corridors. One non-invasive method is the muon tomography. By placing a detector which allows to register the muon direction after the structure, it is possible to have an idea of its composition based on the attenuation of the muon flux, which depends on the material length and density that muons have crossed. This technique, alone or together with other exploration techniques as seismic tomography or electrical resistivity tomography, can provide useful information about the internal structure of the archaeological form that can not be obtained by conventional archaeological methods. In this work, the time measurement necessary to obtain a significant result about the composition of an archaeological structure is estimated. To do that, a Monte Carlo simulation framework based on the MUSIC software, properly tuned for this study, has been developed. The particular case of the Kastas Amfipoli Macedonian tumulus has been considered to perform the simulations.

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

Dichotomous noise models of gene switches

Molecular noise in gene regulatory networks has two intrinsic components, one part being due to fluctuations caused by the birth and death of protein or mRNA molecules which are often present in small numbers and the other part arising from gene state switching, a single molecule event. Stochastic dynamics of gene regulatory circuits appears to be largely responsible for bifurcations into a set of multi-attractor states that encode different cell phenotypes. The interplay of dichotomous single molecule gene noise with the nonlinear architecture of genetic networks generates rich and complex phenomena. In this paper, we elaborate on an approximate framework that leads to simple hybrid multi-scale schemes well suited for the quantitative exploration of the steady state properties of large-scale cellular genetic circuits. Through a path sum based analysis of trajectory statistics, we elucidate the connection of these hybrid schemes to the underlying master equation and provide a rigorous justification for using dichotomous noise based models to study genetic networks. Numerical simulations of circuit models reveal that the contribution of the genetic noise of single molecule origin to the total noise is significant for a wide range of kinetic regimes.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

Indium tin oxide refractometer in the visible and near infrared via lossy mode and surface plasmon resonances with Kretschmann configuration

An indium tin oxide (ITO) refractometer based on the generation of lossy mode resonances (LMRs) and surface plasmon resonances (SPRs) is presented. Both LMRs and SPRs are excited, in a single setup, under grazing angle incidence with Kretschmann configuration in an ITO thin-film deposited on a glass slide. The sensing capabilities of the device are demonstrated using several solutions of glycerin and water with refractive indices ranging from 1.33 to 1.47. LMRs are excited in the visible range, from 617 nm to 682 nm under TE polarization and from 533 nm to 637 nm under TM polarization, with a maximum sensitivity of 700 nm/RIU and 1200 nm/RIU, respectively. For the SPRs, a sensing range between 1375 nm and 2494 nm with a maximum sensitivity of 8300 nm/RIU is measured under TM polarization. Experimental results are supported with numerical simulations based on a modification of the plane-wave method for a one-dimensional multilayer waveguide.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Extending the density functional embedding theory to finite temperature and an efficient iterative method for solving for embedding potentials

A key element in the density functional embedding theory (DFET) is the embedding potential. We discuss two major issues related to the embedding potential: (1) its non-uniqueness and (2) the numerical difficulty for solving for it, especially for the spin-polarized systems. To resolve the first issue, we extend DFET to finite temperature: all quantities, such as the subsystem densities and the total system’s density, are calculated at a finite temperature. This is a physical extension since materials work at finite temperatures. We show that the embedding potential is strictly unique at T > 0. To resolve the second issue, we introduce an efficient iterative embedding potential solver. We discuss how to relax the magnetic moments in subsystems and how to equilibrate the chemical potentials across subsystems. The solver is robust and efficient for several non-trivial examples, in all of which good quality spin-polarized embedding potentials were obtained. We also demonstrate the solver on an extended periodic system: iron body-centered cubic (110) surface, which is related to the modeling of the heterogeneous catalysis involving iron, such as the Fischer-Tropsch and the Haber processes. This work would make it efficient and accurate to perform embedding simulations of some challenging material problems, such as the heterogeneous catalysis and the defects of complicated spin configurations in electronic materials.

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

Density functional theory based simulations of silicon nanowire field effect transistors

First-principles density functional theory (DFT) based, atomistic, self-consistent device simulations are performed for realistically sized Si nanowire field effect transistors (NW FETs) having tens of thousands of atoms. Through mode space transformation, DFT Hamiltonian and overlap matrices are reduced in size from a few thousands to around one hundred. Ultra-efficient quantum-mechanical transport calculations in the non-equilibrium Green's function formalism in a non-orthogonal basis are therefore made possible. The n-type and p-type Si NW FETs are simulated and found to exhibit similar device performance in the nanoscale regime.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Ab initio-informed maximum entropy modeling of rovibrational relaxation and state-specific dissociation with application to the O{sub 2} + O system

Quasi-classical trajectory (QCT) calculations are used to study state-specific ro-vibrational energy exchange and dissociation in the O{sub 2} + O system. Atom-diatom collisions with energy between 0.1 and 20 eV are calculated with a double many body expansion potential energy surface by Varandas and Pais [Mol. Phys. 65, 843 (1988)]. Inelastic collisions favor mono-quantum vibrational transitions at translational energies above 1.3 eV although multi-quantum transitions are also important. Post-collision vibrational favoring decreases first exponentially and then linearly as Δv increases. Vibrationally elastic collisions (Δv = 0) favor small ΔJ transitions while vibrationally inelastic collisions have equilibrium post-collision rotational distributions. Dissociation exhibits both vibrational and rotational favoring. New vibrational-translational (VT), vibrational-rotational-translational (VRT) energy exchange, and dissociation models are developed based on QCT observations and maximum entropy considerations. Full set of parameters for state-to-state modeling of oxygen is presented. The VT energy exchange model describes 22 000 state-to-state vibrational cross sections using 11 parameters and reproduces vibrational relaxation rates within 30% in the 2500–20 000 K temperature range. The VRT model captures 80 × 10{sup 6} state-to-state ro-vibrational cross sections using 19 parameters and reproduces vibrational relaxation rates within 60% in the 5000–15 000 K temperature range. The developed dissociation model reproduces state-specific and equilibrium dissociation rates within 25% using just 48 parameters. The maximum entropy framework makes it feasible to upscale ab initio simulation to full nonequilibrium flow calculations.

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

Energy deposition evaluation for ultra-low energy electron beam irradiation systems using calibrated thin radiochromic film and Monte Carlo simulations

For evaluation of on-site dosimetry and process design in industrial use of ultra-low energy electron beam (ULEB) processes, we evaluate the energy deposition using a thin radiochromic film and a Monte Carlo simulation. The response of film dosimeter was calibrated using a high energy electron beam with an acceleration voltage of 2 MV and alanine dosimeters with uncertainty of 11% at coverage factor 2. Using this response function, the results of absorbed dose measurements for ULEB were evaluated from 10 kGy to 100 kGy as a relative dose. The deviation between the responses of deposit energy on the films and Monte Carlo simulations was within 15%. As far as this limitation, relative dose estimation using thin film dosimeters with response function obtained by high energy electron irradiation and simulation results is effective for ULEB irradiation processes management.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2012

RADIATION PRESSURE AND MASS EJECTION IN {rho}-LIKE STATES OF GRS 1915+105

We present a unifying scenario to address the physical origin of the diversity of X-ray light curves within the {rho} variability class of the microquasar GRS 1915+105. This 'heartbeat' state is characterized by a bright flare that recurs every {approx}50-100 s, but the profile and duration of the flares vary significantly from observation to observation. Based on a comprehensive, phase-resolved study of heartbeats in the Rossi X-ray Timing Explorer archive, we demonstrate that very different X-ray light curves do not require origins in different accretion processes. Indeed, our detailed comparison of the phase-resolved spectra of a double-peaked oscillation and a single-peaked oscillation shows that different cycles can have basically similar X-ray spectral evolution. We argue that all heartbeat oscillations can be understood as the result of a combination of a thermal-viscous radiation pressure instability, a local Eddington limit in the disk, and a sudden, radiation-pressure-driven evaporation or ejection event in the inner accretion disk. This ejection appears to be a universal, fundamental part of the {rho} state, and is largely responsible for a hard X-ray pulse seen in the light curve of all cycles. We suggest that the detailed shape of oscillations in the mass accretion rate through the disk is responsible for the phenomenological differences between different {rho}-type light curves, and we discuss how future time-dependent simulations of disk instabilities may provide new insights into the role of radiation pressure in the accretion flow.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

STAR FORMATION LAWS: THE EFFECTS OF GAS CLOUD SAMPLING

Recent observational results indicate that the functional shape of the spatially resolved star formation-molecular gas density relation depends on the spatial scale considered. These results may indicate a fundamental role of sampling effects on scales that are typically only a few times larger than those of the largest molecular clouds. To investigate the impact of this effect, we construct simple models for the distribution of molecular clouds in a typical star-forming spiral galaxy and, assuming a power-law relation between star formation rate (SFR) and cloud mass, explore a range of input parameters. We confirm that the slope and the scatter of the simulated SFR-molecular gas surface density relation depend on the size of the sub-galactic region considered, due to stochastic sampling of the molecular cloud mass function, and the effect is larger for steeper relations between SFR and molecular gas. There is a general trend for all slope values to tend to {approx}unity for region sizes larger than 1-2 kpc, irrespective of the input SFR-cloud relation. The region size of 1-2 kpc corresponds to the area where the cloud mass function becomes fully sampled. We quantify the effects of selection biases in data tracing the SFR, either as thresholds (i.e., clouds smaller than a given mass value do not form stars) or as backgrounds (e.g., diffuse emission unrelated to current star formation is counted toward the SFR). Apparently discordant observational results are brought into agreement via this simple model, and the comparison of our simulations with data for a few galaxies supports a steep (>1) power-law index between SFR and molecular gas.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

EQUATION OF STATE AND NEUTRON STAR PROPERTIES CONSTRAINED BY NUCLEAR PHYSICS AND OBSERVATION

Microscopic calculations of neutron matter based on nuclear interactions derived from chiral effective field theory, combined with the recent observation of a 1.97 {+-} 0.04 M{sub Sun} neutron star, constrain the equation of state of neutron-rich matter at sub- and supranuclear densities. We discuss in detail the allowed equations of state and the impact of our results on the structure of neutron stars, the crust-core transition density, and the nuclear symmetry energy. In particular, we show that the predicted range for neutron star radii is robust. For use in astrophysical simulations, we provide detailed numerical tables for a representative set of equations of state consistent with these constraints.

79 ASTRONOMY AND ASTROPHYSICS↗