• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제20권3호
• /
• pp.243-259
• /
• 2016

The Richards equation for water movement in unsaturated soil is highly nonlinear partial differential equations which are not solvable analytically unless unrealistic and oversimplifying assumptions are made regarding the attributes, dynamics, and properties of the physical systems. Therefore, conventionally, numerical solutions are the only feasible procedures to model flow in partially saturated porous media. The standard Finite element numerical technique is usually coupled with an Euler time discretizations scheme. Except for the fully explicit forward method, any other Euler time-marching algorithm generates nonlinear algebraic equations which should be solved using iterative procedures such as Newton and Picard iterations. In this study, lumped mass and distributed mass in the frame of Picard and Newton iterative techniques were evaluated to determine the most efficient method to solve the Richards equation with finite element model. The accuracy and computational efficiency of the scheme and of the Picard and Newton models are assessed for three test problems simulating one-dimensional flow processes in unsaturated porous media. Results demonstrated that, the conventional mass distributed finite element method suffers from numerical oscillations at the wetting front, especially for very dry initial conditions. Even though small mesh sizes are applied for all the test problems, it is shown that the traditional mass-distributed scheme can still generate an incorrect response due to the highly nonlinear properties of water flow in unsaturated soil and cause numerical oscillation. On the other hand, non oscillatory solutions are obtained and non-physics solutions for these problems are evaded by using the mass-lumped finite element method.

• Wind and Structures
• /
• 제19권3호
• /
• pp.233-247
• /
• 2014

Vortex-induced oscillation is a type of aeroelastic phenomenon, to which extended structures such as long-span bridges are most susceptible. The vortex-induced vibration (VIV) behaviors of a concerned bridge were investigated conventionally in virtue of wind tunnel tests on string-mounted sectional models. This necessitates the building of a linkage between the response of the sectional model and that of the prototype structure. Although many released literatures have related to this issue and provided suggestions, there is a lack of consistency among them. In this study, some theoretical models describing the vortex-induced structural motion, including the linear empirical model, the nonlinear empirical model and the modified (or generalized) nonlinear empirical model, are firstly reviewed. Then, the concept of equivalent mass density is introduced based on the principle that an equal input of energy should result in identical structural amplitudes. Based on these, the theoretical linkages between the amplitude of a section model and that corresponding to the prototype bridge are discussed with different analytical models. Theoretical derivation indicates that such connections are dependent mainly on two factors, one is the presupposed shape of deformation, and the other is the theoretical VIV model employed. The theoretical analysis in this study shows that, in comparison to the nonlinear empirical models, the linear one can result in obvious larger estimations of the full bridges' responses, especially in cases of cable-stayed bridges.

• Wind and Structures
• /
• 제7권4호
• /
• pp.265-280
• /
• 2004

This paper presents some fundamental results on the dynamics of the periodic Karman wake behind a circular cylinder. The wake is treated like a dynamical system. External forcing is then introduced and its effect investigated. The main result obtained is the following. Perturbation of the wake, by controlled cylinder oscillations in the flow direction at a frequency equal to the Karman vortex shedding frequency, leads to instability of the Karman vortex structure. The resulting wake structure oscillates at half the original Karman vortex shedding frequency. For higher frequency excitation the primary pattern involves symmetry breaking of the initially shed symmetric vortex pairs. The Karman shedding phenomenon can be modeled by a nonlinear oscillator. The symmetrical flow perturbations resulting from the periodic cylinder excitation can also be similarly represented by a nonlinear oscillator. The oscillators represent two flow modes. By considering these two nonlinear oscillators, one having inline shedding symmetry and the other having the Karman wake spatio-temporal symmetry, the possible symmetries of subsequent flow perturbations resulting from the modal interaction are determined. A theoretical analysis based on symmetry (group) theory is presented. The analysis confirms the occurrence of a period-doubling instability, which is responsible for the frequency halving phenomenon observed in the experiments. Finally it is remarked that the present findings have important implications for vortex shedding control. Perturbations in the inflow direction introduce 'control' of the Karman wake by inducing a bifurcation which forces the transfer of energy to a lower frequency which is far from the original Karman frequency.

• Advances in aircraft and spacecraft science
• /
• 제7권3호
• /
• pp.215-228
• /
• 2020

