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,531 records · Page 85

Simulation of effects of direction and air flow speed on temperature distribution in the room covered by various roof materials

The objective of this research is to get a comparison of the distribution of the room temperature by using three materials, namely plastic-rubber composite, clay, and asbestos. The simulation used Ansys Fluent to get the temperature distribution. There were two conditions in this simulations, first the air passing beside the room and second the air passing in front of the room. Each condition will be varied with the air speed of 1 m/s, 2 m/s, 3 m/s, 4 m/s, 5 m/s for each material used. There are three heat transfers in this simulation, namely radiation, convection, and conduction. Based on the ANSI/ ASHRAE Standard 55-2004, the results of the simulation showed that the best temperature distribution was the roof of plastic-rubber composites.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Study of CdTe/CdS solar cell at low power density for low-illumination applications

In this paper, we numerically investigate CdTe/CdS PV cell properties using a simulation program Solar Cell Capacitance Simulator in 1D (SCAPS-1D). A simple structure of CdTe PV cell has been optimized to study the effect of temperature, absorber thickness and work function at very low incident power. Objective of this research paper is to build an efficient and cost effective solar cell for portable electronic devices such as portable computers and cell phones that work at low incident power because most of such devices work at diffused and reflected sunlight. In this report, we simulated a simple CdTe PV cell at very low incident power, which gives good efficiency.

14 SOLAR ENERGY↗

Spin dependent correlations in a homogeneous electron gas at finite temperature

We have studied theoretically the magnetic structure factor of a three-dimensional homogeneous electron gas at finite temperature. The spin density response function has been derived using the Singwi-Tosi-Land-Sjolander (STLS) theory that incorporates the correlation effects through spin anti-symmetric local field correction factor. The numerical results so obtained are compared against the recent path-integral Monte Carlo Simulation data of Brown et al. in the warm-dense regime for various temperature values. We find almost exact agreement at small temperature for different coupling parameter (r{sub s}) values. However, with increasing temperature and decreasing density, there has been observed noticeable disagreement with simulation results. This is attributed to the known failure of the STLS theory in dealing separately with the spin-resolved correlation functions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Evaluating gyro-viscosity in the Kelvin-Helmholtz instability by kinetic simulations

In the present paper, the finite-Larmor-radius (gyro-viscous) term [K. V. Roberts and J. B. Taylor, Phys. Rev. Lett. 8, 197–198 (1962)] is evaluated by using a full kinetic Vlasov simulation result of the Kelvin-Helmholtz instability (KHI). The velocity field and the pressure tensor are calculated from the high-resolution data of the velocity distribution functions obtained by the Vlasov simulation, which are used to approximate the Finite-Larmor-Radius (FLR) term according to Roberts and Taylor [Phys. Rev. Lett. 8, 197–198 (1962)]. The direct comparison between the pressure tensor and the FLR term shows an agreement. It is also shown that the anisotropic pressure gradient enhanced the linear growth of the KHI when the inner product between the vorticity of the primary velocity shear layer and the magnetic field is negative, which is consistent with the previous FLR-magnetohydrodynamic simulation result. This result suggests that it is not sufficient for reproducing the kinetic simulation result by fluid simulations to include the FLR term (or the pressure tensor) only in the equation of motion for fluid.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

THE INFLUENCE OF OUTER SOLAR SYSTEM ARCHITECTURE ON THE STRUCTURE AND EVOLUTION OF THE OORT CLOUD

We study the influence of outer solar system architecture on the structural evolution of the Oort Cloud (OC) and the flux of Earth-crossing comets. In particular, we seek to quantify the role of the giant planets as ''planetary protectors''. To do so, we have run simulations in each of four different planetary mass configurations to understand the significance of each of the giant planets. Because the outer planets modify the structure of the OC throughout its formation, we integrate each simulation over the full age of the solar system. Over this time, we follow the evolution of cometary orbits from their starting point in the protoplanetary disk to their injection into the OC to their possible re-entry into the inner planetary region. We find that the overall structure of the OC, including the location of boundaries and the relative number of comets in the inner and outer parts, does not change significantly between configurations; however, as planetary mass decreases, the trapping efficiency (TE) of comets into the OC and the flux of comets into the observable region increases. We determine that those comets that evolve onto Earth-crossing orbits come primarily from the inner OC but show no preference for initial protoplanetary disk location. We also find that systems that have at least a Saturn-mass object are effective at deflecting possible Earth-crossing comets but the difference in flux between systems with and without such a planet is less than an order of magnitude. We conclude by discussing the individual roles of the planets and the implications of incorporating more realistic planetary accretion and migration scenarios into simulations, particularly on existing discrepancies between low TE and the mass of the protoplanetary disk and on determining the structural boundaries of the OC.

