Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computerized simulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,657 records · Page 92
All sourcesNASA NTRSDOE OSTI
DOE OSTI2012

AGE AND MASS FOR 920 LARGE MAGELLANIC CLOUD CLUSTERS DERIVED FROM 100 MILLION MONTE CARLO SIMULATIONS

We present new age and mass estimates for 920 stellar clusters in the Large Magellanic Cloud (LMC) based on previously published broadband photometry and the stellar cluster analysis package, MASSCLEANage. Expressed in the generic fitting formula, d{sup 2} N/dMdt{proportional_to}M{sup {alpha}} t{sup {beta}}, the distribution of observed clusters is described by {alpha} = -1.5 to -1.6 and {beta} = -2.1 to -2.2. For 288 of these clusters, ages have recently been determined based on stellar photometric color-magnitude diagrams, allowing us to gauge the confidence of our ages. The results look very promising, opening up the possibility that this sample of 920 clusters, with reliable and consistent age, mass, and photometric measures, might be used to constrain important characteristics about the stellar cluster population in the LMC. We also investigate a traditional age determination method that uses a {chi}{sup 2} minimization routine to fit observed cluster colors to standard infinite-mass limit simple stellar population models. This reveals serious defects in the derived cluster age distribution using this method. The traditional {chi}{sup 2} minimization method, due to the variation of U, B, V, R colors, will always produce an overdensity of younger and older clusters, with an underdensity of clusters in the log (age/yr) = [7.0, 7.5] range. Finally, we present a unique simulation aimed at illustrating and constraining the fading limit in observed cluster distributions that includes the complex effects of stochastic variations in the observed properties of stellar clusters.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

THREE-DIMENSIONAL NUMERICAL SIMULATIONS OF FAST-TO-ALFVEN CONVERSION IN SUNSPOTS

The conversion of fast waves to the Alfven mode in a realistic sunspot atmosphere is studied through three-dimensional numerical simulations. An upward propagating fast acoustic wave is excited in the high-{beta} region of the model. The new wave modes generated at the conversion layer are analyzed from the projections of the velocity and magnetic field in their characteristic directions, and the computation of their wave energy and fluxes. The analysis reveals that the maximum efficiency of the conversion to the slow mode is obtained for inclinations of 25 Degree-Sign and low azimuths, while the Alfven wave conversions peak at high inclinations and azimuths between 50 Degree-Sign and 120 Degree-Sign . Downward propagating Alfven waves appear at the regions of the sunspot where the orientation of the magnetic field is in the direction opposite to the wave propagation, since at these locations the Alfven wave couples better with the downgoing fast magnetic wave which is reflected due to the gradients of the Alfven speed. The simulations show that the Alfven energy at the chromosphere is comparable to the acoustic energy of the slow mode, being even higher at high inclined magnetic fields.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

A MATHEMATICAL MODEL FOR AN HOURGLASS MAGNETIC FIELD

Starting with a mathematical boundary value problem for the magnetic vector potential in an axisymmetric cylindrical coordinate system, we derive a general solution for any arbitrary current distribution using the method of Green's functions. We use this to derive an analytic form for an hourglass magnetic field pattern created by electrical currents that are concentrated near (but not confined within) the equatorial plane of a cylindrical coordinate system. Our solution is not characterized by a cusp at the equatorial plane, as in previous solutions based on a current sheet. The pattern we derive provides a very good fit to hourglass magnetic field patterns emerging from three-dimensional numerical simulations of core formation, and can in principle be used for source-fitting of observed magnetic hourglass patterns.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

JOVIAN STRATOSPHERE AS A CHEMICAL TRANSPORT SYSTEM: BENCHMARK ANALYTICAL SOLUTIONS

We systematically investigated the solvable analytical benchmark cases in both one- and two-dimensional (1D and 2D) chemical-advective-diffusive systems. We use the stratosphere of Jupiter as an example but the results can be applied to other planetary atmospheres and exoplanetary atmospheres. In the 1D system, we show that CH{sub 4} and C{sub 2}H{sub 6} are mainly in diffusive equilibrium, and the C{sub 2}H{sub 2} profile can be approximated by modified Bessel functions. In the 2D system in the meridional plane, analytical solutions for two typical circulation patterns are derived. Simple tracer transport modeling demonstrates that the distribution of a short-lived species (such as C{sub 2}H{sub 2}) is dominated by the local chemical sources and sinks, while that of a long-lived species (such as C{sub 2}H{sub 6}) is significantly influenced by the circulation pattern. We find that an equator-to-pole circulation could qualitatively explain the Cassini observations, but a pure diffusive transport process could not. For slowly rotating planets like the close-in extrasolar planets, the interaction between the advection by the zonal wind and chemistry might cause a phase lag between the final tracer distribution and the original source distribution. The numerical simulation results from the 2D Caltech/JPL chemistry-transport model agree well with the analytical solutions for various cases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

