• 제목/요약/키워드: Runge-Kutta fourth order method

검색결과 68건 처리시간 0.026초

Nonlinear low-velocity impact response of graphene platelets reinforced metal foams doubly curved shells

  • Hao-Xuan Ding;Yi-Wen Zhang;Yin-Ping Li;Gui-Lin She
    • Steel and Composite Structures
    • /
    • 제49권3호
    • /
    • pp.281-291
    • /
    • 2023
  • Due to the fact that the nonlinear low-velocity impact response of graphene platelets reinforced metal foams (GPLRMF) doubly curved shells have not been investigated in the existing works, this paper aims to solve this issue. Using Reddy's high-order shear deformation theory (HSDT), the nonlinear governing equations of GPLRMF doubly curved shells are obtained by Euler-Lagrange method, discretized by Galerkin principle, and solved by the fourth-order Runge-Kutta method to obtain the impact force and central deflection. The nonlinear Hertz contact law is applied to determine the contact force. Finally, the impacts of graphene platelets (GPLs) distribution pattern, porosity distribution form, porosity coefficient, damping coefficient, impact parameters (radius and initial velocity), GPLs weight fraction, pre-stressing force and different shell types on the low-velocity impact curves are analyzed. It can be found that, among the four shell structures, the impact resistance of spherical shell is the best, while that of cylindrical shell is the worst.

최적회된 고차-고해상도 집적 유한 차분법을 이용한 초음속 제트 스크리치 톤 수치 해석 (Numerical Analysis of Supersonic Axisymmetric Screech Tone Noise Using Optimized High-Order, High-Resolution Compact Scheme)

  • 이인철;이덕주
    • The Journal of the Acoustical Society of Korea
    • /
    • 제25권1E호
    • /
    • pp.32-35
    • /
    • 2006
  • The screech tone of underexpanded jet is numerically calculated without any specific modeling for the screech tone itself. Fourth-order optimized compact scheme and fourth-order Runge-Kutta method are used to solve the 2D axisymmetric Euler equation. Adaptive nonlinear artificial dissipation model and generalized characteristic boundary condition are also used. The screech tone, generated by a closed loop between instability waves and quasi-periodic shock cells at the near field, is reasonably analyzed with present numerical methods for the underexpanded jet having Mach number 1.13. First of all, the centerline mean pressure distribution is calculated and compared with experimental and other numerical results. The instantaneous density contour plot shows Mach waves due to mixing layer convecting supersonically, which propagate downstream. The pressure signal and its Fourier transform at upstream and downstream shows the directivity pattern of screech tone very clearly. Most of all, we can simulate the axisymmetric mode change of screech tone very precisely with present method. It can be concluded that the basic phenomenon of screech tone including the frequency can be calculated by using high-order and high-resolution schemes without any specific numerical modeling for screech tone feedback loop.

효율적인 상세 반응 기구 해석을 위한 민감도 기반의 부분 음해법 (Partial Preconditioning Approach for the Solution of Detailed Kinetics Problems Based on Sensitivity Analysis)

  • 강기하;문성영;노진현;원수희;최정열
    • 한국연소학회지
    • /
    • 제13권1호
    • /
    • pp.17-22
    • /
    • 2008
  • A partly implicit/quasi-explicit method is introduced for the solution of detailed chemical kinetics with stiff source terms based on the standard fourth-order Runge-Kutta scheme. Present method solves implicitly only the stiff reaction rate equations, whereas the others explicitly. The stiff equations are selected based on the survey of the chemical Jaconian matrix and its Eigenvalues. As an application of the present method constant pressure combustion was analyzed by a detailed mechanism of hydrogen-air combustion with NOx chemistry. The sensitivity analysis reveals that only the 4 species in NOx chemistry has strong stiffness and should be solved implicitly among the 13 species. The implicit solution of the 4 species successfully predicts the entire process with same accuracy and efficiency at half the price.

  • PDF

Comparative study on cracked beam with different types of cracks carrying moving mass

  • Jena, Shakti P.;Parhi, Dayal R.;Mishra, Devasis
    • Structural Engineering and Mechanics
    • /
    • 제56권5호
    • /
    • pp.797-811
    • /
    • 2015
  • An analytical-computational method along with finite element analysis (FEA) has been employed to analyse the dynamic behaviour of deteriorated structures excited by time- varying mass. The present analysis is focused on the comparative study of a double cracked beam with inclined edge cracks and transverse open cracks subjected to traversing mass. The assumed computational method applied is the fourth order Runge-Kutta method. The analysis of the structure has been carried out at constant transit mass and speed. The response of the structure is determined at different crack depth and crack inclination angles. The influence of the parameters like crack depth and crack inclination angles are investigated on the dynamic behaviour of the structure. The results obtained from the assumed computational method are compared with those of the FEA for validation and found good agreements with FEA.

Numerical simulation of tuned liquid tank- structure systems through σ-transformation based fluid-structure coupled solver

  • Eswaran, M.;Reddy, G.R.
    • Wind and Structures
    • /
    • 제23권5호
    • /
    • pp.421-447
    • /
    • 2016
  • Wind-induced and earthquake-induced excitations on tall structures can be effectively controlled by Tuned Liquid Damper (TLD). This work presents a numerical simulation procedure to study the performance of tuned liquid tank- structure system through ${\sigma}$-transformation based fluid-structure coupled solver. For this, a 'C' based computational code is developed. Structural equations are coupled with fluid equations in order to achieve the transfer of sloshing forces to structure for damping. Structural equations are solved by fourth order Runge-Kutta method while fluid equations are solved using finite difference based sigma transformed algorithm. Code is validated with previously published results. The minimum displacement of structure is observed when the resonance condition of the coupled system is satisfied through proper tuning of TLD. Since real-time excitations are random in nature, the performance study of TLD under random excitation is also carried out in which the Bretschneider spectrum is used to generate the random input wave.

