• The Bulletin of The Korean Astronomical Society
• /
• v.39 no.1
• /
• pp.70.2-70.2
• /
• 2014

We developed a three-dimensional (3D) magnetohydrodynamic (MHD) simulation code to reproduce the structure of a solar wind and the propagation of a coronal mass ejection (CME) through it. This code is constructed by a finite volume method based on a total variation diminishing (TVD) scheme using an unstructured grid system (Tanaka 1994). The grid system can avoid the singularity arising in the spherical coordinate system. In this study, we made an improvement of the code focused on the propagation of a CME through a solar wind, which extends a previous work done by Nakamizo et al. (2009). We first reconstructed a solar wind in a steady state from physical values obtained at 50 solar radii away from the Sun via an MHD tomography applied to interplanetary scintillation (IPS) data (Hayashi et al. 2003). We selected CR2057 and inserted a spheromak-type CME (Kataoka et al. 2009) into a reconstructed solar wind. As a result, we found that our simulation well captures the velocity, temperature and density profiles of an observed solar wind. Furthermore, we successfully reproduce the general characteristics of an interplanetary coronal mass ejection (ICME) obtained by the Helios 1/2 spacecraft (R. J. FORSYTH et al. 2006).

• Journal of the Korean Institute of Gas
• /
• v.15 no.2
• /
• pp.47-56
• /
• 2011

Flame propagation velocity is the one ofmainmechanismof the stabilization of triple flame. To quantify the triple flame propagation velocity, Bilger presents the triple flame propagation velocity depending on the mixture fraction gradient, based on the laminar jet flow theory. However, in spite of these many analyses, there was not presented any relation of these variables, triple flame propagation velocity, radius of flame curvature and scalar dissipation rate indirectly. In the present research, we have checked the results of numerical simulation with experiment and numerical analysis and verified the flame propagation velocity with a scalar dissipation rate proposed by Bilger through the numerical simulation. Also we have clarified that flame propagation velocity was depended on the radius of flame curvature and scalar dissipation rate.

• Ocean Systems Engineering
• /
• v.7 no.2
• /
• pp.143-155
• /
• 2017

For complex flexible structures such as nets, the determination of drag forces and its deformation is a challenging task. The accurate prediction of loads on cages is one of the key steps in designing fish farm facilities. The basic physics with a simple cage, can be addressed by the use of experimental studies. However, to design more complex cage system for various environmental conditions, a reliable numerical simulation tool is essential. In this work, the current load on a cage is calculated using a Morison-force model applied at instantaneous positions of equivalent-net modeling. Variations of solidity ratio ($S_n$) of the net and current speed are considered. An equivalent array of cylinders is built to represent the physical netting. Based on the systematic comparisons between the published experimental data for Raschel nets and the current numerical simulations, carried out using the commercial software OrcaFlex, a new formulation for $C_d$ values, used in the equivalent-net model, is presented. The similar approach can also be applied to other netting materials following the same procedure. In case of high solidity ratio and current speed, the hybrid model defines $C_d$ as a function of Re (Reynolds number) and $S_n$ to better represent the corresponding weak diffraction effects. Otherwise, the conventional $C_d$ values depending only on Re can be used with including shielding effects for downstream elements. This new methodology significantly improves the agreement between numerical and experimental data.

• Bulletin of the Korean Space Science Society
• /
• 1993.10a
• /
• pp.15-15
• /
• 1993

Revised Yale Isochrone(Green et al, 1987)과 Lee의 HB simulation data(LEE '||'&'||' Demarque 1990)을 이용하여 Park이 만든 CMD simulation program(Park 1983).에 쌍성을 첨가하면 CMD에 어느 정도의 영향을 나타내게 되는가를 조사하였다. MS, SGB와 RGB는 큰 변화가 없었으나, HB에서는 많은 변화를 나타내었다. 193919 S02401193919 ^x We have observed the emission of HC3N J= 4-3, 5-4 10-9 and 12-11 transitions toward the Sgr B2 central region in an area of 150"XISO" with resolutions of16"-48". The intensities and central velocities of line profiles show significant variations with positions. In contrast to the in tensities of the low J-level transitions which gradually increase from the central source toward the outside region the HC3N emission of the high J-level transition become stronger toward the central radio continuum source MO5. Systematic chainge in the radial velocity of each line profile occurs along north-south direction there are a few peaks in most line profiles, and these indicate that there not multiple velocity components along the line of sight. Distributions of excitation temperah at and column density which were estimated from the excitation calculations show the existence of a small(IxEpc),hot(Tex> SOK) core which contains two tempegatlue peaks at -15" east and north of MDS. The column density of HCaN is (1-3):n1014cm-2. Column density at distant position from MD5 is larger than that in the (:entral region. We have deduced that this hot-core has a mass of 10sR1 which i:s about an order of magnitude larger those obtained by previous studies.previous studies.