THE SPACE MOTION OF LEO I: THE MASS OF THE MILKY WAY'S DARK MATTER HALO

We combine our Hubble Space Telescope measurement of the proper motion of the Leo I dwarf spheroidal galaxy (presented in a companion paper) with the highest resolution numerical simulations of Galaxy-size dark matter halos in existence to constrain the mass of the Milky Way's dark matter halo (M{sub vir,MW}). Despite Leo I's large Galactocentric space velocity (200 km s{sup -1}) and distance (261 kpc), we show that it is extremely unlikely to be unbound if Galactic satellites are associated with dark matter substructure, as 99.9% of subhalos in the simulations are bound to their host. The observed position and velocity of Leo I strongly disfavor a low-mass Milky Way: if we assume that Leo I is the least bound of the Milky Way's classical satellites, then we find that M{sub vir,MW} > 10{sup 12} M{sub Sun} at 95% confidence for a variety of Bayesian priors on M{sub vir,MW}. In lower mass halos, it is vanishingly rare to find subhalos at 261 kpc moving as fast as Leo I. Should an additional classical satellite be found to be less bound than Leo I, this lower limit on M{sub vir,MW} would increase by 30%. Imposing a mass-weighted {Lambda}CDM prior, we find a median Milky Way virial mass of M{sub vir,MW} = 1.6 Multiplication-Sign 10{sup 12} M{sub Sun }, with a 90% confidence interval of [1.0-2.4] Multiplication-Sign 10{sup 12} M{sub Sun }. We also confirm a strong correlation between subhalo infall time and orbital energy in the simulations and show that proper motions can aid significantly in interpreting the infall times and orbital histories of satellites.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

TOWARD A THEORY OF ASTROPHYSICAL PLASMA TURBULENCE AT SUBPROTON SCALES

We present an analytical study of subproton electromagnetic fluctuations in a collisionless plasma with a plasma beta of the order of unity. In the linear limit, a rigorous derivation from the kinetic equation is conducted focusing on the role and physical properties of kinetic-Alfvén and whistler waves. Then, nonlinear fluid-like equations for kinetic-Alfvén waves and whistler modes are derived, with special emphasis on the similarities and differences in the corresponding plasma dynamics. The kinetic-Alfvén modes exist in the lower-frequency region of phase space, ω << k v{sub Ti} , where they are described by the kinetic-Alfvén system. These modes exist both below and above the ion-cyclotron frequency. The whistler modes, which are qualitatively different from the kinetic-Alfvén modes, occupy a different region of phase space, k v{sub Ti} << ω << k{sub z}v{sub Te} , and they are described by the electron magnetohydrodynamics (MHD) system or the reduced electron MHD system if the propagation is oblique. Here, k{sub z} and k are the wavenumbers along and transverse to the background magnetic field, respectively, and v{sub Ti} and v{sub Te} are the ion and electron thermal velocities, respectively. The models of subproton plasma turbulence are discussed and the results of numerical simulations are presented. We also point out possible implications for solar-wind observations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

MAGNETIC FIELD COMPONENTS ANALYSIS OF THE SCUPOL 850 μm POLARIZATION DATA CATALOG

We present an extensive analysis of the 850 μm polarization maps of the SCUBA Polarimeter Legacy (SCUPOL) Catalogue produced by Matthews et al., focusing exclusively on the molecular clouds and star-forming regions. For the sufficiently sampled regions, we characterize the depolarization properties and the turbulent-to-mean magnetic field ratio of each region. Similar sets of parameters are calculated from two-dimensional synthetic maps of dust-emission polarization produced with three-dimensional magnetohydrodynamics (MHD) numerical simulations scaled to the S106, OMC-2/3, W49, and DR21 molecular cloud polarization maps. For these specific regions, the turbulent MHD regimes retrieved from the simulations, as described by the turbulent Alfvén and Sonic Mach numbers, are consistent within a factor one to two with the values of the same turbulent regimes estimated from the analysis of Zeeman measurements data provided by Crutcher. Constraints on the values of the inclination angle α of the mean magnetic field with respect to the line of sight are also given. The values obtained from the comparison of the simulations with the SCUPOL data are consistent with the estimates made by using two observational methods provided by other authors. Our main conclusion is that simple, ideal, isothermal, and non-self-gravitating MHD simulations are sufficient in order to describe the large-scale observed physical properties of the envelopes of this set of regions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Small-scale gradients of charged particles in the heliospheric magnetic field

