• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.12
• /
• pp.1555-1569
• /
• 2000

Numerical simulation of fluid flow with moving free surface has been carried out. For the free surface flow, a VOF(Volume of Fluid)-based algorithm utilizing a fixed grid system has been investigated. In order to reduce numerical smearing at the free surface represented on a fixed grid system, a new free surface tracking algorithm based on the donor-acceptor scheme has been presented. Novel features of the proposed algorithm are characterized as two numerical tools; the orientation vector to represent the free surface orientation in each cell and the baby-cell to determine the fluid volume flux at each cell boundary. The proposed algorithm can be easily implemented in any irregular non-uniform grid systems that are usual in finite element method (FEM). Moreover, the proposed algorithm can be extended and applied to the 3-D free surface flow problem without additional efforts. For computation of unsteady incompressible flow, a finite element approximation based on the explicit fractional step method has been adopted. In addition, the SUPG(streamline upwind/Petrov-Galerkin) method has been implemented to deal with convection dominated flows. Combination of the proposed free surface tracking scheme and explicit fractional step formulation resulted in an efficient solution algorithm. Validity of the present solution algorithm was demonstrated from its application to the broken dam and the solitary wave propagation problems.

• Journal of the Korea Institute of Military Science and Technology
• /
• v.7 no.3 s.18
• /
• pp.21-29
• /
• 2004

Numerical analysis of aerodynamic heating for KPSAM is performed using aerodynamic heating model suitable to KPSAM, which has complex flow field resulting from the spike attached to the dome, such as large separation area and the strong shock/boundary layer interaction region around reattachment point on the dome. The aerodynamic heating model is validated and modified through the comparison between the flight test measurement and the thermal analysis results. TFD temperature sensors are installed on the dome to measure surface temperature during the flight. Computation results, obtained from the heat transfer analysis on the sensors, agree well with flight test data. The aerodynamic heating model provides heat transfer rate into surface as a boundary condition of unsteady 1D/axisymmetric thermal analysis on the missile structure. The axisymmetric thermal analysis using FLUENT is more versatile than the 1D analysis and can be applied to the heating problem related with complex structures and multi-dimensional heat transfer problems such as prediction of temperature rise at contact surface of different materials.

• Structural Engineering and Mechanics
• /
• v.76 no.4
• /
• pp.529-540
• /
• 2020

Concurrent topology optimization of macrostructure and microstructure has attracted significant interest due to its high structural performance. However, most of the existing works are carried out under deterministic conditions, the obtained design may be vulnerable or even cause catastrophic failure when the load position exists uncertainty. Therefore, it is necessary to take load position uncertainty into consideration in structural design. This paper presents a computational method for robust concurrent topology optimization with consideration of load position uncertainty. The weighted sum of the mean and standard deviation of the structural compliance is defined as the objective function with constraints are imposed to both macro- and micro-scale structure volume fractions. The Bivariate Dimension Reduction method and Gauss-type quadrature (BDRGQ) are used to quantify and propagate load uncertainty to calculate the objective function. The effective properties of microstructure are evaluated by the numerical homogenization method. To release the computation burden, the decoupled sensitivity analysis method is proposed for microscale design variables. The bi-directional evolutionary structural optimization (BESO) method is used to obtain the black-and-white designs. Several 2D and 3D examples are presented to validate the effectiveness of the proposed robust concurrent topology optimization method.

• Journal of Mechanical Science and Technology
• /
• v.20 no.12
• /
• pp.2292-2314
• /
• 2006

