Coupled CFD-DEM simulation of interfacial fluid–particle interaction during binder jet 3D printing

Coupled CFD-DEM simulation of interfacial fluid–particle interaction during binder jet 3D printing

바인더 제트 3D 프린팅 중 계면 유체-입자 상호 작용에 대한 CFD-DEM 결합 시뮬레이션

Joshua J. Wagner, C. Fred Higgs III

https://doi.org/10.1016/j.cma.2024.116747

Abstract

The coupled dynamics of interfacial fluid phases and unconstrained solid particles during the binder jet 3D printing process govern the final quality and performance of the resulting components. The present work proposes a computational fluid dynamics (CFD) and discrete element method (DEM) framework capable of simulating the complex interfacial fluid–particle interaction that occurs when binder microdroplets are deposited into a powder bed. The CFD solver uses a volume-of-fluid (VOF) method for capturing liquid–gas multifluid flows and relies on block-structured adaptive mesh refinement (AMR) to localize grid refinement around evolving fluid–fluid interfaces. The DEM module resolves six degrees of freedom particle motion and accounts for particle contact, cohesion, and rolling resistance. Fully-resolved CFD-DEM coupling is achieved through a fictitious domain immersed boundary (IB) approach. An improved method for enforcing three-phase contact lines with a VOF-IB extension technique is introduced. We present several simulations of binder jet primitive formation using realistic process parameters and material properties. The DEM particle systems are experimentally calibrated to reproduce the cohesion behavior of physical nickel alloy powder feedstocks. We demonstrate the proposed model’s ability to resolve the interdependent fluid and particle dynamics underlying the process by directly comparing simulated primitive granules with one-to-one experimental counterparts obtained from an in-house validation apparatus. This computational framework provides unprecedented insight into the fundamental mechanisms of binder jet 3D printing and presents a versatile new approach for process parameter optimization and defect mitigation that avoids the inherent challenges of experiments.

바인더 젯 3D 프린팅 공정 중 계면 유체 상과 구속되지 않은 고체 입자의 결합 역학이 결과 구성 요소의 최종 품질과 성능을 좌우합니다. 본 연구는 바인더 미세액적이 분말층에 증착될 때 발생하는 복잡한 계면 유체-입자 상호작용을 시뮬레이션할 수 있는 전산유체역학(CFD) 및 이산요소법(DEM) 프레임워크를 제안합니다.

CFD 솔버는 액체-가스 다중유체 흐름을 포착하기 위해 VOF(유체량) 방법을 사용하고 블록 구조 적응형 메쉬 세분화(AMR)를 사용하여 진화하는 유체-유체 인터페이스 주위의 그리드 세분화를 국지화합니다. DEM 모듈은 6개의 자유도 입자 운동을 해결하고 입자 접촉, 응집력 및 구름 저항을 설명합니다.

완전 분해된 CFD-DEM 결합은 가상 도메인 침지 경계(IB) 접근 방식을 통해 달성됩니다. VOF-IB 확장 기술을 사용하여 3상 접촉 라인을 강화하는 향상된 방법이 도입되었습니다. 현실적인 공정 매개변수와 재료 특성을 사용하여 바인더 제트 기본 형성에 대한 여러 시뮬레이션을 제시합니다.

DEM 입자 시스템은 물리적 니켈 합금 분말 공급원료의 응집 거동을 재현하기 위해 실험적으로 보정되었습니다. 우리는 시뮬레이션된 기본 과립과 내부 검증 장치에서 얻은 일대일 실험 대응물을 직접 비교하여 프로세스의 기본이 되는 상호 의존적인 유체 및 입자 역학을 해결하는 제안된 모델의 능력을 보여줍니다.

이 계산 프레임워크는 바인더 제트 3D 프린팅의 기본 메커니즘에 대한 전례 없는 통찰력을 제공하고 실험에 내재된 문제를 피하는 공정 매개변수 최적화 및 결함 완화를 위한 다용도의 새로운 접근 방식을 제시합니다.

Introduction

Binder jet 3D printing (BJ3DP) is a powder bed additive manufacturing (AM) technology capable of fabricating geometrically complex components from advanced engineering materials, such as metallic superalloys and ultra-high temperature ceramics [1], [2]. As illustrated in Fig. 1(a), the process is comprised of many repetitive print cycles, each contributing a new cross-sectional layer on top of a preceding one to form a 3D CAD-specified geometry. The feedstock material is first delivered from a hopper to a build plate and then spread into a thin layer by a counter-rotating roller. After powder spreading, a print head containing many individual inkjet nozzles traverses over the powder bed while precisely jetting binder microdroplets onto select regions of the spread layer. Following binder deposition, the build plate lowers by a specified layer thickness, leaving a thin void space at the top of the job box that the subsequent powder layer will occupy. This cycle repeats until the full geometries are formed layer by layer. Powder bed fusion (PBF) methods follow a similar procedure, except they instead use a laser or electron beam to selectively melt and fuse the powder material. Compared to PBF, binder jetting offers several distinct advantages, including faster build rates, enhanced scalability for large production volumes, reduced machine and operational costs, and a wider selection of suitable feedstock materials [2]. However, binder jetted parts generally possess inferior mechanical properties and reduced dimensional accuracy [3]. As a result, widescale adoption of BJ3DP to fabricate high-performance, mission-critical components, such as those common to the aerospace and defense sectors, is contingent on novel process improvements and innovations [4].

A major obstacle hindering the advancement of BJ3DP is our limited understanding of how various printing parameters and material properties collectively influence the underlying physical mechanisms of the process and their effect on the resulting components. To date, the vast majority of research efforts to uncover these relationships have relied mainly on experimental approaches [5], [6], [7], [8], [9], [10], [11], [12], [13], [14], [15], [16], [17], [18], [19], which are often expensive and time-consuming and have inherent physical restrictions on what can be measured and observed. For these reasons, there is a rapidly growing interest in using computational models to circumvent the challenges of experimental investigations and facilitate a deeper understanding of the process’s fundamental phenomena. While significant progress has been made in developing and deploying numerical frameworks aimed at powder spreading [20], [21], [22], [23], [24], [25], [26], [27] and sintering [28], [29], [30], [31], [32], simulating the interfacial fluid–particle interaction (IFPI) in the binder deposition stage is still in its infancy. In their exhaustive review, Mostafaei et al. [2] point out the lack of computational models capable of resolving the coupled fluid and particle dynamics associated with binder jetting and suggest that the development of such tools is critical to further improving the process and enhancing the quality of its end-use components.

We define IFPI as a multiphase flow regime characterized by immiscible fluid phases separated by dynamic interfaces that intersect the surfaces of moving solid particles. As illustrated in Fig. 1(b), an elaborate IFPI occurs when a binder droplet impacts the powder bed in BJ3DP. The momentum transferred from the impacting droplet may cause powder compaction, cratering, and particle ejection. These ballistic disturbances can have deleterious effects on surface texture and lead to the formation of large void spaces inside the part [5], [13]. After impact, the droplet spreads laterally on the bed surface and vertically into the pore network, driven initially by inertial impact forces and then solely by capillary action [33]. Attractive capillary forces exerted on mutually wetted particles tend to draw them inward towards each other, forming a packed cluster of bound particles referred to as a primitive [34]. A single-drop primitive is the most fundamental building element of a BJ3DP part, and the interaction leading to its formation has important implications on the final part characteristics, such as its mechanical properties, resolution, and dimensional accuracy. Generally, binder droplets are deposited successively as the print head traverses over the powder bed. The traversal speed and jetting frequency are set such that consecutive droplets coalesce in the bed, creating a multi-drop primitive line instead of a single-drop primitive granule. The binder must be jetted with sufficient velocity to penetrate the powder bed deep enough to provide adequate interlayer binding; however, a higher impact velocity leads to more pronounced ballistic effects.

A computational framework equipped to simulate the interdependent fluid and particle dynamics in BJ3DP would allow for unprecedented observational and measurement capability at temporal and spatial resolutions not currently achievable by state-of-the-art imaging technology, namely synchrotron X-ray imaging [13], [14], [18], [19]. Unfortunately, BJ3DP presents significant numerical challenges that have slowed the development of suitable modeling frameworks; the most significant of which are as follows:

  • 1.Incorporating dynamic fluid–fluid interfaces with complex topological features remains a nontrivial task for standard mesh-based CFD codes. There are two broad categories encompassing the methods used to handle interfacial flows: interface tracking and interface capturing [35]. Interface capturing techniques, such as the popular volume-of-fluid (VOF) [36] and level-set methods [37], [38], are better suited for problems with interfaces that become heavily distorted or when coalescence and fragmentation occur frequently; however, they are less accurate in resolving surface tension and boundary layer effects compared to interface tracking methods like front-tracking [39], arbitrary Lagrangian–Eulerian [40], and space–time finite element formulations [41]. Since interfacial forces become increasingly dominant at decreasing length scales, inaccurate surface tension calculations can significantly deteriorate the fidelity of IFPI simulations involving <100 μm droplets and particles.
  • 2.Dynamic powder systems are often modeled using the discrete element method (DEM) introduced by Cundall and Strack [42]. For IFPI problems, a CFD-DEM coupling scheme is required to exchange information between the fluid and particle solvers. Fully-resolved CFD-DEM coupling suggests that the flow field around individual particle surfaces is resolved on the CFD mesh [43], [44]. In contrast, unresolved coupling volume averages the effect of the dispersed solid phase on the continuous fluid phases [45], [46], [47], [48]. Comparatively, the former is computationally expensive but provides detailed information about the IFPI in question and is more appropriate when contact line dynamics are significant. However, since the pore structure of a powder bed is convoluted and evolves with time, resolving such solid–fluid interfaces on a computational mesh presents similar challenges as fluid–fluid interfaces discussed in the previous point. Although various algorithms have been developed to deform unstructured meshes to accommodate moving solid surfaces (see Bazilevs et al. [49] for an overview of such methods), they can be prohibitively expensive when frequent topology changes require mesh regeneration rather than just modification through nodal displacement. The pore network in a powder bed undergoes many topology changes as particles come in and out of contact with each other, constantly closing and opening new flow channels. Non-body-conforming structured grid approaches that rely on immersed boundary (IB) methods to embed the particles in the flow field can be better suited for such cases [50]. Nevertheless, accurately representing these complex pore geometries on Cartesian grids requires extremely high mesh resolutions, which can impose significant computational costs.
  • 3.Capillary effects depend on the contact angle at solid–liquid–gas intersections. Since mesh nodes do not coincide with a particle surface when using an IB method on structured grids, imposing contact angle boundary conditions at three-phase contact lines is not straightforward.

While these issues also pertain to PBF process modeling, resolving particle motion is generally less crucial for analyzing melt pool dynamics compared to primitive formation in BJ3DP. Therefore, at present, the vast majority of computational process models of PBF assume static powder beds and avoid many of the complications described above, see, e.g., [51], [52], [53], [54], [55], [56], [57], [58], [59]. Li et al. [60] presented the first 2D fully-resolved CFD-DEM simulations of the interaction between the melt pool, powder particles, surrounding gas, and metal vapor in PBF. Following this work, Yu and Zhao [61], [62] published similar melt pool IFPI simulations in 3D; however, contact line dynamics and capillary forces were not considered. Compared to PBF, relatively little work has been published regarding the computational modeling of binder deposition in BJ3DP. Employing the open-source VOF code Gerris [63], Tan [33] first simulated droplet impact on a powder bed with appropriate binder jet parameters, namely droplet size and impact velocity. However, similar to most PBF melt pool simulations described in the current literature, the powder bed was fixed in place and not allowed to respond to the interacting fluid phases. Furthermore, a simple face-centered cubic packing of non-contacting, monosized particles was considered, which does not provide a realistic pore structure for AM powder beds. Building upon this approach, we presented a framework to simulate droplet impact on static powder beds with more practical particle size distributions and packing arrangements [64]. In a study similar to [33], [64], Deng et al. [65] used the VOF capability in Ansys Fluent to examine the lateral and vertical spreading of a binder droplet impacting a fixed bimodal powder bed with body-centered packing. Li et al. [66] also adopted Fluent to conduct 2D simulations of a 100 μm diameter droplet impacting substrates with spherical roughness patterns meant to represent the surface of a simplified powder bed with monosized particles. The commercial VOF-based software FLOW-3D offers an AM module centered on process modeling of various AM technologies, including BJ3DP. However, like the above studies, particle motion is still not considered in this codebase. Ur Rehman et al. [67] employed FLOW-3D to examine microdroplet impact on a fixed stainless steel powder bed. Using OpenFOAM, Erhard et al. [68] presented simulations of different droplet impact spacings and patterns on static sand particles.

Recently, Fuchs et al. [69] introduced an impressive multipurpose smoothed particle hydrodynamics (SPH) framework capable of resolving IFPI in various AM methods, including both PBF and BJ3DP. In contrast to a combined CFD-DEM approach, this model relies entirely on SPH meshfree discretization of both the fluid and solid governing equations. The authors performed several prototype simulations demonstrating an 80 μm diameter droplet impacting an unconstrained powder bed at different speeds. While the powder bed responds to the hydrodynamic forces imparted by the impacting droplet, the particle motion is inconsistent with experimental time-resolved observations of the process [13]. Specifically, the ballistic effects, such as particle ejection and bed deformation, were drastically subdued, even in simulations using a droplet velocity ∼ 5× that of typical jetting conditions. This behavior could be caused by excessive damping in the inter-particle contact force computations within their SPH framework. Moreover, the wetted particles did not appear to be significantly influenced by the strong capillary forces exerted by the binder as no primitive agglomeration occurred. The authors mention that the objective of these simulations was to demonstrate their codebase’s broad capabilities and that some unrealistic process parameters were used to improve computational efficiency and stability, which could explain the deviations from experimental observations.

In the present paper, we develop a novel 3D CFD-DEM numerical framework for simulating fully-resolved IFPI during binder jetting with realistic material properties and process parameters. The CFD module is based on the VOF method for capturing binder–air interfaces. Surface tension effects are realized through the continuum surface force (CSF) method with height function calculations of interface curvature. Central to our fluid solver is a proprietary block-structured AMR library with hierarchical octree grid nesting to focus enhanced grid resolution near fluid–fluid interfaces. The GPU-accelerated DEM module considers six degrees of freedom particle motion and includes models based on Hertz-Mindlin contact, van der Waals cohesion, and viscoelastic rolling resistance. The CFD and DEM modules are coupled to achieve fully-resolved IFPI using an IB approach in which Lagrangian solid particles are mapped to the underlying Eulerian fluid mesh through a solid volume fraction field. An improved VOF-IB extension algorithm is introduced to enforce the contact angle at three-phase intersections. This provides robust capillary flow behavior and accurate computations of the fluid-induced forces and torques acting on individual wetted particles in densely packed powder beds.

We deploy our integrated codebase for direct numerical simulations of single-drop primitive formation with powder beds whose particle size distributions are generated from corresponding laboratory samples. These simulations use jetting parameters similar to those employed in current BJ3DP machines, fluid properties that match commonly used aqueous polymeric binders, and powder properties specific to nickel alloy feedstocks. The cohesion behavior of the DEM powder is calibrated based on the angle of repose of the laboratory powder systems. The resulting primitive granules are compared with those obtained from one-to-one experiments conducted using a dedicated in-house test apparatus. Finally, we demonstrate how the proposed framework can simulate more complex and realistic printing operations involving multi-drop primitive lines.

Section snippets

Mathematical description of interfacial fluid–particle interaction

This section briefly describes the governing equations of fluid and particle dynamics underlying the CFD and DEM solvers. Our unified framework follows an Eulerian–Lagrangian approach, wherein the Navier–Stokes equations of incompressible flow are discretized on an Eulerian grid to describe the motion of the binder liquid and surrounding gas, and the Newton–Euler equations account for the positions and orientations of the Lagrangian powder particles. The mathematical foundation for

CFD solver for incompressible flow with multifluid interfaces

This section details the numerical methodology used in our CFD module to solve the Navier–Stokes equations of incompressible flow. First, we introduce the VOF method for capturing the interfaces between the binder and air phases. This approach allows us to solve the fluid dynamics equations considering only a single continuum field with spatial and temporal variations in fluid properties. Next, we describe the time integration procedure using a fractional-step projection algorithm for

DEM solver for solid particle dynamics

This section covers the numerical procedure for tracking the motion of individual powder particles with DEM. The Newton–Euler equations (Eqs. (10), (11)) are ordinary differential equations (ODEs) for which many established numerical integrators are available. In general, the most challenging aspects of DEM involve processing particle collisions in a computationally efficient manner and dealing with small time step constraints that result from stiff materials, such as metallic AM powders. The

Unified CFD-DEM solver

The preceding sections have introduced the CFD and DEM solution algorithms separately. Here, we discuss the integrated CFD-DEM solution algorithm and related details.

Binder jet process modeling and validation experiments

In this section, we deploy our CFD-DEM framework to simulate the IFPI occurring during the binder droplet deposition stage of the BJ3DP process. The first simulations attempt to reproduce experimental single-drop primitive granules extracted from four nickel alloy powder samples with varying particle size distributions. The experiments are conducted with a dedicated in-house test apparatus that allows for the precision deposition of individual binder microdroplets into a powder bed sample. The

Conclusions

This paper introduces a coupled CFD-DEM framework capable of fully-resolved simulation of the interfacial fluid–particle interaction occurring in the binder jet 3D printing process. The interfacial flow of binder and surrounding air is captured with the VOF method and surface tension effects are incorporated using the CSF technique augmented by height function curvature calculations. Block-structured AMR is employed to provide localized grid refinement around the evolving liquid–gas interface.

CRediT authorship contribution statement

Joshua J. Wagner: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Software, Visualization, Writing – original draft, Writing – review & editing. C. Fred Higgs III: Conceptualization, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Writing – original draft, Writing – review & editing.

Declaration of competing interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Acknowledgments

This work was supported by a NASA Space Technology Research Fellowship, United States of America, Grant No. 80NSSC19K1171. Partial support was also provided through an AIAA Foundation Orville, USA and Wilbur Wright Graduate Award, USA . The authors would like to gratefully acknowledge Dr. Craig Smith of NASA Glenn Research Center for the valuable input he provided on this project.

References (155)

Fig. 1 Oscillation of a free surface due to the step reduction of gravity acceleration from kzi ≈ 9.81 to kz ≈ 0

Reorientation of Cryogenic Fluids Upon Step Reduction of Gravity

단계적 중력 감소 시 극저온 유체의 방향 전환

Malte Stief∗, Jens Gerstmann∗∗, and Michael E. Dreyer∗∗∗
ZARM, Center of Applied Space Technology and Microgravity, University of Bremen, Am Fallturm, D-28359 Bremen
Experiments to observe the surface oscillation of cryogenic liquids have been performed with liquid nitrogen inside a 50 mm
diameter right circular cylinder. The surface oscillation is driven by the capillary force that becomes dominant after a sudden
reduction of the gravity acceleration acting on the liquid. The experiments show differences from the speculated behavior and
enables one to observe new features.

Introduction and motivation

최근 몇 년 동안 Bremen의 낙하탑에서 중력의 단계적 감소 시 방향 재지향 거동과 표면 진동을 조사하기 위해 수많은 실험이 수행되었습니다[1]. 이 실험의 원리는 그림 1에 나와 있습니다.