As is shown in the paper, the Koltunov-Rzhanitsyn singular kernel of heredity (when constructing mathematical models of the dynamics problem of the hereditary theory of viscoelasticity) adequately describes real mechanical processes, best approximates experimental data for a long period of time. A mathematical model of the problem of the flutter of viscoelastic plates moving in a gas with a high supersonic velocity is given. Using the Bubnov-Galerkin method, discrete models of the problem of the flatter of viscoelastic plates flowed over by supersonic gas flow are obtained. A numerical method is developed to solve nonlinear integro-differential equations (IDE) for the problem of the hereditary theory of viscoelasticity with weakly singular kernels. A general computational algorithm and a system of application programs have been developed, which allow one to investigate the nonlinear dynamic problems of the hereditary theory of viscoelasticity with weakly singular kernels. On the basis of the proposed numerical method and algorithm, nonlinear problems of the flutter of viscoelastic plates flowed over in a gas flow at an arbitrary angle are investigated. In a wide range of changes in various parameters of the plate, the critical velocity of the flutter is determined. It is shown that the singularity parameter α affects not only the oscillations of viscoelastic systems, but the critical velocity of the flutter as well.

• Advances in aircraft and spacecraft science
• /
• 제3권3호
• /
• pp.271-298
• /
• 2016

In this paper, we study a problem of passive suppression of helicopter Ground Resonance (GR) using a single degree freedom Nonlinear Energy Sink (NES), GR is a dynamic instability involving the coupling of the blades motion in the rotational plane (i.e. the lag motion) and the helicopter fuselage motion. A reduced linear system reproducing GR instability is used. It is obtained using successively Coleman transformation and binormal transformation. The analysis of the steadystate responses of this model is performed when a NES is attached on the helicopter fuselage. The NES involves an essential cubic restoring force and a linear damping force. The analysis is achieved applying complexification-averaging method. The resulting slow-flow model is finally analyzed using multiple scale approach. Four steady-state responses corresponding to complete suppression, partial suppression through strongly modulated response, partial suppression through periodic response and no suppression of the GR are highlighted. An algorithm based on simple criterions is developed to predict these steady-state response regimes. Numerical simulations of the complete system confirm this analysis of the slow-flow dynamics. A parametric analysis of the influence of the NES damping coefficient and the rotor speed on the response regime is finally proposed.

• 한국추진공학회지
• /
• 제15권6호
• /
• pp.98-106
• /
• 2011

본 논문에서는 많은 가스터빈 산업체 및 연구기관에서 연소불안정 현상과 관련된 변수들을 예측하기 위해 가장 보편적으로 이루어지고 있는 열음향 해석 모델에 대한 기술 소개 및 최근의 연구 동향을 분석하였다. 선형 시스템 해석을 통하여 연소 불안정이 발생하는 고유 주파수 및 불안정 초기 성장률의 예측이 가능하다. 이를 위하여 정의된 시스템에서의 음향파와 열발생율 섭동간의 선형 관계식을 선형 음향 이론으로부터 유도할 수 있고, 이 관계식의 해를 구하기 위해서 가장 중요한 부분은 화염 전달 함수로부터 n-${\tau}$ 함수를 구하여 열발생율 섭동 결과에 대한 정보를 얻는 것이다. 현재까지의 연구 결과로부터 선형 특성 해석에는 상당한 진보가 이루어져 왔고, 실제 가스터빈 연소기에 적용하는 노력이 있었으나, 한계 진폭과 과도기 현상 예측을 위해서 요구되는 비선형 동적 특성 모델링 기술 개발은 현재 간단한 연소기와 버너의 적용에 머물러 있는 실정이다. 실제 복잡한 가스터빈과 같은 연소 시스템에 적용되기 위해서는 비선형 경계 조건을 고려한 시스템 동적 특성 연구와 화염의 비선형 거동을 더욱 정확히 설명할 수 있는 전달 함수에 대한 예측 기술이 선행되어야 한다.

• 대한조선학회논문집
• /
• 제31권2호
• /
• pp.45-56
• /
• 1994