Efficient numerical method is developed for the prediction of aerodynamic noise generation and propagation in low Mach number flows such as aeolian tone noise. The proposed numerical method is based on acoustic/viscous splitting techniques of which acoustic solvers use simplified linearised Euler equations, full linearised Euler equations and nonlinear perturbation equations as acoustic governing equations. All of acoustic equations are forced with immersed surface dipole model which is developed for the efficient computation of aerodynamic noise generation and propagation in low Mach number flows in which dipole source, originating from unsteady pressure fluctuation on a solid surface, is known to be more efficient than quadrupole sources. Multi-scale overset grid technique is also utilized to resolve the complex geometries. Initially, aeolian tone from single cylinder is considered to examine the effects that the immersed surface dipole models combined with the different acoustic governing equations have on the overall accuracy of the method. Then, the current numerical method is applied to the simulation of the aeolian tones from twin cylinders aligned perpendicularly to the mean flow and separated 3 diameters between their centers. In this configuration, symmetric vortices are shed from twin cylinders, which leads to the anti-phase of the lift dipoles and the in-phase of the drag dipoles. Due to these phase differences, the directivity of the fluctuating pressure from the lift dipoles shows the comparable magnitude with that from the drag dipoles at 10 diameters apart from the origin. However, the directivity at 100 diameters shows that the lift-dipole originated noise has larger magnitude than, but still comparable to, that of the drag-dipole one. Comparison of the numerical results with and without mean flow effects on the acoustic wave emphasizes the effects of the sheared background flows around the cylinders on the propagating acoustic waves, which is not generally considered by the classic acoustic analogy methods. Through the comparison of the results using the immersed surface dipole models with those using point sources, it is demonstrated that the current methods can allow for the complex interactions between the acoustic wave and the solid wall and the effects of the mean flow on the acoustic waves.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.22 no.6
• /
• pp.84-93
• /
• 2018

Three-dimensional (3D) numerical simulations have been carried out to study how coolant injection conditions influence the cooling efficiency and projectile ejection performance in a mixture-gas ejection system (or gas-steam launch system). The 3D single-phase computational model was verified using a 1D model constructed with reference to the previous research and then a two-phase flow computation simulating coolant injection on to hot gas was performed using a DPM (Discrete Phase Model). As a result of varying the coolant flow rate and number of injection holes, cooling efficiency was improved when the number of injection holes were increased. In addition, the change of the coalescence frequency and spatial distribution of coolant droplets caused by the injection condition variation resulted in a change of the droplet diameter, affecting the evaporation rate of coolant. The evaporation was found to be a critical factor in the design optimization of the ejection system by suppressing the pressure drop while the temperature decreases inside the breech.

• Wind and Structures
• /
• v.14 no.5
• /
• pp.413-434
• /
• 2011

Without output covariance estimation, one reference-based Stochastic Subspace Technique (SST) for extracting modal parameters and flutter derivatives of bridge deck is developed and programmed. Compared with the covariance-driven SST and the oscillation signals incurred by oncoming or signature turbulence that adopted by previous investigators, the newly-presented identification scheme is less time-consuming in computation and a more desired accuracy should be contributed to high-quality free oscillated signals excited by specific initial displacement. The reliability and identification precision of this technique are confirmed by a numerical example. For the 3-DOF sectional models of Sutong Bridge deck (streamlined) and Suramadu Bridge deck (bluff) in wind tunnel tests, with different wind velocities, the lateral bending, vertical bending, torsional frequencies and damping ratios as well as 18 flutter derivatives are extracted by using SST. The flutter derivatives of two kinds of typical decks are compared with the pseudo-steady theoretical values, and the performance of $H_1{^*}$, $H_3{^*}$, $A_1{^*}$, $A_3{^*}$ is very stable and well-matched with each other, respectively. The lateral direct flutter derivatives $P_5{^*}$, $P_6{^*}$ are comparatively more accurate than other relevant lateral components. Experimental procedure seems to be more critical than identification technique for refining the estimation precision.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.23 no.3
• /
• pp.28-34
• /
• 2019

The thrust control calculation according to the operation of plug nozzle for cold air test and the analysis of the flow characteristics of the pintle shape and operation pressure are performed. The numerical computation was verified by comparing the flow structure and the coefficient of thrust with the experimental data. It was confirmed that the nozzle throat was formed at the design position on each pintle shape, and thrust control up to 1:8 was achieved only by the stroke change. Finally, although the aerospike nozzle is autonomous, it is unfavorable in the under_expansion condition, if it is designed for a very low nozzle pressure ratio.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.29 no.9 s.240
• /
• pp.997-1005
• /
• 2005

