• 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.

• 대한의용생체공학회:의공학회지
• /
• 제13권3호
• /
• pp.181-188
• /
• 1992

This paper provides an overview of system analysis of immunology. The theoretical research in this area is aimed at an understanding of the precise manner by which the immune system controls Infec pious diseases, cancer, and AIDS. This can provide a systematic plan for immunological experimentation by means of an integrated program of immune system analysis, mathematical modeling and computer simulation. Biochemical reactions and cellular fission are naturally modeled as nonlinear dynamical processes to synthesize the human immune system! as well as the complete organism it is intended to protect. A foundation for the control of tumors is presented, based upon the formulation of a realistic, knowledge based mathematical model of the interaction between tumor cells and the immune system. Ordinary bilinear differential equations which are coupled by such nonlinear term as saturation are derived from the basic physical phenomena of cellular and molecular conservation. The parametric control variables relevant to the latest experimental data are also considered. The model consists of 12 states, each composed of first-order, nonlinear differential equations based on cellular kinetics and each of which can be modeled bilinearly. Finally, tumor control as an application of immunotherapy is analyzed from the basis established.

• Coupled systems mechanics
• /
• 제10권1호
• /
• pp.39-60
• /
• 2021

The present paper is concerned with the study of nonlinear ultrasonic waves in a magneto thermo (MT) elastic armchair single-walled carbon nanotube (ASWCNT) resting on polymer matrix. The analytical formulation is developed based on Eringen's nonlocal elasticity theory to account small scale effect. After developing the formal solution of the mathematical model consisting of partial differential equations, the frequency equations have been analyzed numerically by using the nonlinear foundations supported by Winkler-Pasternak model. The solution is obtained by ultrasonic wave dispersion relations. Parametric work is carried out to scrutinize the influence of the non local scaling, magneto-mechanical loadings, foundation parameters, various boundary condition and length on the dimensionless frequency of nanotube. It is noticed that the boundary conditions, nonlocal parameter, and tube geometrical parameters have significant effects on dimensionless frequency of nano tubes. The results presented in this study can provide mechanism for the study and design of the nano devices like component of nano oscillators, micro wave absorbing, nano-electron technology and nano-electro- magneto-mechanical systems (NEMMS) that make use of the wave propagation properties of armchair single-walled carbon nanotubes embedded on polymer matrix.

• Structural Engineering and Mechanics
• /
• 제90권3호
• /
• pp.219-232
• /
• 2024

In current study, for the first time, Nonlinear Bending of a skew microplate made of a laminated composite strengthened with graphene nanosheets is investigated. A mixture of mechanical and thermal stresses is applied to the plate, and the reaction is analyzed using the First Shear Deformation Theory (FSDT). Since different percentages of graphene sheets are included in the multilayer structure of the composite, the characteristics of the composite are functionally graded throughout its thickness. Halpin-Tsai models are used to characterize mechanical qualities, whereas Schapery models are used to characterize thermal properties. The microplate's non-linear strain is first calculated by calculating the plate shear deformation and using the Green-Lagrange tensor and von Karman assumptions. Then the elements of the Couple and Cauchy stress tensors using the Modified Coupled Stress Theory (MCST) are derived. Next, using the Hamilton Principle, the microplate's governing equations and associated boundary conditions are calculated. The nonlinear differential equations are linearized by utilizing auxiliary variables in the nonlinear solution by applying the Frechet approach. The linearized equations are rectified via an iterative loop to precisely solve the problem. For this, the Differential Quadrature Method (DQM) is utilized, and the outcomes are shown for the basic support boundary condition. To ascertain the maximum values of microplate deflection for a range of circumstances-such as skew angles, volume fractions, configurations, temperatures, and length scales-a parametric analysis is carried out. To shed light on how the microplate behaves in these various circumstances, the resulting results are analyzed.

• 대한기계학회논문집A
• /
• 제27권7호
• /
• pp.1138-1145
• /
• 2003

In this paper, a study on the dynamic behavior and lubrication characteristics of a reciprocating compression mechanism used in small refrigeration compressor is performed. In the problem formulation of the compressor dynamics, the viscous frictional force between piston and cylinder wall is considered in order to determine the coupled dynamic behaviors of piston and crankshaft. The solutions of the equations of motion of the reciprocating mechanism along with the time dependent Reynolds equations for the lubricating film between piston and cylinder wall and oil films of the journal bearings are obtained simultaneously. The hydrodynamic forces of journal bearings are calculated using finite bearing model and Gumbel boundary condition. And, a Newton-Raphson procedure was employed in solving the nonlinear equations of piston and crankshaft. The results explored the effects of design parameters on the stability and lubrication characteristics of the compression mechanism.

