• Journal of the Korean Mathematical Society
• /
• v.59 no.3
• /
• pp.519-548
• /
• 2022

In this paper, a stabilized-penalized collocated finite volume (SPCFV) scheme is developed and studied for the stationary generalized Navier-Stokes equations with mixed Dirichlet-traction boundary conditions modelling an incompressible biological fluid flow. This method is based on the lowest order approximation (piecewise constants) for both velocity and pressure unknowns. The stabilization-penalization is performed by adding discrete pressure terms to the approximate formulation. These simultaneously involve discrete jump pressures through the interior volume-boundaries and discrete pressures of volumes on the domain boundary. Stability, existence and uniqueness of discrete solutions are established. Moreover, a convergence analysis of the nonlinear solver is also provided. Numerical results from model tests are performed to demonstrate the stability, optimal convergence in the usual L² and discrete H¹ norms as well as robustness of the proposed scheme with respect to the choice of the given traction vector.

• The Journal of the Acoustical Society of Korea
• /
• v.16 no.1E
• /
• pp.29-34
• /
• 1997

Numerical prediction of aerodynamic noise radiated by subsonic rotors are carried out. A computer program has been developed which incorporates both the discrete frequency noise as well as the broadband noise arising from the ingestion of turbulence. Acoustic analogy is used in conjunction with Homicz's formulation of turbulence ingestion noise. Formulation 1A of Farassat is used to enhance the numerical analysis performance of Ffowcs-Williams Hawkings equation by eliminating the numericla time differentiation. Homicz's trubulence ingestion noise prediction technique is used to understand the characteristics of broadband noise radiated by isotropic trubulence in gestion. Numerical predictions are carried out for a number of rotor configurations and compared with experimental data. Monopole consideration of transonic rotor agrees well with both the experimental data and the linear theory. Noise radiation characteristics of rotor at lifting hover are investigated utilizing simple blade loading obtained by thin wing section theory. By incorporating discrete noise prediction of steady loading with broadband spectrum, much better agreement with experimental data is obtained in the low frequency region. The contributions from different noise mechanisms can also be analyzed through this method.

• Interaction and multiscale mechanics
• /
• v.6 no.2
• /
• pp.107-125
• /
• 2013

A strain smoothing meshfree formulation with stabilized conforming nodal integration is presented for modeling the consolidation process in saturated porous media. In the present method, nodal strain smoothing is consistently introduced into the meshfree approximation of strain and pore pressure gradient variables associated with the saturated porous media. Meanwhile, in order to achieve a consistent numerical implementation, a smoothing approximation of the meshfree shape function within a nodal representative domain is also proposed in the stiffness construction. The resulting discrete system of equations is all expressed in smoothed nodal measures that are very efficient for numerical evaluation. Subsequently the space-time fully discrete equations are further established by the generalized trapezoidal rule for time integration. The effectiveness of the proposed meshfree consolidation analysis method is systematically illustrated by several benchmark problems.

• Journal of Ship and Ocean Technology
• /
• v.7 no.3
• /
• pp.34-55
• /
• 2003

Explicit Kutta condition approximation, proved useful in existing time-domain solver of the unsteady propeller problem, requires a specified functional behavior of the vorticity in space near the trailing edge. In this paper, the strength of the discrete vortices is controlled to have a specified behavior in space in the frequency domain approach. A new formulation is introduced and is implemented for analysis of a lifting surface of a rectangular planform. Sample computations carried out according to the new formulation compares well with that of existing unsteady lifting problem in the time domain.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.4
• /
• pp.509-517
• /
• 1997

Results of flat plate compressible boundary layer calculation, based on discrete formulation of DSMC method, are presented in low Mach number and low Knudsen number range. The free stream is a uniform flow of pure nitrogen at various Mach numbers in low pressures (i.e. rarefied gas). Complete thermal accommodation and diffuse molecular reflections are used as the wall boundary condition, replacing unreal no-slip condition used in continuum calculations. In the discrete formulation of DSMC method, there is no need to use ad hoc assumptions on transport properties like viscosity and thermal conductivity, instead viscosity is calculated from values of other field variables (velocity and shear stress). Also the results are compared with existing self-similar continuum solutions. In all Mach number cases computed, velocity slip is most pronounced in regions near the leading edge where continuum formulation renders the solution singular. As the boundary layer develops further downstream, velocity slips asymptote to values that are between 10 to 20% of the magnitude of free stream velocity. When the free stream number density is reduced, so the gas more rarefied, the velocity slip increases as expected.

