• Structural Engineering and Mechanics
• /
• 제72권2호
• /
• pp.203-216
• /
• 2019

Reliability analysis of composite structures considering random variation of involved parameters is quite important as composite materials revealed large statistical variations in their mechanical properties. The reliability analysis of such structures by the first order reliability method (FORM) and Monte Carlo Simulation (MCS) based approach involves repetitive evaluations of performance function. The response surface method (RSM) based metamodeling technique has emerged as an effective solution to such problems. In the application of metamodeling for uncertainty quantification and reliability analysis of composite structures; the finite element model is usually formulated by either classical laminate theory or first order shear deformation theory. But such theories show significant error in calculating the structural responses of composite structures. The present study attempted to apply the RSM based MCS for reliability analysis of composite shell structures where the surrogate model is constructed using higher order shear deformation theory (HSDT) of composite structures considering the uncertainties in the material properties, load, ply thickness and radius of curvature of the shell structure. The sensitivity of responses of the shell is also obtained by RSM and finite element method based direct approach to elucidate the advantages of RSM for response sensitivity analysis. The reliability results obtained by the proposed RSM based MCS and FORM are compared with the accurate reliability analysis results obtained by the direct MCS by considering two numerical examples.

• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 2002년도 가을 학술발표회 논문집
• /
• pp.583-590
• /
• 2002

Point collocation method based on the fast derivatives approximation of meshfree shape function is applied to solid mechanics in this study. Enhanced meshfree approximation with approximated derivative of shape function is reviewed, and formulation of linear elastic solid mechanics by point collocation method is presented. It implies that governing equation of solid mechanics with strong form is directly formulated without no numerical integration cells or grid. The regularity of weight function is not required due to a use of approximated derivative, so we propose the exponential type weight function that is discontinuous in first derivative. The convergence and stability of the proposed method is verified by passing the generalized patch test. Also, the efficiency and applicability of the proposed method in solid mechanics is verified by solving types of solid problems. Numerical results show that not only a use of proposed weight function leads lower error and higher convergence rate than that of the conventional weight functions, but also the improved collocation method with derivative approximation enables to compute the derivatives of shape function very fast and accurately enough to replace the classical direct derivative calculation.

• 대한조선학회논문집
• /
• 제55권5호
• /
• pp.394-400
• /
• 2018

Various numerical methods have been adopted for indoor noise assessments of ship plant. Acoustical radiosity method is one of the high frequency approaches for acoustic field analysis, which assumes diffuse reflections by boundaries so that it could be efficiently applied to the acoustically diffused indoor space noise analysis. In this study, an acoustic field analysis program has been developed based on radiosity method, which could apply for acoustically large enclosures such as ship's indoor space. For this purpose, the procedure of the acoustical radiosity method has been summarized and implemented to an acoustic field analysis program using MATLAB. Numerical example for a rectangular indoor space has investigated validity of the implemented program. Steady state sound pressure levels calculated for a continuous acoustic source signal have shown good agreement with those by other solutions such as an analytic solution and a ray tracing method. Instantaneous sound pressure levels calculated for an impulsive acoustic signal have provided the clues of direct/reflected acoustic field and reverberation time.

• 한국철도학회:학술대회논문집
• /
• 한국철도학회 2011년도 정기총회 및 추계학술대회 논문집
• /
• pp.1594-1600
• /
• 2011

본 연구에서는 가장 최근에 NIST(National Institute of Standards and Technology)에서 LES(large eddy simulation)기법을 사용하여 개발한 3D 화재유동 해석용 FDS5의 중요한 기능 중 Werner-Wengle wall law의 성능을 평가하기 위하여 평행 평판에서의 유동을 조사하였다. 격자 형성을 위하여 $y^+$ 값은 11 이상으로 유지하도록 하였으며, 총 사용격자는 $32{\times}32{\times}32$를 사용하였다. 입구와 출구에서는 반복(periodic) 경계조건이 주어졌고, 양측면에서는 대칭(symmetry) 경계조건이 주어졌다. 충분히 발달된 난류에서의 유동조건을 조사하기 위하여 Re=10,700을 사용하였다. 시뮬레이션으로 구하여진 결과는 DNS(Direct Numerical Simulation) 결과 및 이론값과 비교하여 분석하였다. 또한 FDS 결과를 RANS의 난류 유동 해석 결과와도 비교하였다.

• 물과 미래
• /
• 제21권3호
• /
• pp.291-291
• /
• 1988

2차원 흐름이 원주주위를 지날 때 발생하는 흐름의 변화가 기본방정식인 연속방정식과 운동량방정식에 의하여 수치적으로 해석된다. 수치해석 과정은 기본방정식에 유함수, 와도 및 흐름의 특성을 나타내는 무차원 매개변수를 도입하여 무차원 유함수-와소수송식을 유도한 후, successive over relaxation scheme과 alternating direct implicit scheme으로 수행된다. 수치실험은 레이놀즈수 125-275를 기존의 수치해석에서는 주로 수치실험 결과와 비교한다. 원주표면의 압력을 구하는 방법에 있어서 기존의 수치해석에서는 주로 방사 운동량방정식만을 사용하였으나, 본 논문에서는 기존의 방법외에 방사 운동량방정식 및 접선 운동량방정식에 의해 압력을 계한하고, 두 값을 비교하여 레이놀즈수에 따른 압력을 구하는 방법을 제시한다. 또한 와도의 분포를 도시하여 원주에 의한 후류의 영향을 받지 않는 외부경계의 한계를 새로이 설정한다.