79 ASTRONOMY AND ASTROPHYSICS↗

SIMULATIONS OF TURBULENT DYNAMOS DRIVEN BY THE MAGNETIC HELICITY FLUX

We present results of numerical simulations of driven magnetohydrodynamic turbulence in a periodic box with an imposed large-scale sinusoidal shear. We produce strong large-scale magnetic fields through dynamo action. This dynamo action is highly dependent on the spectral properties of the forcing function. In particular, we can produce either dynamo or antidynamo action by choosing particular forcing functions for the small-scale turbulence. Our results show that the large-scale electromotive force is well correlated with the local accumulation of eddy-scale magnetic helicity. Conversely, the electromotive force is very weakly correlated with the large-scale magnetic field strength or its product with kinetic helicity. We argue that in these simulations large-scale dynamo action is driven by an inverse cascade of magnetic helicity, ultimately caused by the tendency of anisotropic turbulence to create a magnetic helicity flux.

79 ASTRONOMY AND ASTROPHYSICS↗

RADIATION-HYDRODYNAMIC SIMULATIONS OF PROTOSTELLAR OUTFLOWS: SYNTHETIC OBSERVATIONS AND DATA COMPARISONS

We present results from three-dimensional, self-gravitating, radiation-hydrodynamic simulations of low-mass protostellar outflows. We construct synthetic observations in {sup 12}CO in order to compare with observed outflows and evaluate the effects of beam resolution and outflow orientation on inferred outflow properties. To facilitate the comparison, we develop a quantitative prescription for measuring outflow opening angles. Using this prescription, we demonstrate that, in both simulations and synthetic observations, outflow opening angles broaden with time similarly to observed outflows. However, the interaction between the outflowing gas and the turbulent core envelope produces significant asymmetry between the redshifted and blueshifted outflow lobes. We find that applying a velocity cutoff may result in outflow masses that are underestimated by a factor five or more, and masses derived from optically thick CO emission further underpredict the mass of the high-velocity gas by a factor of 5-10. Derived excitation temperatures indicate that outflowing gas is hotter than the ambient gas with temperature rising over time, which is in agreement with the simulation gas temperatures. However, excitation temperatures are otherwise not well correlated with the actual gas temperature.

79 ASTRONOMY AND ASTROPHYSICS↗

THE METALLICITY EVOLUTION OF INTERACTING GALAXIES

Nuclear inflows of metal-poor interstellar gas triggered by galaxy interactions can account for the systematically lower central oxygen abundances observed in local interacting galaxies. Here, we investigate the metallicity evolution of a large set of simulations of colliding galaxies. Our models include cooling, star formation, feedback, and a new stochastic method for tracking the mass recycled back to the interstellar medium from stellar winds and supernovae. We study the influence of merger-induced inflows, enrichment, gas consumption, and galactic winds in determining the nuclear metallicity. The central metallicity is primarily a competition between the inflow of low-metallicity gas and enrichment from star formation. An average depression in the nuclear metallicity of {approx}0.07 is found for gas-poor disk-disk interactions. Gas-rich disk-disk interactions, on the other hand, typically have an enhancement in the central metallicity that is positively correlated with the gas content. The simulations fare reasonably well when compared to the observed mass-metallicity and separation-metallicity relationships, but further study is warranted.

79 ASTRONOMY AND ASTROPHYSICS↗

DWARF GALAXY FORMATION WITH H{sub 2}-REGULATED STAR FORMATION

