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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

A BARYONIC EFFECT ON THE MERGER TIMESCALE OF GALAXY CLUSTERS

Accurate estimation of the merger timescales of galaxy clusters is important for understanding the cluster merger process and further understanding the formation and evolution of the large-scale structure of the universe. In this paper, we explore a baryonic effect on the merger timescale of galaxy clusters by using hydrodynamical simulations. We find that the baryons play an important role in accelerating the merger process. The merger timescale decreases upon increasing the gas fraction of galaxy clusters. For example, the merger timescale is shortened by a factor of up to 3 for merging clusters with gas fractions of 0.15, compared with the timescale obtained with 0 gas fractions. The baryonic effect is significant for a wide range of merger parameters and is particularly more significant for nearly head-on mergers and high merging velocities. The baryonic effect on the merger timescale of galaxy clusters is expected to have an impact on the structure formation in the universe, such as the cluster mass function and massive substructures in galaxy clusters, and a bias of “no-gas” may exist in the results obtained from the dark matter-only cosmological simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Time delay analysis of the lensed quasar SDSS J1001+5027

We modify the algorithm we proposed in Aghamousa and Shafieloo for the time delay estimation of strongly lensed systems incorporating the weighted cross-correlation and weighted summation of correlation coefficients. We show the high performance of this algorithm by applying it to Time Delay Challenge (TDC1) simulated data. We apply then our proposed method to the light curves of the lensed quasar SDSS J1001+5027 since this system has been well studied by other groups, to compare our results with their findings. In this work we propose a new estimator, the “mirror” estimator, along with a list of criteria for reliability testing of the estimation. Our mirror estimator results are −117.1{sub −3.7}{sup +7.1} and −117.1{sub −8.8}{sup +7.2} using simple Monte Carlo simulations and simulated light curves provided by Rathna Kumar et al., respectively. Although the TDC1 simulations do not reflect the properties of the SDSS J1001+5027 light curves, using these simulations results in a smaller uncertainty, which shows that the higher quality observations can lead to a substantially more precise time delay estimation. Our time delay estimation is in agreement with the findings of the other groups for this strongly lensed system, and the difference in the size of the error bars reflects the importance of appropriate light curve simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Tuning linear response dynamics near the Dirac points in the bosonic Mott insulator

Optical lattice systems offer the possibility of creating and tuning Dirac points which are present in the tight-binding lattice dispersions. For example, such a behavior can be achieved in the staggered flux lattice or honeycomb type of lattices. Here we focus on the strongly correlated bosonic dynamics in the vicinity of Dirac points. In particular, we investigate bosonic Mott insulator phase in which quasiparticle excitations have a simple particle–hole interpretation. We show that linear response dynamics around Dirac points, can be significantly engineered at least in two ways: by the type of external perturbation or by changing the lattice properties. The key role is played by the interband transitions. Moreover, we explain that the behavior of these transitions is directly connected to different energy scales of the effective hopping amplitudes for particles and holes. Presented in this work theoretical study about tunability of linear response dynamics near the Dirac points can be directly simulated in the optical lattice systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

Numerical stability of relativistic beam multidimensional PIC simulations employing the Esirkepov algorithm

Rapidly growing numerical instabilities routinely occur in multidimensional particle-in-cell computer simulations of plasma-based particle accelerators, astrophysical phenomena, and relativistic charged particle beams. Reducing instability growth to acceptable levels has necessitated higher resolution grids, high-order field solvers, current filtering, etc. except for certain ratios of the time step to the axial cell size, for which numerical growth rates and saturation levels are reduced substantially. This paper derives and solves the cold beam dispersion relation for numerical instabilities in multidimensional, relativistic, electromagnetic particle-in-cell programs employing either the standard or the Cole–Karkkainnen finite difference field solver on a staggered mesh and the common Esirkepov current-gathering algorithm. Good overall agreement is achieved with previously reported results of the WARP code. In particular, the existence of select time steps for which instabilities are minimized is explained. Additionally, an alternative field interpolation algorithm is proposed for which instabilities are almost completely eliminated for a particular time step in ultra-relativistic simulations.

