• Proceedings of the Korean Society of Coastal and Ocean Engineers Conference
• /
• 2003.08a
• /
• pp.208-215
• /
• 2003

Convective-diffusion equations appear in various disciplines such as hydrology, chemical engineering and oceanography dealing with the transport problem of scalar quantities. Since it is nonlinear, numerical methods are generally used to obtain its solution. Very limited number of analytical solutions are available usually in cases when the convective velocity is constant or has a simple functional form (for some collection of the solutions, see Noye, 1987). There is however a continuing need to develop analytical solutions because of its practical importance. Analytical solutions of the convection-diffusion equation are valuable not only for the better understanding on the transport process but the verification of numerical schemes. (omitted)

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.29 no.9 s.240
• /
• pp.1049-1056
• /
• 2005

A conservative pressure-based finite-volume numerical method has been developed for computing flow and heat transfer by using an unstructured grid system. The method admits arbitrary convex polyhedra. Care is taken in the discretization and solution procedures to avoid formulations that are cell-shape-specific. A collocated variable arrangement formulation is developed, i.e. all dependent variables such as pressure and velocity are stored at cell centers. Gradients required for the evaluation of diffusion fluxes and for second-order-accurate convective operators are found by a novel second-order accurate spatial discretization. Momentum interpolation is used to prevent pressure checkerboarding and the SIMPLE algorithm is used for pressure-velocity coupling. The resulting set of coupled nonlinear algebraic equations is solved by employing a segregated approach, leading to a decoupled set of linear algebraic equations fer each dependent variable, with a sparse diagonally dominant coefficient matrix. These equations are solved by an iterative preconditioned conjugate gradient solver which retains the sparsity of the coefficient matrix, thus achieving a very efficient use of computer resources.

• Korea-Australia Rheology Journal
• /
• v.17 no.2
• /
• pp.47-62
• /
• 2005

The present study deals with a mathematical model describing the dynamic response of heat and mass transfer in blood flow through bifurcated arteries under stenotic condition. The geometry of the bifurcated arterial segment possessing constrictions in both the parent and the daughter arterial lumen frequently appearing in the diseased arteries causing malfunction of the cardiovascular system, is formulated mathematically with the introduction of the suitable curvatures at the lateral junction and the flow divider. The blood flowing through the artery is treated to be Newtonian. The nonlinear unsteady flow phenomena is governed by the Navier-Stokes equations while those of heat and mass transfer are controlled by the heat conduction and the convection-diffusion equations respectively. All these equations together with the appropriate boundary conditions describing the present biomechanical problem following the radial coordinate transformation are solved numerically by adopting finite difference technique. The respective profiles of the flow field, the temperature and the concentration and their distributions as well are obtained. The influences of the stenosis, the arterial wall motion and the unsteady behaviour of the system in terms of the heat and mass transfer on the blood stream in the entire arterial segment are high­lighted through several plots presented at the end of the paper in order to illustrate the applicability of the present model under study.

• Structural Engineering and Mechanics
• /
• v.15 no.6
• /
• pp.669-683
• /
• 2003

The stochastic optimal nonlinear control of coupled adjacent building structures is studied based on the stochastic dynamical programming principle and the stochastic averaging method. The coupled structures with control devices under random seismic excitation are first condensed to form a reduced-order structural model for the control analysis. The stochastic averaging method is applied to the reduced model to yield stochastic differential equations for structural modal energies as controlled diffusion processes. Then a dynamical programming equation for the energy processes is established based on the stochastic dynamical programming principle, and solved to determine the optimal nonlinear control law. The seismic response mitigation of the coupled structures is achieved through the structural energy control and the dimension of the optimal control problem is reduced. The seismic excitation spectrum is taken into account according to the stochastic dynamical programming principle. Finally, the nonlinear controlled structural response is predicted by using the stochastic averaging method and compared with the uncontrolled structural response to evaluate the control efficacy. Numerical results are given to demonstrate the response mitigation capabilities of the proposed stochastic optimal control method for coupled adjacent building structures.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.3
• /
• pp.426-434
• /
• 2000

This paper deals with a generalized similarity solution for the one-dimensional solidification of ternary or higher-order multicomponent alloys. The present approach not only retains the existing features of binary systems such as temperature- solute coupling, shrinkage-induced flow, solid-liquid property differences, and finite back diffusion, but also is capable of handling a multicomponent alloy without restrictions on the partition coefficient and microsegregation parameter. For an alloy of N-solute species, governing equations in the mushy region reduce to (N+2) nonlinear ordinary differential equations via similarity transformation, which are to be solved along with the closed-form solutions for the solid and liquid regions. A linearized correction scheme adopted in the solution procedure facilitates to determine the solidus and liquidus positions stably. The result for a sample ternary alloy agrees excellently with the numerical prediction as well as the reported similarity solution. Additional calculations are also presented to show the utility of this study. Finally, it is concluded that the present analysis includes the previous analytical approaches as subsets.