총괄 변수 모델을 이용한 DPF 재생 성능에 관한 연구 (A Study on the Regeneration Performance of DPF using Lumped Parameter Model)

  • 전문수
    • 융복합기술연구소 논문집
    • /
    • 제1권1호
    • /
    • pp.41-47
    • /
    • 2011
  • With the world-wide demand on the emission minimization, the needs on the diesel aftertreatment devices with high efficiency are also increasing. In order to effectively develop or design a high-performance diesel particulate filter, a clear understanding on the deposition and regeneration mechanism is required. In the present study, a theory on the lumped parameter model for wall-flow type diesel particulate filters is described focusing on the deposition efficiency, pressure drop inside the filter. The fourth order explicit Runge-Kutta method is utilized for the mass flow rate computation. Engine operation modes with controlled and uncontrolled regeneration options are selected. The computational lumped parameter model is validated by comparing the computed results with the measured data.

  • PDF

연안지역의 정체수역에서 방류되는 하$\cdot$폐수의 근역거동 수치모의 (Numerical Simulation for Near Field-Behavior of Wastewater Discharged into Stagnant Ambient in Coastal Region)

  • 권석재;서일원
    • 한국해안해양공학회지
    • /
    • 제17권3호
    • /
    • pp.166-177
    • /
    • 2005
  • 본 연구에서는 정체수역에서 유입이론과 일부 형상계수의 조정을 통해 근역에서의 하$\cdot$폐수 혼합거동을 해석할 수 있는 근역제트적분모형을 개발하기위해 총 6개의 상미분 보존방정식에 6개의 미지수를 가지는 문제를 수치적으로 풀기 위한 4차의 Runge-Kutta기법을 사용하였다. 또한, LIF 시스템을 이용하여 검정과정을 통해서 단일수평부력제트의 수리실험을 수행하였다. 그리고 기존의 모형 CORMIX 1, WSJET,그리고 본 모형의 계산결과를 수리실험결과와 서로 비교하였다. VISJET모형에 의해 예측된 중심선 제적이 운동량과 부력이 지배적인 구간에서는 실험결과와 근접한 반면에 본 제트적분모형에 의해 예측된 결과는 천이영역에서 측정된 궤적과 잘 일치하였다. 중심선희석률에 있어서 운동량과 부력이 지배적인 구간에서 CORMIX1 모형의 결과와 잘 일치하는 반면에 초기영역과 천이 영역에서 본 모형의 결과와 대체로 잘 일치하는 경향을 보였다.

A virtual shaker testing experience: Modeling, computational methodology and preliminary results

  • Nali, Pietro;Bettacchioli, Alain;Landi, Guglielmo;Gnoffo, Marco
    • Advances in aircraft and spacecraft science
    • /
    • 제5권2호
    • /
    • pp.251-258
    • /
    • 2018
  • This work illustrates the progress of a TAS activity at exploring the challenges and the benefits of the Virtual Shaker Testing (VST) approach. The definition and the validation of new computational methodologies with respect to the state of the art were encouraged throughout this activity. The shaker Finite Element (FE) model in lateral configuration was built for the purpose and it was merged with the SpaceCraft (S/C) FE model, together with the S/C-Shaker adapter. FE matrices were reduced through the Craig-Bampton method. The VST transient analysis was performed in MATLAB(R) numerical computing environment. The closed-loop vibration control is accounted for and the solution is obtained through the fourth-order Runge Kutta method. The use of pre-existing built-in functions was limited by authors with the aim of tracing the impact of all the problems' parameters in the solution. Assumptions and limitations of the proposed methodology are detailed throughout this paper. Some preliminary results pertaining to the current progress of the activity are thus illustrated before the conclusions.

Flow Characteristics of Gaseous Leak flows in Narrow Cracks

  • Hong, Chung-Pyo
    • 한국유체기계학회 논문집
    • /
    • 제11권4호
    • /
    • pp.14-21
    • /
    • 2008
  • The prediction for gaseous leak flows through a narrow crack is important for a leak-before-break (LBB) analysis. Therefore, the methodology to obtain the flow characteristics of gaseous leak flow in a narrow crack for the wide range by using the product of friction factor and Reynolds number correlations (fRe) for a micro-channel is developed and presented. The correlation applied here was proposed by the previous study. The fourth-order Runge-Kutta method was employed to integrate the nonlinear ordinary differential equation for the pressure and the regular-Falsi method was also employed to find the inlet Mach number. A narrow crack whose opening displacement ranges from 10 to $100{\mu}m$ with a crack length in the range from 2 to 200mm was chosen for sample prediction. The present results are compared with both numerical simulation results and available experimental measurements. The results are in excellent agreement with them. The leak flow rate can be approximately predicted by using proposed methodology.

GENERALIZED EULER PROCESS FOR SYSTEMS OF NONLINEAR DIFFERENTIAL EQUATIONS

  • Yu, Dong-Won
    • Journal of applied mathematics & informatics
    • /
    • 제7권3호
    • /
    • pp.941-958
    • /
    • 2000
  • Euler method is generalized to solve the system of nonlinear differential equations. The generalization is carried out by taking a special constant matrix S so that exp(tS) can be exactly computed. Such a matrix S is extracted from the Jacobian matrix of the given problem. Stability of the generalized Euler process is discussed. It is shown that the generalized Euler process is comparable to the fourth order Runge-Kutta method. We also exemplify that the important qualitative and geometric features of the underlying dynamical system can be recovered by the generalized Euler process.