43 PARTICLE ACCELERATORS↗
DOE OSTI2017

A generalized transport-velocity formulation for smoothed particle hydrodynamics

The standard smoothed particle hydrodynamics (SPH) method suffers from tensile instability. In fluid-dynamics simulations this instability leads to particle clumping and void regions when negative pressure occurs. In solid-dynamics simulations, it results in unphysical structure fragmentation. In this work the transport-velocity formulation of Adami et al. (2013) is generalized for providing a solution of this long-standing problem. Other than imposing a global background pressure, a variable background pressure is used to modify the particle transport velocity and eliminate the tensile instability completely. Furthermore, such a modification is localized by defining a shortened smoothing length. The generalized formulation is suitable for fluid and solid materials with and without free surfaces. The results of extensive numerical tests on both fluid and solid dynamics problems indicate that the new method provides a unified approach for multi-physics SPH simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

A unidirectional subwavelength focusing near-field plate

Near-field plates consist of non-periodically patterned surfaces that can overcome the diffraction limit and confine electromagnetic fields to subwavelength dimensions. Previous near-field plates experimentally demonstrated extreme field tailoring capabilities. However, their performance suffered from radiation/reflection in undesired directions, those other than the subwavelength focus. This issue can limit the practical use of near-field plates. In this paper, we address this issue by designing a unidirectional near-field plate that can form a subwavelength focal pattern, while suppressing the field radiated/reflected in other directions. The design and operation of the proposed unidirectional near-field plate are verified through full-wave simulation. The unidirectional near-field plate may find application in high resolution imaging and probing, high density data storage, and wireless power transfer systems. As an example, its utility as a high resolution probe is demonstrated through full-wave electromagnetic simulation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

Nonlinear physics of shear Alfvén waves

Shear Alfvén waves (SAW) play fundamental roles in thermonuclear plasmas of fusion interest, since they are readily excited by energetic particles in the MeV range as well as by the thermal plasma components. Thus, understanding fluctuation induced transport in burning plasmas requires understanding nonlinear SAW physics. There exist two possible routes to nonlinear SAW physics: (i) wave-wave interactions and the resultant spectral energy transfer; (ii) nonlinear wave-particle interactions of SAW instabilities with energetic particles. Within the first route, it is advantageous to understand and describe nonlinear processes in term of proximity of the system to the Alfvénic state, where wave-wave interactions are minimized due to the cancellation of Reynolds and Maxwell stresses. Here, various wave-wave nonlinear dynamics are elucidated in terms of how they break the Alfvénic state. In particular, we discuss the qualitative and quantitative modification of the SAW parametric decay process due to finite ion compressibility and finite ion Larmor radius. We also show that toroidal geometry plays a crucial role in the nonlinear excitation of zonal structures by Alfvén eigenmodes. Within the second route, the coherent nonlinear dynamics of structures in the energetic particle phase space, by which secular resonant particle transport can occur on meso- and macro-scales, must be addressed and understood. These 'nonlinear equilibria' or 'phase-space zonal structures' dynamically evolve on characteristic (fluctuation induced) turbulent transport time scales, which are generally of the same order of the nonlinear time scale of the underlying fluctuations. In this work, we introduce the general structure of nonlinear Schrödinger equations with complex integro-differential nonlinear terms, which govern these physical processes. To elucidate all these aspects, theoretical analyses are presented together with numerical simulation results.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Optical loss due to diffraction by concentrator Fresnel lenses

Fresnel lenses are widely used in concentrating photovoltaic (CPV) systems as a primary optical element. They focus sunlight on small solar cells or on the entrance apertures of secondary optical elements. A Fresnel lens consists of several prism rings and diffraction by these prism rings is unavoidable. Some of the light that would reach a designated target area according to geometric optics will miss it due to diffraction. This diffraction loss may be of relevant magnitude for CPV applications. The results of published analytical calculations are evaluated, discussed, and compared to computer simulations and measurements.