We describe cosmological galaxy formation simulations with the adaptive mesh refinement code Enzo that incorporate a star formation prescription regulated by the local abundance of molecular hydrogen. We show that this H{sub 2}-regulated prescription leads to a suppression of star formation in low-mass halos (M{sub h} {approx}< 10{sup 10} M{sub Sun }) at z > 4, alleviating some of the dwarf galaxy problems faced by theoretical galaxy formation models. H{sub 2} regulation modifies the efficiency of star formation of cold gas directly, rather than indirectly reducing the cold gas content with 'supernova feedback'. We determine the local H{sub 2} abundance in our most refined grid cells (76 proper parsec in size at z = 4) by applying the model of Krumholz, McKee, and Tumlinson, which is based on idealized one-dimensional radiative transfer calculations of H{sub 2} formation-dissociation balance in {approx}100 pc atomic-molecular complexes. Our H{sub 2}-regulated simulations are able to reproduce the empirical (albeit lower z) Kennicutt-Schmidt relation, including the low {Sigma}{sub gas} cutoff due to the transition from atomic to molecular phase and the metallicity dependence thereof, without the use of an explicit density threshold in our star formation prescription. We compare the evolution of the luminosity function, stellar mass density, and star formation rate density from our simulations to recent observational determinations of the same at z = 4-8 and find reasonable agreement between the two.

79 ASTRONOMY AND ASTROPHYSICS↗

ANOMALOUS TRANSPORT OF HIGH-ENERGY COSMIC RAYS IN GALACTIC SUPERBUBBLES. I. NUMERICAL SIMULATIONS

We present a simple continuous-time random-walk model for the transport of energetic particles accelerated by a collection of supernova explosions in a galactic superbubble, developed to simulate and highlight signatures of anomalous transport on the particles' evolution and their spectra in a multi-shock context. We assume standard diffusive shock acceleration (DSA) theory for each shock encounter. The superbubble (an OB stars association) is idealized as a heterogeneous region of particle sources and sinks bounded by a random surface. The model is based on two coupled stochastic differential equations and is applied for protons and alpha particles. Using characteristic values for a typical bubble, our simulations suggest that acceleration and transport in the bubble may be sub-diffusive. In addition, a spectral break in the particles' evolution and spectra is evident located at Almost-Equal-To 10{sup 15} eV for protons and Almost-Equal-To 3 Multiplication-Sign 10{sup 15} eV for alphas. Our simulations are consistent with a bubble's mean magnetic field strength of Almost-Equal-To 1 {mu}G and a shock separation distance {approx}0.1 Multiplication-Sign the characteristic radius of the bubble. The simulations imply that the diffusion coefficient (for the elementary shock acceleration process) is {approx}< 10{sup 27} cm{sup 2} s{sup -1} at 1 GeV/c. While the sub-diffusive transport is readily attributed to the stochastic nature of the acceleration time according to DSA theory, the spectral break appears to be an artifact of transport in a finite medium. These simulations point to a new and intriguing phenomenon associated with the statistical nature of collective acceleration of high-energy cosmic rays in galactic superbubbles.

79 ASTRONOMY AND ASTROPHYSICS↗

COMPARING THE EFFECT OF RADIATIVE TRANSFER SCHEMES ON CONVECTION SIMULATIONS

We examine the effect of different radiative transfer schemes on the properties of three-dimensional (3D) simulations of near-surface stellar convection in the superadiabatic layer, where energy transport transitions from fully convective to fully radiative. We employ two radiative transfer schemes that fundamentally differ in the way they cover the 3D domain. The first solver approximates domain coverage with moments, while the second solver samples the 3D domain with ray integrations. By comparing simulations that differ only in their respective radiative transfer methods, we are able to isolate the effect that radiative efficiency has on the structure of the superadiabatic layer. We find the simulations to be in good general agreement, but they show distinct differences in the thermal structure in the superadiabatic layer and atmosphere.

79 ASTRONOMY AND ASTROPHYSICS↗

A ROBUST MEASURE OF COSMIC STRUCTURE BEYOND THE POWER SPECTRUM: COSMIC FILAMENTS AND THE TEMPERATURE OF DARK MATTER

