• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2004.03a
• /
• pp.550-559
• /
• 2004

Theoretical studies have been carried out to examine the influence of the grain geometry-dependent driving forces, which control the internal flow pattern of solid rockets. Numerical studies have been executed with the help of a two-dimensional code. This code solves standard k-omega turbulence equations using the coupled second order implicit unsteady formulation. It has been concluded that the grain port divergence angles have significant leverage on the formation of recirculation bubbles leading for pressure oscillations, flow separation and reattachment. In solid rockets flow reattachment will favour secondary ignition and that will add to the complexity of the starting transient prediction.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2004.03a
• /
• pp.528-531
• /
• 2004

For a simple one-dimensional ignition problem a mathematical model is described to investigate the difficulties in numerical simulations. Some computation results are obtained and comparison is made with analytical solution. Discussions are made on topics such as 1) coordinate transformation, 2) gas-phase and solid-phase analysis; (divergence form of the governing system, a finite-volume discretization, implicit time integration, upwind split flux, spatial accuracy improvement are described. Mass, reagent mass, and energy conservations are solved.), and 3) method to determine quantities on the burning surface (matching). Results obtained for small values of the non-dimensional pressure show a steady-combustion and good agreement with the analytical solution. Numerical instability appeared for larger values of the pressure, discussion on the cause of the problem is made. This effort is a part of a study of flame spread phenomena on solid propellant surface.

• Nuclear Engineering and Technology
• /
• v.53 no.12
• /
• pp.3861-3878
• /
• 2021

The steady-state simplified spherical harmonics equations (SP_N equations) are a higher order approximation to the neutron transport equations than the neutron diffusion equation that also have reasonable computational demands. This work extends these results for the analysis of transients by comparing of two formulations of time-dependent SP_N equations considering different treatments for the time derivatives of the field moments. The first is the full system of equations and the second is a diffusive approximation of these equations that neglects the time derivatives of the odd moments. The spatial discretization of these methodologies is made by using a high order finite element method. For the time discretization, a semi-implicit Euler method is used. Numerical results show that the diffusive formulation for the time-dependent simplified spherical harmonics equations does not present a relevant loss of accuracy while being more computationally efficient than the full system.

• Steel and Composite Structures
• /
• v.48 no.1
• /
• pp.73-88
• /
• 2023

This paper presents a novel implicit level set method for topology optimization of functionally graded (FG) structures with pre-existing discontinuities (pre-cracks) using radial basis functions (RBF). The mathematical formulation of the optimization problem is developed by incorporating RBF-based nodal densities as design variables and minimizing compliance as the objective function. To accurately capture crack-tip behavior, crack-tip enrichment functions are introduced, and an eXtended Finite Element Method (X-FEM) is employed for analyzing the mechanical response of FG structures with strong discontinuities. The enforcement of boundary conditions is achieved using the Hamilton-Jacobi method. The study provides detailed mathematical expressions for topology optimization of systems with defects using FG materials. Numerical examples are presented to demonstrate the efficiency and reliability of the proposed methodology.

• 한국전산유체공학회:학술대회논문집
• /
• 2004.03a
• /
• pp.146-153
• /
• 2004

The Numerical study has been carried out to investigate the effects of chemical reaction and thermal radiation on the rocket plume flow-field at various altitudes. The theoretical formulation is based on the Navier-Stokes equations for compressible flows along with the infinitely fast chemistry and thermal radiation. The governing equations were solved by a finite volume fully-implicit TVD(Total Variation Diminishing) code which uses Roe's approximate Riemann solver and MUSCL(Monotone Upstream-centered Schemes for Conservation Laws) scheme. LU-SGS (Lower Upper Symmetric Gauss Seidel) method is used for the implicit solution strategy. An equilibrium chemistry module for hydrocarbon mixture with detailed thermo-chemical properties and a thermal radiation module for optically thin media were incorporated with the fluid dynamics code. In this study, kerosene-fueled rocket was assumed operating at O/F ratio of 2.34 with a nozzle expansion ratio of 6.14. Flight conditions considered were Mach number zero at ground level, Mach number 1.16 at altitude 5.06km and Mach number 2.9 at altitude 17.34km. Numerical results gave the understandings on the detailed plume structures at different altitude conditions. The diffusive effect of the thermal radiation on temperature field and the effect of chemical recombination during the expansion process could be also understood. By comparing the results from frozen flow and infinitely fast chemistry assumptions, the excess temperature of the exhaust gas resulting from the chemical recombination seems to be significant and cannot be neglected in the view point of performance, thermal protection and flow physics.

• Journal of the Korea Concrete Institute
• /
• v.20 no.3
• /
• pp.317-323
• /
• 2008