• 한국전산유체공학회:학술대회논문집
• /
• 1996.05a
• /
• pp.74-84
• /
• 1996

Characteristic boundary conditions are discussed in conjunction with a flux-difference splitting formulation as modified from Roe's linearization. Details of how one can implement the characteristic boundary conditions which are compatible with the discrete formulation at interior points are given for different types of boundaries including subsonic outflow and adiabatic wall. The latter conditions are demonstrated through computation of supersonic ogive-cylinder flow at high angle of attack and the computed wall pressure distribution is compared with experiment.

• Industrial Engineering and Management Systems
• /
• v.8 no.3
• /
• pp.155-161
• /
• 2009

We focus on discrete event systems with a structure of parallel processing, synchronization, and no-concurrency. We use max-plus algebra, which is an effective approach for controller design for this type of system, for modeling and formulation. Since a typical feature of this type of system is that the initial schedule is frequently changed due to unpredictable disturbances, we use a simple model and numerical examples to examine the possibility of applying the concepts of the feeding buffer and the project buffer of critical chain project management (CCPM) on max-plus linear discrete event systems in order to control the occurrence of an undesirable state change. The application of a CCPM-based framework on a max-plus linear discrete event system was proven to be effective.

• Structural Engineering and Mechanics
• /
• v.26 no.2
• /
• pp.191-210
• /
• 2007

The three different approaches (reformulations) of evolution strategies (ESs) have been proposed in the literature as extensions of the technique for solving discrete problems. This study implements an extensive research on application, evaluation and comparison of them in discrete optimum design of steel frames. A unified formulation is first developed to explain these approaches, so that differences and similarities between their inherent search mechanisms can clearly be identified. Two examples from practical design of steel frames are studied next to measure their performances in locating the optimum. Extensive numerical experimentations are performed in both examples to facilitate a statistical analysis of their convergence characteristics. The results obtained are presented in the histograms demonstrating the distribution of the best designs located by each approach. In addition, an average improvement of the best design during the course of evolution is plotted in each case to compare their relative convergence rates.

• Computers and Concrete
• /
• v.12 no.1
• /
• pp.19-36
• /
• 2013

In the finite element analysis of reinforced concrete structures, discrete representation of the steel reinforcing bars is considered advantageous over smeared representation because of the more realistic modelling of their bond-slip behaviour. However, there is up to now limited research on how to simulate the dowel action of discrete reinforcing bars, which is an important component of shear transfer in cracked concrete structures. Herein, a numerical model for the dowel action of discrete reinforcing bars is developed. It features derivation of the dowel stiffness based on the beam-on-elastic-foundation theory and direct assemblage of the dowel stiffness matrix into the stiffness matrices of adjoining concrete elements. The dowel action model is incorporated in a nonlinear finite element program based on secant stiffness formulation and application to deep beams tested by others demonstrates that the incorporation of dowel action can improve the accuracy of the finite element analysis.

• Steel and Composite Structures
• /
• v.47 no.3
• /
• pp.365-374
• /
• 2023

This paper presents a new scheme for constructing locking-free finite elements in thick and thin plates, called Discrete Shear Gap element (DSG), using multiphase material topology optimization for triangular elements of Reissner-Mindlin plates. Besides, common methods are also presented in this article, such as quadrilateral element (Q4) and reduced integration method. Moreover, when the plate gets too thin, the transverse shear-locking problem arises. To avoid that phenomenon, the stabilized discrete shear gap technique is utilized in the DSG3 system stiffness matrix formulation. The accuracy and efficiency of DSG are demonstrated by the numerical examples, and many superior properties are presented, such as being a strong competitor to the common kind of Q4 elements in the static topology optimization and its computed results are confirmed against those derived from the three-node triangular element, and other existing solutions.