• The Bulletin of The Korean Astronomical Society
• /
• v.45 no.1
• /
• pp.55.1-55.1
• /
• 2020

In an attempt to investigate the nonlinear dynamics such as shock, shear, and turbulence associated with ultra-relativistic jets, we develop a new relativistic hydrodynamics (RHD) code based on the weighted essentially non-oscillatory (WENO) scheme. It is a 5th-order accurate, finite-difference scheme, which has been widely used for solving hyperbolic systems of conservation equations. The code is parallelized with MPI and OpenMP. Through an extensive set of tests, the accuracy and efficiency of different WENO reconstructions, and different time discretizations are assessed. Different implementations of the equation of state (EOS) for relativistic fluid are incorporated, As the fiducial setup for simulations of ultra-relativistic jets, we adopt the EOS in Ryu et al. (2006) to treat arbitrary adiabatic index of relativistic fluid, the WENO-Z reconstructions to minimize numerical dissipation without loss of stability, and the strong stability preserving Runge-Kutta (SSPRK) method to achieve stable time stepping with large CFL numbers. In addition, the code includes a high-order flux averaging along the transverse directions for multi-dimensional problems, and the modified eigenvalues for the acoustic modes to effectively control the carbuncle instability. We find that the new code performs satisfactorily simulations of ultra-relativistic jets.

• Journal of The Korean Astronomical Society
• /
• v.38 no.2
• /
• pp.129-133
• /
• 2005

The temperature (T) and entropy (S) fields of baryonic gas, or intergalactic medium (IGM), in the ACDM cosmology are analyzed using simulation samples produced by a hybrid cosmological hydrodynamic/N-body code based on the weighted essentially non-oscillatory scheme. We demonstrate that, in the nonlinear regime, the dynamical similarity between the IGM and dark matter will be broken in the presence of strong shocks in the IGM. The heating and entropy production by the shocks breaks the IGM into multiple phases. The multiphase and non-Gaussianity of the IGM field would be helpful to account for the high-temperature and high-entropy gas observed in groups and clusters with low-temperature IGM observed by Ly$\alpha$ forest lines and the intermittency observed by the spikes of quasi-stellar object's absorption spectrum.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.17 no.1 s.118
• /
• pp.41-47
• /
• 2007

Scattering fields of two dimensional acoustic waves by a circular cylinder are investigated. The present numerical approach for the acoustic scattering problem has difficulties of numerical robustness, long-time stability and suitability of far-field boundary treatments. The time-dependent periodic acoustic source is used to analyze Interference patterns between incident waves and waves reflected by the cylinder. Characteristic boundary algorithms coupled with 4th order Modified-Flux-Approach ENO(essentially non-oscillatory) schemes are employed in generalized coordinates to examine the effect of the wane frequency on the interference patterns. Non-reflecting boundary conditions, which is crustal for accurate computations of aeroacoustic problems, are used not to contaminate scattering fields by reflected waves at the outer boundary. Computed scattering fields show the circumferential acoustic modes generated by interacting between acoustic sources and scattered waves. At a lower frequency, the wave passes almost straight through the cylinder without Interacting with circular cylinder. Simulation results are presented and compared with the analytic solution. Computed RMS-pressure distribution on the cylinder wall is good agreement with exact solution.

• Journal of Korea Water Resources Association
• /
• v.34 no.1
• /
• pp.45-55
• /
• 2001

Transcritical flow is a term intended to denote the existence of both supercritical and subcritical flows within a computational domain. The major problems that need to be addressed while modeling transcritical flows include handling the differing features of signal propagation in subcritical and supercritical flow regions and maintaining conservation. The present study proposes the implicit ENO method as a high-resolution scheme for transcritical flow. This implicit ENO scheme is based on the ENO method, a new class of uniformly high-order-accurate essentially non-oscillatory implicit scheme, which has the advantage of unconditional stability. The implicit ENO scheme has not been used for the transcritical flow in open channel until now. As a result of application to the hypothetical dam-break flow, the implicit ENO scheme was ploved to produce accurate results with good robustness even though in the case of verb strong shock wave.

• Proceedings of the Acoustical Society of Korea Conference
• /
• autumn
• /
• pp.389-392
• /
• 2001

To investigate the generation mechanism of the shock-associated noise, an underexpanded supersonic jet from an axisymmetic nozzle is simulated under the conditions of the Nozzle exit Mach number of 2 and the exit pressure ratio of Pe/Pe =1.5. The present simulation is performed based on the high-order accuracy and high-resolution ENO (Essentially Non-Oscillatory) scheme to capture the time-dependent flow structure representing the sound source. It was found that the shock-associated noise is generated by the weak interaction between the downstream propagating large turbulence structures of the jet flow and the quasi-periodic shock cell structure during the one is passing through the other. The directivity of propagating waves to the upstream is clearly shown in the visualization of pressure field. It is shown that the present calculation of the centerline pressure distribution is in fare agreement with the experimental data at the location of first shock cell.

