• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2000.11a
• /
• pp.527-533
• /
• 2000

It is investigated that the characteristics of rotating cantilevered composite beam with a breathing crack. In the present study, the crack is modeled as a breathing crack which opens and closes with the motion of the unidirectional graphite-fiber reinforced polyimide beam. For the finite element analysis, the cracked element is modelled by the local flexibility matrix calculated on the basis of fracture mechanics using Castiligano theorem. Rotating beam is considered only transverse bending motion so that the element includes two degrees of freedom per node such as the transverse deflection and slope. The time history and frequency response function of the beam with a breathing crack are studied by Newmark direct time integration method and FFT(Fast Fourier Transform)simulation. Effects of various parameters such as the crack depths, crack locations, ply angles, volume fraction ratios, and rotating speeds of the beam are also studied. Numerical results indicate that it is more reliable to be modelled as a breathing crack than an open crack.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2008.03a
• /
• pp.602-605
• /
• 2008

The anode sheath structure in the hollow anode of an anode-layer type Hall thruster was numerically computed using a fully kinetic 2D3V Particle-in-Cell and Direct Simulation Monte Carlo(PIC-DSMC) code. By treating both ions and electrons as particles, anode surface region, which is electrically non-neutral, was analyzed. In order to analyze in detail, the calculation code was parallelized using Message Passing Interface (MPI). The code successfully simulated the discharge current oscillation. In the low magnetic induction case, ion sheath appears in the anode surface because ionization is enough to maintain the plasma occurs in the anode hollow. As the magnetic induction increases, main ionization region move to outside of the anode. At the same time, anode sheath voltage decreases. In the high magnetic induction case, electron sheath appears on the anode surface periodically because the ionization occurs mainly in the discharge channel. This anode sheath condition shift can be explained using the simple sheath model.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2011.04a
• /
• pp.367-370
• /
• 2011

A high-resolution numerical study is carried out to investigate the transient process of the combustion and the shock-train developments in an ethylene-fueled direct-connect dual-mode scramjet combustor. Air-throttling is then applied at the expansion part of the combustor to provide mass addition to block the flow to subsonic speed, hence to enhance the fuel-air mixing and ignition. Present simulation shows the detailed results for the better understanding of transient processes of the operation regimes in the dual-mode scramjet combustor.

• Journal of Ship and Ocean Technology
• /
• v.9 no.4
• /
• pp.1-10
• /
• 2005

The effect of buoyancy orientation on turbulent channel flow has been investigated using DNS (direct numerical simulation). Grashof number is kept at 9.6 $\times 10^{5}$ while changing the orientation of the buoyancy vector to be parallel or perpendicular to the channel walls. Four study cases can be distinguished during this research namely; streamwise, wall-normal unstable stratification, wall-normal stable stratification and spanwise oriented buoyancy. The driving mean pressure gradient used in all cases is adjusted to keep mass flow rate constant while friction Reynolds number is around 150. At this Grashof number, the skin friction shows decrement in the unstable and stable stratification and increment in the other two cases. Analyses of the changes of flow structure for the four cases are presented highlighting on the mean quantities and second order statistics.

• Journal of Ship and Ocean Technology
• /
• v.12 no.2
• /
• pp.1-15
• /
• 2008

The main idea of this paper introduce stochastic structural parameters and random dynamic excitation directly into the dynamic functional variational formulations, and developed the nonlinear dynamic analysis of a stochastic variational principle and the corresponding stochastic finite element method via the weighted residual method and the small parameter perturbation technique. An interpolation method was adopted, which is based on representing the random field in terms of an interpolation rule involving a set of deterministic shape functions. Direct integration Wilson-${\theta}$ Method was adopted to solve finite element equations. Numerical examples are compared with Monte-Carlo simulation method to show that the approaches proposed herein are accurate and effective for the nonlinear dynamic analysis of structures with random parameters.

• Communications for Statistical Applications and Methods
• /
• v.26 no.2
• /
• pp.175-189
• /
• 2019

The problem of the estimation of quantitative sensitive variable using the item sum technique (IST) on successive occasions has been discussed. IST difference, IST regression, and IST general class of estimators have been proposed to estimate quantitative sensitive variable at the current occasion in two occasion successive sampling. The proposed new estimators have been elaborated under Trappmann et al. (Journal of Survey Statistics and Methodology, 2, 58-77, 2014) as well as Perri et al. (Biometrical Journal, 60, 155-173, 2018) allocation designs to allocate long list and short list samples of IST. The properties of all proposed estimators have been derived including optimum replacement policy. The proposed estimators have been mutually compared under the above mentioned allocation designs. The comparison has also been conducted with a direct method. Numerical applications through empirical as well as simplistic simulation has been used to show how the illustrated IST on successive occasions may venture in practical situations.

