• Proceedings of the KSME Conference
• /
• 2008.11a
• /
• pp.442-447
• /
• 2008

This paper presents a new approach to simulate fluid-solid interaction problems involving non-matching interfaces. The coupling between fluid and solid domains with dissimilar finite element meshes consisting of 4-node quadrilateral elements is achieved by using the interface element method (IEM). Conditions of compatibility between fluid and solid meshes are satisfied exactly by introducing the interface elements defined on interfacing regions. Importantly, a consistent transfer of loads through matching interface element meshes guarantees the present method to be an efficient approach of the solution strategy to fluid-solid interaction problems. An arbitrary Lagrangian-Eulerian (ALE) description is adopted for the fluid domain, while for the solid domain an updated Lagrangian formulation is considered to accommodate finite deformations of an elastic structure. The stabilized equal order velocity-pressure elements for incompressible flows are used in the motion of fluids. Fully coupled equations are solved simultaneously in a single computational domain. Numerical results are presented for fluid-solid interaction problems involving nonmatching interfaces to demonstrate the effectiveness of the methodology.

• Journal of Mechanical Science and Technology
• /
• v.20 no.12
• /
• pp.2230-2241
• /
• 2006

This paper presents a sole application of boundary element method to the conjugate heat transfer problem of thermally developing laminar flow in a thick walled pipe when the fluid velocities are fully developed. Due to the coupled mechanism of heat conduction in the solid region and heat convection in the fluid region, two separate solutions in the solid and fluid regions are sought to match the solid-fluid interface continuity condition. In this method, the dual reciprocity boundary element method (DRBEM) with the axial direction marching scheme is used to solve the heat convection problem and the conventional boundary element method (BEM) of axisymmetric model is applied to solve the heat conduction problem. An iterative and numerically stable BEM solution algorithm is presented, which uses the coupled interface conditions explicitly instead of uncoupled conditions. Both the local convective heat transfer coefficient at solid-fluid interface and the local mean fluid temperature are initially guessed and updated as the unknown interface thermal conditions in the iterative solution procedure. Two examples imposing uniform temperature and heat flux boundary conditions are tested in thermally developing region and compared with analytic solutions where available. The benchmark test results are shown to be in good agreement with the analytic solutions for both examples with different boundary conditions.

• The Pure and Applied Mathematics
• /
• v.19 no.2
• /
• pp.147-169
• /
• 2012

We introduce a Heaviside-function formulation of the interaction between incompressible two-phase fluid and a non-deformable solid. Fluid and solid interact in two ways : fluid satises the Dirichlet boundary condition imposed by the velocity field of solid, and solid is accelerated by the surface traction exerted by fluid. The two-way couplings are formulated by the Heaviside function to the interface between solid and fluid. The cumbersome treatment of interface is taken care of by the Heaviside function, and the interaction is discretized in a simple manner. The discretization results in a stable and accurate projection method.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.18 no.5
• /
• pp.451-457
• /
• 2006

Conjugate heat transfer (CHT) is the simultaneous, coupled heat transfer within a fluid and an adjoining solid, and the interface treatment plays an important role in its analysis, particularly when using unstructured grid system. In the present paper a new solid-fluid interface treatment in CHT analysis is presented and applied to two typical CHT problems, i.e. natural convections in both concentric thick-walled cylinders and cavity with a centered solid body. The present interface treatment for unstructured mesh clearly demonstrates the same accuracy and robustness as that for typical structured mesh.

• Interaction and multiscale mechanics
• /
• v.2 no.1
• /
• pp.45-68
• /
• 2009

