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

Spin wave spectra in perpendicularly magnetized permalloy rings

The dynamic behavior of perpendicularly magnetized permalloy circular rings is systematically investigated as a function of film thickness using broadband field modulated ferromagnetic resonance spectroscopy. We observed the splitting of one spin wave mode into a family of dense resonance peaks for the rings, which is markedly different from the single mode observed for continuous films of the same thickness. As the excitation frequency is increased, the mode family observed for the rings gradually converges into one mode. With the increase in the film thickness, a sparser spectrum of modes is observed. Our experimental results are in qualitative agreement with the dynamic micromagnetic simulations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Ultrafast defect dynamics: A new approach to all optical broadband switching employing amorphous selenium thin films

Optical switches offer higher switching speeds than electronics, however, in most cases utilizing the interband transitions of the active medium for switching. As a result, the signal suffers heavy losses. In this article, we demonstrate a simple and yet efficient ultrafast broadband all-optical switching on ps timescale in the sub-bandgap region of the a-Se thin film, where the intrinsic absorption is very weak. The optical switching is attributed to short-lived transient defects that form localized states in the bandgap and possess a large electron-phonon coupling. We model these processes through first principles simulation that are in agreement with the experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Enhancement of the droplet nucleation in a dense supersaturated Lennard-Jones vapor

The vapor–liquid nucleation in a dense Lennard-Jones system is studied analytically and numerically. A solution of the nucleation kinetic equations, which includes the elementary processes of condensation/evaporation involving the lightest clusters, is obtained, and the nucleation rate is calculated. Based on the equation of state for the cluster vapor, the pre-exponential factor is obtained. The latter diverges as a spinodal is reached, which results in the nucleation enhancement. The work of critical cluster formation is calculated using the previously developed two-parameter model (TPM) of small clusters. A simple expression for the nucleation rate is deduced and it is shown that the work of cluster formation is reduced for a dense vapor. This results in the nucleation enhancement as well. To verify the TPM, a simulation is performed that mimics a steady-state nucleation experiments in the thermal diffusion cloud chamber. The nucleating vapor with and without a carrier gas is simulated using two different thermostats for the monomers and clusters. The TPM proves to match the simulation results of this work and of other studies.

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

Validation metrics for turbulent plasma transport

Developing accurate models of plasma dynamics is essential for confident predictive modeling of current and future fusion devices. In modern computer science and engineering, formal verification and validation processes are used to assess model accuracy and establish confidence in the predictive capabilities of a given model. This paper provides an overview of the key guiding principles and best practices for the development of validation metrics, illustrated using examples from investigations of turbulent transport in magnetically confined plasmas. Particular emphasis is given to the importance of uncertainty quantification and its inclusion within the metrics, and the need for utilizing synthetic diagnostics to enable quantitatively meaningful comparisons between simulation and experiment. As a starting point, the structure of commonly used global transport model metrics and their limitations is reviewed. An alternate approach is then presented, which focuses upon comparisons of predicted local fluxes, fluctuations, and equilibrium gradients against observation. The utility of metrics based upon these comparisons is demonstrated by applying them to gyrokinetic predictions of turbulent transport in a variety of discharges performed on the DIII-D tokamak [J. L. Luxon, Nucl. Fusion 42, 614 (2002)], as part of a multi-year transport model validation activity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Features in simulation of crystal growth using the hyperbolic PFC equation and the dependence of the numerical solution on the parameters of the computational grid

We investigate the three-dimensional mathematical model of crystal growth called PFC (Phase Field Crystal) in a hyperbolic modification. This model is also called the modified model PFC (originally PFC model is formulated in parabolic form) and allows to describe both slow and rapid crystallization processes on atomic length scales and on diffusive time scales. Modified PFC model is described by the differential equation in partial derivatives of the sixth order in space and second order in time. The solution of this equation is possible only by numerical methods. Previously, authors created the software package for the solution of the Phase Field Crystal problem, based on the method of isogeometric analysis (IGA) and PetIGA program library. During further investigation it was found that the quality of the solution can strongly depends on the discretization parameters of a numerical method. In this report, we show the features that should be taken into account during constructing the computational grid for the numerical simulation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

