• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.6
• /
• pp.1510-1517
• /
• 1995

A theoretical method is described to simulate the propagation of turbulent premixed flames in a closed vessel. The objective is to develop and test an efficient technique to predict the propagation speed of flame as well as the geometric structure of the flame surfaces. Flame is advected by the statistically generated turbulent flow field and propagates as a wave by solving twodimensional Hamilton-Jacobi equation. In the simulation of the unburned gas flow field, following turbulence properties were satisfied: mean velocity field, turbulence intensities, spatial and temporal correlations of velocity fluctuations. It is assumed that these properties are not affected by the expansion of the burned gas region. Predictions were compared with existing experimental data for flames propagating in a closed vessel charged with hydrogen/air mixture with various turbulence intensities and Reynolds numbers. Comparisons were made in flame radius growth rate, rms flame radius fluctuations, and average perimeter and fractal dimensions of the flame boundaries. Two dimensional time dependent simulation resulted in correct trends of the measured flame data. The reasonable behavior and high efficiency proves the usefulness of this method in difficult problems of flame propagation such as in internal combustion engines.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.14 no.6
• /
• pp.1357-1363
• /
• 1994

A Monte-Carlo simulation model is developed to predict size distribution produced by the coalescence of air bubbles in turbulent shear f1ow. The simulation consists of generating a population of air bubbles into the initial positions at each time step and tracking them by simulating motions and checking collisions. The radial displacement of air bubbles in the simulation model is produced by numerically solving an advective diffusion equation. Longitudinal displacements are generated from the logarithmic flow velovity distribution and the bubble rise velocity. Collision of air bubbles for each time step is detected by a geometric test using their relative positions at the beginning of the time step and relative displacements during the time step. At the end of the time step, the total number of bubbles, their positions, and sizes are updated. The computer program is coded such that minimum storages for sizes and positions of bubbles are required.

• Structural Engineering and Mechanics
• /
• v.27 no.6
• /
• pp.713-739
• /
• 2007

Nowadays computers can perform symbolic computations in addition to mere number crunching operations for which they were originally designed. Symbolic computation opens up exciting possibilities in Structural Mechanics and engineering. Classical areas have been increasingly neglected due to the advent of computers as well as general purpose finite element software. But now, classical analysis has reemerged as an attractive computer option due to the capabilities of symbolic computation. The repetitive cycles of simultaneous - equation sets required by the finite element technique can be eliminated by solving a single set in symbolic form, thus generating a truly closed-form solution. This consequently saves in data preparation, storage and execution time. The power of Symbolic computation is demonstrated by six examples by applying symbolic computation 1) to solve coupled shear wall 2) to generate beam element matrices 3) to find the natural frequency of a shear frame using transfer matrix method 4) to find the stresses of a plate subjected to in-plane loading using Levy's approach 5) to draw the influence surface for deflection of an isotropic plate simply supported on all sides 6) to get dynamic equilibrium equations from Lagrange equation. This paper also presents yet another computationally efficient and accurate numerical method which is based on the concept of derivative of a function expressed as a weighted linear sum of the function values at all the mesh points. Again this method is applied to solve the problems of 1) coupled shear wall 2) lateral buckling of thin-walled beams due to moment gradient 3) buckling of a column and 4) static and buckling analysis of circular plates of uniform or non-uniform thickness. The numerical results obtained are compared with those available in existing literature in order to verify their accuracy.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2007.05a
• /
• pp.859-865
• /
• 2007

This paper presents a method to analyze the unbalance response of a high speed polygon mirror scanner motor supported by sintered bearing and flexible supporting structures by using the finite element method and the mode superposition method. The appropriate finite element equations for polygon mirror are described by rotating annular sector element using Kirchhoff plate theory and von Karman non-linear strain, and its rigid body motion is also considered. The rotating components except for the polygon mirror are modeled by Timoshenko beam element including the gyroscopic effect. The flexible supporting structures are modeled by using a 4-node tetrahedron element and 4-node shell element with rotational degrees of freedom. Finite element equations of each component of the polygon mirror scanner motor and the flexible supporting structures are consistently derived by satisfying the geometric compatibility in the internal boundary between each component. The rigid link constraints are also imposed at the interface area between sleeve and sintered bearing to describe the physical motion at this interface. A global matrix equation obtained by assembling the finite element equations of each substructure is transformed to a state-space matrix-vector equation, and both damped natural frequencies and modal damping ratios are calculated by solving the associated eigenvalue problem by using the restarted Arnoldi iteration method. Unbalance responses in time and frequency domain are performed by superposing the eigenvalues and eigenvectors from the free vibration analysis. The validity of the proposed method is verified by comparing the simulated unbalance response with the experimental results. This research also shows that the flexibility of supporting structures plays an important role in determining the unbalance response of the polygon mirror scanner motor.

• Journal of Digital Contents Society
• /
• v.16 no.6
• /
• pp.909-916
• /
• 2015

This paper presents a method to display the 3D vector fields by analyzing phase space. This method is based on the connections between ordinary differential equations and the topology of vector fields. The phase space analysis should be geometric interpolation of an autonomous system of equation in the form of the phase space. Every solution of it system of equations corresponds not to a curve in a space, but the motion of a point along the curve. This analysis is the basis of this paper. This new method is required to decompose the hexahedral cell into five or six tetrahedral cells for 3D vector fields. The critical points can be easily found by solving a simple linear system for each tetrahedron. The tangent curves can be integrated by finding the intersection points of an integral curve traced out by the general solution of each tetrahedron and plane containing a face of the tetrahedron.