14 SOLAR ENERGY↗
DOE OSTI2014

On the time scale associated with Monte Carlo simulations

Uniform-acceptance force-bias Monte Carlo (fbMC) methods have been shown to be a powerful technique to access longer timescales in atomistic simulations allowing, for example, phase transitions and growth. Recently, a new fbMC method, the time-stamped force-bias Monte Carlo (tfMC) method, was derived with inclusion of an estimated effective timescale; this timescale, however, does not seem able to explain some of the successes the method. In this contribution, we therefore explicitly quantify the effective timescale tfMC is able to access for a variety of systems, namely a simple single-particle, one-dimensional model system, the Lennard-Jones liquid, an adatom on the Cu(100) surface, a silicon crystal with point defects and a highly defected graphene sheet, in order to gain new insights into the mechanisms by which tfMC operates. It is found that considerable boosts, up to three orders of magnitude compared to molecular dynamics, can be achieved for solid state systems by lowering of the apparent activation barrier of occurring processes, while not requiring any system-specific input or modifications of the method. We furthermore address the pitfalls of using the method as a replacement or complement of molecular dynamics simulations, its ability to explicitly describe correct dynamics and reaction mechanisms, and the association of timescales to MC simulations in general.

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

New expression for collisionless magnetic reconnection rate

For 2D, symmetric, anti-parallel, collisionless magnetic reconnection, new expressions for the reconnection rate in the electron diffusion region are introduced. It is shown that these expressions can be derived in just a few simple steps from a physically intuitive starting point; the derivations are given in their entirety, and the validity of each step is confirmed. The predictions of these expressions are compared to the results of several long-duration, open-boundary particle-in-cell reconnection simulations to demonstrate excellent agreement.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Testing the high turbulence level breakdown of low-frequency gyrokinetics against high-frequency cyclokinetic simulations

This paper presents numerical simulations of the nonlinear cyclokinetic equations in the cyclotron harmonic representation [R. E. Waltz and Zhao Deng, Phys. Plasmas 20, 012507 (2013)]. Simulations are done with a local flux-tube geometry and with the parallel motion and variation suppressed using a newly developed rCYCLO code. Cyclokinetic simulations dynamically follow the high-frequency ion gyro-phase motion which is nonlinearly coupled into the low-frequency drift-waves possibly interrupting and suppressing gyro-averaging and increasing the transport over gyrokinetic levels. By comparing the more fundamental cyclokinetic simulations with the corresponding gyrokinetic simulations, the breakdown of gyrokinetics at high turbulence levels is quantitatively tested over a range of relative ion cyclotron frequency 10 < Ω*{sup }< 100 where Ω*{sup }= 1/ρ*, and ρ* is the relative ion gyroradius. The gyrokinetic linear mode rates closely match the cyclokinetic low-frequency rates for Ω*{sup }> 5. Gyrokinetic transport recovers cyclokinetic transport at high relative ion cyclotron frequency (Ω*{sup }≥ 50) and low turbulence level as required. Cyclokinetic transport is found to be lower than gyrokinetic transport at high turbulence levels and low-Ω* values with stable ion cyclotron (IC) modes. The gyrokinetic approximation is found to break down when the density perturbations exceed 20%. For cyclokinetic simulations with sufficiently unstable IC modes and sufficiently low Ω*{sup }∼ 10, the high-frequency component of cyclokinetic transport level can exceed the gyrokinetic transport level. However, the low-frequency component of the cyclokinetic transport and turbulence level does not exceed that of gyrokinetics. At higher and more physically relevant Ω*{sup }≥ 50 values and physically realistic IC driving rates, the low-frequency component of the cyclokinetic transport and turbulence level is still smaller than that of gyrokinetics. Thus, the cyclokinetic simulations do not account for the so-called “L-mode near edge short fall” seen in some low-frequency gyrokinetic transport and turbulence simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

