• Structural Engineering and Mechanics
• /
• v.72 no.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.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2002.10a
• /
• 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.

• Journal of the Society of Naval Architects of Korea
• /
• v.55 no.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.

• Proceedings of the KSR Conference
• /
• 2011.10a
• /
• pp.1594-1600
• /
• 2011

Recently developed FDS5 CFD code has employed a near-wall flow treatment method which is Werner-Wengle wall law provided by NIST(National Institute of Standards and Technology). In this study, the wall law has been verified against DNS(Direct Numerical Simulation) data in the parallel plate. The $y^+$ was kept above 11 to fulfill the near-wall flow requirement in the grid generation. The total grid was $32{\times}32{\times}32$. The boundary condition for inlet and outlet was periodic condition and for both side, symmetric condition was used. The fully developed turbulent flow was generated and Re = 10,700. The simulated results were compared with DNS data. RANS results were also used for verification.

• Water for future
• /
• v.21 no.3
• /
• pp.291-291
• /
• 1988

The two dimensional time dependent flow past a circular cylinder is analyzed numerically. In the analysis, equations of conservation of mass and momentum are transformed to equations of stream function-vorticity and vorticity transport, and nondimensionalized by nondimensional parameters representing flow characteristics, The resulting stream function-vorticity equstion and vorticity transport equation are solved by successive over relaxation scheme and alternating direct implicit scheme. Numerical experments are performed for the flow in the range of Reynolds number 125 to 275. The time dependent streamlines, vorticities, pressure on cylinder surface, separation angle, and drag and lift coefficients are calculated, and the method for estimation of pressure on cylinder surface and the outer boundary limit are developed.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.16 no.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.

• Water for future
• /
• v.31 no.4
• /
• pp.291-297
• /
• 1998

The two dimensional time dependent flow past a circular cylinder is analyzed numerically. In the analysis, equations of conservation of mass and momentum are transformed to equations of stream function-vorticity and vorticity transport, and nondimensionalized by nondimensional parameters representing flow characteristics, The resulting stream function-vorticity equation and vorticity transport equation are solved by successive over relaxation scheme and alternating direct implicit scheme. Numerical experiments are performed for the flow in the range of Reynolds number 125 to 275. The time dependent streamlines, vorticities, pressure on cylinder surface, separation angle, and drag and lift coefficients are calculated, and the method for estimation of pressure on cylinder surface and the outer boundary limit are developed.

• Journal of computational fluids engineering
• /
• v.11 no.2 s.33
• /
• 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.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.32 no.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
• /
• v.59 no.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.