• Advances in aircraft and spacecraft science
• /
• v.6 no.4
• /
• pp.333-348
• /
• 2019

Low fuel regression rate is the main drawback of hybrid rocket which should be overcome. One of the improvement techniques to this problem is usage of a solid fuel grain with a complicated geometry port, which has been promoted owing to the recent development of additive manufacturing technologies. In the design of a hybrid rocket fuel grain with a complicated geometry port, the understanding of fuel regression behavior is very important. Numerical investigations of fuel regression behavior requires a capturing method of solid fuel surface, i.e. gas-solid interface. In this study, level set method is employed as such a method and the preliminary numerical tool for capturing a hybrid rocket solid fuel surface is developed. At first, to test the adequacy of the numerical modeling, the simulation results for circular port are compared to the experimental results in open literature. The regression rates and oxidizer to fuel ratios show good agreements between the simulations and the experiments, after passing enough time. However, during the early period of combustion, there are the discrepancies between the simulations and the experiments, owing to transient phenomena. Second, the simulations of complicated geometry ports are demonstrated. In this preliminary step, a star shape is employed as complicated geometry of port. The slot number effect in star port is investigated. The regression rate decreases with increasing the slot number, except for the star port with many slots (8 slots) in the latter half of combustion. The oxidizer to fuel ratio increases with increasing the slot number.

• Wind and Structures
• /
• v.34 no.6
• /
• pp.499-510
• /
• 2022

Computational Wind Engineering has rapidly grown in the last decades and it is currently reaching a relatively mature state. The prediction of wind loading by means of numerical simulations has been proved effective in many research studies and applications to design practice are rapidly spreading. Despite such success, caution in the use of simulations for wind loading assessment is still advisable and, indeed, required. The computational burden and the know-how needed to run high-fidelity simulations is often unavailable and the possibility to use simplified models extremely attractive. In this paper, the applicability of some well-known 2D unsteady RANS models, particularly the k-ω SST, in the aerodynamic characterization of extruded bodies with bluff sections is investigated. The main focus of this paper is on the drag coefficient prediction. The topic is not new, but, in the authors' opinion, worth a careful revisitation. In fact, despite their great technical relevance, a systematic study focussing on sections which manifest a fully detached flow configuration has been overlooked. It is here shown that the considered 2D RANS exhibit a pathological behaviour, failing to reproduce the transition between reattached and fully detached flow regime.

• Coupled systems mechanics
• /
• v.12 no.5
• /
• pp.409-418
• /
• 2023

A coupled algorithm is proposed which first considers the creation of porous structure of the material and then the simulations of response of mechanical components with porous structure to a variable load history. The simulations are carried out by the Prandtl operator approach in the finite element method (FEM) which enables structural simulations of mechanical components subjected to variable thermomechanical loads. Temperature-dependent material properties and multilinear kinematic hardening of the material can be taken into account by this approach. Several simulations are then performed for a tensile-compressive specimen made of a generic porous structure and mechanical properties of Aluminium alloy AlSi9Cu3. Variable mechanical load history has been applied to the specimens under constant temperature conditions. Comparison of the simulation results shows a considerable elastoplastic stress-strain response in the vicinity of pores whilst the surface of the gauge-length of the specimen remains in the elastic region of the material. Moreover, the distribution of the pore sizes seems more influential to the stress-strain field during the loading than their radial position in the gauge-length.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.28 no.3
• /
• pp.159-164
• /
• 2024

Recent earthquakes in Korea, like Gyeongju and Pohang, have highlighted the need for accurate seismic hazard assessment. The lack of substantial ground motion data necessitates stochastic simulation methods, traditionally used with a simplistic point-source assumption. However, as earthquake magnitude increases, the influence of finite faults grows, demanding the adoption of finite faults in simulations for accurate ground motion estimates. We analyzed variations in simulated ground motions with and without the finite fault method for earthquakes with magnitude (Mw) ranging from 5.0 to 7.0, comparing pseudo-spectral acceleration. We also studied how slip distribution and hypocenter location affect simulations for a virtual earthquake that mimics the Gyeongju earthquake with Mw 5.4. Our findings reveal that finite fault effects become significant at magnitudes above Mw 5.8, particularly at high frequencies. Notably, near the hypocenter, the virtual earthquake's ground motion significantly changes using a finite fault model, especially with heterogeneous slip distribution. Therefore, applying finite fault models is crucial for simulating ground motions of large earthquakes (Mw ≥ 5.8 magnitude). Moreover, for accurate simulations of actual earthquakes with complex rupture processes having strong localized slips, incorporating finite faults is essential even for more minor earthquakes.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.9 no.4
• /
• pp.446-459
• /
• 2017

