• The Journal of the Korea Contents Association
• /
• v.11 no.9
• /
• pp.1-8
• /
• 2011

A fuzzy inference technique for real-time textile animation without integration at textile model based Mass-Spring model is introduced. Until now many techniques have used the Mass-Spring model to describe elastically deformable objects like textile. A textile object is able to represent as a deformable surface composed of spring and masses, the movement of textile surface which is analysed through the numerical integration by the fundamental law of dynamics such as Hooke's law. However, the integration methods have 'instability problems' if the explicit Euler's method is applied or 'large amounts of calculation' if the implicit Euler's method is applied. A simple and fast animation technique for Mass-Spring model of a textile with fuzzy inference is proposed. The stabilized simulation result is obtained the state of each mass-point in real-time for the n of mass-points by a relatively simple calculation.

• Journal of the Society of Naval Architects of Korea
• /
• v.38 no.3
• /
• pp.1-13
• /
• 2001

RANS solver has been developed to solve the flows past a ship with a drift angle. The solver employs a finite volume method for the spatial discretization and Euler implicit method for the time integration. Turbulent flows are simulated by Spalart-Allmaras one-equation model. Developed solver is applied to analyze the hydrodynamic forces and flows of two tankers with a same forebody but different afterbodies. The computed flows and hydrodynamic forces are compared with the measured flows and captive model test data. The computed results show good agreements with experimental data and show clearly the effects of stern hull form on the hydrodynamic forces and the flows.

• Journal of the Korean Geotechnical Society
• /
• v.22 no.9
• /
• pp.45-59
• /
• 2006

The successful performance of any numerical geotechnical simulation depends on the accuracy and efficiency of the numerical implementation of constitutive model used to simulate the stress-strain (constitutive) response of the soil. The corner stone of the numerical implementation of constitutive models is the numerical integration of the incremental form of soil-plasticity constitutive equations over a discrete sequence of time steps. In this paper a well known two-surface soil plasticity model is implemented using a generalized implicit return mapping algorithm to arbitrary convex yield surfaces referred to as the Closest-Point-Projection method (CPPM). The two-surface model describes the nonlinear behavior of coarse-grained materials by incorporating a bounding surface concept together with isotropic and kinematic hardening as well as fabric formulation to account for the effect of fabric formation on the unloading response. In the course of investigating the performance of the CPPM integration method, it is proven that the algorithm is an accurate, robust, and efficient integration technique useful in finite element contexts. It is also shown that the algorithm produces a consistent tangent operator $\frac{d\sigma}{d\varepsilon}$ during the iterative process with quadratic convergence rate of the global iteration process.

• International Journal of Aeronautical and Space Sciences
• /
• v.2 no.2
• /
• pp.28-38
• /
• 2001

A numerical method for the assessment and correction of tunnel wall interference effects on forced-oscillation testing is presented. The method is based on the wall pressure signature method using computed wall pressure distributions. The wall pressure field is computed using unsteady three-dimensional full Navier-Stokes solver for a 70-degree pitching delta wing in a wind tunnel. Approximately-factorized alternate direction implicit (AF-ADI) scheme is advanced in time by solving block tri-diagonal matrices. The algebraic Baldwin-Lomax turbulence, model is included to simulate the turbulent flow effect. Also, dual time sub-iteration with, local, time stepping is implemented to improve the convergence. The computed wall pressure field is then imposed as boundary conditions for Euler re-simulation to obtain the interference flow field. The static computation shows good agreement with experiments. The dynamic computation demonstrates reasonable physical phenomena with a good convergence history. The effects of the tunnel wall in upwash and blockage are analyzed using the computed interference flow field for several reduced frequencies and amplitudes. The corrected results by pressure signature method agree well with the results of free air conditions.

• Journal of the Society of Naval Architects of Korea
• /
• v.41 no.4
• /
• pp.17-29
• /
• 2004