The extinction characteristics of low strain rate normal gravity (1-g) nonpremixed methane-air flames were studied numerically and experimentally. A time-dependent axisymmetric two-dimensional (2D) model considering buoyancy effects and radiative heat transfer was developed to capture the structure and extinction limits of 1-g flames. One-dimensional (1D) computations were also conducted to provide information on 0-g flames. A 3-step global reaction mechanism was used in both the 1D and 2D computations to predict the measured extinction limit and flame temperature. A specific maximum heat release rate was introduced to quantify the local flame strength and to elucidate the extinction mechanism. Overall fractional contribution by each term in the energy equation to the heat release was evaluated to investigate the multi-dimensional structure and radiative extinction of 1-g flames. Images of flames were taken for comparison with the model calculation undergoing extinction. The two-dimensional numerical model was validated by comparing flame temperature profiles and extinction limits with experiments and ID computation results. The 2D computations yielded insight into the extinction mode and flame structure of 1-g flames. Two combustion regimes depending on the extinction mode were identified. Lateral heat loss effects and multi-dimensional flame structure were also found. At low strain rates of 1-g flame ('Regime A'), the flame is extinguished from the weak outer flame edge, which is attributed to multi-dimensional flame structure and flow field. At high strain rates, ('Regime B'), the flame extinction initiates near the flame centerline due to an increased diluent concentration in reaction zone, which is the same as the extinction mode of 1D flame. These two extinction modes could be clearly explained with the specific maximum heat release rate.

• Journal of Welding and Joining
• /
• v.11 no.3
• /
• pp.34-47
• /
• 1993

Finite element models were developed for thermal and residual stress analysis for the specific welding problems. They were used to evaluate the effectiveness of the various welding heat input models, such as ramp heat input function and lumped pass models. Through the parametric studies, thermal-mechanical modeling sensitivity to the ramp function and lumping techniques was determined by comparing the predicted results with experimental data. The kinetics for residual stress formation during welding can be developed by iteration of various proposed mechanisms in the parametric study. A ramp heat input function was developed to gradually apply the heat flux with variable amplitude to the model. This model was used to avoid numerical convergence problems due to an instantaneous increase in temperature near the fusion zone. Additionally, it enables the model to include the effect of a moving arc in a two-dimensional plane. The ramp function takes into account the variation in the out of plane energy flow in a 2-D model as the arc approaches, travels across, and departs from each plane under investigation. A lumped pass model was developed to reduce the computation cost in the analysis of multipass welds. Several weld passes were assumed as one lumped pass in this model. Recommendations were provided about ramp lumping techniques and the optimum number of weld passes that can be combined into a single thermal input.

• Computers and Concrete
• /
• v.30 no.4
• /
• pp.277-287
• /
• 2022

Due to the scarcity of extortionate experimental data, fatigue failure of the reinforced concrete (RC) element might be achieved economically adopting nonlinear finite element (FE) analysis as an alternative approach. However, conventional implicit dynamic analysis is expensive, quasi-static method overlooks interaction effects and inertia, direct cyclic analysis computes stabilized responses. Apart from this, explicit dynamic analysis may provide a numerical operating system for factual long-term responses. The study explores the fatigue behavior based on a simplified explicit dynamic solution employing nonlinear time domain analysis. Among fourteen RC beams, one beam is selected to validate under static loading, one under fatigue with the experimental study and other twelve to check the detail fatigue behavior. The SWOT (Strength, Weakness, Opportunities, Threats) analysis has been carried out to pinpoint the detail scenario in the adoption of numerical approach as an alternative to the experimental study. Excellent agreement of FE and experimental results is seen. The 3D nonlinear RC beam model at service fatigue limits is truthful to be used as an expedient contrivance to envisage the precise fatigue behavior. The simplified analysis approach for RC beam under fatigue offers savings in computation to predict responses providing acceptable accuracy rather than the complicated laboratory investigation. At higher frequency, the flexural failure occurs a bit earlier gradually compared to the repeated loading case of lower frequency. The deflection increases by 6%-10% at the end of first cycle for beams with increasing frequency of cyclic loading. However, at the end of fatigue loading, greater deflection occur earlier for higher load range because of more rapid stiffness degradation. For higher frequency, a slight boost in concrete compressive strains at an initial stage of loading has been seen indicating somewhat stepper increment. Stiffness degradation in larger loading cycle at same duration escalates the upsurge of the rate of strain in case of higher frequency.