High field terahertz emission from relativistic laser-driven plasma wakefields

We propose a method to generate high field terahertz (THz) radiation with peak strength of GV/cm level in the THz frequency gap range of 1–10 THz using a relativistic laser interaction with a gaseous plasma target. Due to the effect of local pump depletion, an initially Gaussian laser pulse undergoes leading edge erosion and eventually evolves to a state with leading edge being step function. Interacting with such a pulse, electrons gain transverse residual momentum and excite net transverse currents modulated by the relativistic plasma frequency. These currents give rise to the low frequency THz emission. We demonstrate this process with one and two dimensional particle-in-cell simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

LUMINOUS SATELLITES VERSUS DARK SUBHALOS: CLUSTERING IN THE MILKY WAY

The observed population of the Milky Way satellite galaxies offers a unique testing ground for galaxy formation theory on small scales. Our novel approach was to investigate the clustering of the known Milky Way satellite galaxies and to quantify the amount of substructure within their distribution using a two-point correlation function statistic in each of three spaces: configuration space, line-of-sight velocity space, and four-dimensional (4D) phase space. These results were compared to those for three sets of subhalos in the Via Lactea II cold dark matter (CDM) simulation defined to represent the luminous dwarfs. We found no evidence at a significance level above 2{sigma} of substructure within the distribution of the Milky Way satellite galaxies in any of the three spaces. The 'luminous' subhalo sets are more strongly clustered than are the Milky Way satellites in all three spaces and over a broader range of scales in 4D phase space. Each of the 'luminous' subhalo sets are clustered as a result of substructure within their line-of-sight velocity space distributions at greater than 3{sigma} significance, whereas the Milky Way satellite galaxies are randomly distributed in line-of-sight velocity space. While our comparison is with only one CDM simulation, the inconsistencies between the Milky Way satellite galaxies and the Via Lactea II subhalo sets for all clustering methods suggest a potential new 'small-scale' tension between CDM theory and the observed Milky Way satellites. Future work will obtain a more robust comparison between the observed Milky Way satellites and CDM theory by studying additional simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Magnetohydrodynamics using path or stream functions

Magnetization in highly conductive plasmas is ubiquitous to astronomical systems. Flows in such media can be described by three path functions Λ{sub α}, or, for a steady flow, by two stream functions λ{sub κ} and an additional field such as mass density ρ, velocity v, or travel time Δt. While typical analyses of a frozen magnetic field B are problem-specific and involve nonlocal gradients of the fluid element position x(t), we derive the general, local (in Λ or λ space) solution B=(∂x/∂Λ{sub α}){sub t} B-tilde {sub α}ρ/ ρ-tilde , where Lagrangian constants denoted by a tilde are directly fixed at a boundary hypersurface H-tilde on which B is known. For a steady flow, ρ-tilde B/ρ=(∂x/∂λ{sub κ}){sub Δt} B-tilde {sub κ}+v B-tilde {sub 3}/ v-tilde ; here the electric field E∼( B-tilde {sub 2}∇λ{sub 1}− B-tilde {sub 1}∇λ{sub 2})/ ρ-tilde depends only on λ{sub κ} and the boundary conditions. Illustrative special cases include compressible axisymmetric flows and incompressible flows around a sphere, showing that viscosity and compressibility enhance the magnetization and lead to thicker boundary layers. Our method is especially useful for directly computing electric fields, and for addressing upstream magnetic fields that vary in spacetime. We thus estimate the electric fields above heliospheres and magnetospheres, compute the draping of magnetic substructure around a planetary body, and demonstrate the resulting inverse polarity reversal layer. Our analysis can be immediately incorporated into existing hydrodynamic codes that are based on stream or path functions, to passively evolve the electromagnetic fields in a simulated flow. Furthermore, in such a prescription the electromagnetic fields are frozen onto the grid, so it may be developed into a fully magnetohydrodynamic, efficient simulation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Analytic properties of force-free jets in the Kerr spacetime—I

