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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 775 records · Page 43

All sources NASA NTRS DOE OSTI

DOE OSTI2018

Exponential surface melting of Cu nanoparticles observed by in-situ TEM

Highlights: • For the first time, surface melting rate of nano-scale Cu particles was investigated by in situ TEM. • The surface melting rates of different morphologies, such as platform or sphere, are different. • Different from ln[Tm/(Tm − T)] or (Tm − T){sup −r} (r is a constant), the actual surface melting of Cu nanoparticles is exponential. - Abstract: In surface melting, a liquid film forms on the crystal surface even below its melting point (T{sub m}). Currently, surface melting theory generally assumes a crystal with an infinite surface, and it is believed that the surface melting layer thickens according to ln [T{sub m}/(T{sub m} − T)] or (T{sub m} − T){sup −r} (r is a constant). In this work, we observed the melting of Cu nanoparticles with diameters of about 30 nm at 890 °C using in-situ heating high-resolution transmission electron microscopy (HRTEM). We found that the thickness of the melting layer on the surface of the Cu nanoparticles increased exponentially. The growth pattern of the melting layer was notably different from the results of previous simulations. This discovery may enhance our understanding of the melting mechanism of Cu nanomaterials.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

DOE OSTI2013

The putative liquid-liquid transition is a liquid-solid transition in atomistic models of water. II