2차원 수중익 주위의 비선형 자유표면 유동을 시간영역에서 시뮬레이션할 수 있는 수치해법이 개발되었다. 본 수치해법은 고차 스펙트럴법과 경계요소법을 조합한 스펙트럴/경계요소법(spectral/boundary-element method)이며, 자유표면은 고차 스펙트럴법에 의해 그리고 수중익과 후류 보오텍스는 경계요소법에 의해 다루어 진다. 본 방법은 자유표면과 수중익의 비선형/비정상 상호작용문제에 폭넓게 적용될 수 있으며, 특히 스펙트럴법을 사용하므로 자유표면 유동을 매우 효율적으로 다룰 수 있다. 적용예로 정지상태에서 출발하여 자유표면 근처에서 균속 전진운동 또는 전진 및 동요운동을 하는 경우가 다루어 졌고, 자유표면파와 동유체력에서의 비선형/비정상 효과들이 보여지고 있다. 비정상해의 특수한 경우로 얻어지는 정상상태의 계산결과들이 다른 이론 또는 실험결과들과 비교되었으며, 좋은 일치를 보이고 있다.

• 제어로봇시스템학회논문지
• /
• 제10권12호
• /
• pp.1155-1163
• /
• 2004

Extracellular signal-regulated kinase (ERK) signaling pathway is one of the mitogen-activated protein kinase (MAPK) signal transduction pathways. This pathway is known as pivotal in many signaling networks that govern proliferation, differentiation and cell survival. The ERK signaling pathway comprises positive and negative feedback loops, depending on whether the terminal kinase stimulates or inhibits the activation of the initial level. In this paper, we attempt to model the ERK pathway by considering both of the positive and negative feedback mechanisms based on Michaelis-Menten kinetics. In addition, we propose a fraction ratio model based on the mass action law. We first develop a mathematical model of the ERK pathway with fraction ratios. Secondly, we analyze the dynamical properties of the fraction ratio model based on simulation studies. Furthermore, we propose a concept of an inhibitor, catalyst, and substrate (ICS) controller which regulates the inhibitor, catalyst, and substrate concentrations of the ERK signal transduction pathway. The ICS controller can be designed through dynamical analysis of the ERK signaling transduction pathway within limited concentration ranges.

• Journal of Magnetics
• /
• 제21권3호
• /
• pp.450-459
• /
• 2016

Experimental results and data of numerical simulations are presented, concerning generation of wideband radio frequency (RF) oscillations in a nonlinear transmission line (NLTL) which contains a pre-magnetized core of ferrite material. Emphasis is made on the means for extracting the RF signal from the line, in order to radiate it into free space. Antennas of two types that can be used for the purpose are considered, both featuring a coaxial design. This permits availing of the principal advantages of coaxial systems, such as operation in the mode of a traveling TEM wave; wide range of the transmitted frequencies, and a reasonably simple design. The antennas studied, specifically a disc-cone dipole (DCD) and an impulse reflector antenna ('Half-IRA', or HIRA type) differ significantly in effective width of the radiated spectrum and in spatial characteristics of the radiated field in far region.

• 대한조선학회지
• /
• 제23권2호
• /
• pp.1-13
• /
• 1986

The nonlinear hydrodynamic forces acting on a two-dimensional circular cylinder, oscillating with large amplitude in the free surface, are calculated by using the Semi-Lagrangian Time-Step-ping Method used by O.M. Faltinsen. In present calculation the position and the potential value of free surface are calculated using the exact kinematic and dynamic free surface boundary condition. At each time step an integral equation is solved to obtain the value of potential and normal velocity along the boundaries, consisting of both the body surface and the free surface. Some effort was devoted to the elimination of instability arising in the range of high frequency. Numerical simulations were performed up to the 3rd or 4th period which seems to be enough for the transient effect to die out. Each harmonic component and time-mean force are obtained by the Fourier transform of forces in time domain. The results are compared with others' experimental and theoretical results. Particularly, the calculation shows the tendency that the acceleration-phase 1st-harmonic component(added mass) increases as the motion amplitude increases and a reverse tendency in the velocity-phase 1st-harmonic component(damping coefficient). The Yamashita's experimental result also shows the same tendency. In general, the present result show relatively good agreement with the Yamashita's experimental result except for the time-mean force.