• Korea-Australia Rheology Journal
• /
• v.14 no.4
• /
• pp.161-174
• /
• 2002

A new technique for numerical calculation of viscoelastic flow based on the combination of Neural Net-works (NN) and Brownian Dynamics simulation or Stochastic Simulation Technique (SST) is presented in this paper. This method uses a "universal approximator" based on neural network methodology in combination with the kinetic theory of polymeric liquid in which the stress is computed from the molecular configuration rather than from closed form constitutive equations. Thus the new method obviates not only the need for a rheological constitutive equation to describe the fluid (as in the original Calculation Of Non-Newtonian Flows: Finite Elements St Stochastic Simulation Techniques (CONNFFESSIT) idea) but also any kind of finite element-type discretisation of the domain and its boundary for numerical solution of the governing PDE's. As an illustration of the method, the time development of the planar Couette flow is studied for two molecular kinetic models with finite extensibility, namely the Finitely Extensible Nonlinear Elastic (FENE) and FENE-Peterlin (FENE-P) models.P) models.

• Nonlinear Functional Analysis and Applications
• /
• v.27 no.3
• /
• pp.679-689
• /
• 2022

In this paper, we establish sufficient conditions for the existence and uniqueness of solution for a class of initial boundary value problems with Dirichlet condition in regard to a category of fractional-order partial differential equations. The results are established by a method based on the theorem of Lax Milligram.

• The Bulletin of The Korean Astronomical Society
• /
• v.46 no.2
• /
• pp.47.3-48
• /
• 2021

We studied the large scale dynamo process in a system forced by helical magnetic field. The dynamo process is basically nonlinear, but can be linearized with 𝛼&𝛽 coefficients and large scale magnetic field $\bar{B}$. This is very useful to the investigation of solar (stellar) dynamo. A coupled semi-analytic equations based on statistical mechanics are used to investigate the exact evolution of 𝛼&𝛽. This equation set needs only magnetic helicity ${\bar{H}}_M({\equiv}{\langle}{\bar{A}}{\cdot}{\bar{B}}{\rangle},\;{\bar{B}}={\nabla}{\times}{\bar{A}})$ and magnetic energy ${\bar{E}}_M({\equiv}{\langle}{\bar{B}}^2{\rangle}/2)$. They are fundamental physics quantities that can be obtained from the dynamo simulation or observation without any artificial modification or assumption. 𝛼 effect is thought to be related to magnetic field amplification. However, in reality the averaged 𝛼 effect decreases very quickly without a significant contribution to ${\bar{B}}$ field amplification. Conversely, 𝛽 effect contributing to the magnetic diffusion maintains a negative value, which plays a key role in the amplification with Laplacian ∇²(= - k²) for the large scale regime. In addition, negative magnetic diffusion accounts for the attenuation of plasma kinetic energy E_V(= 〈 U² 〉/2) (U: plasma velocity) when the system is saturated. The negative magnetic diffusion is from the interaction of advective term - U • ∇ B from magnetic induction equation and the helical velocity field. In more detail, when 'U' is divided into the poloidal component U_pol and toroidal one U_tor in the absence of reflection symmetry, they interact with - B • ∇ U and - U • ∇ B from ∇ × 〈 U × B 〉 leading to 𝛼 effect and (negative) 𝛽 effect, respectively. We discussed this process using the theoretical method and intuitive field structure model supported by the simulation result.

• Journal of Ocean Engineering and Technology
• /
• v.7 no.1
• /
• pp.13-24
• /
• 1993

The dispersion in the oscillatory flow generated by gravitational waves above the spatially periodic repples is studied. The steady parts of equations describing the orbit of the passive particle in a two dimensional field are assumed to be simply trigonometric functions. From the view point of nonlinear dynamics, the motion of the particle is chaotic under externally time-periodic perturbations which come from the wave motion. Two cases considered here are; (i) shallow water, and (ii) deep water approximation.

• Current Optics and Photonics
• /
• v.4 no.4
• /
• pp.380-390
• /
• 2020

In this study, the traditional electron-blocking layer (EBL) in (In,Ga)N/GaN light-emitting diodes is replaced by a circular EBL that is the same size as the n-pad. The three-dimensional (3D) nonlinear Poisson, drift-diffusion, and continuity equations are adopted to simulate current transport in the LED and its characteristics. The results indicate that the local carrier-density distribution obtained for the circular EBL design is more uniform than that for the traditional EBL design. This improves the uniformity of local radiative recombination and local internal quantum efficiency (IQE) at high injection levels, which leads to a higher lumped IQE and lower efficiency droop. With the circular EBL, the lumped IQE is higher in the outer active region and lower in the active region under the n-pad. Since most emissions from the active region under the n-pad are absorbed by the n-pad, obviously, an LED with a circular EBL will have a higher external quantum efficiency (EQE). The results also show that this LED works at lower applied voltages.