We discover that the mass of dark matter particles m {sub DM} is imprinted in phase correlations of the cosmic density field more significantly than in the two-point correlation. In particular, phase correlations trace m {sub DM} out to scales about five times larger than the two-point correlation. This result relies on a new estimator l(r) of pure phase information in Fourier space, which can be interpreted as a parameter-free and scale-invariant tracer of filament-like structure. Based on simulated density fields, we show how m {sub DM} can, in principle, be measured using l(r), given a suitably reconstructed density field.

79 ASTRONOMY AND ASTROPHYSICS↗

COSMOLOGICAL MAGNETOHYDRODYNAMIC SIMULATIONS OF GALAXY CLUSTER RADIO RELICS: INSIGHTS AND WARNINGS FOR OBSERVATIONS

Non-thermal radio emission from cosmic-ray electrons in the vicinity of merging galaxy clusters is an important tracer of cluster merger activity, and is the result of complex physical processes that involve magnetic fields, particle acceleration, gas dynamics, and radiation. In particular, objects known as radio relics are thought to be the result of shock-accelerated electrons that, when embedded in a magnetic field, emit synchrotron radiation in the radio wavelengths. In order to properly model this emission, we utilize the adaptive mesh refinement simulation of the magnetohydrodynamic evolution of a galaxy cluster from cosmological initial conditions. We locate shock fronts and apply models of cosmic-ray electron acceleration that are then input into radio emission models. We have determined the thermodynamic properties of this radio-emitting plasma and constructed synthetic radio observations to compare observed galaxy clusters. We find a significant dependence of the observed morphology and radio relic properties on the viewing angle of the cluster, raising concerns regarding the interpretation of observed radio features in clusters. We also find that a given shock should not be characterized by a single Mach number. We find that the bulk of the radio emission comes from gas with T > 5 Multiplication-Sign 10{sup 7} K, {rho} {approx} 10{sup -28}-10{sup -27} g cm{sup -3}, with magnetic field strengths of 0.1-1.0 {mu}G, and shock Mach numbers of M {approx} 3-6. We present an analysis of the radio spectral index which suggests that the spatial variation of the spectral index can mimic synchrotron aging. Finally, we examine the polarization fraction and position angle of the simulated radio features, and compare to observations.

79 ASTRONOMY AND ASTROPHYSICS↗

THE FORCE-FREE MAGNETOSPHERE OF A ROTATING BLACK HOLE

We revisit the Blandford-Znajek process and solve the fundamental equation that governs the structure of the steady-state force-free magnetosphere around a Kerr black hole. The solution depends on the distributions of the magnetic field angular velocity {omega} and the poloidal electric current I. These are not arbitrary. They are determined self-consistently by requiring that magnetic field lines cross smoothly the two singular surfaces of the problem: the inner ''light surface'' located inside the ergosphere and the outer ''light surface'' which is the generalization of the pulsar light cylinder. We find the solution for the simplest possible magnetic field configuration, the split monopole, through a numerical iterative relaxation method analogous to the one that yields the structure of the steady-state axisymmetric force-free pulsar magnetosphere. We obtain the rate of electromagnetic extraction of energy and confirm the results of Blandford and Znajek and of previous time-dependent simulations. Furthermore, we discuss the physical applicability of magnetic field configurations that do not cross both ''light surfaces''.

79 ASTRONOMY AND ASTROPHYSICS↗

Numerical simulations of active region scale flux emergence: From spot formation to decay

We present numerical simulations of active region scale flux emergence covering a time span of up to 6 days. Flux emergence is driven by a bottom boundary condition that advects a semi-torus of magnetic field with 1.7 × 10{sup 22} Mx flux into the computational domain. The simulations show that, even in the absence of twist, the magnetic flux is able the rise through the upper 15.5 Mm of the convection zone and emerge into the photosphere to form spots. We find that spot formation is sensitive to the persistence of upflows at the bottom boundary footpoints, i.e., a continuing upflow would prevent spot formation. In addition, the presence of a torus-aligned flow (such flow into the retrograde direction is expected from angular momentum conservation during the rise of flux ropes through the convection zone) leads to a significant asymmetry between the pair of spots, with the spot corresponding to the leading spot on the Sun being more axisymmetric and coherent, but also forming with a delay relative to the following spot. The spot formation phase transitions directly into a decay phase. Subsurface flows fragment the magnetic field and lead to intrusions of almost field free plasma underneath the photosphere. When such intrusions reach photospheric layers, the spot fragments. The timescale for spot decay is comparable to the longest convective timescales present in the simulation domain. We find that the dispersal of flux from a simulated spot in the first two days of the decay phase is consistent with self-similar decay by turbulent diffusion.