Internal waves are computed using a ghost fluid method on an unstructured grid. Discontinuities in density and dynamic pressure are captured in one cell without smearing or oscillations along a multimaterial interface. A time-accurate incompressible Navier-Stokes/Euler solver is developed based on a three-point backward difference formula for the physical time marching. Artificial compressibility is introduced with respect to pseudotime and an implicit method is used for the pseudotime iteration. To track evolution of an interface, a level set function is coupled with the governing equations. Roe's flux difference splitting method is used to calculate numerical fluxes of the coupled equations. To get higher order accuracy, dependent variables are reconstructed based on gradients which are calculated using Gauss theorem. For each edge crossing an interface, dynamic pressure is assigned for a ghost node to enforce the continuity of total pressure along the interface. Solitary internal waves are computed and the results are compared with other computational and experimental results.

• Journal of applied mathematics & informatics
• /
• v.40 no.1_2
• /
• pp.25-45
• /
• 2022

In this study, the non-standard finite difference method for the numerical solution of singularly perturbed parabolic reaction-diffusion subject to Robin boundary conditions has presented. To discretize temporal and spatial variables, we use the implicit Euler and non-standard finite difference method on a uniform mesh, respectively. We proved that the proposed scheme shows uniform convergence in time with first-order and in space with second-order irrespective of the perturbation parameter. We compute three numerical examples to confirm the theoretical findings.

• Journal of the Korean Society for Precision Engineering
• /
• v.13 no.2
• /
• pp.94-101
• /
• 1996

The equations of motion for a basic industrial robotic system which has a rigid or a flexible arm are derived by Lagrange's equation, respectively. Especially, for the deflection of the flexible arm, the assumed mode method is employed. These equations are highly nonlinear equations with nonlinear coupling between the variables of motion. In order to design the control law for the rigid-arm robot, Hunt-Su's nonlinear transformation method and Marino's feedback equivalence condition are used with linear quadratic regulator(LQR) theory. The control law for the rigid-arm robot is employed to input the desired path and to provide the required nonlinear transformations for the flexible-arm robot to follow. By using the implicit Euler method to solve the nonlinear equations, the comparison of the motions between the flexible and the rigid robots and the effect of flexibility are examined.

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

• Journal of the Korean Crystal Growth and Crystal Technology
• /
• v.18 no.5
• /
• pp.191-194
• /
• 2008

We will focus on the horizontal Bridgman growth system to analyze the transport phenomena numerically, because the simple furnace system and the confined growth environment allow for the precise understanding of the transport phenomena in solidification process. In conventional melt growth process, the dopant concentration tends to vary significantly along the crystal. In this work, we propose the modification of crucible geometry for improving the productivity of silicon single-crystal growth by controlling axial specific resistivity distribution. Numerical analysis has been performed to study the transport phenomena of dopant impurities in conventional and proposed Bridgman silicon growth using the finite element method and implicit Euler time integration. It has been demonstrated using mathematical models and by numerical analysis that proposed method is useful for obtaining crystals with superior uniformity along the growth direction at a lower cost than can be obtained by the conventional melt growth process.

• 한국전산유체공학회:학술대회논문집
• /
• 1997.10a
• /
• pp.17-22
• /
• 1997

The performances of the CRAYT3E and CRAYC90 were compared in the point of aerodynamics. The CRAYC90 with and without the highest vector option was run, respectively. The CRAYT3E was run with various processors (from 1pe to 32pes). The communication utilities of MPI and SHMEM were used to inform the boundary data to the other processors. The DADI Euler solver, which is implicit scheme and use central difference method, was used. The domain decomposition method was also used. As the result, the CRAYC90 with the highest vector option is 5.7 times faster than the CRAYT3E with 1 processor. However, because of the scalability of the CRAYT3E, the CRAYT3E with more than 6 processors is faster than CRAYC90. In case that 32 processors used, the CRAYT3E is 6 times faster than CRAYC90 with the highest vector option.