• 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 the Earthquake Engineering Society of Korea
• /
• v.27 no.4
• /
• pp.163-170
• /
• 2023

In the fluid-structure interaction analysis, the finite element formulation is performed for the wave equation for dynamic fluid pressure, and the dynamic pressure is defined as a degree of freedom at the fluid nodes. Therefore, to connect the fluid to the structure, it is necessary to connect the degree of freedom of fluid dynamic pressure and the degree of freedom of structure displacement through an interface element derived from the relationship between dynamic pressure and displacement. The previously proposed fluid-structure interface elements use conformal finite element meshes in which the fluid and structure match. However, it is challenging to construct conformal meshes when complex models, such as water purification plants and wastewater treatment facilities, are models. Therefore, to increase modeling convenience, a method is required to model the fluid and structure domains by independent finite element meshes and then connect them. In this study, two fluid-structure interface elements, one based on constraints and the other based on the integration of nonsmooth functions, are proposed in nonconformal finite element meshes for structures and fluids, and their accuracy is verified.

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

• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.307-310
• /
• 2002

A two-dimensional numerical method for inviscid two-fluid flows with evolution of density interface is developed, and an initially stationary two-dimensional fluid sheet surrounded by another fluid is studied. The Interface between two fluids is modeled as a vertex sheet, and the flow field u÷th the evolution of interface is solved by using vortex-in-cell/front-tracking method. The edge of the sheet Is pulled back into the sheet due to surface tension and a blob is formed at the edge. This blob and fluid sheet are connected by a thin neck. In the inviscid limit, such process of the blob and neck formation is examined in detail and their kinematic characteristics are summarized with dimensionless parameters. The edge recedes at $V=1.06({\sigma}/{\rho}h)^{0.5}$ and the capillary wave Propagating into the fluid sheet must be considered for bettor understanding of the edge receding.

• Solar Energy
• /
• v.5 no.2
• /
• pp.11-19
• /
• 1985

In this paper, the interface stability not to occur mixing and entrainment between the adjacent layers has been studied in the case of the selective withdrawal of a stratum and the injection in stratified fluid formed by the density difference in a small solar pond. There are stability parameter, Richardson number, Rayleigh number and Froude number as the parameters governing stability in order to measure the interface stability on the stratified fluid. The model which could measure the interface stability on the stratified fluid was the small solar pond composed by 1 meters wide, 2 meters high, and 5 meters long. In order to measure the interface stability on the stratified fluid at the inlet port, the middle section and the outlet port, Richardson number, Rayleigh number, and Froude number involved in the parameters governing the stability were calculated by means of the data resulted from the test of the study on hydrodynamic stability between the convective and nonconvective layers in that solar pond. Richardson number written by the ratio of inertia force to buoyancy force can be used in order to measure the stability on the stratified fluid related to the buoyancy force generated from the injection of fluid. Rayleigh number written by the product of Grashof number by Prandtl number can be used in order to measure the stability of the fluid related to the heat flux and diffusivity of viscosity. Froude number written by the ratio of gravity force to inertia force can be used in order to measure the stability of the nonhomogeneous fluid related to the density difference. As the result of calculating the parameters governing stability, the interface stability on the stratified fluid couldn't be identified below the 70cm height from the bottom of the solar pond, but it could be identified above the 70cm height from it at the inlet port, the middle section and the outlet port. When compared with such the three parameters as Richardson number, Rayleigh number, Froude number, the calculated result was in accord with them at inlet port, the middle section and the outlet port. Henceforth, it is learned that even though any of the three parameters is used for the purpose of measuring the interface stability on the stratified fluid, the result will be the same with them. It is concluded that all the use of Richardson number, Rayleigh number, and Froude number, is desirable and infallible to measure the interface stability on the stratified fluid in the case of considering the exist of the fluid flow and the heat flux like the model of the solar pond.

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

• Structural Engineering and Mechanics
• /
• v.76 no.5
• /
• pp.619-629
• /
• 2020

In current study, surface/interface effects for pull-in voltage and viscous fluid velocity effects on dimensionless natural frequency (DNF) of fluid-conveying piezoelectric nanosensor (FCPENS) subjected to direct electrostatic voltage DC with nonlinear excitation, harmonic force and also viscoelastic foundation (visco-pasternak medium and structural damping) are investigated using Gurtin-Murdoch surface/interface (GMSIT) theory. For this analysis, Hamilton's principles, the assumed mode method combined with Lagrange-Euler's are used for the governing equations and boundary conditions. The effects of surface/interface parameters of FCPENS such as Lame's constants (λI,S, μI,S), residual stress (τ0I,S), piezoelectric constants (e31psk,e32psk) and mass density (ρI,S) are considered for analysis of dimensionless natural frequency respect to viscous fluid velocity u̅f and pull-in voltage V̅DC.

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

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.11 no.1
• /
• pp.572-583
• /
• 2019

This paper aims to assess the applicability of the Runge Kutta Discontinuous Galerkin-Direct Ghost Fluid Method to the internal explosion inside a water-filled tube, which previously was studied by many researchers in separate works. Once the explosive charge located at the inner center of the water-filled tube explodes, the tube wall is subjected to an extremely high intensity fluid loading and deformed. The deformation causes a modification of the field of fluid flow in the region near the water-structure interface so that has substantial influence on the response of the structure. To connect the structure and the fluid, valid data exchanges along the interface are essential. Classical fluid structure interaction simulations usually employ a matched meshing scheme which discretizes the fluid and structure domains using a single mesh density. The computational cost of fluid structure interaction simulations is usually governed by the structure because the size of time step may be determined by the density of structure mesh. The finer mesh density, the better solution, but more expensive computational cost. To reduce such computational cost, a non-matched meshing scheme which allows for different mesh densities is employed. The coupled numerical approach of this paper has fewer difficulties in the implementation and computation, compared to gas dynamics based approach which requires complicated analytical manipulations. It can also be applied to wider compressible, inviscid fluid flow analyses often found in underwater explosion events.

• Interaction and multiscale mechanics
• /
• v.5 no.3
• /
• pp.287-318
• /
• 2012

We apply a partitioned-solution (iterative-staggered) coupling method based on a fixed Eulerian mesh with the level set function to a large-deformation fluid-structure interaction (FSI) problem where a large-deformable thin structure moves in a high-speed flow field, as an airbag does during deployment. This method combines advanced fluid and structure solvers-specifically, the constrained interpolation profile finite element method (CIP-FEM) for fluid Eulerian mesh and large-deformable structural elements for Lagrangian structural mesh. We express the large-deformable interface as a zero isosurface by the level set function, and introduce virtual nodes with level sets and structural normal velocities to generate the level set function according to the large-deformable interfacial geometry and enforce the kinematic condition at the interface. The virtual nodes are located in the direction normal to the structural mesh. It is confirmed that application of the method to unfolded airbag deployment simulation shows the adequacy of the method.