In the current study, A model for the flow analysis of fresh and highly flowable concrete has been developed using a particle method, the moving particle semi-implicit (MPS) method. The phenomena on the flow of concrete has been considered as a visco-plastic flow problem, and the basic governing equation of concrete particle dynamics has been based on the Navier-Stokes equation in Lagrangian form and the conservation of mass. In order to formulate a visco-plastic flow constitutive law of fresh concrete, concrete is modeled as a highly viscous material in the state of non-flow and as a visco-plastic material in the state of flow after reaching the yield stress of fresh concrete. A flow test of fresh concrete in the L-box was simulated and the predicted flow was well matched with the experimental result. The developed method was well showed the flow motion of concrete particles because it was formulated to be based on the motion of visco-plastic fluid dynamics.

• Computational Structural Engineering
• /
• v.4 no.2
• /
• pp.87-94
• /
• 1991

In this paper, a mixed variational principle applicable to the linear elasticity of inhomogeneous anisotropic materials is presented. For derivation of the general variational principle, a systematic procedure for the variational formulation of linear coupled boundary value problems developed by Sandhu et al. is employed. Consistency condition of the field operators with the boundary operators results in explicit inclusion of boundary conditions in the governing functional. Extensions of admissible state function spaces and specialization to a certain relation in the general governing functional lead to the desired mixed variational principle. In the physical sense, the present variational principle is analogous to the Reissner's recent formulation obtained by applying Lagrange multiplier technique followed by partial Legendre transform to the classical minimum potential energy principle. However, the present one is more advantageous for the application to the general anisotropic materials since Reissner's principle contains an implicit function which is not easily converted to an explicit form.

• Journal of the Society of Naval Architects of Korea
• /
• v.48 no.3
• /
• pp.275-280
• /
• 2011

In the present study, a unified viscoplastic-damage model has been applied in order to describe the mechanical characteristics of fresh water ice such as nonlinear material behavior and volume fraction. The strain softening phenomenon of fresh water ice under quasi-static compressive loading has been evaluated based on unified viscoplastic model. The material degradation such as growth of slip/fraction has quite close relation with material inside damage. The volume fraction phenomenon of fresh water ice has been identified based on volume fraction (nucleation and growth of damage) model. The viscoplastic-damage model has been transformed to the fully implicit formulation and the discretized formulation has been implemented to ABAQUS user defined subroutine (User MATerial: UMAT) for the benefit of application of commercial finite element program. The proposed computational analysis method has been compared to uni-axial compression test of fresh water ice in order to validate the compatibilities, clarities and usefulness.

• Journal of Mechanical Science and Technology
• /
• v.16 no.8
• /
• pp.1102-1108
• /
• 2002

In this study, a topological design sensitivity of the ai. bearing surface (ABS) is suggested by using an adjoint variable method. The discrete form of the generalized lubrication equation based on a control volume formulation is used as a compatible condition. A residual function of the slider is considered as an equality constraint function, which represents the slider in equilibrium. The slider thickness parameters at all grid cells are chosen as design variables since they are the topological parameters determining the ABS shape. Then, a complicated adjoint variable equation is formulated to directly handle the highly nonlinear and asymmetric coefficient matrix and vector in the discrete system equation of air-lubricated slider bearings. An alternating direction implicit (ADI) scheme is utilized for the numerical calculation. This is an efficient iterative solver to solve large-scale problem in special band storage. Then, a computer program is developed and applied to a slider model of a sophisticated shape. The simulation results of design sensitivity analysis (DSA) are directly compared with those of FDM at the randomly selected grid cells to show the effectiveness of the proposed approach. The overall distribution of DSA results are reported, clearly showing the region on the ABS where special attention should be given during the manufacturing process.

• World Technopolis Review
• /
• v.4 no.2
• /
• pp.46-61
• /
• 2015

There is an implicit assumption in most regional innovation policy studies that once a policy has been made the policy will be implemented. This assumption is not valid for regional innovation policies formulated in many countries. This assumption rests upon certain political and organization conditions. It is frequently said that regional innovation system in Korea is one of the successful cases because of properly established institution for the implementation of innovation policies. The components of institution for innovation are defined in this article. For the analysis of institution for regional innovation system in Korea, three aspects such as organization, policies & programs, and governance structure were taken into account. Institution for regional innovation system in Korea is analyzed with three aspects: (1) organization for policy implementation, (2) innovation policy, and (3) governance structure. Firstly, organizations for innovation in Korea are composed of three categories such as organizations for policy formulation, organizations for policy implementation and agencies for coordination. Secondly, there are two categories of policies for innovation: (1) policies for the enhancement of national competitiveness and policies for the regional innovation capacity building, and (2) policies for fostering manpower and policy for regulatory reform. Thirdly, innovation governance in Korea is composed of three layer structure: (1) top level governance which is composed of two committee, three ministries and two agencies, (2) local level governance such as one innovation related offices and one center for regulatory reform, and (3) one category of agency for coordination as a regional platform.