Engineering topics

Kang, Hyesung

Publications and source records attributed to Kang, Hyesung.

13 records

NASA NTRS2000

Cosmic Rays and Their Radiative Processes in Numerical Cosmology

A cosmological hydrodynamic code is described, which includes a routine to compute cosmic ray acceleration and transport in a simplified way. The routine was designed to follow explicitly diffusive, acceleration at shocks, and second-order Fermi acceleration and adiabatic loss in smooth flows. Synchrotron cooling of the electron population can also be followed. The updated code is intended to be used to study the properties of nonthermal synchrotron emission and inverse Compton scattering from electron cosmic rays in clusters of galaxies, in addition to the properties of thermal bremsstrahlung emission from hot gas. The results of a test simulation using a grid of 128 (exp 3) cells are presented, where cosmic rays and magnetic field have been treated passively and synchrotron cooling of cosmic ray electrons has not been included.

Ryu, Dongsu↗

NASA NTRS1994

Cosmic bullets as particle accelerators and radio sources

We have simulated in two dimensions the dynamical evolution of dense gas clouds(`cosmic bullets') moving supersonically through a uniform low-density medium. The diffusive shock acceleration of relativistic protons (cosmic rays) and their dynamical feedback on the background flow are included by the two-fluid model for this process. The acceleration of relativistic electrons is approximated by a test-particle model, and a passive magnetic field is followed by a simple advection scheme. Strong bow shocks, with Mach numbers similar to that of a bullet's motion, are the most important particle accelerators in the flow, while tail shocks and shocks inside the bullets do not play generally significant roles in this regard. For our simulation parameters, approximately greater than 10% of the initial bullet kinetic energy is converted to a combination of internal energy of gas and cosmic-ray protons by the time the bullets begin to be disrupted. Characteristically, the cosmic rays gain several percent of the available kinetic energy. Bullet destruction on timescales only a little larger than the ram pressure bullet crushing time begins in response to Kelvin-Helmholtz and especially to Rayleigh-Taylor instabilities along the forward bullet surface. For dense bullets this happens before the bullet is stopped by ram pressure. According to our simple model for synchrotron emission from relativistic electrons accelerated and transported within the flows, that emission increases rapidly as the bullet begins to fragment, when it is strongly dominated by field enhancement in sheared flows. Synchrotron emission from the acceleration region within the bow shock is, by contrast, much weaker.

Jones, T. W.↗

NASA NTRS1994

A comparison of cosmological hydrodynamic codes

We present a detailed comparison of the simulation results of various hydrodynamic codes. Starting with identical initial conditions based on the cold dark matter scenario for the growth of structure, with parameters h = 0.5 Omega = Omega(sub b) = 1, and sigma(sub 8) = 1, we integrate from redshift z = 20 to z = O to determine the physical state within a representative volume of size L(exp 3) where L = 64 h(exp -1) Mpc. Five indenpendent codes are compared: three of them Eulerian mesh-based and two variants of the smooth particle hydrodynamics 'SPH' Lagrangian approach. The Eulerian codes were run at N(exp 3) = (32(exp 3), 64(exp 3), 128(exp 3), and 256(exp 3)) cells, the SPH codes at N(exp 3) = 32(exp 3) and 64(exp 3) particles. Results were then rebinned to a 16(exp 3) grid with the exception that the rebinned data should converge, by all techniques, to a common and correct result as N approaches infinity. We find that global averages of various physical quantities do, as expected, tend to converge in the rebinned model, but that uncertainites in even primitive quantities such as (T), (rho(exp 2))(exp 1/2) persists at the 3%-17% level achieve comparable and satisfactory accuracy for comparable computer time in their treatment of the high-density, high-temeprature regions as measured in the rebinned data; the variance among the five codes (at highest resolution) for the mean temperature (as weighted by rho(exp 2) is only 4.5%. Examined at high resolution we suspect that the density resolution is better in the SPH codes and the thermal accuracy in low-density regions better in the Eulerian codes. In the low-density, low-temperature regions the SPH codes have poor accuracy due to statiscal effects, and the Jameson code gives the temperatures which are too high, due to overuse of artificial viscosity in these high Mach number regions. Overall the comparison allows us to better estimate errors; it points to ways of improving this current generation ofhydrodynamic codes and of suiting their use to problems which exploit their best individual features.

Kang, Hyesung↗

NASA NTRS1994

Hot gas in the cold dark matter scenario: X-ray clusters from a high-resolution numerical simulation

A new, three-dimensional, shock-capturing hydrodynamic code is utilized to determine the distribution of hot gas in a standard cold dark matter (CDM) model of the universe. Periodic boundary conditions are assumed: a box with size 85 h(exp -1) Mpc having cell size 0.31 h(exp -1) Mpc is followed in a simulation with 270(exp 3) = 10(exp 7.3) cells. Adopting standard parameters determined from COBE and light-element nucleosynthesis, sigma(sub 8) = 1.05, omega(sub b) = 0.06, and assuming h = 0.5, we find the X-ray-emitting clusters and compute the luminosity function at several wavelengths, the temperature distribution, and estimated sizes, as well as the evolution of these quantities with redshift. We find that most of the total X-ray emissivity in our box originates in a relatively small number of identifiable clusters which occupy approximately 10(exp -3) of the box volume. This standard CDM model, normalized to COBE, produces approximately 5 times too much emission from clusters having L(sub x) is greater than 10(exp 43) ergs/s, a not-unexpected result. If all other parameters were unchanged, we would expect adequate agreement for sigma(sub 8) = 0.6. This provides a new and independent argument for lower small-scale power than standard CDM at the 8 h(exp -1) Mpc scale. The background radiation field at 1 keV due to clusters in this model is approximately one-third of the observed background, which, after correction for numerical effects, again indicates approximately 5 times too much emission and the appropriateness of sigma(sub 8) = 0.6. If we have used the observed ratio of gas to total mass in clusters, rather than basing the mean density on light-element nucleosynthesis, then the computed luminosity of each cluster would have increased still further, by a factor of approximately 10. The number density of clusters increases to z approximately 1, but the luminosity per typical cluster decreases, with the result that evolution in the number density of bright clusters is moderate in this redshift range, showing a broad peak near z = 0.7, and then a rapid decline above redshift z = 3. Detailed computations of the luminosity functions in the range L(sub x) = 10(exp 40) - 10(exp 44) ergs/s in various energy bands are presented for both cluster central regions and total luminosities to be used in comparison with ROSAT and other observational data sets. The quantitative results found disagree significantly with those found by other investigators using semianalytic techniques. We find little dependence of core radius on cluster luminosity and a dependence of temperature on luminosity given by log kT(sub x) = A + B log L(sub x), which is slightly steeper (B = 0.38) than is indicated by observations. Computed temperatures are somewhat higher than observed, as expected, in that COBE-normalized CDM has too much power on the relevant scales. A modest average temperature gradient is found, with temperatures dropping to 90% of central values at 0.4 h(exp -1) Mpc and 70% of central values at 0.9 h(exp -1) Mpc. Examining the ratio of gas to total mass in the clusters normalized to Omega(sub B) h(exp 2) = 0.015, and comparing with observations, we conclude, in agreement with White (1991), that the cluster observations argue for an open universe.

Kang, Hyesung↗

NASA NTRS1994

The dynamical and radiative evolution of clumpy supernova ejecta

Numerical simulations describing the dynamical and radiative evolution of clumpy supernova ejecta are compared with observations of optical and radio emission knots in supernova remnant (SNR) Cassiopeia A. Three major phases are identified in the evolution of clumpy ejecta: a bow-shock phase, an instability phase, and a dispersal phase. The phenomenological and radiative signatures of each phase are discussed and compared with multi-epoch measurements of small-scale features in Cas A. Good correspondence is found between theory and observations. Both support the premise that compact radio emission features are controlled more by magnetic field amplification triggered in the instability phase than by in situ acceleration of new relativistic particles.

Anderson, M. C.↗

NASA NTRS1993

A cosmological hydrodynamic code based on the total variation diminishing scheme

An explicitly second-order finite difference code based on a total variation diminishing scheme for self-gravitating cosmological hydrodynamic systems is described. The code was developed to follow correctly the adiabatic changes of extremely supersonic preshock flows with a Mach number larger than 100 as well as very strong shocks. In calculations of the formation of a 1D Zel'dovich pancake, an energy accuracy of 1 percent is obtained for 32 cells per unit wavelength, and the accuracy reaches 0.01 percent as the number of cells approach 1024. It is shown that shocks are well resolved and separate cleanly the hot, dense, collapsed peaks from the cold, low-density, expanding voids.

Ryu, Dongsu↗

NASA NTRS1993

The stability of cosmic-ray-dominated shocks - A secondary instability

Previous studies in one spatial dimension showed that the precursors of cosmic-ray-dominated shocks are unstable against traveling acoustic disturbances. Here we report that a secondary, Rayleigh-Taylor type instability can exist in cosmic-ray-dominated media influenced by the acoustic instability. Using the local WKB approximation, the growth rate of the secondary instability is shown to be comparable to that of the one-dimensional acoustic instability itself in the cases we have considered. The nonlinear development has been followed numerically with a two-dimensional PPM hydrodynamics code that also incorporates the two-fluid cosmic-ray energy equation. We show that the secondary instability may cause the precursor and postshock flows to become highly turbulent.

Ryu, Dongsu↗

NASA NTRS1993

Cosmic-ray acceleration during the impact of shocks on dense clouds

In order to elucidate the properties of diffusive shock acceleration in nonuniform environments, an extensive set of simulations of the dynamical interactions between plane nonradiative shocks and dense gas clouds was carried out initially in static equilibrium with their environments. These time-dependent calculations are based on the two-fluid model for diffusive cosmic ray transport, and include the dynamically active energetic proton component of the cosmic rays as well as passive electron and magnetic field components. Except when the incident shock is itself already dominated by cosmic ray pressure, it is found that the presence of the cloud adds little to the net acceleration efficiency of the original shock and can, in fact, reduce slightly the net amount of energy transferred to cosmic rays after a given time. It is found that, in 2D cloud simulations, the always-weak bow shock and the shock inside the cloud are less important to acceleration during the interaction than the tail shock.

Jones, T. W.↗

NASA NTRS1992

Cosmic-ray-modified supernova remnant shocks

The evolution of cosmic-ray-modified SNRs during the early Sedov phase is investigated using the time-dependent, two-fluid model for diffusive shock acceleration. Consideration is given to the sensitivity of net acceleration efficiency to model assumptions regarding density structure of the external medium as well as time dependence in the diffusion coefficient and the cosmic-ray adiabatic index. Of the model assumptions explored, the greatest sensitivity was found to time variations in the specific heat ratio to cosmic rays. It was found that if tau is greater than or approximately equal to few times 100, dynamically significant cosmic-ray pressures are produced early in the Sedov inefficient cosmic-ray accelerators. For the SNR models under consideration, the total energy channeled into the cosmic rays can be of order 10 percent of the initial blast energy for both the uniform density and inverse-square density interstellar medium models.

Jones, T. W.↗

NASA NTRS1992

Magnetic fields and radiative shocks in protogalaxies and the origin of globular clusters

The paper examines the hypothesis that globular clusters formed from gravitational instability in dense sheets of gas produced behind radiative shocks inside protogalaxies, such as those produced by the collision of subgalactic mass fragments partaking of the virial motions within the protogalaxy, in order to determine the differences which result if a magnetic field is present in the preshock medium. The MHD conservation equations are solved along with rate equations for nonequilibrium ionization, recombination, molecular formation and dissociation, and the equations of radiative transfer for steady-state shocks of velocity 300 km/s in a gas of preshock densities of 0.1-1 cu cm, and magnetic field strengths of 0.1-6 micro-G. The magnetic field is found to limit the degree of postshock compression and, thereby, to reduce the level of external radiation flux required to suppress H2 formation and cooling.

Shapiro, Paul R.↗

NASA NTRS1992

Radiative shocks and hydrogen molecules in pregalactic gas - The effects of postshock radiation

Previous calculations of radiative shocks in pregalactic gas are generalized and substantially improved. The hydrodynamical conservations are solved, along with the rate equations for nonequilibrium ionization, recombination, and molecule formation, and the equation of radiative transfer, for steady-state shocks in a gas of primordial compositions. This calculation self-consistently includes the effects of the diffuse postshock emission as well as a possible external radiation flux on the postshock flow and the preshock ionization levels. A previous result is confirmed that the shocked gas cools faster than it can recombine and, as a result, is able to form an H2 concentration as high as 0.003 via the formation of H(-) and H2(+) intermediaries due to the enhanced nonequilibrium ionization fraction at 10,000 K. It is shown that, as the level of the external radiation flux is increased, the formation of and cooling by H2 molecules can be inhibited and delayed. The detailed radiation spectra and ionizing photon number fluxes emergent from these shocks are presented.

Kang, Hyesung↗

NASA NTRS1992

Numerical simulations of time-dependent cosmic ray mediated shocks

Recent research in unsteady shocks that are influenced by nonlinear aspects of diffusive particle acceleration and that involve multidimensional flows is outlined, and some of the difficulties in making further progress are discussed. It is suggested that in many situations it is probably not appropriate to assume steady states for astrophysical CR shocks. This results both from the need to be concerned about the finite rate at which CR pressure builds and modifies a shock and from the existence of instabilities. Even so, it seems likely that CR-modified shocks can develop to a state with substantial CR pressure. The overall acceleration efficiency and the CR particle momentum distribution produced in time-dependent flows typically depend on several details such as the ratio of various time scales and the initial proportion of relativistic to nonrelativistic particles in the CR populations.

Jones, T. W.↗

NASA NTRS1990

Radiative shocks and nonequilibrium chemistry in the early universe - Galaxy and primordial star formation

The nonequilibrium radiative cooling, recombination, and molecule formation behind steady-state shock waves in primordial composition gas is studied. The hydrodynamical conservation equations are solved as well as equations for nonequilibrium ionization, recombination, and molecule formation and the equation of radiative transfer. The shocked gas is found to cool faster than it can recombine; as a result it is able to form an H2 concentration as high as 0.001 or higher via the formation of H(-) and H2(+) intermediaries due to the enhanced nonequilibrium ionization at 10,000 K. With such an H2 concentration, the gas cools by rotational-vibrational line excitation of H2 molecules to well below the canonical final temperature of 10,000 K for a molecule-free gas without metals. It is shown that, as the level of external ionizing and dissociating radiation flux is increased, the formation of and cooling by H2 molecules can be inhibited and delayed.

Shapiro, Paul R.↗