Nature of the wiggle instability of galactic spiral shocks

Gas in disk galaxies interacts nonlinearly with an underlying stellar spiral potential to form galactic spiral shocks. While numerical simulations typically show that spiral shocks are unstable to wiggle instability (WI) even in the absence of magnetic fields and self-gravity, its physical nature has remained uncertain. To clarify the mechanism behind the WI, we conduct a normal-mode linear stability analysis and nonlinear simulations assuming that the disk is isothermal and infinitesimally thin. We find that the WI is physical, originating from the generation of potential vorticity at a deformed shock front, rather than Kelvin-Helmholtz instabilities as previously thought. Since gas in galaxy rotation periodically passes through the shocks multiple times, the potential vorticity can accumulate successively, setting up a normal mode that grows exponentially with time. Eigenfunctions of the WI decay exponentially downstream from the shock front. Both shock compression of acoustic waves and a discontinuity of shear across the shock stabilize the WI. The wavelength and growth time of the WI depend on the arm strength quite sensitively. When the stellar-arm forcing is moderate at 5%, the wavelength of the most unstable mode is about 0.07 times the arm-to-arm spacing, with the growth rate comparable to the orbital angular frequency, which is found to be in good agreement with the results of numerical simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

PLASMA EMISSION BY NONLINEAR ELECTROMAGNETIC PROCESSES

The plasma emission, or electromagnetic (EM) radiation at the plasma frequency and/or its harmonic(s), is generally accepted as the radiation mechanism responsible for solar type II and III radio bursts. Identification and characterization of these solar radio burst phenomena were done in the 1950s. Despite many decades of theoretical research since then, a rigorous demonstration of the plasma emission process based upon first principles was not available until recently, when, in a recent Letter, Ziebell et al. reported the first complete numerical solution of EM weak turbulence equations; thus, quantitatively analyzing the plasma emission process starting from the initial electron beam and the associated beam-plasma (or Langmuir wave) instability, as well as the subsequent nonlinear conversion of electrostatic Langmuir turbulence into EM radiation. In the present paper, the same problem is revisited in order to elucidate the detailed physical mechanisms that could not be reported in the brief Letter format. Findings from the present paper may be useful for interpreting observations and full-particle numerical simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

I-ball formation with logarithmic potential

A coherently oscillating real scalar field with potential shallower than quadratic one fragments into spherical objects called I-balls. We study the I-ball formation for logarithmic potential which appears in many cosmological models. We perform lattice simulations and find that the I-balls are formed when the potential becomes dominated by the quadratic term. Furthermore, we estimate the I-ball profile assuming that the adiabatic invariant is conserved during formation and obtain the result that agrees to the numerical simulations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2014

THE BLUE STRAGGLER STAR POPULATION IN NGC 1261: EVIDENCE FOR A POST-CORE-COLLAPSE BOUNCE STATE

We present a multi-passband photometric study of the Blue Straggler Star (BSS) population in the Galactic globular cluster (GC) NGC 1261, using available space- and ground-based survey data. The inner BSS population is found to have two distinct sequences in the color-magnitude diagram (CMD), similar to double BSS sequences detected in other GCs. These well defined sequences are presumably linked to single short-lived events such as core collapse, which are expected to boost the formation of BSSs. In agreement with this, we find a BSS sequence in NGC 1261 which can be well reproduced individually by a theoretical model prediction of a 2 Gyr old population of stellar collision products, which are expected to form in the denser inner regions during short-lived core contraction phases. Additionally, we report the occurrence of a group of BSSs with unusually blue colors in the CMD, which are consistent with a corresponding model of a 200 Myr old population of stellar collision products. The properties of the NGC 1261 BSS populations, including their spatial distributions, suggest an advanced dynamical evolutionary state of the cluster, but the core of this GC does not show the classical signatures of core collapse. We argue that these apparent contradictions provide evidence for a post-core-collapse bounce state seen in dynamical simulations of old GCs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Propagation of a Self-Similar Pulse in a Medium with One-Photon Transient Absorption

