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 865 records · Page 48
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Nonplanar positron-acoustic Gardner solitary waves in electron-positron-ion plasmas with superthermal electrons and positrons

Nonplanar (cylindrical and spherical) positron-acoustic (PA) Gardner solitary waves (SWs) in an unmagnetized plasma system consisting of immobile positive ions, mobile cold positrons, and superthermal (kappa distributed) hot positrons and electrons are investigated. The modified Gardner equation is derived by using the reductive perturbation technique. The effects of cylindrical and spherical geometries, superthermal parameter of hot positrons and electrons, relative temperature ratios, and relative number density ratios on the PA Gardner SWs are studied by using the numerical simulations. The implications of our results in various space and laboratory plasma environments are briefly discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Non-linear plasma wake growth of electron holes

We report that an object's wake in a plasma with small Debye length that drifts across the magnetic field is subject to electrostatic electron instabilities. Such situations include, for example, the moon in the solar wind and probes in magnetized laboratory plasmas. The instability drive mechanism can equivalently be considered drift down the potential-energy gradient or drift up the density-gradient. The gradients arise because the plasma wake has a region of depressed density and electrostatic potential into which ions are attracted along the field. The non-linear consequences of the instability are analysed in this paper. At physical ratios of electron to ion mass, neither linear nor quasilinear treatment can explain the observation of large-amplitude perturbations that disrupt the ion streams well before they become ion-ion unstable. Here we show, however, that electron holes, once formed, continue to grow, driven by the drift mechanism, and if they remain in the wake may reach a maximum non-linearly stable size, beyond which their uncontrolled growth disrupts the ions. The hole growth calculations provide a quantitative prediction of hole profile and size evolution. Hole growth appears to explain the observations of recent particle-in-cell simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Experimental and numerical evaluation of freely spacing-tunable multiwavelength fiber laser based on two seeding light signals

A model of multiwavelength erbium-doped fiber laser (MEFL), which takes into account the impact of fiber attenuation on the four-wave-mixing (FWM), is proposed. Using this model, we numerically study the output characteristics of the MEFL based on FWM in a dispersion shift fiber with two seeding light signals (TSLS) and experimentally verify these characteristics. The numerical and experimental results show that the number of output channels can be increased with the increase of the erbium-doped fiber pump power. In addition, by decreasing the spacing of TSLS and increasing the power of TSLS, the number of output channels can be increased. However, when the power of TSLS exceeds a critical value, the number of output channels decreases. The results by numerical simulation are consistent with experimental observations from the MEFL.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Complexation behavior of oppositely charged polyelectrolytes: Effect of charge distribution

Complexation behavior of oppositely charged polyelectrolytes in a solution is investigated using a combination of computer simulations and experiments, focusing on the influence of polyelectrolyte charge distributions along the chains on the structure of the polyelectrolyte complexes. The simulations are performed using Monte Carlo with the replica-exchange algorithm for three model systems where each system is composed of a mixture of two types of oppositely charged model polyelectrolyte chains (EGEG){sub 5}/(KGKG){sub 5}, (EEGG){sub 5}/(KKGG){sub 5}, and (EEGG){sub 5}/(KGKG){sub 5}, in a solution including explicit solvent molecules. Among the three model systems, only the charge distributions along the chains are not identical. Thermodynamic quantities are calculated as a function of temperature (or ionic strength), and the microscopic structures of complexes are examined. It is found that the three systems have different transition temperatures, and form complexes with different sizes, structures, and densities at a given temperature. Complex microscopic structures with an alternating arrangement of one monolayer of E/K monomers and one monolayer of G monomers, with one bilayer of E and K monomers and one bilayer of G monomers, and with a mixture of monolayer and bilayer of E/K monomers in a box shape and a trilayer of G monomers inside the box are obtained for the three mixture systems, respectively. The experiments are carried out for three systems where each is composed of a mixture of two types of oppositely charged peptide chains. Each peptide chain is composed of Lysine (K) and glycine (G) or glutamate (E) and G, in solution, and the chain length and amino acid sequences, and hence the charge distribution, are precisely controlled, and all of them are identical with those for the corresponding model chain. The complexation behavior and complex structures are characterized through laser light scattering and atomic force microscopy measurements. The order of the apparent weight-averaged molar mass and the order of density of complexes observed from the three experimental systems are qualitatively in agreement with those predicted from the simulations.

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

Delay-induced stochastic bifurcations in a bistable system under white noise