그림 1의 왼쪽에 표시된 것처럼 오른쪽 원형 원통형 용기에 테스트 액체를 레벨 h0까지 채웁니다. 처음에 액체는 정지 상태이며 중앙에서 평평한 인터페이스를 형성합니다.

초기 중력 가속도 kzi ≈ 9.81 [m/s2]와 결과적으로 높은 BOND 수(Bo = ρkziR2/σ)로 인해 실린더의 대칭축에서. 낙하탑에서 실험 캡슐의 방출에 의해 확립된 μ-중력 환경 kz ≈ 0 [m/s2]로의 갑작스러운 전환과 함께 자유 표면은 진동 운동으로 새로운 평형 구성을 찾기 시작합니다(그림의 오른쪽) 1). 이러한 움직임은 그림 1의 중앙에 스케치되어 있습니다.

표면 진동의 구동력은 접착력과 결합된 표면 장력이며, 댐핑은 액체의 점도에 의해 제어됩니다. 위치가 zw인 벽에서 접촉선의 이동은 접촉각 γ에 의해 제어됩니다. 접촉각이 작은 액체용 γ ≈ 0◦

In recent years numerous experiments have been carried out to investigate the reorientation behavior and surface oscillations upon step reduction of gravity at the drop tower in Bremen [1]. The principals of these experiments are shown in figure 1. A right circular cylindrical container is filled up to the level h0 with the test liquid, as shown on the left of figure 1. Initially the liquid is quiescent and forms a flat interface at the center, in the symmetry axis of the cylinder, due to the initial gravity acceleration kzi ≈ 9.81 [m/s2] and the resulting high BOND number (Bo = ρkziR2/σ). With the sudden transition to the µ-gravity environment kz ≈ 0 [m/s2], which is established by the release of the experiment capsular in the drop tower, the free surface is initiated to search its new equilibrium configuration (right side of figure 1) with an oscillatory motion. These movements are sketched in the center of figure 1. The driving force for the surface oscillation is the surface tension in combination with the adhesion force where the damping is controlled by the viscosity of the liquid. The movement of the contact line at the wall, with its position zw, is governed by the contact angle γ. For liquids with small contact angle γ ≈ 0◦

Fig. 1 Oscillation of a free surface due to the step reduction of gravity acceleration from kzi ≈ 9.81 to kz ≈ 0
Fig. 1 Oscillation of a free surface due to the step reduction of gravity acceleration from kzi ≈ 9.81 to kz ≈ 0
Fig. 2 Experiment picture-series showing the oscillation of the free surface at different times for a 50 mm diameter cylinder.
Fig. 2 Experiment picture-series showing the oscillation of the free surface at different times for a 50 mm diameter cylinder.

References

[1] M. Michaelis, Kapillarinduzierte Schwingungen freier Fl¨ussigkeitsoberfl¨achen, Dissertation Universit¨at Bremen, Fortschritt-Berichte
Nr. 454 (VDI Verlag, D¨usseldorf, 2003).

Figure 1.2: Left panel: 3D CAD drawing of a printhead prototype showing (a) the melting unit, (b) the filter units, (c) the reservoir, (d) the static pressure hose, (e) the central part, and (f) the electronic driving supply. Image retrieved from [8]. Right panel: A schematic showing a single nozzle uint in the central part (e) of the printhead shown in the left panel.

Lattice Boltzmann method for contact line dynamics

접촉선 역학을 위한 Lattice Boltzmann 방법

ter verkrijging van de graad van doctor aan de
Technische Universiteit Eindhoven, op gezag van de
rector magnificus prof.dr.ir. C.J. van Duijn, voor een
commissie aangewezen door het College voor
Promoties, in het openbaar te verdedigen
op woensdag 7 mei 2014 om 16:00 uur

Introduction

움직이는 접촉선은 본질적으로 어디에나 존재하며, 표면에 미끄러지는 물방울은 우리가 일상에서 만나는 일반적인 예입니다. 유체 역학의 접촉선은 일반적으로 액체, 고체 및 주변 공기/증기 사이의 공통 경계라고합니다.

최근 미세 유체 공학의 발전으로 인해 접촉 라인의 역학을 제어하는 힘과 흐름 조건에 대한 근본적인 이해와 기술에 대한 많은 요구가 제기되었습니다. 이 논문은 접촉선의 물리학, 분석 및 수치 모델링 및 고무적인 산업 기하학과 관련된 측면을 포함합니다.

동기를 부여하는 산업 응용 분야는 이머전 리소그래피 (ASML)와 잉크젯 노즐 (Océ)의 프린트 헤드입니다. 이 두 가지 문제는 몇 가지 특징적인 길이 및 시간 척도, 고도로 구부러진 유체 인터페이스, 다상 흐름 및 복잡한 경계 조건을 포함하므로 분석 및 수치 연구가 어렵습니다.

포토 리소그래피는 서브 마이크론 정확도로 마스크에서 실리콘 웨이퍼로 패턴을 전송할 수 있는 복잡한 절차입니다 [1]. 포토 리소그래피 공정의 핵심 단계 중 하나는 고해상도 광학 시스템을 사용하여 실리콘 웨이퍼에 코팅 된 포토 레지스트를 DUV (심 자외선) 빛으로 노출시키는 것입니다. 광학 시스템을 사용하여 웨이퍼에 마스킹 할 수 있는 가장 작은 특징 또는 임계 치수 CD는 Rayleigh 기준으로 결정됩니다.

여기서 NA는 광학 시스템의 개구 수를 나타내고, λ는 사용 된 빛의 파장이고 k는 공정 종속 상수입니다. 광학 분야에서 광학 시스템의 개구 수 NA = n sin α는 시스템이 빛을 받아들이거나 방출 할 수 있는 각도 범위를 특성화하는 무차원 숫자입니다.

여기서 α는 렌즈의 수용 각도입니다 (0 < α <π / 2) 및 n은 렌즈와 포토 레지스트 사이의 매질의 굴절률입니다. CD의 가치가 감소하면 전자 장치가 더 작고 빨라집니다. 식에 의해 주어진 레일리 기준에 따르면. (1.1), 더 작은 CD 값은 k 또는 λ를 줄이거 나 NA를 늘림으로써 얻을 수 있습니다. 현재 KrF 및 ArF 엑시머 레이저의 경우 빛의 파장은 각각 최대 280nm 및 193nm까지 감소 될 수 있습니다 [1]. k는 분해능 향상 기술을 사용하여 0.4까지 감소 된 공정 의존 상수입니다 [2 ]. 개구 수는 sin α 또는 n을 증가시켜 증가시킬 수 있습니다.

sin α에 대한 실제 한계는 0.93으로, 이론적 한계 | sin α |에 매우 가깝습니다. ≤ 1. n을 늘리는 것이 이머전 리소그래피 사용의 기본 아이디어입니다. Immersion lithography는 렌즈와 포토 레지스트 사이의 에어 갭이 물로 대체되는 포토 리소그래피 기법입니다 (그림 1.1 (왼쪽 패널) 참조). 침지 리소그래피에 사용되는 물은 193nm 파장에 대해 1.44의 굴절률을 가진 고도로 정제 된 탈 이온수입니다 [3]. 이 굴절률 값은 분해 가능한 피처 크기의 해상도를 약 30 % 정도 증가시킵니다 [3].

이 방법은 훨씬 더 비싼 리소그래피 기술 [4]로 큰 변화를 가져 오지 않아도 된다는 장점을 가지고 더 작은 피처 크기를 달성하는 저렴한 방법입니다. 물이 웨이퍼의 포토 레지스트와 직접 접촉하기 때문에 이머전 리소그래피 기술은 주로 렌즈와 포토 레지스트의 오염 가능성과 관련된 몇 가지 문제를 야기합니다.

특히 웨이퍼 플레이트가 렌즈에 비해 Up 속도로 움직일 때 액체-공기-고체 접촉 라인도 움직입니다 (그림 1.1 (오른쪽 패널) 참조). 특정 최소 속도를 넘어 서면 전진 및 후퇴 접촉 선 (그림 1.1, 오른쪽 패널 참조)이 불안정 해지고 각각 공기를 동반하거나 액체 필름을 웨이퍼로 끌 수 있습니다 [5].

공기와 액체 필름은 결국 기포 나 액체 방울로 부서져서 리소그래피 공정에 부정적인 영향을 미칩니다. 이 논문에서 우리는 플레이트의 속도, 웨이퍼의 습윤 특성 및 주변 공기의 점도에 따라 전진 및 후퇴하는 접촉 라인의 안정성 연구에 기여했습니다.

1.1.2 Drop-on-demand inkjet printer

최신 잉크젯 인쇄 기술은 CIJ (연속 잉크젯) 및 DOD (주문형 드롭) 잉크젯의 두 가지 주요 유형으로 나눌 수 있습니다. CIJ 프린터에서 미세 노즐에서 나오는 액체 분사는 RP (Rayleigh-Plateau) 불안정성으로 인해 물방울로 분해됩니다. 이 RP 불안정성은 액체의 흐름을 정확하게 제어 할 수있는 음향 변동을 생성하는 압전 결정에 의해 유발되어 일정한 간격으로 물방울로 분해됩니다 [7].

DOD 잉크젯 프린터는 작동 원리에 따라 두 가지 범주로 더 나눌 수 있습니다 [8]. 여기서는 압전 잉크젯 (PIJ) 프린터에만 중점을 둡니다. PIJ 프린터에서 낙하 형성은 압전 소자에 의해 생성 된 압력 파에 의해 발생합니다. PIJ 프린터의 프린트 헤드 개략도가 그림 1.2에 나와 있습니다.

PIJ 프린터는 CIJ 프린터에 비해 상대적으로 느리지 만 인쇄 품질이 훨씬 더 높습니다 [7]. 프린터의 품질은 일반적으로 평방 인치당 도트 수 (dpi)로 측정되며 최신 응용 프로그램에는 더 작은 물방울 (높은 dpi)과 더 나은 정확도가 필요합니다. 방울의 정확도와 크기에 영향을 미치는 여러 요인 중에서 노즐, 노즐 플레이트의 젖음성 및 방울 형성 ​​빈도 fDOD가 중요한 역할을합니다 [8].

좋은 방울 형성을 위해 접촉 라인의 위치는 노즐 내에서 정밀하게 제어되어야 합니다. 이 논문에서는 PIJ 프린터에서 드롭 형성의 일부 측면에만 중점을 둡니다. 우리의 연구는 노즐 습윤성과 DOD 주파수가 방울 형성 ​​과정에 미치는 영향을 연구 할 수 있는 수치 도구의 개발을 목표로 합니다.