The shape of an optical self-similar pulse that propagates in a medium of finite length with noninstantaneous nonlinear absorption is found analytically. The pulse has an asymmetric shape and, in addition, frequency modulation, which confirms that it belongs to a new type of self-similar propagation modes of laser radiation—chirped pulses with a self-similar shape. The existence of pulses of this kind is also confirmed by results of computer simulation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Processes of diffusion in a limited ensemble of charged particles in a static magnetic field

The dynamics of a limited ensemble of charged particles in a spatially uniform static magnetic field has been studied analytically and numerically. Calculations have been performed for particles with various masses and charges in a wide range of the parameters of the systems under study. It has been shown that the transverse diffusion coefficient in an ensemble of particles in the magnetic field for weakly nonideal systems is described by the Townsend formula. Analytical estimates have been proposed for the analysis of the effective diffusion radius of the ensemble of particles owing to their thermal motion. The proposed relations have been tested by the numerical simulation of the problem.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Effect of Magnetic Field on the Dynamics of Movement of Charged Particles in Cluster Systems

The results of numerical and analytical studies of the influence of a constant magnetic field on the dynamics of the thermal motion of charged particles in an isotropic electrostatic trap are presented. Simulations were performed for cluster systems of particles with Coulomb interaction in a wide range of their parameters. The analysis of the spectral density of displacement of particles and mass transfer processes in simulated ensembles is performed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

X-pinch. Part II

Results of experimental studies of the X-pinch since its invention and implementation in 1982 at the Lebedev Physical Institute are presented. The review consists of two parts. The first part briefly outlines the history of creation and studies of X-pinches, describes the diagnostic techniques and devices developed during these studies, and presents the main results obtained in studying the physical processes occurring in the X-pinch. The second part is devoted to the results of detailed studies of the spatial, temporal, and spectral characteristics of the X-pinch hot spot—the region where the highest plasma parameters are achieved and which is a source of X-ray emission with extreme parameters. Some results of X-pinch simulations are also presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

The sizes of globular clusters as tracers of galactic halo potentials

We present N-body simulations of globular clusters, exploring the effect of different galactic potentials on cluster sizes, r{sub h}. For various galactocentric distances, R{sub G}, we assess how cluster sizes change when we vary the virial mass and concentration of the host galaxy’s dark-matter halo. We show that sizes of GCs are determined by the local galactic mass density rather than the virial mass of the host galaxy. We find that clusters evolving in the inner halos of less concentrated galaxies are significantly more extended than those evolving in more concentrated ones, while the sizes of those orbiting in the outer halo are almost independent of concentration. Adding a baryonic component to our galaxy models does not change these results much, since its effect is only significant in the very inner halo. Our simulations suggest that there is a relation between r{sub h} and R{sub G}, which systematically depends on the physical parameters of the halo. Hence, observing such relations in individual galaxies can put a new observational constraint on dark-matter halo characteristics. However, by varying the halo mass in a wide range of 10{sup 9}⩽M{sub vir}/M{sub ⊙}⩽10{sup 13}, we find that the r{sub h} − R{sub G} relationship will be nearly independent of halo mass, if one assumes M{sub vir} and c{sub vir} as two correlated parameters, as is suggested by cosmological simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Magnetoacoustic Waves in a Stratified Atmosphere with a Magnetic Null Point