In this paper, the effects of noise and time delay on stochastic bifurcations are investigated theoretically and numerically in a time-delayed Duffing-Van der Pol oscillator subjected to white noise. Due to the time delay, the random response is not Markovian. Thereby, approximate methods have been adopted to obtain the Fokker-Planck-Kolmogorov equation and the stationary probability density function for amplitude of the response. Based on the knowledge that stochastic bifurcation is characterized by the qualitative properties of the steady-state probability distribution, it is found that time delay and feedback intensity as well as noise intensity will induce the appearance of stochastic P-bifurcation. Besides, results demonstrated that the effects of the strength of the delayed displacement feedback on stochastic bifurcation are accompanied by the sensitive dependence on time delay. Furthermore, the results from numerical simulations best confirm the effectiveness of the theoretical analyses.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

The effect of voxel size on dose distribution in Varian Clinac iX 6 MV photon beam using Monte Carlo simulation

Recently, Monte Carlo (MC) calculation method has reported as the most accurate method of predicting dose distributions in radiotherapy. The MC code system (especially DOSXYZnrc) has been used to investigate the different voxel (volume elements) sizes effect on the accuracy of dose distributions. To investigate this effect on dosimetry parameters, calculations were made with three different voxel sizes. The effects were investigated with dose distribution calculations for seven voxel sizes: 1 × 1 × 0.1 cm{sup 3}, 1 × 1 × 0.5 cm{sup 3}, and 1 × 1 × 0.8 cm{sup 3}. The 1 × 10{sup 9} histories were simulated in order to get statistical uncertainties of 2%. This simulation takes about 9-10 hours to complete. Measurements are made with field sizes 10 × 10 cm2 for the 6 MV photon beams with Gaussian intensity distribution FWHM 0.1 cm and SSD 100.1 cm. MC simulated and measured dose distributions in a water phantom. The output of this simulation i.e. the percent depth dose and dose profile in d{sub max} from the three sets of calculations are presented and comparisons are made with the experiment data from TTSH (Tan Tock Seng Hospital, Singapore) in 0-5 cm depth. Dose that scored in voxels is a volume averaged estimate of the dose at the center of a voxel. The results in this study show that the difference between Monte Carlo simulation and experiment data depend on the voxel size both for percent depth dose (PDD) and profile dose. PDD scan on Z axis (depth) of water phantom, the big difference obtain in the voxel size 1 × 1 × 0.8 cm{sup 3} about 17%. In this study, the profile dose focused on high gradient dose area. Profile dose scan on Y axis and the big difference get in the voxel size 1 × 1 × 0.1 cm{sup 3} about 12%. This study demonstrated that the arrange voxel in Monte Carlo simulation becomes important.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2016

Study on space charge compensation in negative hydrogen ion beam

Negative hydrogen ion beam can be compensated by the trapping of ions into the beam potential. When the beam propagates through a neutral gas, these ions arise due to gas ionization by the beam ions. However, the high neutral gas pressure may cause serious negative hydrogen ion beam loss, while low neutral gas pressure may lead to ion-ion instability and decompensation. To better understand the space charge compensation processes within a negative hydrogen beam, experimental study and numerical simulation were carried out at Peking University (PKU). The simulation code for negative hydrogen ion beam is improved from a 2D particle-in-cell-Monte Carlo collision code which has been successfully applied to H{sup +} beam compensated with Ar gas. Impacts among ions, electrons, and neutral gases in negative hydrogen beam compensation processes are carefully treated. The results of the beam simulations were compared with current and emittance measurements of an H{sup −} beam from a 2.45 GHz microwave driven H{sup −} ion source in PKU. Compensation gas was injected directly into the beam transport region to modify the space charge compensation degree. The experimental results were in good agreement with the simulation results.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2015

Thermodynamic properties and transport coefficients of air thermal plasmas mixed with ablated vapors of Cu and polytetrafluoroethylene

Because the fault arc in aircraft electrical system often causes a fire, it is particularly important to analyze its energy and transfer for aircraft safety. The calculation of arc energy requires the basic parameters of the arc. This paper is mainly devoted to the calculations of equilibrium composition, thermodynamic properties (density, molar weight, enthalpy, and specific heat at constant pressure) and transport coefficients (thermal conductivity, electrical conductivity, and viscosity) of plasmas produced by a mixture of air, Cu, and polytetrafluoroethylene under the condition of local thermodynamic equilibrium. The equilibrium composition is determined by solving a system of equations around the number densities of each species. The thermodynamic properties are obtained according to the standard thermodynamic relationships. The transport coefficients are calculated using the Chapman-Enskog approximations. Results are presented in the temperature range from 3000 to 30 000 K for pressures of 0.08 and 0.1 MPa, respectively. The results are more accurate and are reliable reference data for theoretical analysis and computational simulation of the behavior of fault arc.

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

Injection locking of an electronic maser in the hard excitation mode