Naval ships are assigned many and varied missions. Their performance is critical for mission success, and depends on the specifications of the components. This is why performance analyses of naval ships are required at the initial design stage. Since the design and construction of naval ships take a very long time and incurs a huge cost, Modeling and Simulation (M & S) is an effective method for performance analyses. Thus in this study, a simulation core is proposed to analyze the performance of naval ships considering their specifications. This simulation core can perform the engineering level of simulations, considering the mathematical models for naval ships, such as maneuvering equations and passive sonar equations. Also, the simulation models of the simulation core follow Discrete EVent system Specification (DEVS) and Discrete Time System Specification (DTSS) formalisms, so that simulations can progress over discrete events and discrete times. In addition, applying DEVS and DTSS formalisms makes the structure of simulation models flexible and reusable. To verify the applicability of this simulation core, such a simulation core was applied to simulations for the performance analyses of a submarine in an Anti-SUrface Warfare (ASUW) mission. These simulations were composed of two scenarios. The first scenario of submarine diving carried out maneuvering performance analysis by analyzing the pitch angle variation and depth variation of the submarine over time. The second scenario of submarine detection carried out detection performance analysis by analyzing how well the sonar of the submarine resolves adjacent targets. The results of these simulations ensure that the simulation core of this study could be applied to the performance analyses of naval ships considering their specifications.

• Wind and Structures
• /
• v.11 no.2
• /
• pp.75-96
• /
• 2008

A better understanding of tornado-induced wind loads is needed to improve the design of typical structures to resist these winds. An accurate understanding of the loads requires knowledge of near-ground tornado winds, but observations in this region are lacking. The first goal of this study was to verify how well a CFD model, when driven by far field radar observations and laboratory measurements, could capture the flow characteristics of both full scale and laboratory-simulated tornadoes. A second goal was to use the model to examine the sensitivity of the simulations to various parameters that might affect the laboratory simulator tornado. An understanding of near-ground winds in tornadoes will require coordinated efforts in both computational and physical simulation. The sensitivity of computational simulations of a tornado to geometric parameters and surface roughness within a domain based on the Iowa State University laboratory tornado simulator was investigated. In this study, CFD simulations of the flow field in a model domain that represents a laboratory tornado simulator were conducted using Doppler radar and laboratory velocity measurements as boundary conditions. The tornado was found to be sensitive to a variety of geometric parameters used in the numerical model. Increased surface roughness was found to reduce the tangential speed in the vortex near the ground and enlarge the core radius of the vortex. The core radius was a function of the swirl ratio while the peak tangential flow was a function of the magnitude of the total inflow velocity. The CFD simulations showed that it is possible to numerically simulate the surface winds of a tornado and control certain parameters of the laboratory simulator to influence the tornado characteristics of interest to engineers and match those of the field.

• Structural Engineering and Mechanics
• /
• v.74 no.6
• /
• pp.771-787
• /
• 2020

The paper provides an inside look into experimental measurements, followed by numerical simulations and their related uncertainties. The goal of the paper is to present findings related to blast loading and the handling of defects that are inherent in experiments. Very often it might seem that experiments are simplified reflections of real-life conditions. In most cases this is true, but there is a good reason for that. The more complex an experiment is, the larger the amount of uncertainties that can be expected. This especially applies when the blast loading of concrete is the subject of research. When simulations fail to reproduce the results of experimental measurements, it does not necessarily mean there is something wrong with the numerical model. The problem could be missing information. Put differently, the numerical simulation may lack information that seemed irrelevant with regard to the experiment. In the presented case, a reference simulation with a proven material model unexpectedly failed to replicate the results of an experiment where concrete slabs were exposed to blast loading. This resulted in a search for possible unknowns. When all of the uncertainties were examined, the missing information turned out to be the orientation of the charge to the concrete slab. Since the experiment was burdened with error, a sensitivity study had to take place so the influence of this factor could be better understood. The findings point to the fact that even the smallest defect during experiments must somehow be taken into account when designing numerical simulations. Otherwise, the simulations are not correlated to the experiments, but merely to some expectations.