We perform nonlinear MHD simulations to study the propagation of magnetoacoustic waves from the photosphere to the low corona. We focus on a 2D system with a gravitationally stratified atmosphere and three photospheric concentrations of magnetic flux that produce a magnetic null point with a magnetic dome topology. We find that a single wavepacket introduced at the lower boundary splits into multiple secondary wavepackets. A portion of the packet refracts toward the null owing to the varying Alfvén speed. Waves incident on the equipartition contour surrounding the null, where the sound and Alfvén speeds coincide, partially transmit, reflect, and mode-convert between branches of the local dispersion relation. Approximately 15.5% of the wavepacket’s initial energy ( E {sub input}) converges on the null, mostly as a fast magnetoacoustic wave. Conversion is very efficient: 70% of the energy incident on the null is converted to slow modes propagating away from the null, 7% leaves as a fast wave, and the remaining 23% (0.036 E {sub input}) is locally dissipated. The acoustic energy leaving the null is strongly concentrated along field lines near each of the null’s four separatrices. The portion of the wavepacket that refracts toward the null, and the amount of current accumulation, depends on the vertical and horizontal wavenumbers and the centroid position of the wavepacket as it crosses the photosphere. Regions that refract toward or away from the null do not simply coincide with regions of open versus closed magnetic field or regions of particular field orientation. We also model wavepacket propagation using a WKB method and find that it agrees qualitatively, though not quantitatively, with the results of the numerical simulation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

The Dynamics of Inhomogeneous Dark Energy

In this paper, by analyzing the dynamics of inhomogeneous quintessence dark energy, we find that the gradient energy of dark energy will oscillate and gradually vanish, which indicates the gradient energy of the scalar field present in the early universe does not affect the current dynamics of the universe. Moreover, with the decaying of gradient energy, there exists a possible mutual transformation between kinetic energy and gradient energy. In the framework of interacting dark energy models, we argue that inhomogeneous dark energy may have a significant effect on the evolution of the cosmic background, the investigation of which still requires fully relativistic N -body numerical simulations in the future.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Magnetic Flux Rope Identification and Characterization from Observationally Driven Solar Coronal Models

Formed through magnetic field shearing and reconnection in the solar corona, magnetic flux ropes are structures of twisted magnetic field, threaded along an axis. Their evolution and potential eruption are of great importance for space weather. Here we describe a new methodology for the automated detection of flux ropes in simulated magnetic fields, utilizing field-line helicity. Our Flux Rope Detection and Organization (FRoDO) code, which measures the magnetic flux and helicity content of pre-erupting flux ropes over time, as well as detecting eruptions, is publicly available. As a first demonstration, the code is applied to the output from a time-dependent magnetofrictional model, spanning 1996 June 15–2014 February 10. Over this period, 1561 erupting and 2099 non-erupting magnetic flux ropes are detected, tracked, and characterized. For this particular model data, erupting flux ropes have a mean net helicity magnitude of 2.66×10{sup 43} Mx{sup 2}, while non-erupting flux ropes have a significantly lower mean of 4.04×10{sup 42} Mx{sup 2}, although there is overlap between the two distributions. Similarly, the mean unsigned magnetic flux for erupting flux ropes is 4.04×10{sup 21} Mx, significantly higher than the mean value of 7.05×10{sup 20} Mx for non-erupting ropes. These values for erupting flux ropes are within the broad range expected from observational and theoretical estimates, although the eruption rate in this particular model is lower than that of observed coronal mass ejections. In the future, the FRoDO code will prove to be a valuable tool for assessing the performance of different non-potential coronal simulations and comparing them with observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

9-D polarized proton transport in the MEIC figure 8 collider ring - first steps

Spin tracking studies in the MEIC figure-8 collider ion ring are presented, based on a very preliminary design of the lattice. They provide numerical illustrations of some of the aspects of the figure-8 concept, including spin-rotator based spin control, and lay out the path towards a complete spin tracking simulation of a figure-8 ring.

43 PARTICLE ACCELERATORS↗