Using numerical simulations of charged-particles propagating in the heliospheric magnetic field, we study small-scale gradients, or 'dropouts,' in the intensity of solar energetic particles seen at 1 AU. We use two turbulence models, the foot-point random motion model and the two-component model, to generate fluctuating magnetic fields similar to spacecraft observations at 1 AU. The turbulence models include a Kolmogorov-like magnetic field power spectrum containing a broad range of spatial scales from those that lead to large-scale field-line random walk to small scales leading to resonant pitch-angle scattering of energetic particles. We release energetic protons (20 keV-10 MeV) from a spatially compact and instantaneous source. The trajectories of energetic charged particles in turbulent magnetic fields are numerically integrated. Spacecraft observations are mimicked by collecting particles in small windows when they pass the windows at a distance of 1 AU. We show that small-scale gradients in the intensity of energetic particles and velocity dispersions observed by spacecraft can be reproduced using the foot-point random motion model. However, no dropouts are seen in simulations using the two-component magnetic turbulence model. We also show that particle scattering in the solar wind magnetic field needs to be infrequent for intensity dropouts to form.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

KINEMATICS OF SUPERBUBBLES AND SUPERSHELLS IN THE IRREGULAR GALAXY, NGC 1569

We present observations in the optical lines of Hα and [S II] (λλ6717, 6731 Å) and in X-rays of the irregular galaxy, NGC 1569. The observations in Hα and [S II] were made with the UNAM scanning Fabry-Perot interferometer (PUMA) and the X-ray data were obtained from the Chandra data archive. We detected several superbubbles, filaments, and supershells in NGC 1569 for which we determined size as well as their kinematic properties. We present a catalog of expansion velocities of 12 superbubbles, listing their positions, diameters, and physical parameters. Likewise, we present a catalog of 15 filaments and 4 supershells. In order to identify possible X-ray emission from the superbubbles in this galaxy, we analyzed the X-ray emission of NGC 1569 in two energy bands: 0.2-2.0 keV (soft X-rays) and 2.0-8.0 keV (hard X-rays). Based on X-ray images, we detected X-ray emission that could possibly be related to some of the superbubbles. The spectrum of the X-ray superbubbles can be described by an optically thin thermal plasma model. In order to identify the possible coexistence of galactic super winds and superbubbles we have performed adiabatic three-dimensional N-body/smoothed particle hydrodynamics simulations to follow the evolution of the most important stellar clusters in this galaxy, SSC A and SSC B, using the GADGET-2 code. Those simulations demonstrate that depending on the specific initial conditions, the formation of superbubbles or a galactic superwind can result in NGC 1569.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Quantifying (dis)agreement between direct detection experiments in a halo-independent way

We propose an improved method to study recent and near-future dark matter direct detection experiments with small numbers of observed events. Our method determines in a quantitative and halo-independent way whether the experiments point towards a consistent dark matter signal and identifies the best-fit dark matter parameters. To achieve true halo independence, we apply a recently developed method based on finding the velocity distribution that best describes a given set of data. For a quantitative global analysis we construct a likelihood function suitable for small numbers of events, which allows us to determine the best-fit particle physics properties of dark matter considering all experiments simultaneously. Based on this likelihood function we propose a new test statistic that quantifies how well the proposed model fits the data and how large the tension between different direct detection experiments is. We perform Monte Carlo simulations in order to determine the probability distribution function of this test statistic and to calculate the p-value for both the dark matter hypothesis and the background-only hypothesis.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Zoomed high-resolution simulations of Multi-coupled Dark Energy: cored galaxy density profiles at high redshift

We perform for the first time high-resolution zoom-in re-simulations of individual halos in the context of the Multi-coupled Dark Energy (McDE) scenario, which is characterised by the existence of two distinct Dark Matter particle species with opposite couplings to a Dark Energy scalar field. We compare the structural properties of the simulated halos to the standard ΛCDM results. The zoomed-in initial conditions are set up using a specifically designed code called ZInCo that we publicly release along with the present paper. Our numerical results allow to investigate in detail and with unprecedented resolution the halo segregation process that characterises McDE cosmologies from its very early stages. In particular, we find that in contrast to what could be inferred from previous numerical analysis at lower resolution, the segregation process is already in place at redshifts as high as z ∼ 7. Most remarkably, we find that the subsequent evolution of the segregation leads to the formation of cored total matter density profiles with a core size that progressively increases in time. The shape of the cored profiles can be accurately predicted as the superposition of two NFW profiles with an increasing offset, thereby confirming the interpretation of the simulations results in terms of the segregation of the two dark matter components of the halo as a consequence of their different coupling to the Dark Energy field.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2011