• Journal for History of Mathematics
• /
• v.20 no.3
• /
• pp.167-180
• /
• 2007

In this study, we explored the role of analysis in the algebra development. For this, we classified ancient geometric analysis into an analysis by reduction and a Pappusian problematic analysis. this shows that both analyses have the function of reduction. Pappus' analysis consists of four steps; transformation, resolution, construction, demonstration. The transformation, by which conditions of given problem is transformed into other conditions which suggest a problem-solving, seems to be a kind of reduction. Mathematicians created new problems as a result of the reductional function of analysis, and became to see mathematics in the different view. An analytical thinking was a background at the birth of symbolic algebra, the reductional function of analysis played an important role in the development of symbolic algebra.

• Journal of the Korean Institute of Telematics and Electronics
• /
• v.18 no.5
• /
• pp.35-45
• /
• 1981

An asymptotic solution of electromagnetic waves scattered by a right-angled dielectric wedge for plane wave incidence is obtained. Scattered fields are constructed by waves reflected and refracted from dielectric interfaces (geometric-optical fields) and a cylindrical wave diffracted from the edge. The edge diffracted field is obtained by adding a correction to the edge diffraction of physical optics approximation, where the correction field is calculated by solving a dual series equation amenable to simple numerical calculation. Validity of this result is assured by two limits of relative dielectric constant $\varepsilon$ of the wedge. The total asymptotic field calculated results in a Rawlins' Neumann series solution for small $\varepsilon$, and the edge diffraction pattern is shown to approach that of a perfectly conducting wedge for large $\varepsilon$. Calculated field patterns are presented and the accuracy of physical optics approximation is discussed.

• The Korean Journal of Applied Statistics
• /
• v.29 no.6
• /
• pp.1129-1145
• /
• 2016

Functional neuro-connectivity is one of the main issues in brain science in the sense that it is closely related to neurodynamics in the brain. In the paper, we choose fMRI as a main form of response data to brain activity due to its high resolution. We review methods for analyzing functional neuro-connectivity assuming that measurements are made on physiological responses to neuron activation. This means that we deal with a state-space and measurement model, where the state-space model is assumed to represent neurodynamics. Analysis methods and their interpretation should vary subject to what was measured. We included analysis results of real fMRI data by applying a high-dimensional autoregressive model, which indicated that different neurodynamics were required for solving different types of geometric problems.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.4
• /
• pp.737-752
• /
• 1992

Analysis of elasto-dynamic deformation of flexible linkage mechanism is conducted using the finite element method. The equations of motion of the system are derived from the static structural problem in which dynamic inertia, gravitational and driving forces are treated as external loads. Linear spring model is included in the formulation of equation of motions to represent the effects of deformation of elastic bearings of revolute joints on the system behavior. A computer program is constructed and applied to analyze a specific crank-lever 4-bar mechanism. The algorithm of the program is as follows. First, the natural frequencies and the mode shapes of the system are calculated by solving the eigenproblem of the mechanism system which can be considered as a static structure by assuming the input shaft (crank shaft) to be fixed at any given configuration of mechanism. And finally, the elasto-dynamic deformation of the whole system is obtained using mode superposition method for the case of constant input speed. The effect of geometric stiffness on the mechamism is included in the program with the axial forces of links obtained through the quasi-static displacement analysis. It is found that the geometric stiffness exerts an important effect upon the elasto-dynamic behavior of the flexible linkage mechanism. Elastic deformation of bearing lowers the natural frequencies of the system, resulting smaller elastic displacement at the mid-point of the links and bigger elestic displacement at the ends of the links than rigid bearing. The above investigation of flexible linkage mechanism shows that the effects of the elastic deformation of bearing on the mechanism should be considered to design the mechanism which satisfies more preciously the purpose and the condition of design.

• Investigative Magnetic Resonance Imaging
• /
• v.23 no.4
• /
• pp.328-340
• /
• 2019

Purpose: To generate phase images with free of motion-induced artifact and susceptibility-induced distortion using 3D radial ultrashort TE (UTE) MRI. Materials and Methods: The field map was theoretically derived by solving Laplace's equation with appropriate boundary conditions, and used to simulate the image distortion in conventional spin-warp MRI. Manufacturer's 3D radial imaging sequence was modified to acquire maximum number of radial spokes in a given time, by removing the spoiler gradient and sampling during both rampup and rampdown gradient. Spoke direction randomly jumps so that a readout gradient acts as a spoiling gradient for the previous spoke. The custom raw data was reconstructed using a homemade image reconstruction software, which is programmed using Python language. The method was applied to a phantom and in-vivo human brain and abdomen. The performance of UTE was compared with 3D GRE for phase mapping. Local phase mapping was compared with T₂* mapping using UTE. Results: The phase map using UTE mimics true field-map, which was theoretically calculated, while that using 3D GRE revealed both motion-induced artifact and geometric distortion. Motion-free imaging is particularly crucial for application of phase mapping for abdomen MRI, which typically requires multiple breathold acquisitions. The air pockets, which are caught within the digestive pathway, induce spatially varying and large background field. T₂* map, that was calculated using UTE data, suffers from non-uniform T₂* value due to this background field, while does not appear in the local phase map of UTE data. Conclusion: Phase map generated using UTE mimicked the true field map even when non-zero susceptibility objects were present. Phase map generated by 3D GRE did not accurately mimic the true field map when non-zero susceptibility objects were present due to the significant field distortion as theoretically calculated. Nonetheless, UTE allows for phase maps to be free of susceptibility-induced distortion without the use of any post-processing protocols.