• International Journal of Aeronautical and Space Sciences
• /
• v.5 no.2
• /
• pp.18-27
• /
• 2004

An axisymmetric supersonic jet is simulated at a Mach number 2.1 and a Reynolds numberof 70000 to identify the mechanism of Mach wave generation and radiation from the jet. In orderto provide the near-field radiated sound directly and resolve the large-scale vortices highly.high-resolution essentially non-oscillatory(ENO) scheme, which is one of the ComputationalAeroAcoustics(CAA) techniques, is newly employed. Perfectly expanded supersonic jet is selectedas a target to see pure shear layer growth and Mach wave radiation without effect of change injet cross section due to expansion or shock wave generated at nozzle exit. The sound field ishighly directional and dominated by Mach waves generated near the end of potential core. Thenear field sound pressure levels as well as the aerodynamic properties of the jet, such asmean-flow parameters are in fare agreement with experimental data.

• The Bulletin of The Korean Astronomical Society
• /
• v.46 no.1
• /
• pp.36.3-36.3
• /
• 2021

We study ultra-relativistic jets on several tens kpc scales through three-dimensional relativistic hydrodynamic (RHD) simulations using a new RHD code based on the weighted essentially non-oscillatory (WENO) scheme. Utilizing the high-resolution and high-accuracy capabilities of the new code, we especially explore the structures and energetics of nonlinear flows, such as shocks, turbulence, velocity shear in different parts of jets. We find that the mildly relativistic shocks which form in the jet backflow are most effective for the shock dissipation of the jet energy, while the turbulent dissipation is largest either in the backflow or in the shocked ICM, depending on the jet parameter. The velocity shear is strongest across the jet flow to the cocoon boundary. Our results should have important implications for the studies of high-energy cosmic-ray production in radio galaxies.

• International Journal of Aeronautical and Space Sciences
• /
• v.16 no.2
• /
• pp.295-310
• /
• 2015

In the present study, an unstructured mixed mesh flow solver was used to conduct a numerical prediction of the aerodynamic performance of the S-76 rotor in hover. For the present mixed mesh methodology, the near-body flow domain was modeled by using body-fitted prismatic/tetrahedral cells while Cartesian mesh cells were filled in the off-body region. A high-order accurate weighted essentially non-oscillatory (WENO) scheme was employed to better resolve the flow characteristics in the off-body flow region. An overset mesh technique was adopted to transfer the flow variables between the two different mesh regions, and computations were carried out for three different blade configurations including swept-taper, rectangular, and swept-taper-anhedral tip shapes. The results of the simulation were compared against experimental data, and the computations were also made to investigate the effect of the blade tip Mach number. The detailed flow characteristics were also examined, including the tip-vortex trajectory, vortex core size, and first-passing tip vortex position that depended on the tip shape.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2004.03a
• /
• pp.62-68
• /
• 2004

A comprehensive numerical analysis has been carried out for both non-reacting and reacting flows in a scramjet engine combustor with and without a cavity. The theoretical formulation treats the complete conservation equations of chemically reacting flows with finite-rate chemistry of hydrogen-air. Turbulence closure is achieved by means of a k-$\omega$ two-equation model. The governing equations are discretized using a MUSCL-type TVD scheme, and temporally integrated by a second-order accurate implicit scheme. Transverse injection of hydrogen is considered over a broad range of injection pressure. The corresponding equivalence ratio of the overall fuel/air mixture ranges from 0.167 to 0.50. The work features detailed resolution of the flow and flame dynamics in the combustor, which was not typically available in most of the previous studies. In particular, the oscillatory flow characteristics are captured at a scale sufficient to identify the .underlying physical mechanisms. Much of the flow unsteadiness is related not only to the cavity, but also to the intrinsic unsteadiness in the flow-field. The interactions between the unsteady flow and flame evolution may cause a large excursion of flow oscillation. The roles of the cavity, injection pressure, and heat release in determining the flow dynamics are examined systematically.

• Nuclear Engineering and Technology
• /
• v.53 no.7
• /
• pp.2084-2094
• /
• 2021

Delta-form-based methods for solving high order spatial discretization schemes are introduced into the reactor S_N transport equation. Due to the nature of the delta-form, the final numerical accuracy only depends on the residuals on the right side of the discrete equations and have nothing to do with the parts on the left side. Therefore, various high order spatial discretization methods can be easily adopted for only the transport term on the right side of the discrete equations. Then the simplest step or other robust schemes can be adopted to discretize the increment on the left hand side to ensure the good iterative convergence. The delta-form framework makes the sweeping and iterative strategies of various high order spatial discretization methods be completely the same with those of the traditional S_N codes, only by adding the residuals into the source terms. In this paper, the flux limiter method and weighted essentially non-oscillatory scheme are used for the verification purpose to only show the advantages of the introduction of delta-form-based solving methods and other high order spatial discretization methods can be also easily extended to solve the S_N transport equations. Numerical solutions indicate the correctness and effectiveness of delta-form-based solving method.