RESIDUAL ENERGY IN MAGNETOHYDRODYNAMIC TURBULENCE

There is mounting evidence in solar wind observations and in numerical simulations that kinetic and magnetic energies are not in equipartition in magnetohydrodynamic (MHD) turbulence. The origin of their mismatch, the residual energy E{sub r} = E{sub v} - E{sub b} , is not understood well. In the present work this effect is studied analytically in the regime of weak MHD turbulence. We find that residual energy is spontaneously generated by turbulent dynamics, and it has a negative sign, in good agreement with the observations. We find that the residual energy condenses around k {sub ||} = 0 with its k {sub ||}-spectrum broadening linearly with k{sub perpendicular}, where k{sub ||} and k{sub perpendicular} are the wavenumbers parallel and perpendicular to the background magnetic field, and the field-perpendicular spectrum of the residual energy has the scaling E{sub r} (k{sub perpendicular}){proportional_to}k {sup -1}{sub perpendicular} in the inertial interval. These results are found to be in agreement with numerical simulations. We propose that residual energy plays a fundamental role in Alfvenic turbulence and it should be taken into account for the correct interpretation of observational and numerical data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

A ROLE OF MAGNETOSONIC PULSES ON VARIATIONS OF VOYAGER-1 MeV ELECTRON INTENSITY IN THE HELIOSHEATH

Voyager 1 (V1) spacecraft observed electrons of 6-14 MeV in the heliosheath which showed several enhancements of significant flux variation. We compare these temporal electron flux variations, from the time when V1 crossed the termination shock (TS) up to mid-2008, with dynamical phenomena in the heliosheath that are obtained from our MHD simulations which are based on Voyager 2 (V2) observed solar-wind data. Our simulations indicate that all electron flux enhancements, except for one, correspond fairly well to the times when a magnetosonic (MS) pulse was driven downstream of the TS due to collision of interplanetary shock (IPS) or shock-driven MS pulse and its reflection in the heliosheath that either passed V1, or collided with the TS or with the plasma sheet near the heliopause (HP). This result suggests that these enhancements in the electron flux should correspond to either direct or indirect effects of MS pulses in the heliosheath driven by IPSs. The scale of the heliosphere is estimated by comparing V1-observed magnetic field intensity with the simulated intensity which suggests that V1 is possibly located near the HP within 4-8 AU at the present time.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Web-based, GPU-accelerated, Monte Carlo simulation and visualization of indirect radiation imaging detector performance

Purpose: Monte Carlo simulations play a vital role in the understanding of the fundamental limitations, design, and optimization of existing and emerging medical imaging systems. Efforts in this area have resulted in the development of a wide variety of open-source software packages. One such package, hybridMANTIS, uses a novel hybrid concept to model indirect scintillator detectors by balancing the computational load using dual CPU and graphics processing unit (GPU) processors, obtaining computational efficiency with reasonable accuracy. In this work, the authors describe two open-source visualization interfaces, webMANTIS and visualMANTIS to facilitate the setup of computational experiments via hybridMANTIS. Methods: The visualization tools visualMANTIS and webMANTIS enable the user to control simulation properties through a user interface. In the case of webMANTIS, control via a web browser allows access through mobile devices such as smartphones or tablets. webMANTIS acts as a server back-end and communicates with an NVIDIA GPU computing cluster that can support multiuser environments where users can execute different experiments in parallel. Results: The output consists of point response and pulse-height spectrum, and optical transport statistics generated by hybridMANTIS. The users can download the output images and statistics through a zip file for future reference. In addition, webMANTIS provides a visualization window that displays a few selected optical photon path as they get transported through the detector columns and allows the user to trace the history of the optical photons. Conclusions: The visualization tools visualMANTIS and webMANTIS provide features such as on the fly generation of pulse-height spectra and response functions for microcolumnar x-ray imagers while allowing users to save simulation parameters and results from prior experiments. The graphical interfaces simplify the simulation setup and allow the user to go directly from specifying input parameters to receiving visual feedback for the model predictions.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2015