The phenomenon of hard excitation is natural for many electronic oscillators. In particular, in a gyrotron, a maximal efficiency is often attained in the hard excitation regime. In this paper, we study the injection-locking phenomena using two models of an electronic maser in the hard excitation mode. First, bifurcation analysis is performed for the quasilinear model described by ordinary differential equations for the slow amplitude and phase. Two main scenarios of transition to the injection-locked mode are described, which are generalizations of the well-known phase-locking and suppression mechanisms. The results obtained for the quasilinear model are confirmed by numerical simulations of a gyrotron with fixed Gaussian structure of the RF field.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Multipole electrodynamic ion trap geometries for microparticle confinement under standard ambient temperature and pressure conditions

Trapping of microparticles and aerosols is of great interest for physics and chemistry. We report microparticle trapping in case of multipole linear Paul trap geometries, operating under standard ambient temperature and pressure conditions. An 8- and 12-electrode linear trap geometries have been designed and tested with an aim to achieve trapping for larger number of particles and to study microparticle dynamical stability in electrodynamic fields. We report emergence of planar and volume ordered structures of microparticles, depending on the a.c. trapping frequency and particle specific charge ratio. The electric potential within the trap is mapped using the electrolytic tank method. Particle dynamics is simulated using a stochastic Langevin equation. We emphasize extended regions of stable trapping with respect to quadrupole traps, as well as good agreement between experiment and numerical simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Evolution of nanodot morphology on polycarbonate (PC) surfaces by 40 keV Ar{sup +}

In the present paper we have discussed the effect of 40 keV Ar{sup +} ions irradiation on nanoscale surface morphology of Polycarbonate (PC) substrate. Specimens were sputtered at off normal incidences of 30°, 40° and 50° with the fluence of 1 × 10{sup 16} Ar{sup +}cm{sup −2}. The topographical behaviour of specimens was studied by using Atomic Force Microscopy (AFM) technique. AFM study demonstrates the evolution of nano dot morphology on PC specimens on irradiating with 1 × 10{sup 16} Ar{sup +}cm{sup −2}. Average size of dots varied from 37-95 nm in this specified range of incidence while density of dots varied from 0.17-3.0 × 107 dotscm{sup −2}. Such variations in morphological features have been supported by estimation of ion range and sputtering yield through SRIM simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Mid-infrared metasurface made of composite right/left-handed transmission-line

We report on the realization of a mid-infrared metasurface based on the concept of composite right/left-handed transmission-line. The metasurface consists of a three-layer metal-insulator-metal structure patterned into transmission-lines by electron-beam lithography. Angle-variable reflection spectroscopy measurements reveal resonant absorption features corresponding to both right- and left-handed propagations in the leaky-wave guided mode region. Material loss is shown to dominate the quality factor of the left-handed modes, while the radiative loss dominates the right-handed ones. The experimental results are in good agreement with full-wave numerical simulations and are explained with an equivalent circuit model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Bond breaking in epoxy systems: A combined QM/MM approach

A novel method to combine quantum mechanics (QM) and molecular mechanics has been developed to accurately and efficiently account for covalent bond breaking in polymer systems under high strain without the use of predetermined break locations. Use of this method will provide a better fundamental understanding of the mechano-chemical origins of fracture in thermosets. Since classical force fields cannot accurately account for bond breaking, and QM is too demanding to simulate large systems, a hybrid approach is required. In the method presented here, strain is applied to the system using a classical force field, and all bond lengths are monitored. When a bond is stretched past a threshold value, a zone surrounding the bond is used in a QM energy minimization to determine which, if any, bonds break. The QM results are then used to reconstitute the system to continue the classical simulation at progressively larger strain until another QM calculation is triggered. In this way, a QM calculation is only computed when and where needed, allowing for efficient simulations. A robust QM method for energy minimization has been determined, as well as appropriate values for the QM zone size and the threshold bond length. Compute times do not differ dramatically from classical molecular mechanical simulations.

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

Atomistically derived cohesive zone model of intergranular fracture in polycrystalline graphene