• 대한기계학회논문집A
• /
• 제26권12호
• /
• pp.2511-2518
• /
• 2002

Dynamic behaviors and stability of an optical disk drive coupled with an automatic ball balancer (ABB) are analyzed by a theoretical approach. The feeding system is modeled a rigid body with six degree-of-freedom. Using Lagrange's equation, we derive the nonlinear equations of motion for a non -autonomous system with respect to the rectangular coordinate. To investigate the dynamic stability of the system in the neighborhood of the equilibrium positions, the monodromy matrix technique is applied to the perturbed equations. On the other hand, time responses are computed by the Runge -Kutta method. We also investigate the effects of the damping coefficient and the position of ABB on the dynamic behaviors of the system.

• 한국소음진동공학회:학술대회논문집
• /
• 한국소음진동공학회 2001년도 추계학술대회논문집 II
• /
• pp.983-988
• /
• 2001

Dynamic behaviors and stability of an optical disk drive coupled with an automatic ball balancer(ABB) are analyzed by a theoretical approach. The feeding system is modeled a rigid body with six degree-of-freedom. Using Lagrange's equation, we derive the nonlinear equations of motion for a non-autonomous system with respect to the rectangular coordinate. To investigate the dynamic stability of the system in the neighborhood of the equilibrium positions, the monodromy matrix technique is applied to the perturbed equations. On the other hand, time responses are computed by the Runge-Kutta method. We also investigate the effects of the damping coefficient and the position of ABB on the dynamic behaviors of the system.

• Kyungpook Mathematical Journal
• /
• 제60권4호
• /
• pp.853-867
• /
• 2020

The aim of this study was to analyze the momentum and heat transfer of a rotating nanofluid with conducting spherical dust particles. The fluid flows over a stretching surface under the influence of an external magnetic field. By applying similarity transformations, the governing partial differential equations were trans-formed into nonlinear coupled ordinary differential equations. These equations were solved with the built-in function bvp4c in MATLAB. Moreover, the effects of the rotation parameter ω, magnetic field parameter M, mass concentration of the dust particles α, and volume fraction of the nano particles 𝜙, on the velocity and temperature profiles of the fluid and dust particles were considered. The results agree well with those in published papers. According to the result the hikes in the rotation parameter ω decrease the local Nusselt number, and the increasing volume fraction of the nano particles 𝜙 increases the local Nusselt number. Moreover the friction factor along the x and y axes increases with increasing volume fraction of the nano particles 𝜙.

• Steel and Composite Structures
• /
• 제42권5호
• /
• pp.599-615
• /
• 2022

In the present paper, the numerical dynamic analysis of a functionally graded nano-scale nonuniform tube was investigated according to the high-order beam theory coupled with the nonlocal gradient strain theory. The supposed cross-section is changed along the pipe length, and the material distribution, which combines both metal and ceramics, is smoothly changed in the pipe length direction, which is called axially functionally graded (AFG) pipe. Moreover, the porosity voids are dispersed in the cross-section and the radial pattern that the existence of both material distribution along the tube length and porosity voids make a two-dimensional functionally graded (2D-FG) truncated conical pipe. On the basis of the Hamilton principle, the governing equations and the associated boundary conditions equations are derived, and then a numerical approach is applied to solve the obtained equations.

• International Journal of Aeronautical and Space Sciences
• /
• 제4권1호
• /
• pp.53-62
• /
• 2003

In this study, nonlinear aeroelastic characteristics of an supersonic missile wing with strong shock interferences are investigated. The missile wing model has a freeplay structural nonlinearity at its pitch axis. To practically consider the effects of freeplay structural nonlinearity, the fictitious mass method is applied to structural vibration analysis based on finite element method. Nonlinear aerodynamic flows with unsteady shock waves are also considered in supersonic flow regions. To solve the nonlinear aeroelastic governing equations including the freeplay effect, a modal-based coupled time-marching technique based on the fictitious mass method is used in the time-domain. Various aeroelastic computations have been performed for the nonlinear wing structure model. Linear and nonlinear aeroelastic analyses have been conducted and compared with each other in supersonic flow regions. Typical nonlinear limit cycle oscillations and phase plots are presented to show the complex vibration phenomena with simultaneous fluid-structure nonlinearities.