Monte Carlo simulation of PET and SPECT imaging of {sup 90}Y

Purpose: Yittrium-90 ({sup 90}Y) is traditionally thought of as a pure beta emitter, and is used in targeted radionuclide therapy, with imaging performed using bremsstrahlung single-photon emission computed tomography (SPECT). However, because {sup 90}Y also emits positrons through internal pair production with a very small branching ratio, positron emission tomography (PET) imaging is also available. Because of the insufficient image quality of {sup 90}Y bremsstrahlung SPECT, PET imaging has been suggested as an alternative. In this paper, the authors present the Monte Carlo-based simulation–reconstruction framework for {sup 90}Y to comprehensively analyze the PET and SPECT imaging techniques and to quantitatively consider the disadvantages associated with them. Methods: Our PET and SPECT simulation modules were developed using Monte Carlo simulation of Electrons and Photons (MCEP), developed by Dr. S. Uehara. PET code (MCEP-PET) generates a sinogram, and reconstructs the tomography image using a time-of-flight ordered subset expectation maximization (TOF-OSEM) algorithm with attenuation compensation. To evaluate MCEP-PET, simulated results of {sup 18}F PET imaging were compared with the experimental results. The results confirmed that MCEP-PET can simulate the experimental results very well. The SPECT code (MCEP-SPECT) models the collimator and NaI detector system, and generates the projection images and projection data. To save the computational time, the authors adopt the prerecorded {sup 90}Y bremsstrahlung photon data calculated by MCEP. The projection data are also reconstructed using the OSEM algorithm. The authors simulated PET and SPECT images of a water phantom containing six hot spheres filled with different concentrations of {sup 90}Y without background activity. The amount of activity was 163 MBq, with an acquisition time of 40 min. Results: The simulated {sup 90}Y-PET image accurately simulated the experimental results. PET image is visually superior to SPECT image because of the low background noise. The simulation reveals that the detected photon number in SPECT is comparable to that of PET, but the large fraction (approximately 75%) of scattered and penetration photons contaminates SPECT image. The lower limit of {sup 90}Y detection in SPECT image was approximately 200 kBq/ml, while that in PET image was approximately 100 kBq/ml. Conclusions: By comparing the background noise level and the image concentration profile of both the techniques, PET image quality was determined to be superior to that of bremsstrahlung SPECT. The developed simulation codes will be very useful in the future investigations of PET and bremsstrahlung SPECT imaging of {sup 90}Y.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2018

A Fiber-Optic System Generating Pulses of High Spectral Density

A cascade fiber-optic system that generates pulses of high spectral density by using the effect of nonlinear spectral compression is proposed. It is demonstrated that the shape of the pulse envelope substantially influences the degree of compression of its spectrum. In so doing, maximum compression is achieved for parabolic pulses. The cascade system includes an optical fiber exhibiting normal dispersion that decreases along the fiber length, thereby ensuring that the pulse envelope evolves toward a parabolic shape, along with diffraction gratings and a fiber spectral compressor. Based on computer simulation, we determined parameters of cascade elements leading to maximum spectral density of radiation originating from a subpicosecond laser pulse of medium energy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2012

Melting and lattice dynamics of sodium at high pressures. Ab initio quantum molecular-dynamics analysis

The melting and lattice dynamics of sodium are studied by quantum molecular dynamics simulation, i.e., with allowance for anharmonicity, at pressures up to 1 Mbar and temperatures up to 1000 K. The simulation results agree well with the experimental data and our earlier calculation performed ab initio in the quasi-harmonic approximation. The simulation results demonstrate that anharmonic interactions weakly affect the melting curve and the phonon frequencies of Na up to near-melting temperatures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2013

Transition micromagnetic structures in asymmetric vortexlike domain walls (Static solutions and dynamic reconstructions)

The possible types of transition structures with a three-dimensional magnetization distribution over regions in the vortex asymmetric domain walls that exist in magnetically uniaxial soft magnetic films with in-plane anisotropy are studied by computer simulation in terms of a micromagnetic approach. It is shown that the possible structure types include both the type of vertical Bloch lines that was discussed earlier in other works and new types, namely, singular (Bloch) points and clusters consisting of vertical Bloch lines and Bloch points. The spatial configurations of the transition structures are calculated and their topological properties are found. The numerical simulation of the dynamics of closely spaced substructure regions reveals various scenarios of their interaction, including annihilation accompanied by energy release and the excitation of nonlinear waves.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