• Journal of the Korean earth science society
• /
• v.30 no.5
• /
• pp.630-645
• /
• 2009

The tropopause fold event that took place on January 1, 1997 over mid-region on the Korean Peninsula is examined by means of a numerical simulation based on a Mesoscale Model (MM5). The purpose of this paper is to investigate the effects of moisture in reinforcing a tropopause fold linked to an explosive cyclone. Two types of simulations were carried out; 1) simulations for moist conditions in which full physical and dynamic processes are considered and 2) simulations for dry conditions in which cumulus parameterization and cloud microphysics process are excluded. The results of the moist condition simulations demonstrate that the intensity of the central pressure of the cyclone was overestimated compared with the observed values and that the location of the center and the pressure deepening rates (-17 hPa/12 hr) complied with the observed values. The potential vorticity (PV) anomaly on the isentropic surface at 305 K continued to move in a southeast direction on January 1, 1997 and thus created a single tube of tropopause fold covering the northern and the middle area of the Korean Peninsula and reaching the ground surface at 0300 UTC and 0600 UTC. The results of the dry condition simulations show that the tropopause descended to 500 and 670 hPa in 0300 and 0600 UTC, respectively at the same location for the moist condition simulation; however, there was no deep tropopause fold observed. A comparison of the simulated data between the moist and the dry conditions suggests that a deep tropopause fold should happen when there is sufficient moist in the atmosphere and significantly large PV in the lower atmosphere pulls down the upper atmosphere rather than when the tropopause descends itself due to dynamic causes. Thus, it is estimated that moisture in the atmosphere should have played a crucial role in a deep tropopause fold process.

• Atmosphere
• /
• v.29 no.5
• /
• pp.551-566
• /
• 2019

Quantitative understanding of a random error that is associated with Lagrangian particle dispersion modeling is a prerequisite for backward-in-time mode simulations. This study aims to quantify the random error of the WRF-FLEXPART model and suggest an optimum number of the Lagrangian particles for backward-in-time simulations over the Seoul metropolitan area. A series of backward-in-time simulations of the WRF-FLEXPART model has conducted at two receptor points by changing the number of Lagrangian particles and the relative error, as a quantitative indicator of random error, is analyzed to determine the optimum number of the release particles. The results show that in the Seoul metropolitan area a 1-day Lagrangian transport contributes 80~90% in residence time and ~100% in atmospheric enhancement of carbon monoxide. The relative errors in both the residence time and the atmospheric concentration enhancement are larger when the particles release in the daytime than in the nighttime, and in the inland area than in the coastal area. The sensitivity simulations reveal that the relative errors decrease with increasing the number of Lagrangian particles. The use of small number of Lagrangian particles caused significant random errors, which is attributed to the random number sampling process. For the particle number of 6000, the relative error in the atmospheric concentration enhancement is estimated as -6% ± 10% with reduction of computational time to 21% ± 7% on average. This study emphasizes the importance of quantitative analyses of the random errors in interpreting backward-in-time simulations of the WRF-FLEXPART model and in determining the number of Lagrangian particles as well.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.16 no.4
• /
• pp.345-351
• /
• 2003

Due to increasing legal and market demands for safety in the automotive design process, the design of integrated seat is important more and mote because it should satisfy the conflict between stronger and lower weight for safety and environmental demands. In this study for crash simulations, the numerical simulations have been carried out using the explicit finite element program LS-Dyna according to the FMVSS 210 standard for safety test of seat. Since crash simulations are very time-consuming and a series of simulations that does not lead to a better result is very costly, the optimization method must be both efficient and reliable. As a result of that, statistical approaches such as design of experiments and response surface model have been successfully implemented to reduce time-consuming LS-Dyna simulations and optimize the safety and environmental demands together with nonlinear optimization algorithm. Design of experiments is used lot exploring the design space of maximum displacement and total weight and for building response surface models in order to minimize the maximum displacement and total weight of integrated seat.