The Blandford–Znajek (BZ) mechanism describes a process for extracting the rotation energy from a spinning black hole (BH) via magnetic field lines penetrating the event horizon of the central BH. In this paper, we present a perturbation approach to study force-free jets launched by the BZ mechanism, and then discuss its two immediate applications: (1) we present a high-order split-monopole perturbation solution to the BZ mechanism which accurately pins down the energy extraction rate E-dot and well describes the structure of the BH magnetosphere for the entire range of BH spins (0 ≤ a ≤ 1); (2) the approach yields an exact constraint for the monopole field configuration in Kerr spacetime, I=Ω(1−A{sub ϕ}{sup 2}), where A{sub ϕ} is the ϕ component of the vector potential of electromagnetic field, Ω is the angular velocity of the magnetic field lines, and I is the poloidal electric current. The constraint is of particular importance to benchmark the accuracy of numerical simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

MERIDIONAL CIRCULATION DYNAMICS FROM 3D MAGNETOHYDRODYNAMIC GLOBAL SIMULATIONS OF SOLAR CONVECTION

The form of solar meridional circulation is a very important ingredient for mean field flux transport dynamo models. However, a shroud of mystery still surrounds this large-scale flow, given that its measurement using current helioseismic techniques is challenging. In this work, we use results from three-dimensional global simulations of solar convection to infer the dynamical behavior of the established meridional circulation. We make a direct comparison between the meridional circulation that arises in these simulations and the latest observations. Based on our results, we argue that there should be an equatorward flow at the base of the convection zone at mid-latitudes, below the current maximum depth helioseismic measures can probe (0.75 R{sub ⊙}). We also provide physical arguments to justify this behavior. The simulations indicate that the meridional circulation undergoes substantial changes in morphology as the magnetic cycle unfolds. We close by discussing the importance of these dynamical changes for current methods of observation which involve long averaging periods of helioseismic data. Also noteworthy is the fact that these topological changes indicate a rich interaction between magnetic fields and plasma flows, which challenges the ubiquitous kinematic approach used in the vast majority of mean field dynamo simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

University of Arizona High Energy Physics Program at the Cosmic Frontier 2014-2016

This is the final technical report from the University of Arizona High Energy Physics program at the Cosmic Frontier covering the period 2014-2016. The work aims to advance the understanding of dark energy using the Large Synoptic Survey Telescope (LSST). Progress on the engineering design of the power supplies for the LSST camera is discussed. A variety of contributions to photometric redshift measurement uncertainties were studied. The effect of the intergalactic medium on the photometric redshift of very distant galaxies was evaluated. Computer code was developed realizing the full chain of calculations needed to accurately and efficiently run large-scale simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

DISENTANGLING THE ICL WITH THE CHEFs: ABELL 2744 AS A CASE STUDY

Measurements of the intracluster light (ICL) are still prone to methodological ambiguities, and there are multiple techniques in the literature to address them, mostly based on the binding energy, the local density distribution, or the surface brightness. A common issue with these methods is the a priori assumption of a number of hypotheses on either the ICL morphology, its surface brightness level, or some properties of the brightest cluster galaxy (BCG). The discrepancy in the results is high, and numerical simulations just place a boundary on the ICL fraction in present-day galaxy clusters in the range 10%–50%. We developed a new algorithm based on the Chebyshev–Fourier functions to estimate the ICL fraction without relying on any a priori assumption about the physical or geometrical characteristics of the ICL. We are able to not only disentangle the ICL from the galactic luminosity but mark out the limits of the BCG from the ICL in a natural way. We test our technique with the recently released data of the cluster Abell 2744, observed by the Frontier Fields program. The complexity of this multiple merging cluster system and the formidable depth of these images make it a challenging test case to prove the efficiency of our algorithm. We found a final ICL fraction of 19.17 ± 2.87%, which is very consistent with numerical simulations.

79 ASTRONOMY AND ASTROPHYSICS↗