An aim of the study is to develop an efficient numerical simulation technique that can handle the two-scale analysis of fluid permeation filters fabricated by the partial sintering technique of small spherical ceramics. A solid-fluid mixture homogenization method is introduced to predict the mechanical characters such as rigidity and permeability of the porous ceramic filters from the micro-scale geometry and configuration of partially-sintered particles. An extended finite element (X-FE) discretization technique based on the enriched interpolations of respective characteristic functions at fluid-solid interfaces is proposed for the non-interface-fitted mesh solution of the micro-scale analysis that needs non-slip condition at the interface between solid and fluid phases of the unit cell. The homogenization and localization performances of the proposed method are shown in a typical two-dimensional benchmark problem whose model has a hole in center. Three-dimensional applications to the body-centered cubic (BCC) and face-centered cubic (FCC) unit cell models are also shown in the paper. The 3D application is prepared toward the computer-aided optimal design of ceramic filters. The accuracy and stability of the X-FEM based method are comparable to those of the standard interface-fitted FEM, and are superior to those of the voxel type FEM that is often used in such complex micro geometry cases.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.03a
• /
• pp.221-226
• /
• 2008

Numerical simulation is performed for microdroplet deposition on a pre-patterned micro-structure. The level-set method for tracking the liquid-gas interface is extended to treat the immersed (or irregular-shaped) solid surface. The no-slip condition at the fluid-solid interface as well as the matching conditions at the liquid-gas interface is accurately imposed by incorporating the ghost fluid approach based on a sharp-interface representation. The method is further extended to treat the contact angle condition at an immersed solid surface. The present computation of a patterning process using microdroplet ejection demonstrates that the multiphase characteristics between the liquid-gas-solid phases can be used to improve the patterning accuracy.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.10a
• /
• pp.221-226
• /
• 2008

Numerical simulation is performed for microdroplet deposition on a pre-patterned micro-structure. The level-set method for tracking the liquid-gas interface is extended to treat the immersed (or irregular-shaped) solid surface. The no-slip condition at the fluid-solid interface as well as the matching conditions at the liquid-gas interface is accurately imposed by incorporating the ghost fluid approach based on a sharp-interface representation. The method is further extended to treat the contact angle condition at an immersed solid surface. The present computation of a patterning process using microdroplet ejection demonstrates that the multiphase characteristics between the liquid-gas-solid phases can be used to improve the patterning accuracy.

• The KSFM Journal of Fluid Machinery
• /
• v.13 no.2
• /
• pp.5-11
• /
• 2010

The nanofluidics is characterized by a large surface-to-volume ratio, so that the surface properties strongly affect the flow resistance. We present here the results showing that the effect of wetting properties and the surface roughness may considerably reduce the friction of fluid past the boundaries. For a simple fluid flowing over hydrophilic and hydrophobic surfaces, the influences of surface roughness are investigated by the nonequilibrium molecular dynamics (NEMD) simulations. The fluid slip at near a solid surface highly depends on the wall-fluid interaction. For hydrophobic surfaces, apparent fluid slips are observed on smooth and rough surfaces. The solid wall is modeled as a rough atomic sinusoidal wall. The effects on the boundary condition of the roughness characteristics are given by the period and amplitude of the sinusoidal wall. It was found that the slip velocity for wetting conditions at interface decreases with increasing effects of surface roughness. The results show the surface rougheness and wettability determines the slip or no-slip boundary conditions. The surface roughness geometry shows significant effects on the boundary conditions at the interface.

• Coupled systems mechanics
• /
• v.1 no.1
• /
• pp.19-37
• /
• 2012

In this work, the iterative coupling of finite element and boundary element methods for the investigation of coupled fluid-fluid, solid-solid and fluid-solid wave propagation models is reviewed. In order to perform the coupling of the two numerical methods, a successive renewal of the variables on the common interface between the two sub-domains is performed through an iterative procedure until convergence is achieved. In the case of local nonlinearities within the finite element sub-domain, it is straightforward to perform the iterative coupling together with the iterations needed to solve the nonlinear system. In particular, a more efficient and stable performance of the coupling procedure is achieved by a special formulation that allows to use different time steps in each sub-domain. Optimized relaxation parameters are also considered in the analyses, in order to speed up and/or to ensure the convergence of the iterative process.

• Proceedings of the KSME Conference
• /
• 2004.11a
• /
• pp.416-418
• /
• 2004

Tracking of evolving solid-fluid interfaces is treated using level set method and MLS-based finite element with variable nodes. Several applications will be illustrated to demonstrate the effectiveness of the present scheme