79 ASTRONOMY AND ASTROPHYSICS↗

Uncorrelated estimates of the primordial power spectrum

We use the localized principle component analysis to detect deviations from scale invariance of the primordial power spectrum of curvature perturbations. With the technique we make uncorrelated estimates of the primordial power spectrum with five wavenumber bins. In the framework of a minimal ΛCDM model, using the latest cosmic microwave background data from the WMAP and ACT experiments we find that more than 95% of the preferred models are incompatible with the assumption of scale-invariance, but still compatible with a power-law primordial spectrum. We also forecast the sensitivity and constraints achievable by the Planck experiment by performing Monte Carlo studies on simulated data. Planck could significantly improve the constraints on the primordial power spectrum, especially at small scales by roughly a factor of 4.

79 ASTRONOMY AND ASTROPHYSICS↗

THRESHING IN ACTION: THE TIDAL DISRUPTION OF A DWARF GALAXY BY THE HYDRA I CLUSTER

We report on the discovery of strong tidal features around a dwarf spheroidal galaxy in the Hydra I galaxy cluster, indicating its ongoing tidal disruption. This very low surface brightness object, HCC-087, was originally classified as an early-type dwarf in the Hydra Cluster Catalogue (HCC), but our re-analysis of the ESO-VLT/FORS images of the HCC unearthed a clear indication of an S-shaped morphology and a large spatial extent. Its shape, luminosity (M{sub V} = -11.6 mag), and physical size (at a half-light radius of 3.1 kpc and a full length of {approx}5.9 kpc) are comparable to the recently discovered NGC 4449B and the Sagittarius dwarf spheroidal, all of which are undergoing clear tidal disruption. Aided by N-body simulations we argue that HCC-087 is currently at its first apocenter, at 150 kpc, around the cluster center and that it is being tidally disrupted by the galaxy cluster's potential itself. An interaction with the nearby (50 kpc) S0 cluster galaxy HCC-005, at M{sub *} {approx} 3 Multiplication-Sign 10{sup 10} M{sub Sun} is rather unlikely, as this constellation requires a significant amount of dynamical friction and thus low relative velocities. The S-shaped morphology and large spatial extent of the satellite would, however, also appear if HCC-087 would orbit the cluster center. These features appear to be characteristic properties of satellites that are seen in the process of being tidally disrupted, independent of the environment of the destruction. An important finding of our simulations is an orientation of the tidal tails perpendicular to the orbit.

79 ASTRONOMY AND ASTROPHYSICS↗

COUNTERROTATION IN MAGNETOCENTRIFUGALLY DRIVEN JETS AND OTHER WINDS

Rotation measurement in jets from T Tauri stars is a rather difficult task. Some jets seem to be rotating in a direction opposite to that of the underlying disk, although it is not yet clear if this affects the totality or part of the outflows. On the other hand, Ulysses data also suggest that the solar wind may rotate in two opposite ways between the northern and southern hemispheres. We show that this result is not as surprising as it may seem and that it emerges naturally from the ideal MHD equations. Specifically, counterrotating jets neither contradict the magnetocentrifugal driving of the flow nor prevent extraction of angular momentum from the disk. The demonstration of this result is shown by combining the ideal MHD equations for steady axisymmetric flows. Provided that the jet is decelerated below some given threshold beyond the Alfven surface, the flow will change its direction of rotation locally or globally. Counterrotation is also possible for only some layers of the outflow at specific altitudes along the jet axis. We conclude that the counterrotation of winds or jets with respect to the source, star or disk, is not in contradiction with the magnetocentrifugal driving paradigm. This phenomenon may affect part of the outflow, either in one hemisphere, or only in some of the outflow layers. From a time-dependent simulation, we illustrate this effect and show that it may not be permanent.

79 ASTRONOMY AND ASTROPHYSICS↗