This paper extends our earlier studies of free energy functions of density and crystalline order parameters for models of supercooled water, which allows us to examine the possibility of two distinct metastable liquid phases [D. T. Limmer and D. Chandler, J. Chem. Phys.135, 134503 (2011) and preprint http://arxiv.org/abs/arXiv:1107.0337 (2011)]. Low-temperature reversible free energy surfaces of several different atomistic models are computed: mW water, TIP4P/2005 water, Stillinger-Weber silicon, and ST2 water, the last of these comparing three different treatments of long-ranged forces. In each case, we show that there is one stable or metastable liquid phase, and there is an ice-like crystal phase. The time scales for crystallization in these systems far exceed those of structural relaxation in the supercooled metastable liquid. We show how this wide separation in time scales produces an illusion of a low-temperature liquid-liquid transition. The phenomenon suggesting metastability of two distinct liquid phases is actually coarsening of the ordered ice-like phase, which we elucidate using both analytical theory and computer simulation. For the latter, we describe robust methods for computing reversible free energy surfaces, and we consider effects of electrostatic boundary conditions. We show that sensible alterations of models and boundary conditions produce no qualitative changes in low-temperature phase behaviors of these systems, only marginal changes in equations of state. On the other hand, we show that altering sampling time scales can produce large and qualitative non-equilibrium effects. Recent reports of evidence of a liquid-liquid critical point in computer simulations of supercooled water are considered in this light.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2013

Higher-order harmonic effect of a three-dimensional helical wiggler on the Larmor rotation of the equilibrium electrons in a free-electron laser with a positive or a reversed guide magnetic field

Analytical formulas of the Larmor rotation are derived in detail for the equilibrium electrons motion in a free-electron laser with combination of a three-dimensional (3-D) helical wiggler and a positive or a reversed guide magnetic field. Generally, the Larmor radius in the configuration of a reversed guide field is much smaller than that in a positive guide field. At non-resonance, a helical orbit governed by the zero-order component of a 3-D wiggler field could hold; meanwhile, the higher-harmonic effect definitely influences those electrons with off-axis guiding centers and induces the electron-beam spreads. At resonance, the Larmor radius in the configuration of a positive guide field has a singularity with a limit tending to infinite, which causes all the electrons to hit the waveguide wall before the exit of the wiggler. Although Larmor-radius singularity does not exist in the configuration of a reversed guide field, at anti-resonance, the first-order harmonic of a 3-D wiggler field induces a transverse displacement which rapidly grows in proportion to a square of time, and leads part of the electron beam to hit the waveguide wall before reaching the wiggler exit, which depends on the specific parameters of the individual electrons. The analytical conclusions derived in the present paper are examined by the nonlinear simulations and the experimental observation. Disagreement with the previous literatures is discussed in detail.

43 PARTICLE ACCELERATORS↗

DOE OSTI2013

Analysis of local warm forming of high strength steel using near infrared ray energy

The automotive industry has been pressed to satisfy more rigorous fuel efficiency requirements to promote energy conservation, safety features and cost containment. To satisfy this need, high strength steel has been developed and used for many different vehicle parts. The use of high strength steels, however, requires careful analysis and creativity in order to accommodate its relatively high springback behavior. An innovative method, called local warm forming with near infrared ray, has been developed to help promote the use of high strength steels in sheet metal forming. For this method, local regions of the work piece are heated using infrared ray energy, thereby promoting the reduction of springback behavior. In this research, a V-bend test is conducted with DP980. After springback, the bend angles for specimens without local heating are compared to those with local heating. Numerical analysis has been performed using the commercial program, DEFORM-2D. This analysis is carried out with the purpose of understanding how changes to the local stress distribution will affect the springback during the unloading process. The results between experimental and computational approaches are evaluated to assure the accuracy of the simulation. Subsequent numerical simulation studies are performed to explore best practices with respect to thermal boundary conditions, timing, and applicability to the production environment.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2014

Development of a general method for obtaining the geometry of microfluidic networks

In the present study, a general method for geometry of fluidic networks is developed with emphasis on pressure-driven flows in the microfluidic applications. The design method is based on general features of network's geometry such as cross-sectional area and length of channels. Also, the method is applicable to various cross-sectional shapes such as circular, rectangular, triangular, and trapezoidal cross sections. Using constructal theory, the flow resistance, energy loss and performance of the network are optimized. Also, by this method, practical design strategies for the fabrication of microfluidic networks can be improved. The design method enables rapid prediction of fluid flow in the complex network of channels and is very useful for improving proper miniaturization and integration of microfluidic networks. Minimization of flow resistance of the network of channels leads to universal constants for consecutive cross-sectional areas and lengths. For a Y-shaped network, the optimal ratios of consecutive cross-section areas (A{sub i+1}/A{sub i}) and lengths (L{sub i+1}/L{sub i}) are obtained as A{sub i+1}/A{sub i} = 2{sup −2/3} and L{sub i+1}/L{sub i} = 2{sup −1/3}, respectively. It is shown that energy loss in the network is proportional to the volume of network. It is also seen when the number of channels is increased both the hydraulic resistance and the volume occupied by the network are increased in a similar manner. Furthermore, the method offers that fabrication of multi-depth and multi-width microchannels should be considered as an integral part of designing procedures. Finally, numerical simulations for the fluid flow in the network have been performed and results show very good agreement with analytic results.

36 MATERIALS SCIENCE↗

DOE OSTI2014

Magnonic band structure, complete bandgap, and collective spin wave excitation in nanoscale two-dimensional magnonic crystals

We present the observation of a complete bandgap and collective spin wave excitation in two-dimensional magnonic crystals comprised of arrays of nanoscale antidots and nanodots, respectively. Considering that the frequencies dealt with here fall in the microwave band, these findings can be used for the development of suitable magnonic metamaterials and spin wave based signal processing. We also present the application of a numerical procedure, to compute the dispersion relations of spin waves for any high symmetry direction in the first Brillouin zone. The results obtained from this procedure have been reproduced and verified by the well established plane wave method for an antidot lattice, when magnetization dynamics at antidot boundaries are pinned. The micromagnetic simulation based method can also be used to obtain iso–frequency contours of spin waves. Iso–frequency contours are analogous of the Fermi surfaces and hence, they have the potential to radicalize our understanding of spin wave dynamics. The physical origin of bands, partial and full magnonic bandgaps have been explained by plotting the spatial distribution of spin wave energy spectral density. Although, unfettered by rigid assumptions and approximations, which afflict most analytical methods used in the study of spin wave dynamics, micromagnetic simulations tend to be computationally demanding. Thus, the observation of collective spin wave excitation in the case of nanodot arrays, which can obviate the need to perform simulations, may also prove to be valuable.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2014

Charge trapping in aligned single-walled carbon nanotube arrays induced by ionizing radiation exposure

The effects of near-interfacial trapping induced by ionizing radiation exposure of aligned single-walled carbon nanotube (SWCNT) arrays are investigated via measurements of gate hysteresis in the transfer characteristics of aligned SWCNT field-effect transistors. Gate hysteresis is attributed to charge injection (i.e., trapping) from the SWCNTs into radiation-induced traps in regions near the SWCNT/dielectric interface. Self-consistent calculations of surface-potential, carrier density, and trapped charge are used to describe hysteresis as a function of ionizing radiation exposure. Hysteresis width (h) and its dependence on gate sweep range are investigated analytically. The effects of non-uniform trap energy distributions on the relationship between hysteresis, gate sweep range, and total ionizing dose are demonstrated with simulations and verified experimentally.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2014

Synchronization in node of complex networks consist of complex chaotic system

A new synchronization method is investigated for node of complex networks consists of complex chaotic system. When complex networks realize synchronization, different component of complex state variable synchronize up to different scaling complex function by a designed complex feedback controller. This paper change synchronization scaling function from real field to complex field for synchronization in node of complex networks with complex chaotic system. Synchronization in constant delay and time-varying coupling delay complex networks are investigated, respectively. Numerical simulations are provided to show the effectiveness of the proposed method.

36 MATERIALS SCIENCE↗

DOE OSTI2015

Solvation of fluoro-acetonitrile in water by 2D-IR spectroscopy: A combined experimental-computational study

The solvent dynamics around fluorinated acetonitrile is characterized by 2-dimensional infrared spectroscopy and atomistic simulations. The lineshape of the linear infrared spectrum is better captured by semiempirical (density functional tight binding) mixed quantum mechanical/molecular mechanics simulations, whereas force field simulations with multipolar interactions yield lineshapes that are significantly too narrow. For the solvent dynamics, a relatively slow time scale of 2 ps is found from the experiments and supported by the mixed quantum mechanical/molecular mechanics simulations. With multipolar force fields fitted to the available thermodynamical data, the time scale is considerably faster—on the 0.5 ps time scale. The simulations provide evidence for a well established CF–HOH hydrogen bond (population of 25%) which is found from the radial distribution function g(r) from both, force field and quantum mechanics/molecular mechanics simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2015

Heat flux viscosity in collisional magnetized plasmas

Momentum transport in collisional magnetized plasmas due to gradients in the heat flux, a “heat flux viscosity,” is demonstrated. Even though no net particle flux is associated with a heat flux, in a plasma there can still be momentum transport owing to the velocity dependence of the Coulomb collision frequency, analogous to the thermal force. This heat-flux viscosity may play an important role in numerous plasma environments, in particular, in strongly driven high-energy-density plasma, where strong heat flux can dominate over ordinary plasma flows. The heat flux viscosity can influence the dynamics of the magnetic field in plasmas through the generalized Ohm's law and may therefore play an important role as a dissipation mechanism allowing magnetic field line reconnection. The heat flux viscosity is calculated directly using the finite-difference method of Epperlein and Haines [Phys. Fluids 29, 1029 (1986)], which is shown to be more accurate than Braginskii's method [S. I. Braginskii, Rev. Plasma Phys. 1, 205 (1965)], and confirmed with one-dimensional collisional particle-in-cell simulations. The resulting transport coefficients are tabulated for ease of application.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2015

A unified heteronuclear decoupling strategy for magic-angle-spinning solid-state NMR spectroscopy

A unified strategy of two-pulse based heteronuclear decoupling for solid-state magic-angle spinning nuclear magnetic resonance is presented. The analysis presented here shows that different decoupling sequences like two-pulse phase-modulation (TPPM), X-inverse-X (XiX), and finite pulse refocused continuous wave (rCW{sup A}) are basically specific solutions of a more generalized decoupling scheme which incorporates the concept of time-modulation along with phase-modulation. A plethora of other good decoupling conditions apart from the standard, TPPM, XiX, and rCW{sup A} decoupling conditions are available from the unified decoupling approach. The importance of combined time- and phase-modulation in order to achieve the best decoupling conditions is delineated. The consequences of different indirect dipolar interactions arising from cross terms comprising of heteronuclear and homonuclear dipolar coupling terms and also those between heteronuclear dipolar coupling and chemical-shift anisotropy terms are presented in order to unfold the effects of anisotropic interactions under different decoupling conditions. Extensive numerical simulation results are corroborated with experiments on standard amino acids.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2015

Cavitation inception by the backscattering of pressure waves from a bubble interface

The secondary cavitation that occurs by the backscattering of focused ultrasound from a primary cavitation bubble caused by the negative pressure part of the ultrasound (Maxwell, et al., 2011) might be useful for the energy exchange due to bubble oscillations in High Intensity Focused Ultrasound (HIFU). The present study is concerned with the cavitation inception by the backscattering of ultrasound from a bubble. In the present experiment, a laser-induced bubble which is generated by a pulsed focused laser beam with high intensity is utilized as a primary cavitation bubble. After generating the bubble, focused ultrasound is emitted to the bubble. The acoustic field and the bubble motion are observed with a high-speed video camera. It is confirmed that the secondary cavitation bubble clouds are generated by the backscattering from the laser-induced bubble. The growth of cavitation bubble clouds is analyzed with the image processing method. The experimental results show that the height and width of the bubble clouds grow in stepwise during their evolution. The direct numerical simulations are also conducted for the backscattering of incident pressure waves from a bubble in order to evaluate a pressure field near the bubble. It is shown that the ratio of a bubble collapse time t{sub 0} to a characteristic time of wave propagation t{sub S}, η = t{sub 0}/t{sub s}, is an important determinant for generating negative pressure region by backscattering. The minimum pressure location by the backscattering in simulations is in good agreement with the experiment.

42 ENGINEERING↗

DOE OSTI2015

Complete synchronization of uncertain chaotic systems via a single proportional adaptive controller: A comparative study

This paper addresses a comparative computational study on the synchronization quality, cost and converging speed for two pairs of identical chaotic and hyperchaotic systems with unknown time-varying parameters. It is assumed that the unknown time-varying parameters are bounded. Based on the Lyapunov stability theory and using the adaptive control method, a single proportional controller is proposed to achieve the goal of complete synchronizations. Accordingly, appropriate adaptive laws are designed to identify the unknown time-varying parameters. The designed control strategy is easy to implement in practice. Numerical simulations results are provided to verify the effectiveness of the proposed synchronization scheme.

42 ENGINEERING↗

DOE OSTI2016

Particle rings and astrophysical accretion discs

Norman Rostoker had a wide range of interests and significant impact on the plasma physics research at Cornell during the time he was a Cornell professor. His interests ranged from the theory of energetic electron and ion beams and strong particle rings to the related topics of astrophysical accretion discs. We outline some of the topics related to rings and discs including the Rossby wave instability which leads to formation of anticyclonic vortices in astrophysical discs. These vorticies are regions of high pressure and act to trap dust particles which in turn may facilitate planetesimals growth in proto-planetary disks and could be important for planet formation. Analytical methods and global 3D magneto-hydrodynamic simulations have led to rapid advances in our understanding of discs in recent years.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2016

A non-additive repulsive contribution in an equation of state: The development for homonuclear square well chains equation of state validated against Monte Carlo simulation

This work consists of the adaptation of a non-additive hard sphere theory inspired by Malakhov and Volkov [Polym. Sci., Ser. A 49(6), 745–756 (2007)] to a square-well chain. Using the thermodynamic perturbation theory, an additional term is proposed that describes the effect of perturbing the chain of square well spheres by a non-additive parameter. In order to validate this development, NPT Monte Carlo simulations of thermodynamic and structural properties of the non-additive square well for a pure chain and a binary mixture of chains are performed. Good agreements are observed between the compressibility factors originating from the theory and those from molecular simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2016

Waveguide transition with vacuum window for multiband dynamic nuclear polarization systems

A low loss waveguide transition section and oversized microwave vacuum window covering several frequency bands (94 GHz, 140 GHz, 188 GHz) is presented. The transition is compact and was optimized for multiband Dynamic Nuclear Polarization (DNP) systems in a full-wave simulator. The window is more broadband than commercially available windows, which are usually optimized for single band operation. It is demonstrated that high-density polyethylene with urethane adhesive can be used as a low loss microwave vacuum window in multiband DNP systems. The overall assembly performance and dimensions are found using full-wave simulations. The practical aspects of the window implementation in the waveguide are discussed. To verify the design and simulation results, the window is tested experimentally at the three frequencies of interest.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2013

THE FRAGMENTATION OF MAGNETIZED, MASSIVE STAR-FORMING CORES WITH RADIATIVE FEEDBACK

We present a set of three-dimensional, radiation-magnetohydrodynamic calculations of the gravitational collapse of massive (300 M{sub Sun }), star-forming molecular cloud cores. We show that the combined effects of magnetic fields and radiative feedback strongly suppress core fragmentation, leading to the production of single-star systems rather than small clusters. We find that the two processes are efficient at suppressing fragmentation in different regimes, with the feedback most effective in the dense, central region and the magnetic field most effective in more diffuse, outer regions. Thus, the combination of the two is much more effective at suppressing fragmentation than either one considered in isolation. Our work suggests that typical massive cores, which have mass-to-flux ratios of about 2 relative to critical, likely form a single-star system, but that cores with weaker fields may form a small star cluster. This result helps us understand why the observed relationship between the core mass function and the stellar initial mass function holds even for {approx}100 M{sub Sun} cores with many thermal Jeans masses of material. We also demonstrate that a {approx}40 AU Keplerian disk is able to form in our simulations, despite the braking effect caused by the strong magnetic field.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2013

ANALYTICAL THEORY FOR THE INITIAL MASS FUNCTION. III. TIME DEPENDENCE AND STAR FORMATION RATE

The present paper extends our previous theory of the stellar initial mass function (IMF) by including time dependence and by including the impact of the magnetic field. The predicted mass spectra are similar to the time-independent ones with slightly shallower slopes at large masses and peak locations shifted toward smaller masses by a factor of a few. Assuming that star-forming clumps follow Larson-type relations, we obtain core mass functions in good agreement with the observationally derived IMF, in particular, when taking into account the thermodynamics of the gas. The time-dependent theory directly yields an analytical expression for the star formation rate (SFR) at cloud scales. The SFR values agree well with the observational determinations of various Galactic molecular clouds. Furthermore, we show that the SFR does not simply depend linearly on density, as is sometimes claimed in the literature, but also depends strongly on the clump mass/size, which yields the observed scatter. We stress, however, that any SFR theory depends, explicitly or implicitly, on very uncertain assumptions like clump boundaries or the mass of the most massive stars that can form in a given clump, making the final determinations uncertain by a factor of a few. Finally, we derive a fully time dependent model for the IMF by considering a clump, or a distribution of clumps accreting at a constant rate and thus whose physical properties evolve with time. In spite of its simplicity, this model reproduces reasonably well various features observed in numerical simulations of converging flows. Based on this general theory, we present a paradigm for star formation and the IMF.

79 ASTRONOMY AND ASTROPHYSICS↗