Pristine single crystal graphene is the strongest known two-dimensional material, and its nonlinear anisotropic mechanical properties are well understood from the atomic length scale up to a continuum description. However, experiments indicate that grain boundaries in the polycrystalline form reduce the mechanical behavior of polycrystalline graphene. Herein, we perform atomistic-scale molecular dynamics simulations of the deformation and fracture of graphene grain boundaries and express the results as continuum cohesive zone models (CZMs) that embed notions of the grain boundary ultimate strength and fracture toughness. To facilitate energy balance, we employ a new methodology that simulates a quasi-static controlled crack propagation which renders the kinetic energy contribution to the total energy negligible. We verify good agreement between Griffith's critical energy release rate and the work of separation of the CZM, and we note that the energy of crack edges and fracture toughness differs by about 35%, which is attributed to the phenomenon of bond trapping. This justifies the implementation of the CZM within the context of the finite element method (FEM). To enhance computational efficiency in the FEM implementation, we discuss the use of scaled traction-separation laws (TSLs) for larger element sizes. As a final result, we have established that the failure characteristics of pristine graphene and high tilt angle bicrystals differ by less than 10%. This result suggests that one could use a unique or a few typical TSLs as a good approximation for the CZMs associated with the mechanical simulations of the polycrystalline graphene.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Bifurcations, chaos and adaptive backstepping sliding mode control of a power system with excitation limitation

The bifurcation and Lyapunov exponent for a single-machine-infinite bus system with excitation model are carried out by varying the mechanical power, generator damping factor and the exciter gain, from which periodic motions, chaos and the divergence of system are observed respectively. From given parameters and different initial conditions, the coexisting motions are developed in power system. The dynamic behaviors in power system may switch freely between the coexisting motions, which will bring huge security menace to protection operation. Especially, the angle divergences due to the break of stable chaotic oscillation are found which causes the instability of power system. Finally, a new adaptive backstepping sliding mode controller is designed which aims to eliminate the angle divergences and make the power system run in stable orbits. Numerical simulations are illustrated to verify the effectivity of the proposed method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Displacement and localisation of a transparent nanosphere by light-pressure forces in the field of a focused laser beam

We report the results of a numerical simulation of the Langevin equation describing the motion of a transparent nanosphere under the action of a resulting light-pressure force in the field of a continuous focused Gaussian laser beam. The conditions for the localisation of the nanosphere near the focal waist of the beam focused by the lens are determined. An analytical solution of the approximate (truncated) equation of motion is found, which almost exactly coincides with the results of the numerical simulation of the initial equation. (paper)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

FIRST INVESTIGATION OF THE COMBINED IMPACT OF IONIZING RADIATION AND MOMENTUM WINDS FROM A MASSIVE STAR ON A SELF-GRAVITATING CORE

Massive stars shape the surrounding interstellar matter (ISM) by emitting ionizing photons and ejecting material through stellar winds. To study the impact of the momentum from the wind of a massive star on the surrounding neutral or ionized material, we implemented a new HEALPix-based momentum-conserving wind scheme in the smoothed particle hydrodynamics (SPH) code SEREN. A qualitative study of the impact of the feedback from an O7.5-like star on a self-gravitating sphere shows that on its own, the transfer of momentum from a wind onto cold surrounding gas has both a compressing and dispersing effect. It mostly affects gas at low and intermediate densities. When combined with a stellar source's ionizing ultraviolet (UV) radiation, we find the momentum-driven wind to have little direct effect on the gas. We conclude that during a massive star's main sequence, the UV ionizing radiation is the main feedback mechanism shaping and compressing the cold gas. Overall, the wind's effects on the dense gas dynamics and on the triggering of star formation are very modest. The structures formed in the ionization-only simulation and in the combined feedback simulation are remarkably similar. However, in the combined feedback case, different SPH particles end up being compressed. This indicates that the microphysics of gas mixing differ between the two feedback simulations and that the winds can contribute to the localized redistribution and reshuffling of gas.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

PROBING THE ROLE OF DYNAMICAL FRICTION IN SHAPING THE BSS RADIAL DISTRIBUTION. I. SEMI-ANALYTICAL MODELS AND PRELIMINARY N-BODY SIMULATIONS

We present semi-analytical models and simplified N-body simulations with 10{sup 4} particles aimed at probing the role of dynamical friction (DF) in determining the radial distribution of blue straggler stars (BSSs) in globular clusters. The semi-analytical models show that DF (which is the only evolutionary mechanism at work) is responsible for the formation of a bimodal distribution with a dip progressively moving toward the external regions of the cluster. However, these models fail to reproduce the formation of the long-lived central peak observed in all dynamically evolved clusters. The results of N-body simulations confirm the formation of a sharp central peak, which remains as a stable feature over time regardless of the initial concentration of the system. In spite of noisy behavior, a bimodal distribution forms in many cases, with the size of the dip increasing as a function of time. In the most advanced stages, the distribution becomes monotonic. These results are in agreement with the observations. Also, the shape of the peak and the location of the minimum (which, in most of cases, is within 10 core radii) turn out to be consistent with observational results. For a more detailed and close comparison with observations, including a proper calibration of the timescales of the dynamical processes driving the evolution of the BSS spatial distribution, more realistic simulations will be necessary.

79 ASTRONOMY AND ASTROPHYSICS↗