Engineering topics

Ryu, Dongsu

Publications and source records attributed to Ryu, Dongsu.

At least 19 records

NASA NTRS2000

Three-Dimensional Simulations of the Parker Instability in a Uniformly-Rotating Disk

We investigate the, nonlinear effects of uniform rotation on the Parker instability in an exponentially-stratified disk through high-resolution simulations. During the linear stage, the speed of gas motion is subsonic and the evolution with the rotation is not much different from that without the rotation. This is because the Coriolis force is small. During the nonlinear stage, oppositely-directed supersonic flows near a magnetic valley are under the influence of the Coriolis force with different directions, resulting in twisted magnetic field lines near the valley. Sheet-like structures, which are tilted with respect to the initial field direction, are formed with an 1.5 enhancement of column density with respect to its initial value. Even though uniform rotation doesn't give much impact on density enhancement, it generates helically twisted field lines, which may become an additional support mechanism of clouds.

Kim, Jongsoo↗

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 NTRS1995

Numerical magetohydrodynamics in astronphysics: Algorithm and tests for one-dimensional flow`

We describe a numerical code to solve the equations for ideal magnetohydrodynamics (MHD). It is based on an explicit finite difference scheme on an Eulerian grid, called the total variation diminishing (TVD) scheme, which is a second-order-accurate extension of the Roe-type upwind scheme. We also describe a non-linear Riemann solver for ideal MHD, which includes rarefractions as well as shocks. The numerical code and the Riemann solver have been used to test each other. Extensive tests encompassing all the possible ideal MHD structures with planar symmetries (i.e., one-dimensional flows) are presented. These include those for which the field structure is two dimensional (i.e., those flows often called '1 + 1/2 dimensional') as well as those for which the magnetic field plane rotates (i.e., those flows often called '1 + 1/2 + 1/2 dimensional'). Results indicate that the code can resolve strong fast, slow, and magnetosonic shocks within two to four cells, but more cells are required if shocks become weak. With proper steepening, we could resolve rotational discontinuities within three to five cells. However, without successful implementation of steepening, contact discontinuities are resolved with approximately 10 cells and tangential discountinuities are resolved with approximately 15 cells. Out tests confirm that slow compound structures with tow-dimensional magnetic fields are composed of intermediate shocks (so-called 2-4 intermediate shocks) followed by slow rarefaction waves. Finally, tests demostrate that in two-dimensional magnetohydrodynamics, fast compound structures, which are composed of intermediate shocks (so-called 1-3 intermediate shocks) preceeded by fast rarefaction waves, are also possible.

Ryu, Dongsu↗

NASA NTRS1995

Time-dependent simulation of oblique MHD cosmic-ray shocks using the two-fluid model

Using a new, second-order accurate numerical method we present dynamical simulations of oblique MHD cosmic-ray (CR)-modified plane shock evolution. Most of the calculations are done with a two-fluid model for diffusive shock acceleration, but we provide also comparisons between a typical shock computed that way against calculations carried out using the more complete, momentum-dependent, diffusion-advection equation. We also illustrate a test showing that these simulations evolve to dynamical equilibria consistent with previously published steady state analytic calculations for such shocks. In order to improve understanding of the dynamical role of magnetic fields in shocks modified by CR pressure we have explored for time asymptotic states the parameter space of upstream fast mode Mach number, M(sub f), and plasma beta. We compile the results into maps of dynamical steady state CR acceleration efficiency, epsilon(sub c). We have run simulations using constant, and nonisotropic, obliquity (and hence spatially) dependent forms of the diffusion coefficient kappa. Comparison of the results shows that while the final steady states achieved are the same in each case, the history of CR-MHD shocks can be strongly modified by variations in kappa and, therefore, in the acceleration timescale. Also, the coupling of CR and MHD in low beta, oblique shocks substantially influences the transient density spike that forms in strongly CR-modified shocks. We find that inside the density spike a MHD slow mode wave can be generated that eventually steepens into a shock. A strong layer develops within the density spike, driven by MHD stresses. We conjecture that currents in the shear layer could, in nonplanar flows, results in enhanced particle accretion through drift acceleration.

Frank, Adam↗

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

Nonlinear evolution of tidally distorted accretion disks: Two-dimensional simulations

According to a previously published linear analysis, the tidal distortion of accretion disks in binary star systems produces a local hydrodynamic instability to m = 1 internal waves, which may have arbitrarily small wavelengths in the absence of viscosity. The instability is three-dimensional and approximately incrompressible. To explore the nonlinear outcome of this instability, we develop a shearing-sheet approximation on scales comparable to the disk thickness. The large-scale azimuthal variation of the disk is represented by varying the local metric with time (local orbital phase). The hydrodynamic equations can then be posed two-dimensionally on local meridional planes. We solve these equations with a second-order gasdynamical code based on the Total-Variation-Diminishing scheme. Our simulations confirm the predicted linear growth rate. The modes saturate chaotically at velocities scaling as the product of the linear growth rate and the wavelength. If the wavelength is small compared with the disk thickness, the modes remain nearly incompressible even when nonlinear. The two-dimensional power spectrum of velocities after saturation is roughly isotropic and extends over a broad range of scales in an approximately power-law fashion. We measure the heating rate associated with the nonlinear dissipation of the modes. The dissipation implies a secular torque on the disk and a return of angular momentum to the secondary star via the tidal potential. The estimated torque is somewhat larger than the tidal torque produced by maximal disk viscosity (alpha approximately 1). At least in these two-dimensional simulations, however, there is no significant angular momentum flux within the disk.

Ryu, Dongsu↗

NASA NTRS1994

Oblique magnetohydrodynamic cosmic-ray-modified shocks: Two-fluid numerical simulations

We present the first results of time-dependent two-fluid cosmic-ray (CR) modified magnetohydrodynamic (MHD) shock simulations. The calculations were carried out with a new numerical code for one-dimensional ideal MHD. By coupling this code with the CR energy transport equation we can simulate the time-dependent evolution of MHD shocks, including the acceleration of the CR and their feedback on the shock structures. We report tests of the combined numerical method including comparisons with analytical steady state results published earlier by Webb, as well as internal consistency checks for more general MHD CR shock structures after they apppear to have converged to dynamical steady states. We also present results from an initial time-dependent simulation which extends the parameter space domain of previous analytical models. These new results support Webb's suggestion that equilibrium oblique shocks are less effective than parallel shocks in the acceleration of CR. However, for realistic models of anisotropic CR diffusion, oblique shocks may achieve dynamical equilibrium on shorter timescales than parallel shocks.

Frank, Adam↗

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 NTRS1992

Convective instability in differentially rotating disks

A normal mode analysis for nonaxisymmetric perturbations in a thin, differentially rotating disk with a vertical structure that is isothermal and convectively unstable is performed. The vertical gravity is assumed to be external and constant. The perturbation scale is assumed to be much shorter than the radius of the disk but comparable to or less than the thickness. The initial value problem is formulated in shearing coordinates. Dispersion relations are obtained for the three limiting cases of zero shear, axisymmetric perturbations, and small radial wavelengths. The full effects of shear are studied by integrating numerically the initial value problem. Nonaxisymmetric local Fourier modes are found to have a radial wavenumber that increases linearly with time in proportion to the shear times the azimuthal wavenumber. While Coriolis forces exert stabilizing effects on the convective modes, reducing their growth rate and the range of unstable wavelengths, shear has destablizing effects inasmuch as it reduces the epicyclic frequency at a given angular velocity. In a Keplerian disk, perturbations with azimuthal wavelengths about 2 times smaller than vertical wavelengths grow exponentially.

Ryu, Dongsu↗

NASA NTRS1991

The dynamic instability of adiabatic blast waves

Adiabatic blastwaves, which have a total energy injected from the center E varies as t(sup q) and propagate through a preshock medium with a density rho(sub E) varies as r(sup -omega) are described by a family of similarity solutions. Previous work has shown that adiabatic blastwaves with increasing or constant postshock entropy behind the shock front are susceptible to an oscillatory instability, caused by the difference between the nature of the forces on the two sides of the dense shell behind the shock front. This instability sets in if the dense postshock layer is sufficiently thin. The stability of adiabatic blastwaves with a decreasing postshock entropy is considered. Such blastwaves, if they are decelerating, always have a region behind the shock front which is subject to convection. Some accelerating blastwaves also have such region, depending on the values of q, omega, and gamma where gamma is the adiabatic index. However, since the shock interface stabilizes dynamically induced perturbations, blastwaves become convectively unstable only if the convective zone is localized around the origin or a contact discontinuity far from the shock front. On the other hand, the contact discontinuity of accelerating blastwaves is subject to a strong Rayleigh-Taylor instability. The frequency spectra of the nonradial, normal modes of adiabatic blastwaves have been calculated. The results have been applied to the shocks propagating through supernovae envelopes. It is shown that the metal/He and He/H interfaces are strongly unstable against the Rayleigh-Taylor instability. This instability will induce mixing in supernovae envelopes. In addition the implications of this work for the evolution of planetary nebulae is discussed.

Ryu, Dongsu↗

NASA NTRS1991

Oblique MHD cosmic-ray modified shocks: Two-fluid numerical simulations

We present the first results of time dependent, two-fluid, cosmic-ray (CR) modified, MHD shock simulations. The calculations were carried out with a new numerical code for 1-D ideal MHD. By coupling this code with the CR energy transport equation we can simulate the time-dependent evolution of MHD shocks including the acceleration of the CR and their feedback on the shock structures. We report tests of the combined numerical method including comparisons with analytical steady state results published earlier by Webb, as well as internal consistency checks for more general MHD CR shock structures after they appear to have converged to dynamical steady states. We also present results from an initial time dependent simulation which extend the parameter space domain of previous analytical models. These new results support Webb's suggestion that equilibrium oblique shocks are less effective than parallel shocks in the acceleration of CR. However, for realistic models of anisotropic CR diffusion, oblique shocks may achieve dynamical equilibrium on shorter timescale than parallel shocks.

Frank, Adam↗

NASA NTRS1991

Phase space distribution of halo particles and detection of WIMPS

We calculate the velocity distribution of halo particles using realistic models for the density profiles of the Galactic disk, spheroid, and halo components. We get the resulting velocity distribution, which is different from the Maxwellian distribution of isothermal halo. We investigate the consequences concerning the time averaged count rate and the amplitude of the annual modulation in WIMPS detection experiments. Compared to the previous estimate based on the Maxwellian distribution, the time averaged count rate is up to 40 percent larger. On the other hand, the amplitude of the annual modulation is up to a factor of two larger.

Ryu, Dongsu↗

NASA NTRS1990

A model for the distribution of the intergalactic medium

The evolution and distribution of the intergalactic medium (IGM) in a universe dominated by cold dark matter with Omega(0) = 1 and h = 0.5 are investigated. Galaxies form and eject energy into the IGM from z about 20 up to the present, and the distribution of the IGM is dominated by large connected structures. The power spectrum and two-point correlation function of the IGM show a suppressed growth due to the energy injected from galaxies and the mass subtraction to form galaxies. The high-temperature regions of the IGM correspond to the low-density regions and the low-temperature regions correspond to the high-density regions. The temperature of the IGM increases from z = 1 to z = 0, while the prsssure decreases. The present temperature distribution shows a peak at about 10 million K. The mass fraction of the IGM with temperature below 100,000 K is negligible, indicating almost all the hydrogen is ionized.

Ryu, Dongsu↗

NASA NTRS1990

Neutron stars and white dwarfs in galactic halos?

The possibility that galactic halos are composed of stellar remnants such as neutron stars and white dwarfs is discussed. On the basis of a simple model for the evolution of galactic halos, researchers follow the history of halo matter, luminosity, and metal and helium abundances. They assume conventional yields for helium and the heavier elements. By comparing with the observational constraints, which may be considered as fairly conservative, it is found that, for an exponentially decreasing star formation rate (SFR) with e-folding time tau, only values between 6 x 10(8) less than similar to tau less than similar to 2 x 10(9) years are allowed together with a very limited range of masses for the initial mass function (IMF). Star formation is allowed for 2 solar mass less than similar to m less than similar to 8 solar mass if tau = 2 x 10(9) years, and for 4 solar mass less than similar to m less than similar to 6 solar mass if tau = 10(9) years. For tau = 6 x 10(8) years, the lower and upper mass limits merge to similar to 5 solar mass. Researchers conclude that, even though the possibility of neutron stars as halo matter may be ruled out, that of white dwarfs may still be a viable hypothesis, though with very stringent constraints on allowed parameters, that merits further consideration.

Ryu, Dongsu↗

NASA NTRS1989

Neutron stars and white dwarfs in galactic halos

Ryu, Dongsu↗

NASA NTRS1989

A model for the distribution of dark matter, galaxies, and the intergalactic medium in a cold dark matter-dominated universe

The spatial distribution of the cold-dark-matter (CDM) and baryonic components of CDM-dominated cosmological models are characterized, summarizing the results of recent theoretical investigations. The evolution and distribution of matter in an Einstein-de Sitter universe on length scales small enough so that the Newtonian approximation is valid is followed chronologically, assuming (1) that the galaxies, CDM, and the intergalactic medium (IGM) are coupled by gravity, (2) that galaxies form by taking mass and momentum from the IGM, and (3) that the IGM responds to the energy input from the galaxies. The results of the numerical computations are presented in extensive graphs and discussed in detail.

Ryu, Dongsu↗