• Journal of Korea Multimedia Society
• /
• v.11 no.4
• /
• pp.516-524
• /
• 2008

This paper presents a new modeling technique for the simulation of cloth specific characteristics with a set of hard constraints using an implicit constraint enforcement scheme. A conventional explicit Baumgarte constraint stabilization method has several defects. It requires users to pick problem-dependent coefficients to achieve fast convergence and has inherent stabilization limits. The proposed implicit constraint enforcement method is stable with large time steps, does not require problem dependent feed-back parameters, and guarantees the natural physics-based motion of an object. In addition, its computational complexity is the same as the explicit Baumgarte method. This paper describes a formulation of implicit constraint enforcement and provides a constraint error analysis. The modeling technique for complex components of cloth such as seams, buttons, sharp creases, wrinkles, and prevention of excessive elongation are explained. Combined with an adaptive constraint activation scheme, the results using the proposed method show the substantial enhancement of the realism of cloth simulations with a corresponding savings in computational cost.

• Geomechanics and Engineering
• /
• v.11 no.1
• /
• pp.1-39
• /
• 2016

The finite element method (FEM), discrete element method (DEM), and Discontinuous deformation analysis (DDA) are among the standard numerical techniques applied in computational geo-mechanics. However, in some cases there no possibility for modelling by traditional finite analytical techniques or other mesh-based techniques. The solution presented in the current study as a completely Lagrangian and mesh-free technique is smoothed particle hydrodynamics (SPH). This method was basically applied for simulation of fluid flow by dividing the fluid into several particles. However, several researchers attempted to simulate soil-water interaction, landslides, and failure of soil by SPH method. In fact, this method is able to deal with behavior and interaction of different states of materials (liquid and solid) and multiphase soil models and their large deformations. Soil indicates different behaviors when interacting with water, structure, instrumentations, or different layers. Thus, study into these interactions using the mesh based grids has been facilitated by mesh-less SPH technique in this work. It has been revealed that the fast development, computational sophistication, and emerge of mesh-less particle modeling techniques offer solutions for problems which are not modeled by the traditional mesh-based techniques. Also it has been found that the smoothed particle hydrodynamic provides advanced techniques for simulation of soil materials as compared to the current traditional numerical methods. Besides, findings indicate that the advantages of applying this method are its high power, simplicity of concept, relative simplicity in combination of modern physics, and particularly its potential in study of large deformations and failures.

• Journal of the Korea Society of Computer and Information
• /
• v.26 no.4
• /
• pp.47-53
• /
• 2021

In this paper, we propose a quadtree-based optimization technique that enables fast Super-resolution(SR) computation by efficiently classifying and dividing physics-based simulation data required to calculate SR. The proposed method reduces the time required for quadtree computation by downscaling the smoke simulation data used as input data. By binarizing the density of the smoke in this process, a quadtree is constructed while mitigating the problem of numerical loss of density in the downscaling process. The data used for training is the COCO 2017 Dataset, and the artificial neural network uses a VGG19-based network. In order to prevent data loss when passing through the convolutional layer, similar to the residual method, the output value of the previous layer is added and learned. In the case of smoke, the proposed method achieved a speed improvement of about 15 to 18 times compared to the previous approach.

• Nuclear Engineering and Technology
• /
• v.54 no.4
• /
• pp.1421-1430
• /
• 2022