• 한국자동차공학회논문집
• /
• 제16권2호
• /
• pp.166-174
• /
• 2008

Numerical study on the impingement process and the fuel film formation of the hollow-cone fuel spray was conducted under vaporization condition, and the effect of the wall cavity angle on spray-wall impingement structure was investigated. A detailed understanding of this phenomena will help in designing injection systems and controlling the strategies to improve engine performance and exhaust emissions of the Gasoline Direct Injection (GDI) engine. The improved Abramzon model was used to model the spray vaporization process and the Gosman model was adopted for modeling of spray-wall impingement process. The calculated results of the spray-wall impingement process were compared with experimental results. The velocity field of the ambient gas, the Sauter Mean Diameter (SMD) and the generated fuel film on the wall, which are difficult to obtain by the experimental method, were also calculated and discussed. It was found that the radial distance after the wall impingement and the SMD decreased with increasing the cavity angle and the temperature.

• 물과 미래
• /
• 제31권4호
• /
• pp.291-297
• /
• 1998

2차원 흐름이 원주주위를 지날 때 발생하는 흐름의 변화가 기본방정식인 연속방정식과 운동량방정식에 의하여 수치적으로 해석된다. 수채해석 과정은 기본방정식에 유함수, 와도 및 흐름의 특성을 나타내는 무차원 매개변수를 도입하여 무차원 유함수-와소수송식을 유도한후, successive over relaxation scheme과 alternating direct implicit scheme으로 수행된다. 수치실험은 레이놀즈수 125-275를 갖는 흐름에 대하여 수행되었으며, 시간에 따른 유선, 와도, 원주표면의 압력을 구하는 방법에 있어서 기존의 수치해석에서는 주로 방사 운동량방정식만을 사용하였으나, 본 논문에서는 기존의 방법외에 방사 운동량방정식 및 접선 운동량방정식에 의해 압력을 계한하고, 두값을 비교하여 레이놀즈수에 따른 압력을 구하는 방법을 제시한다. 또한 와도의 분포를 도시하여 원주에 의한 후류의 영향을 받지 않는 외부경계의한계를 새로이 설정한다.

• 한국전산유체공학회지
• /
• 제11권2호
• /
• pp.8-18
• /
• 2006

A CFD benchmark calculation for a steam blowdown test was performed for 30 seconds to develop the methodology of numerical analysis for the thermal mixing between steam and subcooled water. In the CFD analysis, the grid model simulating the sparger and the IRWST pool were developed by the axisymmetric condition and then the steam condensation phenomena by a direct contact was modelled by the so-called condensation region model. Thermal mixing phenomenon in the subcooled water tank was treated as an incompressible flow, a free surface flow between the air and the water, a turbulent flow, and a buoyancy flow. The comparison of the CFD results with the test data showed a good agreement as a whole, but a small temperature difference was locally found at some locations. The commercial CFD code of CFX4.4 together with the condensation region model can simulate the thermal mixing behavior reasonably well when a sufficient number of mesh distribution and a proper numerical method are adopted.

• 대한기계학회논문집B
• /
• 제32권6호
• /
• pp.429-445
• /
• 2008

Two-phase compressible flow fields of air-water are investigated numerically in the fixed Eulerian grid framework. The phase interface is captured via volume fractions of each phase. A way to model two phase compressible flows as a single phase one is found based on an equivalent equation of states of Tait's type for a multiphase cell. The equivalent single phase field is discretized using the Roe‘s approximate Riemann solver. Two approaches are tried to suppress the pressure oscillation phenomena at the phase interface, a passive advection of volume fraction and a direct pressure relaxation with the compressible form of volume fraction equation. The direct pressure equalizing method suppresses pressure oscillation successfully and generates sharp discontinuities, transmitting and reflecting acoustic waves naturally at the phase interface. In discretizing the compressible form of volume fraction equation, phase interfaces are geometrically reconstructed to minimize the numerical diffusion of volume fraction and relevant variables. The motion of a projectile in a water-filled tube which is fired by the release of highly pressurized air is simulated presuming the flow field as a two dimensional one, and several design factors affecting the projectile movement are investigated.

• Korean Chemical Engineering Research
• /
• 제59권1호
• /
• pp.138-151
• /
• 2021

The effect of a reactant ratio on the growth of a buoyancy-driven instability in an irreversible A+B→C reaction system is analyzed theoretically and numerically. Taking a non-stoichiometric reactant ratio into account, new linear stability equations are derived without the quasi-steady state assumption (QSSA) and solved analytically. It is found that the main parameters to explain the present system are the Damköhler number, the dimensionless density difference of chemical species and the ratio of reactants. The present initial grow rate analysis without QSSA shows that the system is initially unconditionally stable regardless of the parameter values; however, the previous initial growth rate analysis based on the QSSA predicted the system is unstable if the system is physically unstable. For time evolving cases, the present growth rates obtained from the spectral analysis and pseudo-spectral method support each other, but quite differently from that obtained under the conventional QSSA. Adopting the result of the linear stability analysis as an initial condition, fully nonlinear direct numerical simulations are conducted. Both the linear analysis and the nonlinear simulation show that the reactant ratio plays an important role in the onset and the growth of the instability motion.