Figure 1.2: Left panel: 3D CAD drawing of a printhead prototype showing (a) the melting unit, (b) the filter units, (c) the reservoir, (d) the static pressure hose, (e) the central part, and (f) the electronic driving supply. Image retrieved from [8]. Right panel: A schematic showing a single nozzle uint in the central part (e) of the printhead shown in the left panel.
Figure 1.2: Left panel: 3D CAD drawing of a printhead prototype showing (a) the melting unit, (b) the filter units, (c) the reservoir, (d) the static pressure hose, (e) the central part, and (f) the electronic driving supply. Image retrieved from [8]. Right panel: A schematic showing a single nozzle uint in the central part (e) of the printhead shown in the left panel.
Figure 2.2: The liquid-vapor interface at the microscopic length scale obtained from a molecular dynamics (MD) simulation using Lennard-Jones potential
Figure 2.2: The liquid-vapor interface at the microscopic length scale obtained from a molecular dynamics (MD) simulation using Lennard-Jones potential. The vertical axis is in units of the molecular diameter σ and the stress shown in panel (c) is measured in /σ3 . Here,  is the energy scale corresponding to the intermolecular forces. (a) Snapshot of the liquid-vapor interface in the MD simulation. The red dotted line divides the system in two parts: Left and right. (b) Time-averaged normalized density profile ρ ∗ (z) across the interface. (c) Tangential force per unit area exerted by the left part on the right part of the system. The plot shows the difference between the normal and the tangential components of stress tensor: Π(z) = σ n − σ t . Images reproduced from [16].
Figure 2.3: Left panel: Water drops on a glass substrate
Figure 2.3: Left panel: Water drops on a glass substrate (Image source: http: // way2science. com/ molecular-theory-of-surface-tension).The red dotted line in the figure shows the position of the contact line. The shape of the big drops is affected by the force due to gravity. Right panel: Schematics of a liquid drop on a smooth non-deformable solid surface. The figure shows the contact angle, θe, in thermodynamic equilibrium.
Figure 6.1: Left panel: schematic of a single nozzle unit in the printhead
Figure 6.1: Left panel: schematic of a single nozzle unit in the printhead. Right panel: schematic of the channel-nozzle section of the printhead. The axisymmetric channel-nozzle section (right panel) is the simulation domain for our LB simulation (R = Rc).
Figure 1.1: A water droplet with a radius of 1 mm resting on a glass substrate. The surface of the droplet takes on a spherical cap shape. The contact angle θ is defined by the balance of the interfacial forces.

Effect of substrate cooling and droplet shape and composition on the droplet evaporation and the deposition of particles

기판 냉각 및 액적 모양 및 조성이 액적 증발 및 입자 증착에 미치는 영향

by Vahid Bazargan
M.A.Sc., Mechanical Engineering, The University of British Columbia, 2008
B.Sc., Mechanical Engineering, Sharif University of Technology, 2006
B.Sc., Chemical & Petroleum Engineering, Sharif University of Technology, 2006

고착 방울은 평평한 기판에 놓인 액체 방울입니다. 작은 고정 액적이 증발하는 동안 액적의 접촉선은 고정된 접촉 영역이 있는 고정된 단계와 고정된 접촉각이 있는 고정 해제된 단계의 두 가지 단계를 거칩니다. 고정된 접촉 라인이 있는 증발은 액적 내부에서 접촉 라인을 향한 흐름을 생성합니다.

이 흐름은 입자를 운반하고 접촉 선 근처에 침전시킵니다. 이로 인해 일반적으로 관찰되는 “커피 링”현상이 발생합니다. 이 논문은 증발 과정과 고착성 액적의 증발 유도 흐름에 대한 연구를 제공하고 콜로이드 현탁액에서 입자의 침착에 대한 통찰력을 제공합니다. 여기서 우리는 먼저 작은 고착 방울의 증발을 연구하고 증발 과정에서 기판의 열전도도의 중요성에 대해 논의합니다.

현재 증발 모델이 500µm 미만의 액적 크기에 대해 심각한 오류를 생성하는 방법을 보여줍니다. 우리의 모델에는 열 효과가 포함되어 있으며, 특히 증발 잠열의 균형을 맞추기 위해 액적에 열을 제공하는 기판의 열전도도를 포함합니다. 실험 결과를 바탕으로 접촉각의 진화와 관련된 접촉 선의 가상 움직임을 정의하여 고정 및 고정 해제 단계의 전체 증발 시간을 고려합니다.

우리의 모델은 2 % 미만의 오차로 500 µm보다 작은 물방울에 대한 실험 결과와 일치합니다. 또한 유한한 크기의 라인 액적의 증발을 연구하고 증발 중 접촉 라인의 복잡한 동작에 대해 논의합니다. 에너지 공식을 적용하고 접촉 선이 구형 방울의 후퇴 접촉각보다 높은 접촉각을 가진 선 방울의 두 끝에서 후퇴하기 시작 함을 보여줍니다. 그리고 라인 방울 내부의 증발 유도 흐름을 보여줍니다.

마지막으로, 계면 활성제 존재 하에서 접촉 라인의 거동을 논의하고 입자 증착에 대한 Marangoni 흐름 효과에 대해 논의합니다. 열 Marangoni 효과는 접촉 선 근처에 증착 된 입자의 양에 영향을 미치며, 기판 온도가 낮을수록 접촉 선 근처에 증착되는 입자의 양이 많다는 것을 알 수 있습니다.

Figure 1.1: A water droplet with a radius of 1 mm resting on a glass substrate. The surface of the droplet takes on a spherical cap shape. The contact angle θ is defined by the balance of the interfacial forces.
Figure 1.1: A water droplet with a radius of 1 mm resting on a glass substrate. The surface of the droplet takes on a spherical cap shape. The contact angle θ is defined by the balance of the interfacial forces.
Figure 2.1: Evaporation modes of sessile droplets on a substrate: (a) evaporation at constant contact angle (de-pinned stage) and (b) evaporation at constant contact area (pinned stage)
Figure 2.1: Evaporation modes of sessile droplets on a substrate: (a) evaporation at constant contact angle (de-pinned stage) and (b) evaporation at constant contact area (pinned stage)
Figure 2.2: A sessil droplet with its image can be profiled as the equiconvex lens formed by two intersecting spheres with radius of a.
Figure 2.2: A sessil droplet with its image can be profiled as the equiconvex lens formed by two intersecting spheres with radius of a.
Figure 2.3: The droplet life time for both evaporation modes derived from Equation 2.2.
Figure 2.3: The droplet life time for both evaporation modes derived from Equation 2.2.
Figure 2.4: A probability of escape for vapor molecules at two different sites of the surface of the droplet for diffusion controlled evaporation. The random walk path initiated from a vapor molecule is more likely to result in a return to the surface if the starting point is further away from the edge of the droplet.
Figure 2.4: A probability of escape for vapor molecules at two different sites of the surface of the droplet for diffusion controlled evaporation. The random walk path initiated from a vapor molecule is more likely to result in a return to the surface if the starting point is further away from the edge of the droplet.
Figure 2.5: Schematic of the sessile droplet on a substrate
Figure 2.5: Schematic of the sessile droplet on a substrate. The evaporation rate at the surface of the droplet is enhanced toward the edge of the droplet.
Figure 2.6: The domain mesh (a) and the solution of the Laplace equation for diffusion of the water vapor molecule with the concentration of Cv = 1.9×10−8 g/mm3 at the surface of the droplet into the ambient air with the relative humidity of 55%, i.e. φ = 0.55 (b).
Figure 2.6: The domain mesh (a) and the solution of the Laplace equation for diffusion of the water vapor molecule with the concentration of Cv = 1.9×10−8 g/mm3 at the surface of the droplet into the ambient air with the relative humidity of 55%, i.e. φ = 0.55 (b).
Figure 3.1: The portable micro printing setup. A motorized linear stage from Zaber Technologies Inc. was used to control the place and speed of the micro nozzle.
Figure 3.1: The portable micro printing setup. A motorized linear stage from Zaber Technologies Inc. was used to control the place and speed of the micro nozzle.
Figure 4.6: Temperature contours inside the substrate adjacent to the droplet
Figure 4.6: Temperature contours inside the substrate adjacent to the droplet
Figure 4.7: The effect of substrate cooling on the evaporation rate, the basic model shows the same value for all substrates.
Figure 4.7: The effect of substrate cooling on the evaporation rate, the basic model shows the same value for all substrates.

Bibliography

[1] R. G. Picknett and R. Bexon, “The evaporation of sessile or pendant drops in still air,” Journal of Colloid and Interface Science, vol. 61, pp. 336–350, Sept. 1977. → pages viii, 8, 9, 18, 42
[2] H. Y. Erbil, “Evaporation of pure liquid sessile and spherical suspended drops: A review,” Advances in Colloid and Interface Science, vol. 170, pp. 67–86, Jan. 2012. → pages 1
[3] R. Sharma, C. Y. Lee, J. H. Choi, K. Chen, and M. S. Strano, “Nanometer positioning, parallel alignment, and placement of single anisotropic nanoparticles using hydrodynamic forces in cylindrical droplets,” Nano Lett., vol. 7, no. 9, pp. 2693–2700, 2007. → pages 1, 54, 71
[4] S. Tokonami, H. Shiigi, and T. Nagaoka, “Review: Micro- and nanosized molecularly imprinted polymers for high-throughput analytical applications,” Analytica Chimica Acta, vol. 641, pp. 7–13, May 2009. →pages 71
[5] A. A. Sagade and R. Sharma, “Copper sulphide (CuxS) as an ammonia gas sensor working at room temperature,” Sensors and Actuators B: Chemical, vol. 133, pp. 135–143, July 2008. → pages
[6] W. R. Small, C. D. Walton, J. Loos, and M. in het Panhuis, “Carbon nanotube network formation from evaporating sessile drops,” The Journal of Physical Chemistry B, vol. 110, pp. 13029–13036, July 2006. → pages 71
[7] S. H. Ko, H. Lee, and K. H. Kang, “Hydrodynamic flows in electrowetting,” Langmuir, vol. 24, pp. 1094–1101, Feb. 2008. → pages 42
[8] T. T. Nellimoottil, P. N. Rao, S. S. Ghosh, and A. Chattopadhyay, “Evaporation-induced patterns from droplets containing motile and nonmotile bacteria,” Langmuir, vol. 23, pp. 8655–8658, Aug. 2007. → pages 1
[9] R. Sharma and M. S. Strano, “Centerline placement and alignment of anisotropic nanotubes in high aspect ratio cylindrical droplets of nanometer diameter,” Advanced Materials, vol. 21, no. 1, p. 6065, 2009. → pages 1, 54, 71
[10] V. Dugas, J. Broutin, and E. Souteyrand, “Droplet evaporation study applied to DNA chip manufacturing,” Langmuir, vol. 21, pp. 9130–9136, Sept. → pages 2, 71
[11] Y.-C. Hu, Q. Zhou, Y.-F. Wang, Y.-Y. Song, and L.-S. Cui, “Formation mechanism of micro-flows in aqueous poly(ethylene oxide) droplets on a substrate at different temperatures,” Petroleum Science, vol. 10, pp. 262–268, June 2013. → pages 2, 34, 54
[12] T.-S. Wong, T.-H. Chen, X. Shen, and C.-M. Ho, “Nanochromatography driven by the coffee ring effect,” Analytical Chemistry, vol. 83, pp. 1871–1873, Mar. 2011. → pages 71
[13] J.-H. Kim, S.-B. Park, J. H. Kim, and W.-C. Zin, “Polymer transports inside evaporating water droplets at various substrate temperatures,” The Journal of Physical Chemistry C, vol. 115, pp. 15375–15383, Aug. 2011. → pages 54
[14] S. Choi, S. Stassi, A. P. Pisano, and T. I. Zohdi, “Coffee-ring effect-based three dimensional patterning of Micro/Nanoparticle assembly with a single droplet,” Langmuir, vol. 26, pp. 11690–11698, July 2010. → pages
[15] D. Wang, S. Liu, B. J. Trummer, C. Deng, and A. Wang, “Carbohydrate microarrays for the recognition of cross-reactive molecular markers of microbes and host cells,” Nature biotechnology, vol. 20, pp. 275–281, Mar. PMID: 11875429. → pages 2, 54, 71
[16] H. K. Cammenga, “Evaporation mechanisms of liquids,” Current topics in materials science, vol. 5, pp. 335–446, 1980. → pages 3
[17] C. Snow, “Potential problems and capacitance for a conductor bounded by two intersecting spheres,” Journal of Research of the National Bureau of Standards, vol. 43, p. 337, 1949. → pages 9
[18] R. D. Deegan, O. Bakajin, T. F. Dupont, G. Huber, S. R. Nagel, and T. A. Witten, “Contact line deposits in an evaporating drop,” Physical Review E, vol. 62, p. 756, July 2000. → pages 10, 14, 18, 27, 53, 54, 71, 84
[19] H. Hu and R. G. Larson, “Evaporation of a sessile droplet on a substrate,” The Journal of Physical Chemistry B, vol. 106, pp. 1334–1344, Feb. 2002. → pages 12, 18, 29, 43, 44, 48, 49, 53, 61, 71, 84
[20] Y. O. Popov, “Evaporative deposition patterns: Spatial dimensions of the deposit,” Physical Review E, vol. 71, p. 036313, Mar. 2005. → pages 14, 27, 43, 44, 45, 54
[21] H. Gelderblom, A. G. Marin, H. Nair, A. van Houselt, L. Lefferts, J. H. Snoeijer, and D. Lohse, “How water droplets evaporate on a superhydrophobic substrate,” Physical Review E, vol. 83, no. 2, p. 026306,→ pages
[22] F. Girard, M. Antoni, S. Faure, and A. Steinchen, “Influence of heating temperature and relative humidity in the evaporation of pinned droplets,” Colloids and Surfaces A: Physicochemical and Engineering Aspects, vol. 323, pp. 36–49, June 2008. → pages 18
[23] Y. Y. Tarasevich, “Simple analytical model of capillary flow in an evaporating sessile drop,” Physical Review E, vol. 71, p. 027301, Feb. 2005. → pages 19, 54, 62, 72
[24] A. J. Petsi and V. N. Burganos, “Potential flow inside an evaporating cylindrical line,” Physical Review E, vol. 72, p. 047301, Oct. 2005. → pages 22, 55, 62, 68, 71
[25] A. J. Petsi and V. N. Burganos, “Evaporation-induced flow in an inviscid liquid line at any contact angle,” Physical Review E, vol. 73, p. 041201, Apr.→ pages 23, 53, 55, 72
[26] H. Masoud and J. D. Felske, “Analytical solution for stokes flow inside an evaporating sessile drop: Spherical and cylindrical cap shapes,” Physics of Fluids, vol. 21, pp. 042102–042102–11, Apr. 2009. → pages 23, 55, 62, 71, 72
[27] H. Hu and R. G. Larson, “Analysis of the effects of marangoni stresses on the microflow in an evaporating sessile droplet,” Langmuir, vol. 21, pp. 3972–3980, Apr. 2005. → pages 24, 28, 53, 54, 56, 62, 68, 71, 72, 74, 84
[28] R. Bhardwaj, X. Fang, and D. Attinger, “Pattern formation during the evaporation of a colloidal nanoliter drop: a numerical and experimental study,” New Journal of Physics, vol. 11, p. 075020, July 2009. → pages 28
[29] A. Petsi, A. Kalarakis, and V. Burganos, “Deposition of brownian particles during evaporation of two-dimensional sessile droplets,” Chemical Engineering Science, vol. 65, pp. 2978–2989, May 2010. → pages 28
[30] J. Park and J. Moon, “Control of colloidal particle deposit patterns within picoliter droplets ejected by ink-jet printing,” Langmuir, vol. 22, pp. 3506–3513, Apr. 2006. → pages 28
[31] H. Hu and R. G. Larson, “Marangoni effect reverses coffee-ring depositions,” The Journal of Physical Chemistry B, vol. 110, pp. 7090–7094, Apr. 2006. → pages 29, 74
[32] K. H. Kang, S. J. Lee, C. M. Lee, and I. S. Kang, “Quantitative visualization of flow inside an evaporating droplet using the ray tracing method,” Measurement Science and Technology, vol. 15, pp. 1104–1112, June 2004. → pages 34
[33] S. T. Beyer and K. Walus, “Controlled orientation and alignment in films of single-walled carbon nanotubes using inkjet printing,” Langmuir, vol. 28, pp. 8753–8759, June 2012. → pages 42, 71
[34] G. McHale, “Surface free energy and microarray deposition technology,” Analyst, vol. 132, pp. 192–195, Feb. 2007. → pages 42
[35] R. Bhardwaj, X. Fang, P. Somasundaran, and D. Attinger, “Self-assembly of colloidal particles from evaporating droplets: Role of DLVO interactions and proposition of a phase diagram,” Langmuir, vol. 26, pp. 7833–7842, June→ pages 42
[36] G. J. Dunn, S. K. Wilson, B. R. Duffy, S. David, and K. Sefiane, “The strong influence of substrate conductivity on droplet evaporation,” Journal of Fluid Mechanics, vol. 623, no. 1, p. 329351, 2009. → pages 44
[37] M. S. Plesset and A. Prosperetti, “Flow of vapour in a liquid enclosure,” Journal of Fluid Mechanics, vol. 78, pp. 433–444, 1976. → pages 44
[38] S. Das, P. R. Waghmare, M. Fan, N. S. K. Gunda, S. S. Roy, and S. K. Mitra, “Dynamics of liquid droplets in an evaporating drop: liquid droplet coffee stain? effect,” RSC Advances, vol. 2, pp. 8390–8401, Aug. 2012. → pages 53
[39] B. J. Fischer, “Particle convection in an evaporating colloidal droplet,” Langmuir, vol. 18, pp. 60–67, Jan. 2002. → pages 54
[40] J. L. Wilbur, A. Kumar, H. A. Biebuyck, E. Kim, and G. M. Whitesides, “Microcontact printing of self-assembled monolayers: applications in microfabrication,” Nanotechnology, vol. 7, p. 452, Dec. 1996. → pages 54
[41] T. Kawase, H. Sirringhaus, R. H. Friend, and T. Shimoda, “Inkjet printed via-hole interconnections and resistors for all-polymer transistor circuits,” Advanced Materials, vol. 13, no. 21, p. 16011605, 2001. → pages 71
[42] B.-J. de Gans, P. C. Duineveld, and U. S. Schubert, “Inkjet printing of polymers: State of the art and future developments,” Advanced Materials, vol. 16, no. 3, p. 203213, 2004. → pages 71
[43] H. Sirringhaus, T. Kawase, R. H. Friend, T. Shimoda, M. Inbasekaran, W. Wu, and E. P. Woo, “High-resolution inkjet printing of all-polymer transistor circuits,” Science, vol. 290, pp. 2123–2126, Dec. 2000. PMID:→ pages
[44] D. Soltman and V. Subramanian, “Inkjet-printed line morphologies and temperature control of the coffee ring effect,” Langmuir, vol. 24, pp. 2224–2231, Mar. 2008. → pages 54
[45] R. Tadmor and P. S. Yadav, “As-placed contact angles for sessile drops,” Journal of Colloid and Interface Science, vol. 317, pp. 241–246, Jan. 2008. → pages 56
[46] J. Drelich, “The significance and magnitude of the line tension in three-phase (solid-liquid-fluid) systems,” Colloids and Surfaces A: Physicochemical and Engineering Aspects, vol. 116, pp. 43–54, Sept. 1996. → pages 56
[47] R. Tadmor, “Line energy, line tension and drop size,” Surface Science, vol. 602, pp. L108–L111, July 2008. → pages 69
[48] C.-H. Choi and C.-J. C. Kim, “Droplet evaporation of pure water and protein solution on nanostructured superhydrophobic surfaces of varying heights,” Langmuir, vol. 25, pp. 7561–7567, July 2009. → pages 71
[49] K. F. Baughman, R. M. Maier, T. A. Norris, B. M. Beam, A. Mudalige, J. E. Pemberton, and J. E. Curry, “Evaporative deposition patterns of bacteria from a sessile drop: Effect of changes in surface wettability due to exposure to a laboratory atmosphere,” Langmuir, vol. 26, pp. 7293–7298, May 2010.
[50] D. Brutin, B. Sobac, and C. Nicloux, “Influence of substrate nature on the evaporation of a sessile drop of blood,” Journal of Heat Transfer, vol. 134, pp. 061101–061101, May 2012. → pages 71
[51] D. Pech, M. Brunet, P.-L. Taberna, P. Simon, N. Fabre, F. Mesnilgrente, V. Condra, and H. Durou, “Elaboration of a microstructured inkjet-printed carbon electrochemical capacitor,” Journal of Power Sources, vol. 195, pp. 1266–1269, Feb. 2010. → pages 71
[52] J. Bachmann, A. Ellies, and K. Hartge, “Development and application of a new sessile drop contact angle method to assess soil water repellency,” Journal of Hydrology, vol. 231232, pp. 66–75, May 2000. → pages 71
[53] H. Y. Erbil, G. McHale, and M. I. Newton, “Drop evaporation on solid surfaces: constant contact angle mode,” Langmuir, vol. 18, no. 7, pp. 2636–2641, 2002. → pages
[54] X. Fang, B. Li, J. C. Sokolov, M. H. Rafailovich, and D. Gewaily, “Hildebrand solubility parameters measurement via sessile drops evaporation,” Applied Physics Letters, vol. 87, pp. 094103–094103–3, Aug.→ pages
[55] Y. C. Jung and B. Bhushan, “Wetting behaviour during evaporation and condensation of water microdroplets on superhydrophobic patterned surfaces,” Journal of Microscopy, vol. 229, no. 1, p. 127140, 2008. → pages 71
[56] J. Drelich, J. D. Miller, and R. J. Good, “The effect of drop (bubble) size on advancing and receding contact angles for heterogeneous and rough solid surfaces as observed with sessile-drop and captive-bubble techniques,”
Journal of Colloid and Interface Science, vol. 179, pp. 37–50, Apr. 1996. →pages 72, 75
[57] D. Bargeman and F. Van Voorst Vader, “Effect of surfactants on contact angles at nonpolar solids,” Journal of Colloid and Interface Science, vol. 42, pp. 467–472, Mar. 1973. → pages 73
[58] J. Menezes, J. Yan, and M. Sharma, “The mechanism of alteration of macroscopic contact angles by the adsorption of surfactants,” Colloids and Surfaces, vol. 38, no. 2, pp. 365–390, 1989. → pages
[59] T. Okubo, “Surface tension of structured colloidal suspensions of polystyrene and silica spheres at the air-water interface,” Journal of Colloid and Interface Science, vol. 171, pp. 55–62, Apr. 1995. → pages 73, 76
[60] R. Pyter, G. Zografi, and P. Mukerjee, “Wetting of solids by surface-active agents: The effects of unequal adsorption to vapor-liquid and solid-liquid interfaces,” Journal of Colloid and Interface Science, vol. 89, pp. 144–153, Sept. 1982. → pages 73
[61] T. Mitsui, S. Nakamura, F. Harusawa, and Y. Machida, “Changes in the interfacial tension with temperature and their effects on the particle size and stability of emulsions,” Kolloid-Zeitschrift und Zeitschrift fr Polymere, vol. 250, pp. 227–230, Mar. 1972. → pages 73
[62] S. Phongikaroon, R. Hoffmaster, K. P. Judd, G. B. Smith, and R. A. Handler, “Effect of temperature on the surface tension of soluble and insoluble surfactants of hydrodynamical importance,” Journal of Chemical & Engineering Data, vol. 50, pp. 1602–1607, Sept. 2005. → pages 73, 80
[63] V. S. Vesselovsky and V. N. Pertzov, “Adhesion of air bubbles to the solid surface,” Zh. Fiz. Khim, vol. 8, pp. 245–259, 1936. → pages 75
[64] Hideo Nakae, Ryuichi Inui, Yosuke Hirata, and Hiroyuki Saito, “Effects of surface roughness on wettability,” Acta Materialia, vol. 46, pp. 2313–2318, Apr. 1998. → pages
[65] R. J. Good and M. Koo, “The effect of drop size on contact angle,” Journal of Colloid and Interface Science, vol. 71, pp. 283–292, Sept. 1979. → pages

Figure 2. Simulation of droplet separation by EWOD

Non-Linear Electrohydrodynamics in Microfluidic Devices

미세 유체 장치의 비선형 전기 유체 역학

by Jun ZengHewlett-Packard Laboratories, Hewlett-Packard Company, 1501 Page Mill Road, Palo Alto, CA 94304, USAInt. J. Mol. Sci.201112(3), 1633-1649; https://doi.org/10.3390/ijms12031633Received: 24 January 2011 / Revised: 10 February 2011 / Accepted: 24 February 2011 / Published: 3 March 2011

Abstract

Since the inception of microfluidics, the electric force has been exploited as one of the leading mechanisms for driving and controlling the movement of the operating fluid and the charged suspensions. Electric force has an intrinsic advantage in miniaturized devices. Because the electrodes are placed over a small distance, from sub-millimeter to a few microns, a very high electric field is easy to obtain. The electric force can be highly localized as its strength rapidly decays away from the peak. This makes the electric force an ideal candidate for precise spatial control. The geometry and placement of the electrodes can be used to design electric fields of varying distributions, which can be readily realized by Micro-Electro-Mechanical Systems (MEMS) fabrication methods. In this paper, we examine several electrically driven liquid handling operations. The emphasis is given to non-linear electrohydrodynamic effects. We discuss the theoretical treatment and related numerical methods. Modeling and simulations are used to unveil the associated electrohydrodynamic phenomena. The modeling based investigation is interwoven with examples of microfluidic devices to illustrate the applications. 

Keywords: dielectrophoresiselectrohydrodynamicselectrowettinglab-on-a-chipmicrofluidicsmodelingnumerical simulationreflective display

요약

미세 유체학이 시작된 이래로 전기력은 작동 유체와 충전 된 서스펜션의 움직임을 제어하고 제어하는 ​​주요 메커니즘 중 하나로 활용되어 왔습니다. 전기력은 소형 장치에서 본질적인 이점이 있습니다. 전극이 밀리미터 미만에서 수 미크론까지 작은 거리에 배치되기 때문에 매우 높은 전기장을 쉽게 얻을 수 있습니다. 

전기력은 강도가 피크에서 멀어지면서 빠르게 감소하기 때문에 고도로 국부화 될 수 있습니다. 이것은 전기력을 정밀한 공간 제어를 위한 이상적인 후보로 만듭니다.

전극의 기하학적 구조와 배치는 다양한 분포의 전기장을 설계하는 데 사용될 수 있으며, 이는 MEMS (Micro-Electro-Mechanical Systems) 제조 방법으로 쉽게 실현할 수 있습니다. 

이 논문에서 우리는 몇 가지 전기 구동 액체 처리 작업을 검토합니다. 비선형 전기 유체 역학적 효과에 중점을 둡니다. 이론적 처리 및 관련 수치 방법에 대해 논의합니다. 모델링과 시뮬레이션은 관련된 전기 유체 역학 현상을 밝히는 데 사용됩니다. 모델링 기반 조사는 응용 분야를 설명하기 위해 미세 유체 장치의 예와 결합됩니다. 

키워드 : 유전 영동 ; 전기 유체 역학 ; 전기 습윤 ; 랩 온어 칩 ; 미세 유체 ; 모델링 ; 수치 시뮬레이션 ; 반사 디스플레이

Droplet processing array Droplet based BioFlip
igure 1. Example of droplet-based digital microfluidics architecture. Above is an elevation view showing the layered structure of the chip. Below is a diagram illustrating the system (Adapted from [4]).
Figure 2. Simulation of droplet separation by EWOD
Figure 2. Simulation of droplet separation by EWOD. The top two figures illustrate the device configuration. Electric voltages are applied to all four electrodes embedded in the insulating material. The bottom left figure shows transient simulation solution. It illustrates the process of separating one droplet into two via EWOD. The bottom right figure shows the electric potential distribution inside the device. The color indicates the electric potential; the iso-potential surfaces are also drawn. The image shows the electric field is absent within the droplet body indicating the droplet is either conductive or highly polarizable.
Figure 4. Transient sequence of the Taylor cone formation
Figure 4. Transient sequence of the Taylor cone formation: simulation and experiment comparison. Experimental images are shown in the top row. Simulation results are shown in the bottom row. Their correspondence is indicated by the vertical alignment (Adapted from [4]).
Figure 6. Simulation of charge screening effect using a parallel-plate cell
Figure 6. Simulation of charge screening effect using a parallel-plate cell. Top-left image shows the electric current as function of time and driving voltage, top-right image shows the evolution of the species concentration as function of time and space, the bottom image shows the electric current readout after switching the applied voltage.
Figure 7. Transient simulation of electrohydrodynamic instability and the development of the cellular convective flow pattern.
Figure 7. Transient simulation of electrohydrodynamic instability and the development of the cellular convective flow pattern.
Figure 3. Simulation of dielectrophoresis driven axon migration
Figure 3. Simulation of dielectrophoresis driven axon migration. The set of small images on the left shows a transient simulation of single axon migration under an electric field generated by a pin electrode. The image on the right is a snapshot of a simulation where two axons are fused by dielectrophoresis using a pin electrode. Axons are outlined in white. Also shown are the iso-potential curves.

References

  1. Muller, RS. MEMS: Quo vadis in century XXI. Microelectron. Eng 200053(1–4), 47–54. [Google Scholar]
  2. Reyes, DR; Iossifidis, D; Auroux, PA; Manz, A. Micro total analysis systems. 1. Introduction, theory, and technology. Anal.Chem 200274, 2623–2636. [Google Scholar]
  3. Levy, U; Shamai, R. Tunable optofluidic devices. Microfluid. Nanofluid 20084, 97–105. [Google Scholar]
  4. Zeng, J; Korsmeyer, FT. Principles of droplet electrohydrodynamics for lab-on-a-chip. Lab Chip 20044, 265–277. [Google Scholar]
  5. Fair, RB. Digital microfluidics: Is a true lab-on-a-chip possible? Microfluid. Nanofluid 20073, 245–281. [Google Scholar]
  6. Pollack, MG; Fair, RB; Shenderov, AD. Electrowetting-based actuation of liquid droplets for microfluidic applications. Appl. Phys. Lett 200077(11), 1725–1726. [Google Scholar]
  7. Peykov, V; Quinn, A; Ralston, J. Electrowetting: A model for contact-angle saturation. Colloid Polym. Sci 2000278, 789–793. [Google Scholar]
  8. Verheijen, HJJ; Prins, MWJ. Reversible electrowetting and trapping of charge: Model and experiments. Langmuir 199915, 6616–6620. [Google Scholar]
  9. Mugele, F; Baret, J. Electrowetting: From basics to applications. J. Phys. Condens. Matter 200517, R705–R774. [Google Scholar]
  10. Quilliet, C; Berge, B. Electrowetting: A recent outbreak. Curr. Opin. Colloid Interface Sci 20016, 34–39. [Google Scholar]
  11. Probstein, RF. Physicochemical Hydrodynamics; Wiley: New York, NY, USA, 1994. [Google Scholar]
  12. Koo, J; Kleinstreuer, C. Liquid flow in microchannels: Experimental observations and computational analyses of microfluidics effects. J. Micromech. Microeng 200313, 568–579. [Google Scholar]
  13. Hu, G; Li, D. Multiscale phenomena in microfluidics and nanofluidics. Chem. Eng. Sci 200762, 3443–3454. [Google Scholar]
  14. Haus, HA; Melcher, JR. Electromagnetic Fields and Energy; Prentice-Hall: Englewood Cliffs, NJ, USA, 1989. [Google Scholar]
  15. Leal, LG. Laminar Flow and Convective Transport Processes: Scaling Principles and Asymptotic Analysis; Butterworth-Heinemann: Oxford, UK, 1992. [Google Scholar]
  16. Collins, RT; Harris, MT; Basaran, OA. Breakup of electrified jets. J. Fluid Mech 2007588, 75–129. [Google Scholar]
  17. Sista, R; Hua, Z; Thwar, P; Sudarsan, A; Srinivasan, V; Eckhardt, A; Pollack, M; Pamula, V. Development of a digital microfluidic platform for point of care testing. Lab Chip 20088, 2091–2104. [Google Scholar]
  18. Zeng, J. Modeling and simulation of electrified droplets and its application to computer-aided design of digital microfluidics. IEEE Trans. Comput. Aid. Des. Integr. Circ. Syst 200625(2), 224–233. [Google Scholar]
  19. Walker, SW; Bonito, A; Nochetto, RH. Mixed finite element method for electrowetting on dielectric with contact line pinning. Interface. Free Bound 201012, 85–119. [Google Scholar]
  20. Eck, C; Fontelos, M; Grün, G; Klingbeil, F; Vantzos, O. On a phase-field model for electrowetting. Interface. Free Bound 200911, 259–290. [Google Scholar]
  21. Gascoyne, PRC; Vykoukal, JV. Dielectrophoresis-based sample handling in general-purpose programmable diagnostic instruments. Proc. IEEE 200492(1), 22–42. [Google Scholar]
  22. Jones, TB; Gunji, M; Washizu, M. Dielectrophoretic liquid actuation and nanodroplet formation. J. Appl. Phys 200189(3), 1441–1448. [Google Scholar]
  23. Sretavan, D; Chang, W; Keller, C; Kliot, M. Microscale surgery on single axons. Neurosurgery 200557(4), 635–646. [Google Scholar]
  24. Pohl, HA; Crane, JS. Dielectrophoresis of cells. Biophys. J 197111, 711–727. [Google Scholar]
  25. Melcher, JR; Taylor, GI. Electrohydrodynamics: A review of the role of interfacial shear stresses. Annu. Rev. Fluid Mech 19691, 111–146. [Google Scholar]
  26. Saville, DA. Electrohydrodynamics: The taylor-melcher leaky-dielectric model. Annu. Rev. Fluid Mech 199729, 27–64. [Google Scholar]
  27. Schultz, GA; Corso, TN; Prosser, SJ; Zhang, S. A fully integrated monolithic microchip electrospray device for mass spectrometry. Anal. Chem 200072(17), 4058–4063. [Google Scholar]
  28. Killeen, K; Yin, H; Udiavar, S; Brennen, R; Juanitas, M; Poon, E; Sobek, D; van de Goor, T. Chip-MS: A polymeric microfluidic device with integrated mass-spectrometer interface. Micro Total Anal. Syst 2001, 331–332. [Google Scholar]
  29. Dukhin, SS. Electrokinetic phenomena of the second kind and their applications. Adv. Colloid Interface Sci 199135, 173–196. [Google Scholar]
  30. Wang, Y-C; Stevens, AL; Han, J. Million-fold preconcentration of proteins and peptides by nanofluidic filter. Anal. Chem 200577(14), 4293–4299. [Google Scholar]
  31. Kim, SJ; Wang, Y-C; Han, J. Nonlinear electrokinetic flow pattern near nanofluidic channel. Micro Total Anal. Syst 20061, 522–524. [Google Scholar]
  32. Comiskey, B; Albert, JD; Yoshizawa, H; Jacobson, J. An electrophoretic ink for all-printed reflective electronic displays. Nature 1998394(6690), 253–255. [Google Scholar]
  33. Beunis, F; Strubbe, F; Neyts, K; Bert, T; De Smet, H; Verschueren, A; Schlangen, L. P-39: Electric field compensation in electrophoretic ink display. In Proceedings of the Twenty-fifth International Display Research Conference—Eurodisplay 2005; Edinburgh, UK, 19–22 2005; pp. 344–345. [Google Scholar]
  34. Strubbe, F; Verschueren, ARM; Schlangen, LJM; Beunis, F; Neyts, K. Generation current of charged micelles in nonaqueous liquids: Measurements and simulations. J. Colloid Interface Sci 2006300, 396–403. [Google Scholar]
  35. Hsu, MF; Dufresne, ER; Weitz, DA. Charge stabilization in nonpolar solvents. Langmuir 200521, 4881–4887. [Google Scholar]
  36. Hayes, RA; Feenstra, BJ. Video-speed electronic paper based on electrowetting. Nature 2003425, 383–385. [Google Scholar]
  37. Chakrabarty, K; Su, F. Digital Microfluidic Biochips: Synthesis, Testing, and Reconfiguration Techniques; CRC Press: Boca Raton, FL, USA, 2006. [Google Scholar]
  38. Chakrabarty, K; Fair, RB; Zeng, J. Design tools for digital microfluidic biochips: Towards functional diversification and more than Moore. IEEE Trans.CAD Integr. Circ. Syst 201029(7), 1001–1017. [Google Scholar]
Deep 코팅 검증계산

The Coating Application Using the Excellent Flow Modeling Software FLOW-3D

우수한 플로우 모델링 소프트웨어 FLOW-3D를 이용한 코팅 적용 연구

FLOW-3D는 미국 Flow Science Inc.에 의해 개발된 고유한 계산 유체 동적 프로그램입니다. FORE-3D는 FORDR(장애물 표현의 단편 영역 볼륨) 유한 차이 체계를 기반으로 Navier-Stokes 전체 솔버를 가지고 있습니다.

실제 VOF(Volume of Fluid) 알고리즘은 FLOW-3D에 통합되어 신뢰할 수 있는 자유 표면 흐름 분석을 제공합니다. FLOW-3D에는 다양한 물리적 모델이 있습니다. 따라서 FLOW-3D는 잉크젯 또는 코팅 등 광범위한 산업 영역에 사용됩니다.

본 논문에서는 FLOW-3D의 특징과 동적 접촉선의 직접 연산, 코팅 적용 예제를 설명합니다.

확대한 구형 방울
확대한 구형 방울
표면 파를 수반하는 세류의 시트 모양 흐름/세류가 축퇴하여 액적을 형성하는 예
표면 파를 수반하는 세류의 시트 모양 흐름/세류가 축퇴하여 액적을 형성하는 예
Deep 코팅 검증계산
Deep 코팅 검증계산
롤 코팅 검증계산
롤 코팅 검증계산
2層 コー テ ィング計算 例/ゆ っ くりした ウェ ブ接 近
2層 コー テ ィング計算 例/ゆ っ くりした ウェ ブ接 近
カー テ ン塗 布 のエ ッジ近 傍 にお け る塗 液流 れ解 析
カー テ ン塗 布 のエ ッジ近 傍 にお け る塗 液流 れ解 析
Pinned contact line resulting in coffee ring deposits (a). Constant contact angle and mixed mode resulting in moderately more uniform deposits (b).

Inkjet Printability of Electronic Materials Important to the Manufactur Manufacture of Fully Printed O ully Printed OTFTs

Sooman Lim
Western Michigan University, sooman.lim@gmail.com

초록

본 연구에서는 OTFT(Printed Organic Thin Film Transistors) 제작에 중요한 재료의 잉크젯 인쇄성이 조사되었습니다. 잉크젯 인쇄 잉크의 분사 진화를 이해하기 위해 나노 구리 및 나노 입자 은 잉크로 시뮬레이션이 수행되었습니다. 나노 구리 잉크의 잉크젯 적합성을 예측하기 위해 온도 차이가 있는 Z와 Oh 수를 측정했습니다. FLOW-3D를 이용한 시뮬레이션 연구의 결과를 Dimatix 잉크젯 프린터를 사용하여 얻은 실험 결과와 비교했습니다.

반도체 잉크의 경우, 두 유기 반도체의 잉크젯 인쇄성 P2TDC17FT4(poly[(3,7-dipdecdecyltheno[3,2-b]theno[2′,3′:4,5]theno[2,3-diopneo] 티오페인-2,6-diopeo[2,6-diotyl)]입니다.HT(poly-3 hexylthiophene)를 비교하여 낙하 속도, 낙하 볼륨 및 점화 전압 간의 관계를 확인하고, 낙하 간격 및 기판 온도가 인쇄 품질에 미치는 영향을 확인했습니다.

이러한 연구를 통해 인쇄 가능성과 인쇄 품질은 잉크젯으로 인쇄된 상단 게이트 OTFT를 완벽하게 구현하기에 충분했습니다. 주변 조건에서 인쇄되는 P2TDC17FT4의 성능은 저비용 완전 인쇄 OTFT의 실현에 중요한 영향을 미칩니다.

후처리 연구로 은색 잉크의 유망한 대체품인 나노 구리 잉크를 IPL(Incensive Pulsed Light)로 소결시키는 것이 연구되었습니다. 잉크 필름 두께와 소결 시 필요한 에너지 사이의 관계가 확인되었습니다. 잉크 필름 두께와 관련하여 유리와 PET에 소결하는데 필요한 에너지 수준을 비교한 결과, 이 잉크의 처리 요구 사항에 대한 기판의 열적 기여도가 밝혀졌습니다. 이 조사 결과는 자재 특성 요구 사항에 대한 현재의 이해와 완전히 잉크젯으로 인쇄된 OTFT를 달성하기 위한 과제를 진전시킵니다.

Schematic design showing the principles of operation of a continuous inkjet (CIJ) printer.
Schematic design showing the principles of operation of a continuous inkjet (CIJ) printer.
Illustration of the piezo movement under an applied voltage.
Illustration of the piezo movement under an applied voltage.
Construction of a traditional piezoelectric squeeze type print head.
Construction of a traditional piezoelectric squeeze type print head.
Pinned contact line resulting in coffee ring deposits (a). Constant contact angle and mixed mode resulting in moderately more uniform deposits (b).
Pinned contact line resulting in coffee ring deposits (a). Constant contact angle and mixed mode resulting in moderately more uniform deposits (b).
Marangoni effect, where Tc is the CT line temperature, Ta is the drop apex temperature,  and a is the drop apex surface tension
Marangoni effect, where Tc is the CT line temperature, Ta is the drop apex temperature, and a is the drop apex surface tension
Comparison of drop evolution and drop ejection pictures droplet obtained experimentally and using CFD software for the nano copper ink
Comparison of drop evolution and drop ejection pictures droplet obtained experimentally and using CFD software for the nano copper ink

접촉선의 고정(Contact Line Pinning)

접촉선의 고정(Contact Line Pinning)

증발하는 빗방울에서 남은 잔류의 물은 새로 씻은 자동차에서 좋지 못할 수 있습니다. 그러나, 동일한 증발 공정은, 예를 들어, 드롭 잔류 물이 인쇄 된 이미지 또는 텍스트의 일부가되는 잉크젯 인쇄에서 유리할 수있다. 그러나 동일한 증발 과정이 어떤 경우엔 도움이 될 수 있습니다 예를 들면, 잉크 찌꺼기가 인쇄 된 이미지나 텍스트의 일부가 되는 잉크젯 인쇄가 그렇습니다.

액체 방울의 증발로 인한 잔류의 물이 예상치 못한 방식으로 나타날 수 있습니다. 커피 링 얼룩이 잘 알려진 예이며, 커피의 잔류의 물이 물방울의 바깥 쪽 가장자리에 모여 얇은 원형 링 얼룩이 남습니다. 이 현상은 흥미로운 유체역학적인 과정의 결과입니다. 커피 링 얼룩이 형성 되려면 액체가 증착 된 고체 표면에 고정 된 접촉선이 있어야합니다. 고정 된 접촉선은 액체 방울이 고체 기판과 교차하는 액체 방울의 외부의 가장자리가 방울이 증발함에 따라 정지 상태를 유지함을 의미합니다. 증발은 기판의 열에 의해 발생하며 방울의 얇은 외부의 가장자리에서 가장 크게 생깁니다. 표면 장력은 액체가 증발하면서 손실 된 액체를 대체하기 위해 가장자리를 향해 발생하게 됩니다. 이는 결국 더 많은 용질을 가장자리로 운반하며 모든 액체가 증발 한 후, 결과적으로 커피 링 얼룩을 형성하게하는 더 높은 농도의 용질 잔류 물을 생성합니다.

모델링 접근법

FLOW-3D v12.0의 최신 업데이트로 인해 ‘접촉선의 고정’ 모델이 개발되었으며, 소프트웨어의 기능이 표면 장력 중심의 애플리케이션으로도 광범위하게 확장되었습니다. 표면 접촉의 고정 및 비고정 특성은 잉크젯 인쇄, 코팅 및 스프레이 냉각에서 중요한 역할을 합니다. 습윤 특성에 대한 표면 공법은 미세 유체 장치에서 액체 샘플의 이동을 제어하는 ​​데 사용될 수 있습니다. 모델의 주요 특징은 방울의 가장자리를 고정 위치에 고정하는 수단을 제공하는 것입니다. 형상 구성 요소 및 하위 구성 요소중에 표면에 ‘고정’ 속성을 지정할 수 있습니다. 유체의 접촉선은 처음 표면과 접촉하는 곳에 고정됩니다. 전방 속도를 0으로 유지하면 고정이 적용됩니다. 유체는 접촉선과 표면을 따라 이동하는 것이 아니라 롤오버하여 접촉점을 지나야만 이동할 수 있습니다.

커피 링 얼룩 검증

그림 1은 평평한 수평 표면에 놓인 원형 물방울의 결과를 보여줍니다. 표면은 30 ℃의 일정한 온도로 유지됩니다. 초기 유체 온도는 20 ℃이고 주변 공극의 온도는 일정한 20 ℃입니다. 유체는 밀도 0.967 g/cm3, 점도 0.02022 poise, 비열 1.645e+07 cm2/s/K, 열전도도 1.2964e+4 g*cm/s3/K, 표면 장력 계수 33.15 g/cm2의 일반적인 잉크를 나타냅니다.

그림 1. 고정 된 접촉선을 사용하여 건조 공정 중의 물방울 모양의 변화.

액적 표면의 초기 곡률 반경은 7.5e-03 cm이고, 차지하는 공간은 반경 4.5e-03 cm의 원이며, 겉보기의 초기 접촉각은 37.87 도입니다. 그림 1-a를 참조하시기 바랍니다. 지정된 정적 접촉각은 0 도입니다.

정압에 의한 상변화 모델이 활성화됩니다. 공극 내의 증기 분압은 0이고 상변화 수용 계수는 Rsize = 0.01 입니다.

잉크가 건조될 때 기판 상에 고체가 잔류하는 물이 형성되는 것을 포착하기 위해 잔류 물 모델도 켜집니다. 유체에 용해 된 안료의 농도는 초기 농도 0.01 g/cm3 이고 최대 농도 rmax = 1.1625 g/cm3 에서 운반이 가능한 스칼라로 표시됩니다. 용해 된 안료는 질량 평균을 기준으로 안료의 단위질량당 0.05 poise의 속도로 유체의 순 점도를 향상시킵니다.

이 공정은 3.0 도의 방위 방향으로 하나의 셀에 걸쳐있는 축 대칭 원통형 메쉬로 모델링됩니다. (x 간격 = 6e-05 cm, z 간격 = 4e-05 cm.)

그림 1은 유체가 증발함에 따라 접촉선이 고정 된 상태를 유지하고 있음을 보여줍니다. 0 도의 정적 접촉각 조건은 액적의 중심을 향한 압력 구배를 가져오고, 이는 접촉선 방향으로의 유동을 생성합니다. 용해 된 안료의 농도는 증발로 인해 자유 표면 근처에서 증가하며, 흐름을 따라 농도는 접촉선을 향해 더욱 재분배합니다. (그림 2). 액체가 계속 증발함에 따라, 남아있는 액체의 안료 농도는 증가합니다. 농도가 최대 rmax에 도달하면, 과잉된 안료는 고체가 잔류하는 물로 전환됩니다.

그림2. g / cm3 단위의 안료 농도 및 t = 2.0ms에서의 흐름 패턴. 흐름은 고정 된 접촉선을 향하여 안료 농도가 증가합니다.

접촉선 근처의 유체가 먼저 건조되어 고체가 잔류하는 물이 남습니다. 해당 영역의 유체에 안료 농도가 높기 때문에 고체가 잔류하는 물의 특징인 ‘커피 링’ 패턴이 기판 표면에 생성됩니다. (그림 3 및 4). 안료의 총 질량(용해 + 건조 잔류 물)은 초기 질량의 0.025 % 이내로 보존됩니다.

그림 3. 모든 유체가 증발 된 후 기판 표면에 건조된 잔류 물의 분포 (단위 : g / cm3) .
가장 높은 농도는 고정 된 접촉선의 위치에 있으며, 이는 ‘커피 링’ 효과를 만들어냅니다.
그림 4. 유체가 완전히 증발 한 후 초기 액적의 반경을 따라 건조된 잔류 물의 예상 분포.

물방울 벽의 검증

그림5. 수직 벽에 고정 된 물방울의 변형 : t = 0 ms (파란색), t = 4e-02 ms (연한 파랑) t = 0.2 ms (빨간색).
해당 이미지는 “Effects of microscale topography”, Y.V.Kalinin, V.Berejnov and R. E. Thorne, Langmuir 25, 5391-5397. (2009). 에서의 이미지입니다.

접촉선 고정 응용의 두 번째 예는 수직의 벽에 고정 된 한 방울의 액체 알루미늄의 거동입니다. 유체 밀도는 2.7 g / cm3, 표면 장력 계수 200 g / cm2 및 점도 0.27 poise입니다. 정적 접촉각은 0 도입니다.

초기의 겉보기의 접촉각이 90도가 되도록 반경 0.5cm의 물방울을 수직 벽에 놓습니다 (그림 5). 7e+06 cm/s2의 중력 크기는 표면 장력의 복원 작용을 없애고 액적이 눈에 띄도록 변형시키기 위하여 인위적으로 향상되었습니다. 결과들은 비슷한 크기의 물방울에 대한 실험 결과와의 질적 비교를 포함하여 그림 5에서 보여줍니다.

요약

FLOW-3D의 접촉선 고정 모델은 표면 장력 및 벽의 접착 기능을 확장하여 표면 공법에서 복잡한 상호 작용을 모델링합니다. 접촉선 고정이 실제로 응용되는 분야에 관하여 더 많은 예시와 추가적인 참조를 찾으신다면 여기에서 찾을 수 있습니다.

접촉선의 이해(Contact Line Insights)

접촉선의 이해(Contact Line Insights)

FLOW-3D는 코팅 성능 향상에 관심이있는 엔지니어에게 이상적인 수치 모델링 기능을 많이 갖추고 있습니다. 전산 시뮬레이션은 코팅 흐름에 영향을 미치는 여러 물리적 과정의 상대적 중요성과 효과를 연구 할 수있는 훌륭한 방법입니다. 물리적인 테스트에서 항상 프로세스를 분리하거나 해당 프로세스의 크기를 임의로 조정할 수있는 것은 아닙니다. 여기에서는 리 볼렛 형성(rivulet formation), 핑거링(fingering), 증발, 거친 표면에서의 접촉선 이동 및 유체 흡수와  관련하여 정적 및 동적 접촉각에 대하여 FLOW-3D의 처리에 대해 설명합니다.

 

정적 및 동적 접촉각(Static and Dynamic Contact Angles)

FLOW-3D는 정적 접촉각의 함수로 동적 접촉각을 정확하게 계산하고 입력으로 설정하며 자유 표면 인터페이스에서 작용하는 관련된 힘을 정확하게 계산하여 유체의 소수성을 캡처 할 수 있습니다. 아래 시뮬레이션은 물방울이 경사를 따라 내려갈 때 정적 접촉각이 동적 접촉각에 미치는 영향을 보여줍니다.

 

흡수(Absorption)

종이 기판에 액 적의 충격 및 흡수는 전산 유체 역학 소프트웨어를 사용하여 연구 할 수 있습니다. 여기서 FLOW-3D는 섬유층에서 물방울 충돌을 시뮬레이션하는데 사용되며 표면 장력, 접촉각 및 점도와 관련된 유체 전면의 전파를 살펴 봅니다.

 

 

아래의 FLOW-3D 시뮬레이션에서, 낙하는 직경이 40 미크론이며 초기 하향 속도는 300 cm / s입니다. 기재는 종이이고, 기공률이 30 % 인 20 미크론 두께입니다.

 

 

액체 필름의 핑거링(Fingering in Liquid Films)

FLOW-3D에서 동적 접촉선은 동적 접촉각이나 접촉선의 위치를 ​​지정할 필요없이 직접 모델링됩니다. 이는 소량의 유체에서 유체에 영향을 미치는 모든 동적 힘을 포함하는 수치 모델을 사용하여 수행됩니다. 정적 접촉각은 액체-고체 접착력을 특성화 하는데 사용됩니다.

액체 시트의 핑거링. 왼쪽은 0 °, 오른쪽은 70 °

여기서, 이러한 접근법의 힘의 적용은 경사 표면 아래로 흐르는 액체 필름에서 관찰 된 핑거링에 의해 제공됩니다. 실험적 관찰에 따르면 두 가지 뚜렷한 핑거링 패턴이 발생합니다. 첫 번째 패턴은 작은 정적 접촉각(즉, 습윤 조건)이며 상하한이 모두 하향으로 움직이는 쐐기형 핑거를 나타냅니다. 두 번째 패턴은 큰 정적 접촉각(즉, 습윤 조건이 열악함)이며 가장 균일한 폭을 가진 긴 핑거이고 가장 큰 한계점은 하향으로 움직이지 않는 것이 특징입니다.

 

 

증발 효과(Evaporative Effects)

퇴적(Deposit)

분산 된 고체 물질을 함유하는 액 적은 고체 표면에서 건조 될 때, 함유하고 있는 고체 물질을 침전물로서 남깁니다. 이 침전물의 형상이 많은 인쇄 공정, 청소 및 코팅 공정에 중요한 영향을 미칩니다. 한 종류의 퇴적물의 전형적인 예는 위의 이미지와 같이 엎질러 진 커피 패치의 둘레를 따라 링 얼룩이 형성되는 “커피 링” 문제입니다. 이 유형의 링 침전물은 액체의 증발로 인한 표면 장력 구동 흐름의 결과로, 특히 낙하 둘레에서 발생합니다.

 

건조(Drying)

FLOW-3D의 증발 잔류 액체 모델은 건조 후 톨루엔으로 형성된 잔류된 물의 3D형상을 시뮬레이션합니다. (30 배 확대)

건조는 코팅 공정의 중요한 부분입니다. 하지만 건조의 결함으로 잘 도포 된 코팅을 완전히 취소 할 수도 있습니다. 건조 중에 온도 및 용질 구배는 밀도 및 표면 장력 구배로 인해 코팅 내 유동을 유도 할 수 있으며, 이는 코팅 품질을 잠재적으로 파괴 할 수 있습니다. FLOW-3D의 증발 잔류 물 모델을 사용하면 건조로 인한 흐름을 시뮬레이션하고 값 비싼 물리적 실험에 소요되는 시간을 줄일 수 있습니다.

 

모델링 링 형성(Modeling Ring Formation)

증발에 의해 접촉 라인에서 생성 된 흐름 시뮬레이션

윗쪽 그림에서 FLOW-3D는 증발이 가장 큰 접촉선에서의 증착으로 인해 에지 피닝(edge pinning)이 발생함을 보여줍니다. 증발은 증발로 인한 열 손실로 인해 액체를 냉각시킵니다 (색상은 온도를 나타냄). 동시에 고체 표면은 전도에 의해 액체를 가열합니다. 접촉선 주변에서 증발이 가장 커서, 액체가 접촉선을 향해 흘러 정적 조건을 재설정합니다. 최종 결과는 액체가 완전히 증발하는 액체 가장자리에 현탁 된 고체의 증착입니다.

 

 

참고
[1] Deegan, R., Bakajin, O., Dupont, T. et al. Capillary flow as the cause of ring stains from dried liquid drops, Nature 389, 827–829 (1997).

 

Computational Analysis of Drop Formation and Detachment

Computational Analysis of Drop Formation and Detachment

Introduction and Problem Statement

신속, 반복, 작은 물방울의 생성 및 증착, 작은 형상의 프린팅 또는 패터닝 (예 : l = 10-3-1 mm), 스프레이로  균일한 두께의 박막 형성은 다양한 산업에 매우 중요합니다(1-5). 액체 이동과 액적 형성 / 증착 공정은 복잡한 자유 표면 흐름, 자연적인 모세관운동 형성, thinning, pinch-off를 수반한다 (1-5). 단순한 뉴턴 및 비탄성 유체에 대해 액적 생성 및 액적 이동을 분석하기위한 실험적, 이론적 및 1 차원 시뮬레이션 연구가 진행되었지만 프린팅 또는 패터닝에 대한 기계론적인 이해는 여전히 과제로 남아 있습니다. 현재의 계산에 대한 주된 목표는 뉴턴 유체의 pinch-off에 대한 기계론적 이해를 얻기 위해 FLOW-3D에 내장된 VOF(volume-of-fluid) 접근법으로 시험하는 것입니다. 전산해석은 모세관, 관성, 점성 응력의 복잡한 상호 작용을 포착하여 자기유사 모세관의 thinning and pinch-off를 결정합니다. 뉴턴 유체의 물방울 형성 ​​및 분리현상은  전산해석으로부터 얻어진 자기유사 모세관현상 이론, 보편적인 축소화 기법인 1D 시뮬레이션 (1-7)과 실험 (1, 2, 8-12)을 이용하여 설명될 수 있음을 보여준다. 이러한 우리가 진행한 원형흐름 시뮬레이션은 유한한 시간의 비선형 역학, 위성 낙하현상, 복잡한 형상의 프린팅과 같이 어려운 전산해석의 기반이 될 것 입니다.

방울 형성의 전산 분석
그림 1 : FLOW-3D를 사용하여 시뮬레이션 한 저점도 유체의 드롭 형성 및 분리에 대한 전산해석 : (a) 5개의 저점도 유체에 대한 물방울의 necking에 대한 반경이 시간변화에 따라 표시됩니다. 물방울 necking의 반지름이 오른쪽에서 왼쪽으로 시간에 따른 전개를 보여줍니다. 마찬가지로 스냅 샷은 necking의 반경이 오른쪽에서 왼쪽으로 줄어듭니다. 속도의 크기 (단위 : cm/s) 와 화살표의 방향에 대한 컬러 맵을 사용하면 변형장을 결정할 수 있으며 Fluid 5 (표 1 참조)의 경우에는 순식간에 신장이됩니다. 이미지 II에 캡처 된 pinch-off 하기 전에 형성된 원추형 necking은 실험을 통해 얻은 necking 모양과 유사합니다.

Modeling Approach and Parameter Space

표면 장력 및 중력 모델을 적용한 FLOW-3D 에서 균일한 메쉬 크기를 사용하여 노즐에서 드롭 형성 및 분리에 대한 시뮬레이션을 수행하였습니다. 유한 체적의 유체를 떨어뜨리거나 분리하는 일은 물방울의 성장과 드롭, 노즐에 연결되는 모세관 현상, 관성, 점도 및 중력에 대한 상호 작용을 수반합니다. 시뮬레이션에서 스테인레스 강 노즐 ( {{D} _ {0}} = 2 {{R} _ {0}} = 1.7 \, \ text {mm}) 에서 유한 체적의 뉴턴 유체가 발생합니다. 표면 장력이 중력을 겪으면 새로 형성된 액적 분리가 발생합니다 (mg> 2 \ pi \ sigma {{R} _ {0}}). 시뮬레이션은 유체점도의 영향을 설명하기 위해 두 그룹으로 나누어져 있습니다: 저점도 유체 (글리세롤 함량이 40 % 미만인 물과 글리세롤/물 혼합물) 및 점도가 높은 유체 (예 : 글리세롤과 글리세롤/물 혼합물 점도 > 100x 물 점도). 두 그룹의 유체 특성은 각각 표 1과 2에 나와 있습니다.

계산 분석 드롭 형성 저점도

그림 2 : FLOW-3D를 사용하여 시뮬레이션 한 저점도 유체의 드롭형성 및 분리에 대한 전산 해석 : 반경 플롯에서 4개의 고점도 뉴톤유체에 대해 necking 반경을 시간변화에 따라 표시합니다. 낙하 분리 중 모세관 현상이 스냅 샷으로 표시됩니다. 컬러 맵은 Fluid 8의 속도 크기 (단위 : cm/s)의 변화를 포착합니다 (표2 참조). 화살표는 성장하는 물방울과 얇아지는 물방울내에서 흐름방향을 나타냅니다. FLOW-3D 시뮬레이션으로 얻은 necking 모양은 고점도의 뉴턴유체에 대한 특징인 원통형 유체요소로 이어집니다.

 

<표 1 : FLOW-3D를 사용하여 시뮬레이션 된 저점도 유체의 특성>
Fluid PropertyFluid 1Fluid 2Fluid 3Fluid 4Fluid 5
Viscosity [Pa · s]0.050.020.010.00750.005
Surface Tension  [mN / m]6868686868
Density [g / cm 3 ]11111
Ohnesorge Number0.210.080.040.030.021
 저점도 유체 (표 1의 유체 2) 가 노즐에서 떨어지는 것을 시뮬레이션 합니다. 색상변수는 속도크기 (단위 : cm / s)이며 속도벡터가 표시됩니다.

 

<표 2 : FLOW-3D를 사용하여 시뮬레이션 된 고점도 유체의 특성>
Fluid PropertyFluid 6Fluid 7Fluid 8Fluid 9
Viscosity [Pa · s]1.50.80.50.25
Surface Tension  [mN / m ]68686868
Density [g / cm 3 ]1111
Ohnesorge Number6.243.332.081.04

고점도 유체 (표 2의 유체 8) 가 노즐에서 떨어지는 것을 시뮬레이션 합니다. 색상변수는 속도크기 (단위 : cm / s) 이며 속도 벡터가 표시됩니다.

Discussion of the Simulation Results

드롭 형성 및 분리는 표1과 표2에 열거 된 유체에 대해 FLOW-3D 를 사용하여 시뮬레이션 하였고, 시간 경과에 따른 necking 모양, 반경을 분석하였습니다. 물방울의 necking 모양과 저점도에서의 necking에 대한 역학(그림 1 참조)은 실험, 흐름 이론, 1D 시뮬레이션, 자기유사 관성에 대한 모세현상의 특성을 나타냅니다 (1, 2, 6, 7, 13) :

(1)  \ displaystyle \ frac {{R (t)}} {{{{R} _ {0}}}} \ approx 0.8 R {{{{왼쪽} {R} {0} 3}}} 오른쪽}) ^ {{{{frac {1} {3}}} {{왼쪽 {{{{왼쪽}}} {2} {3}}}}

여기서 R (t)가  necking의 순간 반경이고, R0는 노즐의 외부반경이며,  \ displaystyle \ sigma 는 표면 장력,  \ displaystyle \ rho 는 유체의 밀도 tC 는 pinch-off 시간이다. 마찬가지로, 이러한 더 높은 점도의 뉴턴유체에 대한 반경 변화데이터는 시간에 따른 반경의 감소를 나타내는 것이며,  Papageorgiou’s visco-capillary scaling (8, 9)은 아래의 식으로 표현된다.

(2)  \ {0 \} {} {} {} {} {} {} {} {} {} {} {} {} { } ({{t} _ {p}} - t)

모세관 속도(표면 장력과 점도의 비)의 측정 값은 McKinley와 Tripathi (8)에 의해 Capillary Break-Up Extensional Rheometer (CaBER)라고 불리는 상업적으로 이용 가능한 장비를 사용하여 얻은 값과 모세관 속도는 공칭 표면 장력과 점도를 사용하여 계산됩니다.

FLOW-3D 는 물방울의 necking부분을 속도 벡터로 시각화하여 유체의 흐름을 나타낼 수 있습니다. 또한, 이는 그림 1과 같이 전단, 확장을 겪은 후 얇아지는 물방울이 흐르는 과정의 순간을 결정할 수 있는 가능성을 줍니다. 추가로, 낮은 점도의 뉴턴유체는 높은 점도의 뉴턴 유체에 비해 질적으로 다른 거동을 보여준다(그림 2참조). 낮은 점도의 뉴턴 유체에 대한 necking 프로파일은 이론(6,13)에 따라 자기 유사성이 됩니다.

Conclusions, Outlook and Ongoing work

우리의 예비결과는 FLOW-3D 기반의 전산해석이 액적 형성과 탈착의 기초가 되는 프로토타입의 자유 표면흐름을 시뮬레이션하는데 사용될 수 있음을 보여줍니다 . 시뮬레이션된 반경변화 프로파일이 실험적으로 관찰된 높은 유체 및 이론적으로 예측된 유체인 스케일링 법칙 및 pinch-off dynamics과 일치하는 것을 발견하였습니다.

자주 사용되는 1D 또는 2D 모델과 달리 FLOW-3D 는 기본 응력 및 확장 유동장 (균일도 및 크기)의 강도와 얇은 액체 필라멘트 내 흐름에 대한 시각화를 나타낼 수 있습니다(그림1과 2 참조). 확장 유동장과 연관된 흐름 방향 속도 구배는 모세관현상이 나타나는 물방울의 얇은 부분 내에서 발생합니다. 유동학적으로 복잡한 유체에서 non Newtonian shear 및 신장, 점도뿐만 아니라 그외의 탄성 응력이 nonlinear pinch-off dynamics을 급격하게 변화시킵니다(2, 10-12). 우리는 현재 점탄성과 non-Newtonian 유동학을 사용하여 FLow-3D에 복합 유체의 처리 성능평가를 위한 강력한 연산 프로토콜을 개발하고 있습니다.

References

  1. J. Eggers, Nonlinear dynamics and breakup of free-surface flows. Rev. Mod. Phys. 69, 865-929 (1997).
  2. G. H. McKinley, Visco-elasto-capillary thinning and break-up of complex fluids. Rheology Reviews, 1-48 (2005).
  3. B. Derby, Inkjet Printing of Functional and Structural Materials: Fluid Property Requirements, Feature Stability, and Resolution. Annual Review of Materials Research 40, 395-414 (2010).
  4. O. A. Basaran, H. Gao, P. P. Bhat, Nonstandard Inkjets. Annual Review of Fluid Mechanics 45, 85-113 (2013).
  5. S. Kumar, Liquid Transfer in Printing Processes: Liquid Bridges with Moving Contact Lines. Annual Review of Fluid Mechanics 47, 67-94 (2014).
  6. R. F. Day, E. J. Hinch, J. R. Lister, Self-similar capillary pinchoff of an inviscid fluid. Phys. Rev. Lett. 80, 704-707 (1998).
  7. J. Eggers, M. A. Fontelos, Singularities: Formation, Structure, and Propagation. (Cambridge University Press, Cambridge, UK, 2015), vol. 53.
  8. G. H. McKinley, A. Tripathi, How to extract the Newtonian viscosity from capillary breakup measurements in a filament rheometer. J. Rheol. 44, 653-670 (2000).
  9. D. T. Papageorgiou, On the breakup of viscous liquid threads. Phys. Fluids 7, 1529-1544 (1995).
  10. J. Dinic, L. N. Jimenez, V. Sharma, Pinch-off dynamics and dripping-onto-substrate (DoS) rheometry of complex fluids. Lab on a Chip 17, 460-473 (2017).
  11. J. Dinic, Y. Zhang, L. N. Jimenez, V. Sharma, Extensional Relaxation Times of Dilute, Aqueous Polymer Solutions. ACS Macro Letters 4, 804-808 (2015).
  12. V. Sharma et al., The rheology of aqueous solutions of Ethyl Hydroxy-Ethyl Cellulose (EHEC) and its hydrophobically modified Analogue (hmEHEC): Extensional flow response in capillary break-up, jetting (ROJER) and in a cross-slot extensional rheometer. Soft Matter 11, 3251-3270 (2015).
  13. J. R. Castrejón-Pita et al., Plethora of transitions during breakup of liquid filaments. Proc. Natl. Acad. Sci. U.S.A. 112, 4582-4587 (2015).

Simulating the Wetting and Drying of Shallow Flows [얕은 흐름의 습윤 및 건조 시뮬레이션]

Introduction
Shallow flows, characterized by having a thickness much smaller than their lateral extent, can often be modeled by a depth-averaged (shallow-water or 2.5 dimensional) approximation.
Average fluid velocities are computed in the layer and the top fluid surface is free to move, which leads to a changing fluid-layer thickness. The advantages of this approach are its speed
and simplicity over full three-dimensional simulations.
One complication, however, is how to efficiently account for dynamic contact-line effects at lateral boundaries of the fluid. These boundaries are free to move over the underlying solid
surface. Furthermore, the fluid contact angle at these boundaries depends on the local dynamic flow conditions.
In this paper we present a new shallow-flow computational method based on the Volume-of-Fluid (VOF) technique, which conserves fluid mass, while allowing for general wetting and
drying behavior. Non-uniform surface tension and fluid-substrate interactions, defined by a static contact angle, are included in the model. No special prescriptions are needed to locate
contact line locations or define dynamic contact angles.

Microfluidics Bibliography

Microfluidics Bibliography

다음은 Microfluidics Bibliography의 기술 문서 모음입니다.
이 모든 논문은 FLOW-3D  결과를 특징으로  합니다. 미세 유체 공정 및 장치 를 성공적으로 시뮬레이션하기 위해 FLOW-3D 를 사용 하는 방법에 대해 자세히 알아보십시오  .

2024년 11월 20일 Update

109-24 Dileep Karnam, Yu-Lung Lo, Chia-Hua Yang, Spray mist-assisted drilling of through silicon vias (TSV) using nanosecond laser: Influence of CNT nanofluid, Journal of Materials Research and Technology, 31; pp. 679-688, 2024. doi.org/10.1016/j.jmrt.2024.06.109

22-24   Bin-Jie Lai, Li-Tao Zhu, Zhe Chen, Bo Ouyang, Zheng-Hong Luo, Review on blood flow dynamics in lab-on-a-chip systems: an engineering perspective, Chem & Bio Engineering, 1.1; pp. 26-43, 2024. doi.org/10.1021/cbe.3c00014

196-23 Daicong Zhang, Chunhui Jing, Wei Guo, Yuan Xiao, Jun Luo, Lehua Qi, Microchannels formed using metal microdroplets, Micromachines, 14.10; 1922, 2023. doi.org/10.3390/mi14101922

121-23 Feng Lin Ng, Zhanhong Cen, Yi-Chin Toh, Lay Poh Tan, A 3D-printed micro-perfused culture device with embedded 3D fibrous scaffold for enhanced biomimicry, International Journal of Bioprinting, 2023. doi.org/10.36922/ijb.0226

104-23 Cristina González-Fernández, Jenifer Gómez-Pastora, Eugenio Bringas, Inmaculada Ortiz, Computer-aided design of magnetophoretic microfluidic systems for enhanced recovery of target products, 33rd European Symposium on Computer-Aided Engineering (ESCAPE), 2023.

64-23   Tihomir Tjankov, Dimitar Trifonov, Conceptual design and 3D modeling of a microfluidic device for liver cells investigation, Industry 4.0, 8.2; pp. 39-41, 2023.

34-23   Chao Kang, Ikki Ikeda, Motoki Sakaguchi, Recoil and solidification of a paraffin droplet impacted on a metal substrate: Numerical study and experimental verification, Journal of Fluids and Structures, 118; 103839, 2023. doi.org/10.1016/j.jfluidstructs.2023.103839

64-22   Babatunde Aramide, Computational modelling of electrohydrodynamic jetting (Taylor cone formation, dripping & jet evolution): Case study of electrospinning, Thesis, University College London, 2022.

42-22   Islam Hassan, P. Ravi Selvaganapathy, Microfluidic printheads for highly switchable multimaterial 3D printing of soft materials, Advanced Materials Technologies, 2101709, 2022. doi.org/10.1002/admt.202101709

138-21   Enver Guler, Mine Eti, Aydin Cihanoglu, Esra Altiok, Kadriye Ozlem Hamaloglu, Burcu Gokcal, Ali Tuncel, Nalan Kabay, Ion exchange membranes with enhanced antifouling properties to produce energy from renewable sources, Proceedings of the 6th International Symposium on Green and Smart Technologies for a Sustainable Society, Santander, Cantabria, Spain, December 9-10, 2021.

45-21   Navid Tonekaboni, Mahdi Feizbahr, Nima Tonekaboni, Guang-Jun Jiang, Hong-Xia Chen, Optimization of solar CCHP systems with collector enhanced by porous media and nanofluid, Mathematical Problems in Engineering, 2021; 9984940, 2021. doi.org/10.1155/2021/9984840

40-21   B. Hayes, G.L. Whiting, R. MacCurdy, Modeling of contactless bubble–bubble interactions in microchannels with integrated inertial pumps, Physics of Fluids, 33.4; 042002, 2021. doi.org/10.1063/5.0041924

Below is a collection of technical papers in our Microfluidics Bibliography. All of these papers feature FLOW-3D results. Learn more about how FLOW-3D can be used to successfully simulate microfluidic processes and devices.

14-21   Jian-Chiun Liou, Chih-Wei Peng, Philippe Basset, Zhen-Xi Chen, DNA printing integrated multiplexer driver microelectronic mechanical system head (IDMH) and microfluidic flow estimation, Micromachines, 12.1; 25, 2021. doi.org/10.3390/mi12010025

08-20   Li Yong-Qiang, Dong Jun-Yan and Rui Wei, Numerical simulation for capillary driven flow in capsule-type vane tank with clearances under microgravity, Microgravity Science and Technology, 2020. doi.org/10.1007/s12217-019-09773-z

89-19   Tim Dreckmann, Julien Boeuf, Imke-Sonja Ludwig, Jorg Lumkemann, and Jorg Huwyler, Low volume aseptic filling: impact of pump systems on shear stress, European Journal of Pharmeceutics and Biopharmeceutics, in press, 2019. doi:10.1016/j.ejpb.2019.12.006

88-19   V. Amiri Roodan, J. Gomez-Pastora, C. Gonzalez-Fernandez, I.H. Karampelas, E. Bringas, E.P. Furlani, and I. Ortiz, CFD analysis of the generation and manipulation of ferrofluid droplets, TechConnect Briefs, pp. 182-185, 2019. TechConnect World Innovation Conference & Expo, Boston, Massachussetts, USA, June 17-19, 2019.

55-19     Julio Aleman, Sunil K. George, Samuel Herberg, Mahesh Devarasetty, Christopher D. Porada, Aleksander Skardal, and Graça Almeida‐Porada, Deconstructed microfluidic bone marrow on‐a‐chip to study normal and malignant hemopoietic cell–niche interactions, Small, 2019. doi: 10.1002/smll.201902971

37-19     Feng Lin Ng, Miniaturized 3D fibrous scaffold on stereolithography-printed microfluidic perfusion culture, Doctoral Thesis, Nanyang Technological University, Singapore, 2019.

32-19     Jenifer Gómez-Pastora, Ioannis H. Karampelas, Eugenio Bringas, Edward P. Furlani, and Inmaculada Ortiz, Numerical analysis of bead magnetophoresis from flowing blood in a continuous-flow microchannel: Implications to the bead-fluid interactions, Nature: Scientific Reports, Vol. 9, No. 7265, 2019. doi: 10.1038/s41598-019-43827-x

01-19  Jelena Dinic and Vivek Sharma, Computational analysis of self-similar capillary-driven thinning and pinch-off dynamics during dripping using the volume-of-fluid method, Physics of Fluids, Vol. 31, 2019. doi: 10.1063/1.5061715

75-18   Tobias Ladner, Sebastian Odenwald, Kevin Kerls, Gerald Zieres, Adeline Boillon and Julien Bœuf, CFD supported investigation of shear induced by bottom-mounted magnetic stirrer in monoclonal antibody formulation, Pharmaceutical Research, Vol. 35, 2018. doi: 10.1007/s11095-018-2492-4

53-18   Venoos Amiri Roodan, Jenifer Gómez-Pastora, Aditi Verma, Eugenio Bringas, Inmaculada Ortiz and Edward P. Furlani, Computational analysis of magnetic droplet generation and manipulation in microfluidic devices, Proceedings of the 5th International Conference of Fluid Flow, Heat and Mass Transfer, Niagara Falls, Canada, June 7 – 9, 2018; Paper no. 154, 2018.  doi: 10.11159/ffhmt18.154

35-18   Jenifer Gómez-Pastora, Cristina González Fernández, Marcos Fallanza, Eugenio Bringas and Inmaculada Ortiz, Flow patterns and mass transfer performance of miscible liquid-liquid flows in various microchannels: Numerical and experimental studies, Chemical Engineering Journal, vol. 344, pp. 487-497, 2018. doi: 10.1016/j.cej.2018.03.110

16-18   P. Schneider, V. Sukhotskiy, T. Siskar, L. Christie and I.H. Karampelas, Additive Manufacturing of Microfluidic Components via Wax Extrusion, Biotech, Biomaterials and Biomedical TechConnect Briefs, vol. 3, pp. 162 – 165, 2018.

15-18   J. Gómez-Pastora, I.H. Karampelas, A.Q. Alorabi, M.D. Tarn, E. Bringas, A. Iles, V.N. Paunov, N. Pamme, E.P. Furlani, I. Ortiz, CFD analysis and experimental validation of magnetic droplet generation and deflection across multilaminar flow streams, Biotech, Biomaterials and Biomedical TechConnect Briefs, vol. 3, pp. 182-185, 2018.

14-18   J. Gómez-Pastora, C. González-Fernández, I.H. Karampelas, E. Bringas, E.P. Furlani, and I. Ortiz, Design of Magnetic Blood Cleansing Microdevices through Experimentally Validated CFD Modeling, Biotech, Biomaterials and Biomedical TechConnect Briefs, vol. 3, pp. 170-173, 2018.

10-18   A. Gupta, I.H. Karampelas, J. Kitting, Numerical modeling of the formation of dynamically configurable L2 lens in a microchannel, Biotech, Biomaterials and Biomedical TechConnect Briefs, Vol. 3, pp. 186 – 189, 2018.

17-17   I.H. Karampelas, J. Gómez-Pastora, M.J. Cowan, E. Bringas, I. Ortiz and E.P. Furlani, Numerical Analysis of Acoustophoretic Discrete Particle Focusing in Microchannels, Biotech, Biomaterials and Biomedical TechConnect Briefs 2017, Vol. 3

16-17   J. Gómez-Pastora, I.H. Karampelas, E. Bringas, E.P. Furlani and I. Ortiz, CFD analysis of particle magnetophoresis in multiphase continuous-flow bioseparators, Biotech, Biomaterials and Biomedical TechConnect Briefs 2017, Vol. 3

15-17   I.H. Karampelas, S. Vader, Z. Vader, V. Sukhotskiy, A. Verma, G. Garg, M. Tong and E.P. Furlani, Drop-on-Demand 3D Metal Printing, Informatics, Electronics and Microsystems TechConnect Briefs 2017, Vol. 4

102-16   J. Brindha, RA.G. Privita Edwina, P.K. Rajesh and P.Rani, “Influence of rheological properties of protein bio-inks on printability: A simulation and validation study,” Materials Today: Proceedings, vol. 3, no.10, pp. 3285-3295, 2016. doi: 10.1016/j.matpr.2016.10.010

99-16   Ioannis H. Karampelas, Kai Liu, Fatema Alali, and Edward P. Furlani, Plasmonic Nanoframes for Photothermal Energy Conversion, J. Phys. Chem. C, 2016, 120 (13), pp 7256–7264

98-16   Jelena Dinic and Vivek Sharma, Drop formation, pinch-off dynamics and liquid transfer of simple and complex fluidshttp://meetings.aps.org/link/BAPS.2016.MAR.B53.12, APS March Meeting 2016, Volume 61, Number 2, March 14–18, 2016, Baltimore, Maryland

67-16  Vahid Bazargan and Boris Stoeber, Effect of substrate conductivity on the evaporation of small sessile droplets, PHYSICAL REVIEW E 94, 033103 (2016), doi: 10.1103/PhysRevE.94.033103

57-16   Ioannis Karampelas, Computational analysis of pulsed-laser plasmon-enhanced photothermal energy conversion and nanobubble generation in the nanoscale, PhD Dissertation: Department of Chemical and Biological Engineering, University at Buffalo, State University of New York, July 2016

44-16   Takeshi Sawada et al., Prognostic impact of circulating tumor cell detected using a novel fluidic cell microarray chip system in patients with breast cancer, EBioMedicine, Available online 27 July 2016, doi: 10.1016/j.ebiom.2016.07.027.

39-16   Chien-Hsun Wang, Ho-Lin Tsai, Yu-Che Wu and Weng-Sing Hwang, Investigation of molten metal droplet deposition and solidification for 3D printing techniques, IOP Publishing, J. Micromech. Microeng. 26 (2016) 095012 (14pp), doi: 10.1088/0960-1317/26/9/095012, July 8, 2016

30-16   Ioannis H. Karampelas, Kai Liu and Edward P. Furlani, Plasmonic Nanocages as Photothermal Transducers for Nanobubble Cancer Therapy, Nanotech 2016 Conference & Expo, May 22-25, Washington, DC.

29-16   Scott Vader, Zachary Vader, Ioannis H. Karampelas and Edward P. Furlani, Advances in Magnetohydrodynamic Liquid Metal Jet Printing, Nanotech 2016 Conference & Expo, May 22-25, Washington, DC.

02-16  Stephen D. Hoath (Editor), Fundamentals of Inkjet Printing: The Science of Inkjet and Droplets, ISBN: 978-3-527-33785-9, 472 pages, February 2016 (see chapters 2 and 3 for FLOW-3D results)

125-15   J. Berthier, K.A. Brakke, E.P. Furlani, I.H. Karampelas, V. Poher, D. Gosselin, M. Cubinzolles and P. Pouteau, Whole blood spontaneous capillary flow in narrow V-groove microchannels, Sensors and Actuators B: Chemical, 206, pp. 258-267, 2015.

86-15   Yousub Lee and Dave F. Farson, Simulation of transport phenomena and melt pool shape for multiple layer additive manufacturing, J. Laser Appl. 28, 012006 (2016). doi: 10.2351/1.4935711, published online 2015.

77-15   Ho-Lin Tsai, Weng-Sing Hwang, Jhih-Kai Wang, Wen-Chih Peng and Shin-Hau Chen, Fabrication of Microdots Using Piezoelectric Dispensing Technique for Viscous Fluids, Materials 2015, 8(10), 7006-7016. doi: 10.3390/ma8105355

63-15   Scott Vader, Zachary Vader, Ioannis H. Karampelas and Edward P. Furlani, Magnetohydrodynamic Liquid Metal Jet Printing, TechConnect World Innovation Conference & Expo, Washington, D.C., June 14-17, 2015

46-15   Adwaith Gupta, 3D Printing Multi-Material, Single Printhead Simulation, Advanced Qualification of Additive Manufacturing Materials Workshop, July 20 – 21, 2015, Santa Fe, NM

28-15   Yongqiang Li, Mingzhu Hu, Ling Liu, Yin-Yin Su, Li Duan, and Qi Kang, Study of Capillary Driven Flow in an Interior Corner of Rounded Wall Under MicrogravityMicrogravity Science and Technology, June 2015

20-15   Pamela J. Waterman, Diversity in Medical Simulation Applications, Desktop Engineering, May 2015, pp 22-26,

16-15   Saurabh Singh, Ann Junghans, Erik Watkins, Yash Kapoor, Ryan Toomey, and Jaroslaw Majewski, Effects of Fluid Shear Stress on Polyelectrolyte Multilayers by Neutron Scattering Studies, © 2015 American Chemical Society, DOI: 10.1021/acs.langmuir.5b00037, Langmuir 2015, 31, 2870−2878, February 17, 2015

11-15   Cheng-Han Wu and Weng-Sing Hwang, The effect of process condition of the ink-jet printing process on the molten metallic droplet formation through the analysis of fluid propagation direction, Canadian Journal of Physics, 2015. doi: 10.1139/cjp-2014-0259

03-15 Hanchul Cho, Sivasubramanian Somu, Jin Young Lee, Hobin Jeong and Ahmed Busnaina, High-Rate Nanoscale Offset Printing Process Using Directed Assembly and Transfer of Nanomaterials, Adv. Materials, doi: 10.1002/adma.201404769, February 2015

122-14  Albert Chi, Sebastian Curi, Kevin Clayton, David Luciano, Kameron Klauber, Alfredo Alexander-Katz, Sebastián D’hers and Noel M Elman, Rapid Reconstitution Packages (RRPs) implemented by integration of computational fluid dynamics (CFD) and 3D printed microfluidics, Research Gate, doi: 10.1007/s13346-014-0198-7, July 2014

113-14 Cihan Yilmaz, Arif E. Cetin, Georgia Goutzamanidis, Jun Huang, Sivasubramanian Somu, Hatice Altug, Dongguang Wei and Ahmed Busnaina, Three-Dimensional Crystalline and Homogeneous Metallic Nanostructures Using Directed Assembly of Nanoparticles, 10.1021/nn500084g, © 2014 American Chemical Society, April 2014

110-14 Koushik Ponnuru, Jincheng Wu, Preeti Ashok, Emmanuel S. Tzanakakis and Edward P. Furlani, Analysis of Stem Cell Culture Performance in a Microcarrier Bioreactor System, Nanotech, Washington, D.C., June 15-18, 2014

109-14   Ioannis H. Karampelas, Young Hwa Kim and Edward P. Furlani, Numerical Analysis of Laser Induced Photothermal Effects using Colloidal Plasmonic Nanostructures, Nanotech, Washington, D.C., June 15-18, 2014

108-14   Chenxu Liu, Xiaozheng Xue and Edward P. Furlani, Numerical Analysis of Fully-Coupled Particle-Fluid Transport and Free-Flow Magnetophoretic Sorting in Microfluidic Systems, Nanotech, Washington, D.C., June 15-18, 2014

95-14   Cheng-Han Wu, Weng-Sing Hwang, The effect of the echo-time of a bipolar pulse waveform on molten metallic droplet formation by squeeze mode piezoelectric inkjet printing, Accepted November 2014, Microelectronics Reliability (2014) , © 2014 Elsevier Ltd. All rights reserved.

85-14   Sudhir Srivastava, Lattice Boltzmann method for contact line dynamics, ISBN: 978-90-386-3608-5, Copyright © 2014 S. Srivastava

61-14   Chenxu Liu, A Computational Model for Predicting Fully-Coupled Particle-Fluid Dynamics and Self-Assembly for Magnetic Particle Applications, Master’s Thesis: State University of New York at Buffalo, 2014, 75 pages; 1561583, http://gradworks.umi.com/15/61/1561583.html

41-14 Albert Chi, Sebastian Curi, Kevin Clayton, David Luciano, Kameron Klauber, Alfredo Alexander-Katz, Sebastian D’hers, and Noel M. Elman, Rapid Reconstitution Packages (RRPs) implemented by integration of computational fluid dynamics (CFD) and 3D printed microfluidics, Drug Deliv. and Transl. Res., DOI 10.1007/s13346-014-0198-7, # Controlled Release Society 2014. Available for purchase online at SpringerLink.

21-14  Suk-Hee Park, Ung Hyun Koh, Mina Kim, Dong-Yol Yang, Kahp-Yang Suh and Jennifer Hyunjong Shin, Hierarchical multilayer assembly of an ordered nanofibrous scaffold via thermal fusion bonding, Biofabrication 6 (2014) 024107 (10pp), doi:10.1088/1758-5082/6/2/024107, IOP Publishing, 2014. Available for purchase online at IOP.

17-14   Vahid Bazargan, Effect of substrate cooling and droplet shape and composition on the droplet evaporation and the deposition of particles, Ph.D. Thesis: Department of Mechanical Engineering, The University of British Columbia, March 2014, © Vahid Bazargan, 2014

73-13  Oliver G. Harlen, J. Rafael Castrejón-Pita, and Arturo Castrejon-Pita, Asymmetric Detachment from Angled Nozzles Plates in Drop-on Demand Inkjet Printing, NIP & Digital Fabrication Conference, 2013 International Conference on Digital Printing Technologies. Pages 253-549, pp. 277-280(4)

63-13  Fatema Alali, Ioannis H. Karampelas, Young Hwa Kim, and Edward P. Furlani, Photonic and Thermofluidic Analysis of Colloidal Plasmonic Nanorings and Nanotori for Pulsed-Laser Photothermal ApplicationsJ. Phys. Chem. C, Article ASAP, DOI: 10.1021/jp406986y, Copyright © 2013 American Chemical Society, September 2013.

25-13  Sudhir Srivastava, Theo Driessen, Roger Jeurissen, Herma Wijshoff, and Federico Toschi, Lattice Boltzmann Method to Study the Contraction of a Viscous Ligament, International Journal of Modern Physics © World Scientific Publishing Company, May 2013.

11-13  Li-Chieh Hsu, Yong-Jhih Chen, Jia-Huang Liou, Numerical Investigation in the Factors on the Pool Boiling, Applied Mechanics and Materials Vol. 311 (2013) pp 456-461, © (2013) Trans Tech Publications, Switzerland, doi:10.4028/www.scientific.net/AMM.311.456. Available for purchase online at Scientific.Net.

10-13 Pamela J. Waterman, CFD: Shaping the Medical World, Desktop Engineering, April 2013. Full article available online at Desktop Engineering.

90-12 Charles R. Ortloff and Martin Vogel, Spray Cooling Heat Transfer- Test and CFD Analysis, Electronics Cooling, June 2012. Available online at Electronics Cooling.

79-12    Daniel Parsaoran Siregar, Numerical simulation of evaporation and absorption of inkjet printed droplets, Ph.D. Thesis: Technische Universiteit Eindhoven, September 18, 2012, Copyright 2012 by D.P. Siregar, ISBN: 978-90-386-3190-5.

71-12   Jong-hyeon Chang, Kyu-Dong Jung, Eunsung Lee, Minseog Choi, Seungwan Lee, and Woonbae Kim, Varifocal liquid lens based on microelectrofluidic technology, Optics Letters, Vol. 37, Issue 21, pp. 4377-4379 (2012) http://dx.doi.org/10.1364/OL.37.004377

70-12   Jong-hyeon Chang, Kyu-Dong Jung, Eunsung Lee, Minseog Choi, and Seunwan Lee, Microelectrofluidic Iris for Variable ApertureProc. SPIE 8252, MOEMS and Miniaturized Systems XI, 82520O (February 9, 2012); doi:10.1117/12.906587

69-12   Jong-hyeon Chang, Eunsung Lee, Kyu-Dong Jung, Seungwan Lee, Minseog Choi, and  Woonbae Kim, Microelectrofluidic Lens for Variable CurvatureProc. SPIE 8486, Current Developments in Lens Design and Optical Engineering XIII, 84860X (October 11, 2012); doi:10.1117/12.925852.

61-12  Biddut Bhattacharjee, Study of Droplet Splitting in an Electrowetting Based Digital Microfluidic System, Thesis: Doctor of Philosophy in the College of Graduate Studies (Applied Sciences), The University of British Columbia, September 2012, © Biddut Bhattacharjee.

55-12 Hejun Li, Pengyun Wang, Lehua Qi, Hansong Zuo, Songyi Zhong, Xianghui Hou, 3D numerical simulation of successive deposition of uniform molten Al droplets on a moving substrate and experimental validation, Computational Materials Science, Volume 65, December 2012, Pages 291–301. Available for purchase online at SciVerse.

54-12   Edward P. Furlani, Anthony Nunez, Gianmarco Vizzeri, Modeling Fluid Structure-Interactions for Biomechanical Analysis of the Human Eye, Nanotech Conference & Expo, June 18-21, 2012, Santa Clara, CA.

53-12   Xinyun Wu, Richard D. Oleschuk and Natalie M. Cann, Characterization of microstructured fibre emitters in pursuit of improved nano electrospray ionization performance, The Royal Society of Chemistry 2012, http://pubs.rsc.org, DOI: 10.1039/c2an35249d, May 2012

25-12    Edward P. Furlani, Ioannis H. Karampelas and Qian Xie, Analysis of Pulsed Laser Plasmon-assisted Photothermal Heating and Bubble Generation at the Nanoscale, Lab on a Chip, 10.1039/C2LC40495H, Received 01 May 2012, Accepted 07 Jun 2012. First published on the web 13 Jun 2012.

22-12  R.A. Sultanov, D. Guster, Numerical Modeling and Simulations of Pulsatile Human Blood Flow in Different 3D-Geometries, Book chapter #21 in Fluid Dynamics, Computational Modeling and Applications (2012), ISBN: 978-953-51-0052-2, p. 475 [18 pages]. Available online at INTECH.

21-12  Guo-Wei Huang, Tzu-Yi Hung, and Chin-Tai Chen, Design, Simulation, and Verification of Fluidic Light-Guide Chips with Various Geometries of Micro Polymer Channels, NEMS 2012, Kyoto, Japan, March 5-8, 2012. Available for purchase online at IEEE.

103-11   Suk-Hee Park, Development of Three-Dimensional Scaffolds containing Electrospun Nanofibers and their Applications to Tissue Regeneration, Ph.D. Thesis: School of Mechanical, Aersospace and Systems Engineering, Division of Mechanical Engineering, KAIST, 2011.

81-11   Xinyun Wu, Modeling and Characterization of Microfabricated Emitters-In Pursuit of Improved ESI-MS Performance, thesis: Department of Chemistry, Queen’s University, December 2011, Copyright © Xinyun Wu, 2011

79-11  Cong Lu, A Cell Preparation Stage for Automatic Cell Injection, thesis: Graduate Department of Mechanical and Industrial Engineering, University of Toronto, Copyright © Cong Lu, 2011

77-11 Ge Bai, W. Thomas Leach, Computational fluid dynamics (CFD) insights into agitation stress methods in biopharmaceutical development, International Journal of Pharmaceutics, Available online 8 December 2011, ISSN 0378-5173, 10.1016/j.ijpharm.2011.11.044. Available online at SciVerse.

72-11  M.R. Barkhudarov, C.W. Hirt, D. Milano, and G. Wei, Comments on a Comparison of CFD Software for Microfluidic Applications, Flow Science Technical Note #93, FSI-11-TN93, December 2011

45-11  Chang-Wei Kang, Jiak Kwang Tan, Lunsheng Pan, Cheng Yee Low and Ahmed Jaffar, Numerical and experimental investigations of splat geometric characteristics during oblique impact of plasma spraying, Applied Surface Science, In Press, Corrected Proof, Available online 20 July 2011, ISSN 0169-4332, DOI: 10.1016/j.apsusc.2011.06.081. Available to purchase online at SciVers

33-11  Edward P. Furlani, Mark T. Swihart, Natalia Litchinitser, Christopher N. Delametter and Melissa Carter, Modeling Nanoscale Plasmon-assisted Bubble Nucleation and Applications, Nanotech Conference and Expo 2011, Boston, MA, June 13-16, 2011

32-11  Lu, Cong and Mills, James K., Three cell separation design for realizing automatic cell injection, Complex Medical Engineering (CME), 2011 IEEE/ICME, pp: 599 – 603, Harbin, China, 10.1109/ICCME.2011.5876811, June 2011. Available online at IEEEXplore.

25-11 Issam M. Bahadur, James K. Mills, Fluidic vacuum-based biological cell holding device with piezoelectrically induced vibration, Complex Medical Engineering (CME), 2011 IEEE/ICME International Conference on, 22-25 May 2011, pp: 85 – 90, Harbin, China. Available online at: IEEE Xplore.

14-11  Edward P. Furlani, Roshni Biswas, Alexander N. Cartwright and Natalia M. Litchinitser, Antiresonant guiding optofluidic biosensor, doi:10.1016/j.optcom.2011.04.014, Optics Communication, April 2011

05-11 Hyeju Eom and Keun Park, Integrated numerical analysis to evaluate replication characteristics of micro channels in a locally heated mold by selective induction, International Journal of Precision Engineering and Manufacturing, Volume 12, Number 1, 53-60, DOI: 10.1007/s12541-011-0007-x, 2011. Available online at: SpringerLink.

70-10  I.N. Volnov, V.S. Nagornyi, Modeling Processes for Generation of Streams of Monodispersed Fluid Droplets in Electro-inkjet Applications, Science and Technology News, St. Petersburg State Polytechnic University, 4, pp 294-300, 2010. In Russian.

62-10  F. Mobadersani, M. Eskandarzade, S. Azizi and S. Abbasnezhad, Effect of Ambient Pressure on Bubble Growth in Micro-Channel and Its Pumping Effect, ESDA2010-24436, pp. 577-584, doi:10.1115/ESDA2010-24436, ASME 2010 10th Biennial Conference on Engineering Systems Design and Analysis (ESDA2010), Istanbul, Turkey, July 12–14, 2010. Available online at the ASME Digital Library.

58-10 Tsung-Yi Ho, Jun Zeng, and Chakrabarty, K, Digital microfluidic biochips: A vision for functional diversity and more than moore, Computer-Aided Design (ICCAD), 2010 IEEE/ACM International Conference on, DOI: 10.1109/ICCAD.2010.5654199, © IEEE, November 2010. Available online at IEEE Explore.

51-10  Regina Bleul, Marion Ritzi-Lehnert, Julian Höth, Nico Scharpfenecker, Ines Frese, Dominik Düchs, Sabine Brunklaus, Thomas E. Hansen-Hagge, Franz-Josef Meyer-Almes, Klaus S. Drese, Compact, cost-efficient microfluidics-based stopped-flow device, Anal Bioanal Chem, DOI 10.1007/s00216-010-4446-5, Available online at Springer, November 2010

22-10    Krishendu Chakrabarty, Richard B. Fair and Jun Zeng, Design Tools for Digital Microfluidic Biochips Toward Functional Diversification and More than Moore, IEEE Transactions on Computer-Aided Design of Integrated Circuits and Systems, Vol. 29, No. 7, July 2010

14-10 E. P. Furlani and M. S. Hanchak, Nonlinear analysis of the deformation and breakup of viscous microjets using the method of lines, International Journal for Numerical Methods in Fluids (2010), © 2010 John Wiley & Sons, Ltd., Published online in Wiley InterScience. DOI: 10.1002/fld.2205

55-09 R.A. Sultanov, and D. Guster, Computer simulations of  pulsatile human blood flow through 3D models of the human aortic arch, vessels of simple geometry and a bifurcated artery, Proceedings of the 31st Annual International Conference of the IEEE EMBS (Engineering in Medicine and Biology Society), Minneapolis, September 2-6, 2009, p.p. 4704-4710.

30-09 Anurag Chandorkar and Shayan Palit, Simulation of Droplet Dynamics and Mixing in Microfluidic Devices using a VOF-Based Method, Sensors & Transducers journal, ISSN 1726-5479 © 2009 by IFSA, Vol.7, Special Issue “MEMS: From Micro Devices to Wireless Systems,” October 2009, pp. 136-149.

13-09 E.P. Furlani, M.C. Carter, Analysis of an Electrostatically Actuated MEMS Drop Ejector, Presented at Nanotech Conference & Expo 2009, Houston, Texas, USA, May 3-7, 2009

12-09 A. Chandorkar, S. Palit, Simulation of Droplet-Based Microfluidics Devices Using a Volume-of-Fluid Approach, Presented at Nanotech Conference & Expo 2009, Houston, Texas, USA, May 3-7, 2009

3-09 Christopher N. Delametter, FLOW-3D Speeds MEMS Inkjet Development, Desktop Engineering, January 2009

42-08  Tien-Li Chang, Jung-Chang Wang, Chun-Chi Chen, Ya-Wei Lee, Ta-Hsin Chou, A non-fluorine mold release agent for Ni stamp in nanoimprint process, Microelectronic Engineering 85 (2008) 1608–1612

26-08 Pamela J. Waterman, First-Pass CFD Analyses – Part 2, Desktop Engineering, November 2008

09-08 M. Ren and H. Wijshoff, Thermal effect on the penetration of an ink droplet onto a porous medium, Proc. Eurotherm2008 MNH, 1 (2008)

04-08 Delametter, Christopher N., MEMS development in less than half the time, Small Times, Online Edition, May 2008

02-08 Renat A. Sultanov, Dennis Guster, Brent Engelbrekt and Richard Blankenbecler, 3D Computer Simulations of Pulsatile Human Blood Flows in Vessels and in the Aortic Arch – Investigation of Non-Newtonian Characteristics of Human Blood, The Journal of Computational Physics, arXiv:0802.2362v1 [physics.comp-ph], February 2008

01-08 Herman Wijshoff, thesis: University of Twente, Structure- and fluid dynamics in piezo inkjet printheads, ISBN 978-90-365-2582-4, Venlo, The Netherlands January 2008.

30-07 A. K. Sen, J. Darabi, and D. R. Knapp, Simulation and parametric study of a novel multi-spray emitter for ESI–MS applications, Microfluidics and Nanofluidics, Volume 3, Number 3, June 2007, pp. 283-298(16)

28-07 Dan Soltman and Vivek Subramanian, Inkjet-Printed Line Morphologies and Temperature Control of the Coffee Ring Effect, Langmuir; 2008; ASAP Web Release Date: 16-Jan-2008; (Research Article) DOI: 10.1021/la7026847

23-07 A K Sen and J Darabi, Droplet ejection performance of a monolithic thermal inkjet print head, Journal of Micromechanical and Microengineering,vol.17, pp.1420-1427 (2007) doi:10.1088/0960-1317/17/8/002; Abstract only.

18-07 Herman Wisjhoff, Better Printheads Via Simulation, Desktop Engineering, October 2007, Vol. 13, Issue 2

17-07 Jos de Jong, Ph.D. Thesis: University of Twente, Air entrapment in piezo inkjet printing, ISBN 978-90-365-2483-4, April 2007

15-07 Krishnendu Chakrabarty and Jun Zeng, (Ed.), Design Automation Methods and Tools for Microfluidics-Based Biochips, Springer, September 2006.

14-07 Fei Su and Jun Zeng, Computer-aided design and test for digital microfluidics, IEEE Design & Test of Computers, 24(1), 2007, 60-70.

13-07 Jun Zeng, Modeling and simulation of electrified droplets and its application to computer-aided design of digital microfluidics, IEEE Transactions on Computer-Aided Design of Integrated Circuits and Systems, 25(2), 2006, 224-233.

12-07 Krishnendu Chakrabarty and Jun Zeng, (2005), Automated top-down design for microfluidic biochips, ACM Journal on Emerging Technologies in Computing Systems, 1(3), 2005, 186–223.

01-07 Wijshoff, Herman, Drop formation mechanisms in piezo-acoustic inkjet, NSTI-Nanotech 2007, ISBN 1420061844 Vol. 3, 2007)

23-06 John J. Uebbing, Stephan Hengstler, Dale Schroeder, Shalini Venkatesh, and Rick Haven, Heat and Fluid Flow in an Optical Switch Bubble, Journal of Microelectromechanical Systems, Vol. 15, No. 6, December 2006

21-06 Wijshoff, Herman, Manipulating Drop Formation in Piezo Acoustic Inkjet, Proc. IS&T’s NIP22, 79 (2006)

20-06 J. de Jong, H. Reinten, M. van den Berg, H. Wijshoff, M. Versluis, G. de Bruin, A. Prosperetti and D. Lohse, Air entrapment in piezo-driven inkjet printheads, J. Acoust. Soc. Am. 120(3), 1257 (2006)

11-06 A. K. Sen, J. Darabi, D. R. Knapp and J. Liu, Modeling and Characterization of a Carbon Fiber Emitter for Electrospray Ionization, 1 MEMS and Microsystems Laboratory, Department of Mechanical Engineering, University of South Carolina, 300 Main Street, Columbia, SC 29208, USA, 2 Department of Pharmacology, Medical University of South Carolina, Charleston, SC

5-06 E. P. Furlani, B. G. Price, G. Hawkins, and A. G. Lopez, Thermally Induced Marangoni Instability of Liquid Microjets with Application to Continuous Inkjet Printing, Proceedings of NSTI Nanotech Conference 2006, Vol. 2, pp 534-537.

28-05 O B Fawehinmi, P H Gaskell, P K Jimack, N Kapur, and H M Thompson, A combined experimental and computational fluid dynamics analysis of the dynamics of drop formation, May 2005. DOI: 10.1243/095440605X31788

5-05 E. P. Furlani, Thermal Modulation and Instability of Newtonian Liquid Microjets, presented at Nanotech 2005, Anaheim, CA, May 8-12, 2005.

1-05 C.W. Hirt, Electro-Hydrodynamics of Semi-Conductive Fluids: With Application to Electro-Spraying, Flow Science Technical Note #70, FSI-05-TN70

19-04 G. F. Yao, Modeling of Electroosmosis Without Resolving Physics Inside a Electric Double Layer, Flow Science Technical Note (FSI-04-TN69)

12-04 Jun Zeng and Tom Korsmeyer, Principles of Droplet Electrohydrodynamics for Lab-on-a-Chip, Lab. Chip. Journal, 2004, 4(4), 265-277

9-04 Constantine N. Anagnostopoulos, James M. Chwalek, Christopher N. Delametter, Gilbert A. Hawkins, David L. Jeanmaire, John A. Lebens, Ali Lopez, and David P. Trauernicht, Micro-Jet Nozzle Array for Precise Droplet Metering and Steering Having Increased Droplet Deflection, Proceedings of the 12th International Conference on Solid State Sensors, Actuators and Microsystems, sponsored by IEEE, Boston, June 8-12, 2003, pp. 368-71

8-04 Christopher N. Delametter, David P. Trauernicht, James M. Chwalek, Novel Microfluidic Jet Deflection – Significant Modeling Challenge with Great Application Potential, Technical Proceedings of the 2002 International Conference on Modeling and Simulation of Microsystems sponsored by NSTI, San Juan, Puerto Rico, April 21-25, 2002, pp. 44-47

6-04 D. Vadillo*, G. Desie**, A Soucemarianadin*, Spreading Behavior of Single and Multiple Drops, *Laboratoire des Ecoulements Geophysiques et Industriels (LEGI), and **AGFA-Gevaert Group N.V., XXI ICTAM, 15-21 August 2004, Warsaw, Poland

2-04 Herman Wijshoff, Free Surface Flow and Acousto-Elastic Interaction in Piezo Inkjet, Nanotech 2004, sponsored by the Nano Science & Technology Institute, Boston, MA, March 2004

30-03 D Souders, I Khan and GF Yao, Alessandro Incognito, and Matteo Corrado, A Numerical Model for Simulation of Combined Electroosmotic and Pressure Driven Flow in Microdevices, 7th International Symposium on Fluid Control, Measurement and Visualization

27-03 Jun Zeng, Daniel Sobek and Tom Korsmeyer, Electro-Hydrodynamic Modeling of Electrospray Ionization – CAD for a µFluidic Device-Mass Spectrometer Interface, Agilent Technologies Inc, paper presented at Transducers 2003, June 03 Boston (note: Reference #10 is to FLOW-3D)

17-03 John Uebbing, Switching Fiber-optic Circuits with Microscopic Bubbles, Sensors Magazine, May 2003, Vol 20, No 5, p 36-42

16-03 CFD Speeds Development of MEMS-based Printing Technology, MicroNano Magazine, June 2003, Vol 8, No 6, p 16

3-03 Simulation Speeds Design of Microfluidic Medical Devices, R&D Magazine, March 2003, pp 18-19

1-03 Simulations Help Microscopic Bubbles Switch Fiber-Optic Circuits, Agilent Technologies, Fiberoptic Product News, January 2003, pp 22-23

27-02 Feng, James Q., A General Fluid Dynamic Analysis of Drop Ejection in Drop-on-Demand Ink Jet Devices, Journal of Imaging Science and Technology®, Volume 46, Number 5, September/October 2002

1-02 Feixia Pan, Joel Kubby, and Jingkuang Chen, Numerical Simulation of Fluid Structure Interaction in a MEMS Diaphragm Drop Ejector, Xerox Wilson Research Center, Institute of Physics Publishing, Journal of Micromechanics and Microengineering, 12 (2002), PII: SO960-1317(02)27439-2, pp. 70-76

48-01   Rainer Gruber, Radial Mass Transfer Enhancement in Bubble-Train Flow, PhD thesis in Engineering Sciences, Rheinisch- Westf alischen Technische Hochschule Aachen, December 2001.

34-01 Furlani, E.P., Delametter, C.N., Chwalek, J.M., and Trauernicht, D., Surface Tension Induced Instability of Viscous Liquid Jets, Fourth International Conference on Modeling and Simulation of Microsystems, April 2001

12-01 C. N. Delametter, Eastman Kodak Company, Micro Resolution, Mechanical Engineering, Col 123/No 7, July 2001, pp 70-72

11-01 C. N. Delametter, Eastman Kodak Company, Surface Tension Induced Instability of Viscous Liquid Jets, Technical Proceeding of the Fourth International Conference on Modeling and Simulation of Microsystems, April 2001

9-01 Aman Khan, Unipath Limited Research and Development, Effects of Reynolds Number on Surface Rolling in Small Drops, PVP-Col 431, Emerging Technologies for Fluids, Structures and Fluids, Structures and Fluid Structure Interaction — 2001

2-00 Narayan V. Deshpande, Significance of Inertance and Resistance in Fluidics of Thermal Ink-Jet Transducers, Journal of Imaging Science and Technology, Volume 40, Number 5, Sept./Oct. 1996, pp.457-461

4-98 D. Deitz, Connecting the Dots with CFD, Mechanical Engineering Magazine, pp. 90-91, March 1998

14-94 M. P. O’Hare, N. V. Deshpande, and D. J. Drake, Drop Generation Processes in TIJ Printheads, Xerox Corporation, Adv. Imaging Business Unit, IS&T’s Tenth International Congress on Advances in Non-Impact Printing, Tech. 1994

14-92 Asai, A.,Three-Dimensional Calculation of Bubble Growth and Drop Ejection in a Bubble Jet Printer, Journal of Fluids Engineering Vol. 114 December 1992:638-641

Coating Bibliography

아래는 코팅 참고 문헌의 기술 문서 모음입니다. 
이 모든 논문은 FLOW-3D  결과를 포함하고 있습니다. FLOW-3D를 사용하여 코팅 공정을 성공적으로 시뮬레이션  하는 방법에 대해 자세히 알아보십시오.

Coating Bibliography

2024년 11월 20일 Update

98-24 Fabiano I. Indicatti, Bo Cheng, Michael Rädler, Elisabeth Stammen, Klaus Dilger, Experimental and numerical investigation of the squeegee process during stencil printing of thick adhesive sealings, The Journal of Adhesion, 2024. doi.org/10.1080/00218464.2024.2356105

130-22   Md Didarul Islam, Himendra Perera, Benjamin Black, Matthew Phillips, Muh-Jang Chen, Greyson Hodges, Allyce Jackman, Yuxuan Liu, Chang-Jin Kim, Mohammed Zikry, Saad Khan, Yong Zhu, Mark Pankow, Jong Eun Ryu, Template-free scalable fabrication of linearly periodic microstructures by controlling ribbing defects phenomenon in forward roll coating for multifunctional applications, Advanced Materials Interfaces, 9.27; 2201237, 2022. doi.org/10.1002/admi.202201237

03-21   Delong Jia, Peng Yi, Yancong Liu, Jiawei Sun, Shengbo Yue, Qi Zhao, Effect of laser­ textured groove wall interface on molybdenum coating diffusion and metallurgical bonding, Surface and Coatings Technology, 405; 126561, 2021. doi.org/10.1016/j.surfcoat.2020.126561

50-19     Peng Yi, Delong Jia, Xianghua Zhan, Pengun Xu, and Javad Mostaghimi, Coating solidification mechanism during plasma-sprayed filling the laser textured grooves, International Journal of Heat and Mass Transfer, Vol. 142, 2019. doi:10.1016/j.ijheatmasstransfer.2019.118451

01-19   Jelena Dinic and Vivek Sharma, Computational analysis of self-similar capillary-driven thinning and pinch-off dynamics during dripping using the volume-of-fluid method, Physics of Fluids, Vol. 31, 2019. doi: 10.1063/1.5061715

85-18   Zia Jang, Oliver Litfin and Antonio Delgado, A semi-analytical approach for prediction of volume flow rate in nip-fed reverse roll coating process, Proceedings in Applied Mathematics and Mechanics, Vol. 18, no. 1, Special Issue: 89th Annual Meeting of the International Association of Applied Mathematics and Mechanics, 2018. doi: 10.1002/pamm.201800317

80-14   Hiroaki Koyama, Kazuhiro Fukada, Yoshitaka Murakami, Satoshi Inoue, and Tatsuya Shimoda, Investigation of Roll-to-Sheet Imprinting for the Fabrication of Thin-film Transistor Electrodes, IEICE TRAN, ELECTRON, VOL.E97-C, NO.11, November 2014

46-14   Isabell Vogeler, Andreas Olbers, Bettina Willinger and Antonio Delgado, Numerical investigation of the onset of air entrainment in forward roll coating, 17th International Coating Science and Technology Symposium September 7-10, 2014 San Diego, CA, USA

17-12  Chi-Feng Lin, Bo-Kai Wang, Carlos Tiu and Ta-Jo Liu, On the Pinning of Downstream Meniscus for Slot Die Coating, Advances in Polymer Technology, Vol. 00, No. 0, 1-9 (2012) © 2012 Wiley Periodicals, Inc. Available online at Wiley.

01-11  Reid Chesterfield, Andrew Johnson, Charlie Lang, Matthew Stainer, and Jonathan Ziebarth, Solution-Coating Technology for AMOLED Displays, Information Display Magazine, 1/11 0362-0972/01/2011-024 © SID 2011.

61-09 Yi-Rong Chang, Chi-Feng Lin and Ta-Jo Liu, Start-up of slot die coating, Polymer Engineering and Science, Vol. 49, pp. 1158-1167, 2009. doi:10.1002/pen.21360

26-06  James M. Brethour, 3-D transient simulation of viscoelastic coating flows, 13th International Coating Science and Technology Symposium, September 2006, Denver, Colorado

19-06  Ivosevic, M., Cairncross, R. A., and Knight, R., 3D Predictions of Thermally Sprayed Polymer Splats Modeling Particle Acceleration, Heating and Deformation on Impact with a Flat Substrate, Int. J. of Heat and Mass Transfer, 49, pp. 3285 – 3297, 2006

9-06  M. Ivosevic, R. A. Cairncross, R. Knight, T. E. Twardowski, V. Gupta, Drexel University, Philadelphia, PA; J. A. Baldoni, Duke University, Durham, NC, Effect of Substrate Roughness on Splatting Behavior of HVOF Sprayed Polymer Particles Modeling and Experiments, International Thermal Spray Conference, Seattle, WA, May 2006.

26-05  Ivosevic, M., Cairncross, R. A., Knight, R., Impact Modeling of Thermally Sprayed Polymer Particles, Proc. International Thermal Spray Conference [ITSC-2005], Eds., DVS/IIW/ASM-TSS, Basel, Switzerland, May 2005.

11-05  Brethour, J., Simulation of Viscoelastic Coating Flows with a Volume-of-fluid Technique, in Proceedings of the 6th European Coating Symposium, Bradford, UK, 2005

1-05 C.W. Hirt, Electro-Hydrodynamics of Semi-Conductive Fluids: With Application to Electro-Spraying, Flow Science Technical Note #70, FSI-05-TN70

38-04 K.H. Ho and Y.Y. Zhao, Modelling thermal development of liquid metal flow on rotating disc in centrifugal atomisation, Materials Science and Engineering, A365, pp. 336-340, 2004. doi:10.1016/j.msea.2003.09.044

30-04  M. Ivosevic, R.A. Cairncross, and R. Knight, Impact Modeling of HVOF Sprayed Polymer Particles, Presented at the 12th International Coating Science and Technology Symposium, Rochester, New York, September 23-25, 2004

29-04  J.M. Brethour and C.W. Hirt, Stains Arising from Dried Liquid Drops, Presented at the 12th International Coating Science and Technology Symposium, Rochester, New York, September 23-25, 2004

20-03  James Brethour, Filling and Emptying of Gravure Cells–A CFD Analysis, Convertech Pacific October 2002, Vol. 10, No 4, p 34-37

4-03   M. Toivakka, Numerical Investigation of Droplet Impact Spreading in Spray Coating of Paper, In Proceedings of 2003 TAPPI 8th Advanced Coating Fundamentals Symposium, TAPPI Press, Atlanta, 2003

28-02  J.M. Brethour and H. Benkreira, Filling and Emptying of Gravure Cells—Experiment and CFD Comparison, 11th International Coating Science and Technology Symposium, September 23-25, 2002, Minneapolis, Minnesota

22-02  Hirt, C.W., and Brethour, J.M., Contact Line on Rough Surfaces with Application to Air Entrainment, Presented at the 11th International Coating Science and Technology Symposium, September 23-25, 2002, Minneapolis, Minnesota. Unpublished.

17-01  J. M. Brethour, C. W. Hirt, Moving Contact Lines on Rough Surfaces, 4th European Coating Symposium, 2001, Belgium

16-01  J. M. Brethour, Filling and Emptying of Gravure Cells–-A CFD Analysis, proceedings of the 4th European Coating Symposium 2001, October 1-4, 2001, Brussels, Belgium

26-00 Ronald H. Miller and Gary S. Strumolo, A Self-Consistent Transient Paint Simulation, Proceedings of IMEC2000: 2000 ASME International Mechanical Engineering Congress and Exposition, November 2000, Orlando, Florida

6-99  C. W. Hirt, Direct Computation of Dynamic Contact Angles and Contact Lines, ECC99 Coating Conference, Erlangen, Germany (FSI-99-00-2), Sept. 1999

7-98 J. E. Richardson and Y. Becker, Three-Dimensional Simulation of Slot Coating Edge Effects, Flow Science Inc, and Polaroid Corporation, presented at the 9th International Coating Science and Technology Symposium, Newark, DE, May 18-20, 1998

6-98  C. W. Hirt and E. Choinski, Simulation of the Wet-Start Process in Slot Coating, Flow Science Inc, and Polaroid Corporation, presented at the 9th International Coating Science and Technology Symposium, Newark, DE, May 18-20, 1998

3-97  C. W. Hirt and J. E. Richardson of Flow Science Inc, and K.S. Chen, Sandia National Laboratory, Simulation of Transient and Three-Dimensional Coating Flows Using a Volume-of-Fluid Technique, presented at the 50th Annual Conference of the Society for Imaging and Science Technology, Boston, MA 18-23 May 1997

2-96 C. W. Hirt, K. S. Chen, Simulation of Slide-Coating Flows Using a Fixed Grid and a Volume-of-Fluid Front-Tracking Technique, presented a the 8th International Coating Process Science & Technology Symposium, February 25-29, 1996, New Orleans, LA

Contact Line Insights

Contact Line Insights

FLOW-3D의 수치 모델링 기능은 코팅 성능 향상에 관심이 있는 엔지니어에게 이상적입니다. 계산 시뮬레이션은 코팅 흐름에 영향을 미치는 다양한 물리적 공정의 상대적 중요성과 효과를 연구하는 훌륭한 방법입니다. 물리적 테스트에서 프로세스를 분리하거나 해당 프로세스의 규모를 임의로 조정하는 것이 항상 가능한 것은 아닙니다. 이 섹션에서는 리 블릿 형성(rivulet formation), 핑거링(fingering), 증발, 거친 표면 위의 접촉선 이동 및 유체 흡수와 관련하여 FLOW-3D의 정적 및 동적 접촉각 처리에 대해 설명합니다.

Static and Dynamic Contact Angles

FLOW-3D는 입력으로 설정된 정적 접촉각의 함수로 동적 접촉각과 자유 표면 인터페이스에서 작용하는 관련 힘을 정확하게 계산하여 유체의 소수성을 캡처 할 수 있습니다. 아래 시뮬레이션은 물방울이 경사면 아래로 이동함에 따라 정적 접촉각이 동적 접촉각에 미치는 영향을 보여줍니다.

L.M. Hocking 박사는 그의 저서 [“A moving fluid interface on a rough surface,” J. Fluid Mech., 76, 801, (1976)]에서 표면에 미세한 요철이 흐름 구조를 유도하기 때문에 Contact line이 고체 표면을 통해 이동할 수 있으며 이는 거시적 관점에서 “velocity slip”로 해석 될 수 있다고 했습니다.

이 가설에 대한 전산 해석은 FLOW-3D를 이용하여 쉽게 수행됩니다. 선택된 테스트는 가로, 규칙적으로 이격 된 직사각형 슬롯 패턴 이차원 고체 표면 구성됩니다. 슬롯은 2mm 깊이 10mm 폭, 그리고 그들 사이 폭 10mm 고체 조각을 갖고 이격 됩니다. 이 크기는 전형적으로 상대적으로 부드러운 표면에 긁힌 모양입니다. 액체와 고체 사이의 정적인 접촉각이 60 °가 되도록 선택 하였습니다. 작동 유체는 물로 선정되었고 시험은 채널을 통해 속도30cm / s의 평균 물높이 15mm의 채널의 바닥에 있는 거친 표면을 두고 구동 이루어져 있습니다. 채널의 상단은 free-slip boundary로 정해집니다.

Hocking의 주장대로 micro-scale 교란이 Large scale 관점에서 보았을 때 계산된 속도장으로 보면 velocity slip의 한 종류로서 해석 될 수 있습니다. 아래는 계산된 수평 속도 분포를 나타내고 있습니다. 이것은 표면 바로 위에 제어 볼륨 층의 계산 된 수평 속도 분포를 제공하는 X-Y 플롯에 그래픽으로 보여 주고 있습니다. 격자 미세화에 의해 표면의 고체 부분의 윗쪽 속도가 영이 되는 경향이 있지만, 슬롯들 위에 있는 속도는 영이 안되게 유지됩니다. 많은 요철 위의 이러한 속도의 평균은 효과적인 슬립으로 해석 될 수 있는 non-zero 수평 이송 속도를 일으킵니다.

Evaporative Effects

분산된 고체 물질을 포함하는 액체 방울이 고체 표면에서 건조되면 고체 물질이 침전물로 남습니다. 이 퇴적물의 패턴은 많은 인쇄, 청소 및 코팅 공정에 중요한 의미를 갖습니다. 한 가지 유형의 침전물의 전형적인 예는 왼쪽 이미지와 같이 유출 된 커피 조각의 둘레를 따라 링 얼룩이 형성되는 “커피 링”문제입니다. 이러한 유형의 링 침전물은 액체의 증발로 인한 표면 장력 구동 흐름의 결과로 발생하며, 특히 방울 주변에서 발생합니다 [1].

Drying

건조는 코팅 공정의 중요한 부분입니다. 잘 도포된 코팅은 건조 결함으로 인해 완전히 손상될 수 있습니다. 건조 중에 온도 및 용질 구배는 밀도 및 표면 장력 구배로 인해 코팅 내 흐름을 유도 할 수 있으며, 이로 인해 잠재적으로 코팅 품질이 손상 될 수 있습니다. FLOW-3D의 증발 잔류물 모델을 통해 사용자는 건조로 인한 흐름을 시뮬레이션하고 값 비싼 물리적 실험에 소요되는 시간을 줄일 수 있습니다.

FLOW-3D’s evaporation residue model simulates a 3D view of residue formed from toluene after drying (magnified 30x)

Modeling Ring Formation

FLOW-3D는 증발이 가장 큰 접촉 라인에서의 증착으로 인해 에지 고정이 발생 함을 보여줍니다.

링 형성 모델링
증발에 의해 접촉 라인에서 생성 된 흐름 시뮬레이션
증발은 증발로 인한 열 손실로 인해 액체를 냉각시킵니다 (색상은 온도를 나타냄). 동시에 고체 표면은 전도에 의해 액체를 가열합니다. 증발은 접촉 라인 근처에서 가장 크므로 액체가 접촉 라인을 향해 흐르게하여 정적 상태를 다시 설정합니다. 최종 결과는 액체가 완전히 증발하는 액체 가장자리에 부유 고체가 증착됩니다.

FLOW-3D의 접촉 선 고정 모델에 대해 자세히 알아보십시오.

Simulation of flow generated at a contact line by evaporation