The computation of the solid angle and the detector efficiency is considering to be one of the most important factors during the measuring process for the radioactivity, especially the cylindrical γ-ray NaI(Tl) detectors nowadays have applications in several fields such as industry, hazardous for health, the gamma-ray radiation detectors grow to be the main essential instruments in radiation protection sector. In the present work, a generic numerical simulation method (NSM) for calculating the efficiency of the γ-ray spectrometry setup is established. The formulas are suitable for any type of source-to-detector shape and can be valuable to determine the full-energy peak and the total efficiencies and P/T ratio of cylindrical γ-ray NaI(Tl) detector setup concerning the truncated conical radioactive source. This methodology is based on estimate the path length of γ-ray radiation inside the detector active medium, inside the source itself, and the self-attenuation correction factors, which typically use to correct the sample attenuation of the original geometry source. The calculations can be completed in general by using extra reasonable and complicate analytical and numerical techniques than the standard models; especially the effective solid angle, and the detector efficiency have to be calculated in case of the truncated conical radioactive source studied condition. Moreover, the (NSM) can be used for the straight calculations of the γ-ray detector efficiency after the computation of improvement that need in the case of γ-γ coincidence summing (CS). The (NSM) confirmation of the development created by the efficiency transfer method has been achieved by comparing the results of the measuring truncated conical radioactive source with certified nuclide activities with the γ-ray NaI(Tl) detector, and a good agreement was obtained after corrections of (CS). The methodology can be unlimited to find the theoretical efficiencies and modifications equivalent to any geometry by essential sufficiently the physical selective considered situation.

• The Bulletin of The Korean Astronomical Society
• /
• v.46 no.2
• /
• pp.47.3-48
• /
• 2021

We studied the large scale dynamo process in a system forced by helical magnetic field. The dynamo process is basically nonlinear, but can be linearized with 𝛼&𝛽 coefficients and large scale magnetic field $\bar{B}$. This is very useful to the investigation of solar (stellar) dynamo. A coupled semi-analytic equations based on statistical mechanics are used to investigate the exact evolution of 𝛼&𝛽. This equation set needs only magnetic helicity ${\bar{H}}_M({\equiv}{\langle}{\bar{A}}{\cdot}{\bar{B}}{\rangle},\;{\bar{B}}={\nabla}{\times}{\bar{A}})$ and magnetic energy ${\bar{E}}_M({\equiv}{\langle}{\bar{B}}^2{\rangle}/2)$. They are fundamental physics quantities that can be obtained from the dynamo simulation or observation without any artificial modification or assumption. 𝛼 effect is thought to be related to magnetic field amplification. However, in reality the averaged 𝛼 effect decreases very quickly without a significant contribution to ${\bar{B}}$ field amplification. Conversely, 𝛽 effect contributing to the magnetic diffusion maintains a negative value, which plays a key role in the amplification with Laplacian ∇²(= - k²) for the large scale regime. In addition, negative magnetic diffusion accounts for the attenuation of plasma kinetic energy E_V(= 〈 U² 〉/2) (U: plasma velocity) when the system is saturated. The negative magnetic diffusion is from the interaction of advective term - U • ∇ B from magnetic induction equation and the helical velocity field. In more detail, when 'U' is divided into the poloidal component U_pol and toroidal one U_tor in the absence of reflection symmetry, they interact with - B • ∇ U and - U • ∇ B from ∇ × 〈 U × B 〉 leading to 𝛼 effect and (negative) 𝛽 effect, respectively. We discussed this process using the theoretical method and intuitive field structure model supported by the simulation result.

• Steel and Composite Structures
• /
• v.44 no.5
• /
• pp.721-740
• /
• 2022

In this paper, an accurate kinematic model has been developed to study the mechanical response of functionally graded (FG) sandwich beams, mainly covering the bending, buckling and free vibration problems. The studied structure with homogeneous hardcore and softcore is considered to be simply supported in the edges. The present model uses a new refined shear deformation beam theory (RSDBT) in which the displacement field is improved over the other existing high-order shear deformation beam theories (HSDBTs). The present model provides good accuracy and considers a nonlinear transverse shear deformation shape function, since it is constructed with only two unknown variables as the Euler-Bernoulli beam theory but complies with the shear stress-free boundary conditions on the upper and lower surfaces of the beam without employing shear correction factors. The sandwich beams are composed of two FG skins and a homogeneous core wherein the material properties of the skins are assumed to vary gradually and continuously in the thickness direction according to the power-law distribution of volume fraction of the constituents. The governing equations are drawn by implementing Hamilton's principle and solved by means of the Navier's technique. Numerical computations in the non-dimensional terms of transverse displacement, stresses, critical buckling load and natural frequencies obtained by using the proposed model are compared with those predicted by other beam theories to confirm the performance of the proposed theory and to verify the accuracy of the kinematic model.