• Advances in aircraft and spacecraft science
• /
• v.9 no.3
• /
• pp.217-242
• /
• 2022

This paper proposes a novel time-domain homogenization model combining the viscoelastic constitutive law with Eshelby's inclusion theory-based micromechanics model to predict the mechanical behavior of the particle reinforced composite material. The proposed model is intuitive and straightforward capable of predicting composites' viscoelastic behavior in the time domain. The isotropization technique for non-uniform stress-strain fields and incremental Mori-Tanaka schemes for high volume fraction are adopted in this study. Effects of the imperfectly bonded interphase layer on the viscoelastic behavior on the dynamic mechanical behavior are also investigated. The proposed model is verified by the direct numerical simulation and DMA (dynamic mechanical analysis) experimental results. The proposed model is useful for multiscale analysis of viscoelastic composite materials, and it can also be extended to predict the nonlinear viscoelastic response of composite materials.

• 제어로봇시스템학회:학술대회논문집
• /
• 1998.10a
• /
• pp.100-104
• /
• 1998

This paper describes the real-time implementation of an adaptive controller fur the robotic manipulator. Digital signal processors(DSPs) are special purpose micro-processors that are particularly powerful for intensive numerical computations involving sums and products of variables. TMS320C50 chips are used in implementing real time adaptive control algorithms to provide an enhanced motion for robotic manipulators. In the proposed scheme, adaptation laws are derived from the improved Lyapunov second stability analysis based on the direct adaptive control theory. The adaptive controller consists of an adaptive feedforward controller and feedback controller. The proposed control scheme is simple in structure, fast in computation, and suitable for real-time control. Moreover, this scheme does not require any accurate dynamic modeling, nor values of manipulator parameters and payload. Performance of the adaptive controller is illustrated by simulation and experimental results for a assembling robot.

• Advances in aircraft and spacecraft science
• /
• v.10 no.6
• /
• pp.521-543
• /
• 2023

Reynolds-averaged Navier-Stokes (RANS) models are extensively employed in industrial settings for the purpose of simulating intricate fluid flows. However, these models are subject to certain limitations. Notably, disparities persist in the Reynolds stresses when comparing the RANS model with high-fidelity data obtained from Direct Numerical Simulation (DNS) or experimental measurements. In this work we propose an approach to mitigate these discrepancies while retaining the favorable attributes of the Menter Shear Stress Transport (SST) model, such as its significantly lower computational expense compared to DNS simulations. This strategy entails incorporating an explicit algebraic model and employing a neural network to correct the turbulent characteristic time. The imposition of realizability constraints is investigated through the introduction of penalization terms. The assimilated Reynolds stress model demonstrates good predictive performance in both in-sample and out-of-sample flow configurations. This suggests that the model can effectively capture the turbulent characteristics of the flow and produce physically realistic predictions.

• Structural Engineering and Mechanics
• /
• v.86 no.5
• /
• pp.607-619
• /
• 2023

This study, it was tried to evaluate the asphalt behavior under tensile loading conditions through indirect Brazilian and direct tensile tests, experimentally and numerically. This paper is important from two points of view. The first one, a new test method was developed for the determination of the direct tensile strength of asphalt and its difference was obtained from the indirect test method. The second one, the effects of particle size and loading rate have been cleared on the tensile fracture mechanism. The experimental direct tensile strength of the asphalt specimens was measured in the laboratory using the compression-to-tensile load converting (CTLC) device. Some special types of asphalt specimens were prepared in the form of slabs with a central hole. The CTLC device is then equipped with this specimen and placed in the universal testing machine. Then, the direct tensile strength of asphalt specimens with different sizes of ingredients can be measured at different loading rates in the laboratory. The particle flow code (PFC) was used to numerically simulate the direct tensile strength test of asphalt samples. This numerical modeling technique is based on the versatile discrete element method (DEM). Three different particle diameters were chosen and were tested under three different loading rates. The results show that when the loading rate was 0.016 mm/sec, two tensile cracks were initiated from the left and right of the hole and propagated perpendicular to the loading axis till coalescence to the model boundary. When the loading rate was 0.032 mm/sec, two tensile cracks were initiated from the left and right of the hole and propagated perpendicular to the loading axis. The branching occurs in these cracks. This shows that the crack propagation is under quasi-static conditions. When the loading rate was 0.064 mm/sec, mixed tensile and shear cracks were initiated below the loading walls and branching occurred in these cracks. This shows that the crack propagation is under dynamic conditions. The loading rate increases and the tensile strength increases. Because all defects mobilized under a low loading rate and this led to decreasing the tensile strength. The experimental results for the direct tensile strengths of asphalt specimens of different ingredients were in good accordance with their corresponding results approximated by DEM software.