KR20220012662A - Fluid analysis simulation method and fluid simulation apparatus - Google Patents
Fluid analysis simulation method and fluid simulation apparatus Download PDFInfo
- Publication number
- KR20220012662A KR20220012662A KR1020200091710A KR20200091710A KR20220012662A KR 20220012662 A KR20220012662 A KR 20220012662A KR 1020200091710 A KR1020200091710 A KR 1020200091710A KR 20200091710 A KR20200091710 A KR 20200091710A KR 20220012662 A KR20220012662 A KR 20220012662A
- Authority
- KR
- South Korea
- Prior art keywords
- virtual
- fluid
- particle
- boundary
- simulation
- Prior art date
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/25—Design optimisation, verification or simulation using particle-based methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Fluid Mechanics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Physics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computing Systems (AREA)
- Algebra (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
Description
본 발명은 입자 기반의 유체 해석 시뮬레이션 방법 및 유체 해석 시뮬레이션 장치에 관한 것이다.The present invention relates to a particle-based fluid analysis simulation method and a fluid analysis simulation device.
전산유체역학(Computational Fluid Dynamics, CFD)은 유체 역학의 한 분야로서, 컴퓨터를 이용하여 유체의 동적인 움직임을 수치해석적 방법으로 계산해 내는 것이다. 전산유체역학은 편미분방정식인 나비에-스토크스 방정식(Naiver-Stokes Equation)을 FDM(Finite Difference Method), FEM(Finite Element Method), FVM(Finite Volume Method) 및 SPH(Smoothed Particle Hydrodynamics) 등의 방법을 통해 이산화함으로써 유체의 유동을 연산한다.Computational Fluid Dynamics (CFD) is a field of fluid mechanics that calculates the dynamic motion of a fluid using a computer in a numerical way. Computational fluid dynamics is a partial differential equation, Naiver-Stokes equation, FDM (Finite Difference Method), FEM (Finite Element Method), FVM (Finite Volume Method) and SPH (Smoothed Particle Hydrodynamics) methods such as Calculate the flow of the fluid by discretizing it through
나비에-스토크스 방정식을 계산하는 방법에는 공간 도메인을 작은 공간 격자(Mesh or Grid)로 이산화하여 연산하는 격자 기반 방법과 유체를 다수의 입자의 집합으로 표현하는 입자 기반 방법이 있다.There are two methods for calculating the Navier-Stokes equation: a grid-based method for discretizing a spatial domain into a small mesh or grid and a particle-based method for expressing a fluid as a set of multiple particles.
입자 기반 방법은 해석 대상을 격자로 표현하는 대신 입자로 표현함으로써, 자연 현상이나 물리 현상에 대한 보다 자연스러운 시뮬레이션이 가능하다. 입자 기반 방법에는 SPH(Smoothed Particle Hydrodynamics), MPS(Moving Particle Semi-implicit), LBM(Lattice Boltzmann Method) 등이 있다.In the particle-based method, a more natural simulation of natural or physical phenomena is possible by expressing the analysis target as particles instead of using a grid. Particle-based methods include Smoothed Particle Hydrodynamics (SPH), Moving Particle Semi-implicit (MPS), and Lattice Boltzmann Method (LBM).
입자 기반 방법 중 하나인 SPH(Smoothed Particle Hydrodynamics) 기반의 유체 해석은, 격자 기반 방법과 달리 격자를 생성하는 단계가 생략되므로 해석의 결과물을 비교적 빠르게 시뮬레이션할 수 있다.In fluid analysis based on smoothed particle hydrodynamics (SPH), which is one of the particle-based methods, unlike the grid-based method, the step of generating a grid is omitted, so the results of the analysis can be simulated relatively quickly.
또한, SPH 기반의 유체 해석은 격자를 생성하지 않고 입자를 이용하여 해석하므로, 액체와 기체의 계면과 같은 자유 표면의 해석을 비교적 용이하게 수행할 수 있다.In addition, since the SPH-based fluid analysis uses particles without generating a lattice, analysis of free surfaces such as liquid-gas interfaces can be performed relatively easily.
또한, SPH 기반의 유체 해석은 기체, 액체, 고체 중 두 개 이상을 포함하는 다상 유동의 해석을 비교적 정확하게 수행할 수 있다.In addition, the SPH-based fluid analysis can perform the analysis of multiphase flow including two or more of gas, liquid, and solid relatively accurately.
이러한 장점들로 인해 유체의 유동을 시뮬레이션하는데 있어 최근 SPH가 많이 이용되고 있다.Due to these advantages, recently, SPH has been widely used in simulating the flow of a fluid.
종래의 입자 기반 방법을 이용하여 복수의 입자의 유동을 계산하는 경우, 시뮬레이션 영역의 내부 영역에 위치한 하나의 기준 입자로부터 일정 반경 내에 위치한 복수의 인접 입자를 탐색함으로써, 하나의 입자 및 복수의 인접 입자 간의 유동을 계산한다. When the flow of a plurality of particles is calculated using a conventional particle-based method, one particle and a plurality of adjacent particles are searched for by searching for a plurality of adjacent particles located within a predetermined radius from one reference particle located in the inner region of the simulation area. Calculate the flow between
일반적으로, 입자 기반 방식의 유체 해석은 기준 입자가 이의 인접 입자로부터 받는 물리적 힘(즉, 밀도에 의한 힘, 압력에 의한 힘 및 점성에 의한 힘 등)을 고려함으로써 해당 입자의 유동을 계산한다. In general, particle-based fluid analysis calculates the flow of a reference particle by considering the physical forces it receives from its neighboring particles (ie, density-induced forces, pressure-induced forces, and viscous-induced forces).
일반적으로, 입자 기반 방식의 유체 해석은 기준 입자가 이의 인접 입자로부터 받는 물리적 힘(즉, 압력에 의한 힘, 점성에 의한 힘 및 외력에 의한 힘 등)을 고려함으로써 해당 입자의 유동을 계산한다. In general, particle-based fluid analysis calculates the flow of a reference particle by considering the physical force (ie, pressure force, viscous force, external force, etc.) that a reference particle receives from its adjacent particle.
한편, 입자 기반 방식의 유체 해석에서 시뮬레이션 영역의 경계 부근에 위치한 입자는 경계 외부에 배치된 가상 입자를 이용하여 유동 데이터를 계산한다. 이 때, 가상 입자는 크기에 따라 일정한 간격으로 배치되어야 한다. 이는, 각 유체 입자가 커널 함수를 통해 주변 인접 입자의 물리량을 참조하여 해당 유체 입자의 물리량을 계산하므로, 만약, 가상 입자가 불균일하게 배치될 경우, 해당 경계의 수치 오차가 증가하게 된다. 따라서, 유체 해석의 정확도는 가상 입자의 배치 형태 및 가상 입자의 물리량에 의해 좌우된다. Meanwhile, in the particle-based fluid analysis, for particles located near the boundary of the simulation region, flow data is calculated using virtual particles placed outside the boundary. In this case, the virtual particles should be arranged at regular intervals according to their size. In this case, since each fluid particle calculates the physical quantity of the corresponding fluid particle by referring to the physical quantity of the neighboring particles through the kernel function, if the virtual particles are non-uniformly arranged, the numerical error of the boundary increases. Therefore, the accuracy of fluid analysis depends on the arrangement shape of the virtual particles and the physical quantity of the virtual particles.
최근, 입자법에서는 격자법에서 사용되는 경계 조건 중 하나인 노이만(Neumann) 경계 조건이 이용되기 시작하였다. 노이만 경계 조건이란, 미리 알고 있는 물리량을 이용하여 미지의 값을 유추하는 외삽법(extrapolation)으로, 입자법에서는 유체 입자의 물리량을 가상 마커를 통해 가상 입자에게 전달하는 방식으로 구현된다. Recently, in the particle method, the Neumann boundary condition, which is one of the boundary conditions used in the lattice method, has started to be used. The Neumann boundary condition is an extrapolation method for inferring an unknown value using a known physical quantity. In the particle method, the physical quantity of a fluid particle is transmitted to the virtual particle through a virtual marker.
기본적으로, 노이만 경계 조건은 하기의 수학식 1과 같은 압력 구배를 통해 경계 내의 유체 입자의 물리량을 정확히 계산함과 동시에 유체 입자가 경계 밖으로 나가게 되는 벽면 침투 현상을 방지한다.Basically, the Neumann boundary condition accurately calculates the physical quantity of the fluid particles in the boundary through the pressure gradient as in
수학식 1에 따르면, 노이만 경계 조건에서는 법선 벡터를 구하기 위해 마커의 법선 방향을 결정할 필요가 있다.According to
이하에서는 도 1a 및 도 1b를 참조하여, 입자법에서 노이만 경계 조건을 이용하여 유체 해석 시뮬레이션하는 과정을 설명하도록 한다. Hereinafter, a process of fluid analysis simulation using the Neumann boundary condition in the particle method will be described with reference to FIGS. 1A and 1B .
도 1a 및 도 1b는 종래의 유체 해석 시뮬레이션 과정을 설명하기 위한 도면이다. 구체적으로, 도 1a는 종래의 시뮬레이션 영역에서 경계(120)가 수직으로 접하는 경우를 도시한 도면이고, 도 1b는 종래의 시뮬레이션 영역에서 경계(120)가 예각으로 접하는 경우를 도시한 도면이다. 여기서, 경계(120)란 구조물의 벽 근처를 의미하는 것일 수 있다. 1A and 1B are diagrams for explaining a conventional fluid analysis simulation process. Specifically, FIG. 1A is a diagram illustrating a case in which the
도 1a 및 도 1b를 참조하면, 입자법에서 노이만 경계 조건을 이용하는 경우, 유체 입자(100)의 물리량이 산출되면, 유체 입자(100)의 물리량이 가상 마커(110)를 통해 가상 입자에게 전달된다. 이를 위해, 가상 마커(110)는 유체 입자(100)가 맞닿는 경계 면에서 각 유체 입자(100)에 대응하여 배치된다. 1A and 1B , when the Neumann boundary condition is used in the particle method, when the physical quantity of the
또한, 가상 마커(110)는 가상 마커(110)가 경계(120)의 어느 면에 소속되었는지에 따라 가상 마커(110) 각각의 법선 방향이 결정된다. 이 때, 경계(120)의 여러 면이 접하는 곳에 배치된 가상 마커(111, 112, 113)의 경우, 경계(120)의 임의의 한 면에 대하여 법선 방향이 결정된다. Also, in the
즉, 시뮬레이션 영역의 모서리에 배치된 가상 마커(111, 112) 또는 시뮬레이션 영역의 꼭지점에 배치된 가상 마커(113)의 경우, 법선 방향을 구하기 위해 별도의 작업이 필요하거나 법선 방향을 구하기 어려웠다. That is, in the case of the
특히, 도 1b와 같이, 경계(120)가 예각으로 접하는 경우, 경계(120)에서 가상 입자가 등간격으로 배치될 수 없어 가상 마커(113)를 각 유체 입자(100)와 대응하여 배치하기 어렵다는 단점이 있었다. In particular, as shown in FIG. 1B , when the
선행기술문헌: 일본 등록특허 제6009075호Prior art document: Japanese Patent Registration No. 6009075
본 발명은 상술한 문제점을 해결하기 위한 것으로서, 유체에 대한 복수의 유체 입자를 시뮬레이션 영역 내에 배치하고, 경계 조건을 계산하기 위한 복수의 가상 입자를 시뮬레이션 영역 외에 배치하여 시뮬레이션 영역을 모델링하는 유체 해석 시뮬레이션 방법 및 장치를 제공하고자 한다. The present invention is to solve the above problems, and fluid analysis simulation in which a plurality of fluid particles for a fluid are placed in a simulation area and a plurality of virtual particles for calculating a boundary condition are placed outside the simulation area to model the simulation area It is intended to provide a method and apparatus.
복수의 유체 입자의 물리량에 기초하여 시뮬레이션 영역의 경계에 위치한 복수의 가상 마커의 물리량을 산출하고, 복수의 가상 마커의 물리량에 기초하여 경계 조건을 계산하여, 복수의 유체 입자의 물리량 및 경계 조건에 기초하여 유체 해석 시뮬레이션을 수행하는 방법 및 장치를 제공하고자 한다. The physical quantities of a plurality of virtual markers located at the boundary of the simulation area are calculated based on the physical quantities of the plurality of fluid particles, and boundary conditions are calculated based on the physical quantities of the plurality of virtual markers, so that the physical quantities of the plurality of fluid particles and the boundary conditions An object of the present invention is to provide a method and apparatus for performing fluid analysis simulation based on the present invention.
다만, 본 실시예가 이루고자 하는 기술적 과제는 상기된 바와 같은 기술적 과제들로 한정되지 않으며, 또 다른 기술적 과제들이 존재할 수 있다.However, the technical problems to be achieved by the present embodiment are not limited to the technical problems described above, and other technical problems may exist.
상술한 기술적 과제를 달성하기 위한 기술적 수단으로서, 본 발명의 일 실시예는 유체에 대한 복수의 유체 입자를 시뮬레이션 영역 내에 배치하고, 경계 조건을 계산하기 위한 복수의 가상 입자를 상기 시뮬레이션 영역 외에 배치하여 상기 시뮬레이션 영역을 모델링하는 단계, 상기 복수의 유체 입자의 물리량에 기초하여 상기 시뮬레이션 영역의 경계에 위치한 복수의 가상 마커의 물리량을 산출하는 단계, 상기 복수의 가상 마커의 물리량에 기초하여 상기 경계 조건을 계산하는 단계 및 상기 복수의 유체 입자의 물리량 및 상기 경계 조건에 기초하여 유체 해석 시뮬레이션을 수행하는 단계를 포함하고, 상기 복수의 가상 입자의 물리량은 각 가상 입자로부터 기설정된 범위 내의 복수의 마커의 물리량의 평균에 기초한 유체 해석 시뮬레이션 방법을 제공할 수 있다. As a technical means for achieving the above technical problem, an embodiment of the present invention arranges a plurality of fluid particles for a fluid in a simulation area, and arranges a plurality of virtual particles for calculating a boundary condition outside the simulation area, modeling the simulation region; calculating physical quantities of a plurality of virtual markers located at a boundary of the simulation region based on physical quantities of the plurality of fluid particles; determining the boundary condition based on physical quantities of the plurality of virtual markers; calculating and performing a fluid analysis simulation based on the physical quantities of the plurality of fluid particles and the boundary conditions, wherein the physical quantities of the plurality of virtual particles are physical quantities of a plurality of markers within a preset range from each virtual particle It is possible to provide a fluid analysis simulation method based on the average of .
또한, 본 발명의 다른 실시예는 유체에 대한 복수의 유체 입자를 시뮬레이션 영역 내에 배치하고, 경계 조건을 계산하기 위한 복수의 가상 입자를 상기 시뮬레이션 영역 외에 배치하여 상기 시뮬레이션 영역을 모델링하는 모델링부, 상기 복수의 유체 입자의 물리량에 기초하여 상기 시뮬레이션 영역의 경계에 위치한 복수의 가상 마커의 물리량을 산출하는 물리량 산출부, 상기 복수의 가상 마커의 물리량에 기초하여 상기 경계 조건을 계산하는 경계 조건 계산부 및 상기 복수의 유체 입자의 물리량 및 상기 경계 조건에 기초하여 유체 해석 시뮬레이션을 수행하는 유체 시뮬레이션 수행부를 포함하고, 상기 복수의 가상 입자의 물리량은 각 가상 입자로부터 기설정된 범위 내의 복수의 마커의 물리량의 평균에 기초한 것인 유체 해석 시뮬레이션 장치를 제공할 수 있다. In addition, another embodiment of the present invention provides a modeling unit for modeling the simulation region by arranging a plurality of fluid particles for a fluid in a simulation region and arranging a plurality of virtual particles for calculating a boundary condition outside the simulation region; A physical quantity calculator for calculating physical quantities of a plurality of virtual markers located at the boundary of the simulation region based on physical quantities of a plurality of fluid particles, a boundary condition calculator for calculating the boundary conditions based on physical quantities of the plurality of virtual markers, and and a fluid simulation performing unit performing a fluid analysis simulation based on the physical quantities of the plurality of fluid particles and the boundary condition, wherein the physical quantities of the plurality of virtual particles are an average of physical quantities of a plurality of markers within a preset range from each virtual particle. It is possible to provide a fluid analysis simulation device that is based on
상술한 과제 해결 수단은 단지 예시적인 것으로서, 본 발명을 제한하려는 의도로 해석되지 않아야 한다. 상술한 예시적인 실시예 외에도, 도면 및 발명의 상세한 설명에 기재된 추가적인 실시예가 존재할 수 있다.The above-described problem solving means are merely exemplary, and should not be construed as limiting the present invention. In addition to the exemplary embodiments described above, there may be additional embodiments described in the drawings and detailed description.
전술한 본 발명의 과제 해결 수단 중 어느 하나에 의하면, 본 발명은 가상 마커가 구조물의 모서리에 위치한 경우라 하더라도 별도의 작업 없이 손쉽게 법선 방향을 구하도록 하는 유체 해석 시뮬레이션 방법 및 장치를 제공할 수 있다. According to any one of the above-described problem solving means of the present invention, the present invention can provide a fluid analysis simulation method and apparatus for easily obtaining the normal direction without a separate operation even when the virtual marker is located at the edge of the structure. .
또한, 본 발명은 제 1 경계와 제 2 경계가 예각으로 만나는 경우라 하더라도 각 유체 입자에 대응하여 가상 마커를 배치할 수 있는 유체 해석 시뮬레이션 방법 및 장치를 제공할 수 있다.Also, the present invention may provide a fluid analysis simulation method and apparatus capable of disposing virtual markers corresponding to each fluid particle even when the first boundary and the second boundary meet at an acute angle.
도 1a 및 도 1b는 종래의 입자 기반의 유체 해석 시뮬레이션 과정을 설명하기 위한 도면이다.
도 2는 본 발명의 일 실시예에 따른 유체 해석 시뮬레이션 장치의 블록도이다.
도 3은 본 발명의 일 실시예에 따른 시뮬레이션 영역의 경계에 위치한 복수의 가상 마커의 물리량을 산출하는 과정을 설명하기 위한 예시적인 도면이다.
도 4a 및 도 4b는 본 발명의 일 실시예에 따른 가상 입자로부터 소정의 탐색 반경 내에 위치한 인접 가상 마커를 탐색하는 과정을 설명하기 위한 예시적인 도면이다.
도 5a 및 도 5b는 본 발명의 일 실시예에 따른 제 1 가상 마커의 단위 법선 벡터를 도출하는 과정을 설명하기 위한 예시적인 도면이다.
도 6은 본 발명의 일 실시예에 따른 제 1 경계 조건을 계산하는 과정을 설명하기 위한 예시적인 도면이다.
도 7은 본 발명의 일 실시예에 따른 제 2 경계 조건을 계산하는 과정을 설명하기 위한 예시적인 도면이다.
도 8은 본 발명의 일 실시예에 따른 유체 해석 시뮬레이션 장치에서 유체 해석 시뮬레이션 방법의 순서도이다.1A and 1B are diagrams for explaining a conventional particle-based fluid analysis simulation process.
2 is a block diagram of a fluid analysis simulation apparatus according to an embodiment of the present invention.
FIG. 3 is an exemplary diagram for explaining a process of calculating physical quantities of a plurality of virtual markers located at the boundary of a simulation area according to an embodiment of the present invention.
4A and 4B are exemplary views for explaining a process of searching for an adjacent virtual marker located within a predetermined search radius from a virtual particle according to an embodiment of the present invention.
5A and 5B are exemplary views for explaining a process of deriving a unit normal vector of a first virtual marker according to an embodiment of the present invention.
6 is an exemplary diagram for explaining a process of calculating a first boundary condition according to an embodiment of the present invention.
7 is an exemplary diagram for explaining a process of calculating a second boundary condition according to an embodiment of the present invention.
8 is a flowchart of a fluid analysis simulation method in the fluid analysis simulation apparatus according to an embodiment of the present invention.
아래에서는 첨부한 도면을 참조하여 본 발명이 속하는 기술 분야에서 통상의 지식을 가진 자가 용이하게 실시할 수 있도록 본 발명의 실시예를 상세히 설명한다. 그러나 본 발명은 여러 가지 상이한 형태로 구현될 수 있으며 여기에서 설명하는 실시예에 한정되지 않는다. 그리고 도면에서 본 발명을 명확하게 설명하기 위해서 설명과 관계없는 부분은 생략하였으며, 명세서 전체를 통하여 유사한 부분에 대해서는 유사한 도면 부호를 붙였다. DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS Hereinafter, embodiments of the present invention will be described in detail with reference to the accompanying drawings so that those of ordinary skill in the art can easily implement them. However, the present invention may be embodied in several different forms and is not limited to the embodiments described herein. And in order to clearly explain the present invention in the drawings, parts irrelevant to the description are omitted, and similar reference numerals are attached to similar parts throughout the specification.
명세서 전체에서, 어떤 부분이 다른 부분과 "연결"되어 있다고 할 때, 이는 "직접적으로 연결"되어 있는 경우뿐 아니라, 그 중간에 다른 소자를 사이에 두고 "전기적으로 연결"되어 있는 경우도 포함한다. 또한 어떤 부분이 어떤 구성요소를 "포함"한다고 할 때, 이는 특별히 반대되는 기재가 없는 한 다른 구성요소를 제외하는 것이 아니라 다른 구성요소를 더 포함할 수 있는 것을 의미하며, 하나 또는 그 이상의 다른 특징이나 숫자, 단계, 동작, 구성요소, 부분품 또는 이들을 조합한 것들의 존재 또는 부가 가능성을 미리 배제하지 않는 것으로 이해되어야 한다. Throughout the specification, when a part is "connected" with another part, this includes not only the case of being "directly connected" but also the case of being "electrically connected" with another element interposed therebetween. . In addition, when a part "includes" a certain component, it means that other components may be further included, rather than excluding other components, unless otherwise stated, and one or more other features However, it is to be understood that the existence or addition of numbers, steps, operations, components, parts, or combinations thereof is not precluded in advance.
본 명세서에 있어서 '부(部)'란, 하드웨어에 의해 실현되는 유닛(unit), 소프트웨어에 의해 실현되는 유닛, 양방을 이용하여 실현되는 유닛을 포함한다. 또한, 1 개의 유닛이 2 개 이상의 하드웨어를 이용하여 실현되어도 되고, 2 개 이상의 유닛이 1 개의 하드웨어에 의해 실현되어도 된다.In this specification, a "part" includes a unit realized by hardware, a unit realized by software, and a unit realized using both. In addition, one unit may be implemented using two or more hardware, and two or more units may be implemented by one hardware.
본 명세서에 있어서 단말 또는 디바이스가 수행하는 것으로 기술된 동작이나 기능 중 일부는 해당 단말 또는 디바이스와 연결된 서버에서 대신 수행될 수도 있다. 이와 마찬가지로, 서버가 수행하는 것으로 기술된 동작이나 기능 중 일부도 해당 서버와 연결된 단말 또는 디바이스에서 수행될 수도 있다.Some of the operations or functions described as being performed by the terminal or device in this specification may be instead performed by a server connected to the terminal or device. Similarly, some of the operations or functions described as being performed by the server may also be performed in a terminal or device connected to the server.
이하 첨부된 도면을 참고하여 본 발명의 일 실시예를 상세히 설명하기로 한다.Hereinafter, an embodiment of the present invention will be described in detail with reference to the accompanying drawings.
도 2는 본 발명의 일 실시예에 따른 유체 해석 시뮬레이션 장치의 블록도이다. 도 2를 참조하면, 유체 해석 시뮬레이션 장치(1)는 서버, 데스크탑, 노트북, 키오스크(KIOSK) 및 스마트폰(smartphone), 태블릿 PC를 포함할 수 있다. 다만, 유체 해석 시뮬레이션 장치(1)는 앞서 예시된 것들로 한정 해석되는 것은 아니다. 즉, 유체 해석 시뮬레이션 장치(1)는 후술하는 SPH 기반의 유체 해석 시뮬레이션 방법을 수행하는 프로세서를 탑재한 모든 장치를 포함할 수 있다.2 is a block diagram of a fluid analysis simulation apparatus according to an embodiment of the present invention. Referring to FIG. 2 , the fluid
유체 해석 시뮬레이션 장치(1)는 유체의 3차원 유동 해석을 수행한다. 즉, 유체 해석 시뮬레이션 장치(1)는 3차원 시뮬레이션 영역 및 3차원 시뮬레이션 영역에 위치하는 복수의 입자를 모델링하고, 복수의 입자의 3차원 시뮬레이션 영역 내에서의 유동을 해석한다. 다만, 본원에서는 설명의 편의를 위해 시뮬레이션 영역 및 입자를 2차원으로 표현하여 설명한다.The fluid
유체 해석 시뮬레이션 장치(1)는 SPH(Smoothed Particle Hydrodynamics) 기반으로 유체를 해석하기 위한 시뮬레이션을 수행할 수 있다. SPH(Smoothed Particle Hydrodynamics)는 전산유체역학(Computational Fluid Dynamics, CFD)에서 사용될 수 있는 입자 방식의 유체 해석 기법의 하나이다. SPH는 유체의 움직임을 시뮬레이션하기 위하여, 해석 대상인 유체를 하나 이상의 입자로 표현할 수 있다. SPH는 각 입자를 추적하면서 입자가 가지는 물리량을 계산할 수 있고, 계산 결과에 기초하여 유체 해석 시뮬레이션을 수행할 수 있다.The fluid
본 발명에 따른 유체 해석 시뮬레이션 방법은 유체 해석 시뮬레이션이 실시간으로 계산되는 응용 분야를 포함하지만 이에 한정되지 않고, 유체 해석 시뮬레이션을 필요로 하는 다양한 응용 분야에서도 적용된다. The fluid analysis simulation method according to the present invention includes, but is not limited to, an application field in which a fluid analysis simulation is calculated in real time, and is applied to various application fields requiring fluid analysis simulation.
예시적인 응용 분야는 예를 들어 컴퓨터 게임, 의학 시뮬레이션, 과학 응용 및 컴퓨터 애니메이션을 포함한다. Exemplary applications include, for example, computer games, medical simulations, scientific applications, and computer animation.
도 2를 참조하면, 유체 해석 시뮬레이션 장치(1)는 모델링부(210), 물리량 산출부(220), 경계 조건 계산부(230) 및 유체 해석 시뮬레이션 수행부(240)를 포함할 수 있다. Referring to FIG. 2 , the fluid analysis and
모델링부(210)는 유체에 대한 복수의 유체 입자를 시뮬레이션 영역 내에 배치하고, 경계 조건을 계산하기 위한 복수의 가상 입자를 시뮬레이션 영역 외에 배치하여 시뮬레이션 영역을 모델링할 수 있다. 예를 들어, 모델링부(210)는 사용자로부터 키보드, 마우스, 조이스틱, 터치 스크린 및 마이크 등을 이용하여 지형 정보, 구조물 정보, 경계 조건 정보, 입자 물성 정보 및 중력 가속도 정보 중 적어도 하나를 입력받고, 입력받은 적어도 하나에 기초하여 시뮬레이션 영역을 모델링할 수 있다. The
여기서, 구조물 정보는 밀도, 반발 계수 및 마찰 계수 중 적어도 하나를 포함할 수 있다.Here, the structure information may include at least one of a density, a coefficient of restitution, and a coefficient of friction.
또한, 입자 물성 정보는 입자 반경, 밀도, 점성, 음속 및 초기 속도 중 적어도 하나를 포함할 수 있다.In addition, the particle property information may include at least one of a particle radius, a density, a viscosity, a speed of sound, and an initial velocity.
물리량 산출부(220)는 복수의 유체 입자의 물리량에 기초하여 시뮬레이션 영역의 경계에 위치한 복수의 가상 마커의 물리량을 산출할 수 있다. 복수의 가상 마커의 물리량을 산출하는 과정에 대해서는 도 3을 통해 상세히 설명하도록 한다. The
도 3은 본 발명의 일 실시예에 따른 시뮬레이션 영역의 경계에 위치한 복수의 가상 마커의 물리량을 산출하는 과정을 설명하기 위한 예시적인 도면이다. 도 3을 참조하면, 물리량 산출부(220)는 다음의 수학식 2를 이용하여 복수의 가상 마커(302) 및 복수의 유체 입자(300) 간의 위치 벡터의 차를 산출할 수 있다. FIG. 3 is an exemplary view for explaining a process of calculating physical quantities of a plurality of virtual markers located at the boundary of a simulation area according to an embodiment of the present invention. Referring to FIG. 3 , the
물리량 산출부(220)는 다음의 수학식 3을 이용하여 기존의 커널 함수인 'W'로부터 보정된 커널 함수인 ''를 도출할 수 있다. 여기서, 기존의 커널 함수인 'W'는 유체 입자(300)의 물리량 산출을 위해 이용되고, 보정된 커널 함수인 ''는 가상 마커(302)의 물리량 산출을 위해 이용된다. 이는 커널반경내의 입자들이 적을 경우 보정된 커널을 사용하여 부족한 주변입자들을 보완하기 위함이다.The physical
수학식 3을 참조하면, 보정된 커널 함수인 ''는 기존의 커널 함수(W)를 이용한 복수의 가상 마커(302) 및 복수의 유체 입자(300) 간의 위치 벡터의 차(rij) 및 영향 반경(h, 310)에 기초하여 도출될 수 있다. Referring to Equation 3, the corrected kernel function ' ' can be derived based on the difference (r ij ) of the position vectors between the plurality of virtual markers 302 and the plurality of fluid particles 300 and the radius of influence (h, 310) using the existing kernel function (W). .
물리량 산출부(220)는 다음의 수학식 4를 이용하여 복수의 유체 입자(300)에 대한 복수의 가상 마커(302) 각각의 속도를 산출할 수 있다. The
수학식 4를 참조하면, 물리량 산출부(220)는 유체 입자의 밀도(), 유체 입자의 질량(m), 보정된 커널 함수()를 이용한 위치 벡터의 차(rij) 및 영향 반경(h)에 기초하여 복수의 유체 입자(300)에 대한 복수의 가상 마커(302) 각각의 속도를 산출할 수 있다. Referring to Equation 4, the physical
물리량 산출부(220)는 다음의 수학식 5를 이용하여 복수의 유체 입자에 대한 복수의 가상 마커 각각의 압력을 산출할 수 있다. The
수학식 5를 참조하면, 물리량 산출부(220)는 유체 입자의 밀도(), 유체 입자의 질량(m), 보정된 커널 함수()를 이용한 위치 벡터의 차(rij) 및 영향 반경(h)에 기초하여 복수의 유체 입자(300)에 대한 복수의 가상 마커(302) 각각의 압력을 산출할 수 있다. Referring to Equation 5, the physical
다시 도 2로 돌아와서, 경계 조건 계산부(230)는 복수의 가상 마커의 물리량에 기초하여 경계 조건을 계산할 수 있다. 복수의 가상 마커의 물리량에 기초하여 경계 조건을 계산하는 과정에 대해서는 도 4a 내지 도 7b를 통해 상세히 설명하도록 한다. 경계 조건 계산부(230)는 탐색부(235)를 포함할 수 있다. Returning to FIG. 2 , the
도 4a 및 도 4b는 본 발명의 일 실시예에 따른 가상 입자로부터 소정의 탐색 반경 내에 위치한 인접 가상 마커를 탐색하는 과정을 설명하기 위한 예시적인 도면이다. 4A and 4B are exemplary views for explaining a process of searching for an adjacent virtual marker located within a predetermined search radius from a virtual particle according to an embodiment of the present invention.
도 4a를 참조하면, 탐색부(235)는 경계(403)로부터 제 1 가상 입자(400)까지의 거리에 기초하여 탐색 반경(401)을 결정할 수 있다. 예를 들어, 제 1 가상 입자(400)가 경계(403)와 인접한 경우, 탐색 반경(401)은 유효 반경(h)의 2배일 수 있다(즉, 2h). Referring to FIG. 4A , the
이후, 탐색부(235)는 결정된 탐색 반경(401) 내에 위치한 적어도 하나의 인접 가상 마커(402)를 탐색할 수 있다. Thereafter, the
도 4b를 참조하면, 탐색부(235)는 예를 들어, 제 2 가상 입자(410)가 경계(403)와 인접한 있는 제 1 가상 입자(400)보다 경계(403)로부터 멀리 배치된 경우, 유효 반경(h)의 4배일 수 있다(즉, 4h).Referring to FIG. 4B , the
이후, 탐색부(235)는 결정된 탐색 반경(411) 내에 위치한 적어도 하나의 인접 가상 마커(412)를 탐색할 수 있다. Thereafter, the
여기서, 도 4a의 제 1 가상 입자(400)와 도 4b의 제 2 가상 입자(410) 간의 탐색 반경(401, 402)을 각각 다르게 결정하는 이유는, 도 4b에서 제 1 가상 입자(400)와 동일한 탐색 반경으로 제 2 가상 입자(410)로부터 인접한 가상 마커(412)를 탐색하고자 하는 경우, 인접 가상 마커(412)가 탐색되지 않을 수 있기 때문이다. 이와 반대로, 도 4a에서 제 2 가상 입자(410)와 동일한 탐색 반경으로 제 1 가상 입자(400)로부터 인접한 가상 마커(402)를 탐색하는 경우, 인접 가상 마커(402)를 과도하게 많이 탐색하게 되기 때문이다.Here, the reason for differently determining the
즉, 본 발명에 따르면, 경계로부터 가상 입자까지의 거리에 따라 가상 마커를 탐색할 탐색 반경을 결정함으로써, 각 가상 입자마다 유체 해석에 필요한 가상 마커가 적절하게 탐색될 수 있다. That is, according to the present invention, a virtual marker required for fluid analysis can be appropriately searched for each virtual particle by determining a search radius for searching for a virtual marker according to the distance from the boundary to the virtual particle.
물리량 산출부(220)는 각 가상 입자의 탐색 반경 내에 위치한 적어도 하나의 인접 가상 마커의 물리량을 이용하여 각 가상 입자의 물리량으로 산출할 수 있다. 예를 들어, 각 가상 입자의 물리량은 각 가상 입자의 탐색 반경 내에 위치한 적어도 하나의 인접 가상 마커의 물리량의 평균값을 이용하여 각 가상 입자의 물리량으로 산출할 수 있다. 여기서, 복수의 가상 입자의 물리량은 각 가상 입자로부터 기설정된 범위 내의 복수의 마커의 물리량의 평균에 기초한 것일 수 있다.The
가상 입자의 탐색 반경 내의 복수의 가상 입자의 물리량은 다음의 수학식 6을 통해 산출될 수 있다. The physical quantity of the plurality of virtual particles within the search radius of the virtual particle may be calculated through the following Equation (6).
수학식 6을 참조하면, n은 탐색된 적어도 하나의 가상 마커의 수를 의미하고, X는 임의의 물리량을 의미할 수 있다. Referring to Equation 6, n may mean the number of at least one searched virtual marker, and X may mean an arbitrary physical quantity.
도 5a 및 도 5b는 본 발명의 일 실시예에 따른 제 1 가상 마커의 단위 법선 벡터를 도출하는 과정을 설명하기 위한 예시적인 도면이다. 도 5a는 시뮬레이션 영역에서 경계가 수직으로 접하는 경우를 도시한 도면이고, 도 5b는 시뮬레이션 영역에서 경계가 예각으로 접하는 경우를 도시한 도면이다.5A and 5B are exemplary views for explaining a process of deriving a unit normal vector of a first virtual marker according to an embodiment of the present invention. FIG. 5A is a diagram illustrating a case in which a boundary is vertically tangent in a simulation region, and FIG. 5B is a diagram illustrating a case in which a boundary is in contact with an acute angle in the simulation region.
종래의 경우, 가상 마커는 각 유체 입자에 대응하여 배치되어야 했지만, 본 발명에 따르면, 가상 마커는 경계 상에서 등간격으로 배치된다.In the conventional case, the virtual markers had to be arranged corresponding to each fluid particle, but according to the present invention, the virtual markers are arranged at equal intervals on the boundary.
경계 조건 계산부(230)는 탐색된 적어도 하나의 인접 가상 마커의 경계에 대한 법선 벡터를 도출할 수 있다. The
경계 조건 계산부(230)는 도출된 적어도 하나의 인접 가상 마커의 경계에 대한 법선 벡터의 평균을 산출하고, 산출된 평균에 기초하여 제 1 가상 입자의 경계에 대한 단위 법선 벡터를 도출할 수 있다. 예를 들어, 제 1 가상 입자의 탐색 반경 내에 위치한 적어도 하나의 가상 마커의 법선 벡터의 산술 평균에 기초하여 제 1 가상 입자의 경계에 대한 단위 법선 벡터가 도출될 수 있다. The
도 1a 및 도 5a를 참조하면, 구조물의 경계가 수직으로 접하고, 제 1 가상 마커(501)가 구조물의 모서리에 위치하는 경우, 종래에는 도 1a와 같이 제 1 가상 마커의 경계에 대한 법선 벡터를 같이 제 1 가상 마커에 대응하는 가상 입자의 경계에 대한 단위 법선 벡터로서 도출하였으나, 본 발명은 도 5a와 같이 제 1 가상 입자의 탐색 반경 내에 위치한 적어도 하나의 인접 가상 마커의 법선 벡터의 평균을 제 1 가상 입자의 경계에 대한 단위 법선 벡터(502)로 도출할 수 있다. Referring to FIGS. 1A and 5A , when the boundary of the structure is perpendicular to the boundary and the first
도 1b 및 도 5b를 참조하면, 구조물의 경계면이 예각으로 접하고, 제 1 가상 마커(511)가 구조물의 꼭지점에 위치하는 경우, 종래에는 도 1b와 같이 제 1 가상 마커의 경계에 대한 법선 벡터를 같이 제 1 가상 마커에 대응하는 가상 입자의 경계에 대한 단위 법선 벡터로서 도출하였으나, 본 발명은 도 5b와 같이 제 1 가상 입자의 탐색 반경 내에 위치한 적어도 하나의 가상 마커의 법선 벡터의 평균을 제 1 가상 입자의 경계에 대한 단위 법선 벡터(512)로 도출할 수 있다. Referring to FIGS. 1B and 5B , when the boundary surface of the structure is in contact with an acute angle and the first
도 1a 및 도 1b와 같이 종래의 경우, 가상 마커와 가상 입자가 1:1로 대응되어야 하므로, 가상 마커는 가상 입자의 위치에 종속되어야 했다. 이로 인해, 선과 선이 만나는 모서리 부분에서 가상 마커의 배치가 어렵고, 가상 입자 또는 가상 마커의 배치가 어긋난 경우, 입자의 관통 및 수치오차의 증대 등의 문제가 발생하였다. 또한, 가상 입자에 종속되어 가상 마커를 배치하는 경우, 가상 마커가 일정하게 배치되지 않은 경우가 많이 발생됨으로써, 가상 입자에 잘못된 값의 물리량이 전달될 수 있다는 단점이 발생되었다. In the conventional case, as shown in FIGS. 1A and 1B , since the virtual marker and the virtual particle have to correspond 1:1, the virtual marker has to depend on the location of the virtual particle. For this reason, it is difficult to arrange the virtual markers at the corners where the lines meet, and when the virtual particles or the virtual markers are misaligned, problems such as penetration of particles and an increase in numerical error occurred. In addition, when the virtual markers are arranged depending on the virtual particles, there are many cases in which the virtual markers are not uniformly arranged, so that an erroneous physical quantity may be transmitted to the virtual particles.
그러나 본 발명은 도 5a 및 도 5b와 같이 가상 마커가 가상 입자와 별도의 대응 관계를 가지 않고, 가상 입자의 위치와 관계 없이 구조물의 경계에 등간격으로 배치되므로 가상 입자의 비정렬 배치의 영향을 크게 받지 않음으로써, 적절한 물리량을 가상 마커에게 전달할 수 있다는 장점을 가지고 있다. 또한, 여러 개의 가상 마커를 참조하여 가상 입자의 물리량을 정의함으로써, 종래에 계산하기 어려웠던 모서리 부근의 계산이 용이해진다는 장점을 갖는다. However, in the present invention, since the virtual markers do not have a separate correspondence with the virtual particles as shown in FIGS. 5A and 5B and are arranged at equal intervals on the boundary of the structure regardless of the location of the virtual particles, the influence of the non-aligned arrangement of the virtual particles is reduced. By not receiving a large amount, it has the advantage that an appropriate physical quantity can be transmitted to the virtual marker. In addition, by defining the physical quantity of a virtual particle with reference to several virtual markers, it has an advantage that it becomes easy to calculate the vicinity of an edge, which was difficult to calculate in the prior art.
도 6은 본 발명의 일 실시예에 따른 제 1 경계 조건을 계산하는 과정을 설명하기 위한 예시적인 도면이다. 6 is an exemplary diagram for explaining a process of calculating a first boundary condition according to an embodiment of the present invention.
도 6을 참조하면, 경계 조건 계산부(230)는 적어도 하나의 인접 가상 마커의 속도 및 압력의 평균을 산출하고, 산출된 평균에 기초하여 제 1 가상 입자의 속도 및 압력을 도출할 수 있다. 이 때, 경계 조건 계산부(230)는 제 1 가상 입자의 경계에 대한 단위 법선 벡터 및 제 1 가상 입자의 속도에 기초하여 제 1 경계 조건(600)을 계산할 수 있다. 여기서, 제 1 경계 조건(600)으로 경계 근처에서 프리 슬립(free-slip) 또는 논 슬립(non-slip)을 구현하기 위한 경계 조건으로 벽면 경계 조건(wall boundary condition)이 이용될 수 있다. 이를 위해, 제 1 가상 입자의 속도는 다음의 수학식 7을 이용하여 계산하고, 제 1 경계 조건(600)은 다음의 수학식 8을 이용하여 계산할 수 있다. Referring to FIG. 6 , the
수학식 7을 참조하면, Vm은 제 1 가상 입자의 탐색 반경 내에 위치한 적어도 하나의 가상 마커의 속도의 평균을 의미하는 것으로, n은 제 1 가상 입자의 탐색 반경 내에 위치한 적어도 하나의 가상 마커의 수를 의미하고, 은 탐색 반경 내에 위치한 적어도 하나의 가상 마커의 물리량을 의미할 수 있다. Referring to Equation 7, Vm means the average of the velocities of at least one virtual marker located within the search radius of the first virtual particle, and n is the number of at least one virtual marker located within the search radius of the first virtual particle. means, may mean a physical quantity of at least one virtual marker located within the search radius.
수학식 8을 참조하면, nwall은 경계에 대한 단위 법선 벡터이고, Vmarker는 수학식 7을 통해 계산된 가상 마커의 속도이고, 은 계수로서, 은 논-슬립(non-slip) 조건에서는 값이 '0'이고, 슬립(slip) 조건에서는 값이 '1'일 수 있다. 여기서, 제 1 가상 입자의 물리량은 탐색 반경 내에 위치한 하나 이상의 가상 마커의 평균값일 수 있다. Referring to Equation 8, n wall is the unit normal vector to the boundary, Vmarker is the velocity of the virtual marker calculated through Equation 7, is the coefficient, may have a value of '0' under a non-slip condition, and may have a value of '1' under a slip condition. Here, the physical quantity of the first virtual particle may be an average value of one or more virtual markers located within the search radius.
도 7은 본 발명의 일 실시예에 따른 제 2 경계 조건을 계산하는 과정을 설명하기 위한 예시적인 도면이다. 7 is an exemplary diagram for explaining a process of calculating a second boundary condition according to an embodiment of the present invention.
도 6 및 도 7을 참조하면, 경계 조건 계산부(230)는 제 1 가상 입자의 압력 및 외력에 기초하여 제 2 경계 조건(700)을 계산할 수 있다. 여기서, 제 2 경계 조건(700)으로 압력 노이만 경계 조건(Pressure Neumann boundary condition)이 이용될 수 있다. 압력 노이만 경계 조건은 벽면에서 물리량의 구배(gradient)를 줌으로써, 압력이 연속적으로 나타나게 한다. 기본적으로, 압력 구배()를 줌으로써, 압력이 연속적으로 나타나는데, 이를 통해 도메인 내에 있는 유체 입자의 물리량을 엄밀하게 계산함과 동시에 경계 밖으로 나가는 벽면 침투 현상을 방지할 수 있다. 6 and 7 , the
이를 위해, 제 1 가상 입자의 압력은 다음의 수학식 9를 이용하여 계산하고, 제 1 가상 입자의 외력은 다음의 수학식 10을 이용하여 계산하고, 제 2 경계 조건(700)은 다음의 수학식 11을 이용하여 계산할 수 있다.To this end, the pressure of the first virtual particle is calculated using the following Equation 9, the external force of the first virtual particle is calculated using the following Equation 10, and the
수학식 9를 참조하면, Pm은 제 1 가상 입자의 탐색 반경 내에 위치한 적어도 하나의 가상 마커의 압력의 평균을 의미하는 것으로, n은 제 1 가상 입자의 탐색 반경 내에 위치한 적어도 하나의 가상 마커의 수를 의미하고, 은 탐색 반경 내에 위치한 적어도 하나의 가상 마커의 물리량을 의미할 수 있다. Referring to Equation 9, Pm denotes an average of pressures of at least one virtual marker located within the search radius of the first virtual particle, and n is the number of at least one virtual marker located within the search radius of the first virtual particle. means, may mean a physical quantity of at least one virtual marker located within the search radius.
수학식 10을 참조하면, 가상 입자의 외력인 Fext는 Visc가 유체 입자의 점성에 의한 가속도 항을 의미하고, g가 중력가속도를 의미하고, n은 제 1 가상 입자의 탐색 반경 내에 위치한 적어도 하나의 가상 마커의 수를 의미할 수 있다.Referring to Equation 10, as for the external force of the virtual particle, F ext , Visc means an acceleration term due to the viscosity of the fluid particle, g means the gravitational acceleration, and n is at least one located within the search radius of the first virtual particle. may mean the number of virtual markers.
수학식 11을 참조하면, 경계에 대한 압력인 Pwall는 Pmarker가 가상 입자의 압력을 의미하고, d는 경계와 가상 입자 사이의 최단거리를 의미하고, Fext는 유체 입자의 외력을 의미할 수 있다. Referring to Equation 11, P wall , the pressure on the boundary, P marker means the pressure of the virtual particle, d means the shortest distance between the boundary and the virtual particle, and F ext is the external force of the fluid particle. can
유체 해석 시뮬레이션 수행부(240)는 복수의 유체 입자의 물리량 및 경계 조건에 기초하여 유체 해석 시뮬레이션을 수행할 수 있다. The fluid analysis
구체적으로, 유체 해석 시뮬레이션 수행부(240)는 SPH 알고리즘을 이용하여 각 입자와 이웃 입자 간의 충돌 또는 각 입자와 구조물 모델을 구성하는 다각형 간의 충돌로 인해 발생하는 유동 데이터를 계산하고, 유동 데이터에 기초하여 유체 해석 시뮬레이션을 수행할 수 있다.Specifically, the fluid analysis
SPH 알고리즘은 각 입자의 물성 정보(예컨대, 질량, 속도, 점성 및 가속도)를 이용하여 각 입자의 유동을 연산하는 것으로서, 각 입자의 물성 정보는 각 입자의 위치를 중심으로 한 방사형 베이시스 함수와 같은 커널 함수(Kernel Function) 세트를 사용하여 보간된다.The SPH algorithm calculates the flow of each particle by using the physical property information (eg, mass, velocity, viscosity, and acceleration) of each particle. It is interpolated using a set of kernel functions.
이러한 방식으로 각 입자의 물성 정보를 보간하면 나비에-스토크스 방정식과 같은 표준 방정식을 사용하여 유체의 동역학을 계산하는데 사용할 수 있는 압력 필드 및 점성 필드와 같은 연속 필드가 생성된다. Interpolating the physical property information of each particle in this way produces continuous fields such as pressure and viscosity fields that can be used to calculate the dynamics of a fluid using standard equations such as the Navier-Stokes equation.
예를 들어, 나비에-스토크스 방정식은 유체를 다음과 같이 모델링한다.For example, the Navier-Stokes equation models a fluid as
수학식 12에서 "v"는 입자의 속도, "ρ"는 입자의 밀도, "p"는 입자에 대한 압력, "g"는 중력, "μ"는 유체의 점성계수를 의미한다.In Equation 12, “v” is the velocity of the particles, “ρ” is the density of the particles, “p” is the pressure on the particles, “g” is the gravity, and “μ” is the viscosity coefficient of the fluid.
한편, SPH 알고리즘에 의하면, 각 입자의 밀도는 수학식 13에 의해 도출된다.Meanwhile, according to the SPH algorithm, the density of each particle is derived by Equation 13.
또한, 각 입자의 압력에 의한 힘은 수학식 14에 의해 도출된다.In addition, the force due to the pressure of each particle is derived by Equation (14).
또한, 각 입자의 점성에 의한 힘은 수학식 15에 의해 도출된다.In addition, the force due to the viscosity of each particle is derived by Equation (15).
유체 해석 시뮬레이션 수행부(240)는 SPH 알고리즘을 이용하여 각 입자의 밀도, 압력 및 점성 등의 유동 데이터의 변화값을 연산한다. 예를 들어, 유체 해석 시뮬레이션 수행부(240)는 각 입자의 초기 유동 데이터에 기초하여 다음 시간 스텝(제 1 시간 스텝)에서의 각 입자의 유동 데이터를 연산하고, 이에 기초하여 각 입자의 유동을 연산한다.The fluid analysis
또한, 유체 해석 시뮬레이션 수행부(240)는 제 1 시간 스텝에서의 각 입자의 유동 데이터에 기초하여 그 다음 시간 스텝에서의 각 입자의 유동 데이터를 연산하고, 이에 기초하여 각 입자의 유동을 연산한다.In addition, the fluid analysis
유체 해석 시뮬레이션 수행부(240)는 각 시간 스텝에서의 각 입자의 유동 데이터를 연산하여 각 입자의 유동을 연산함으로써, 유체 해석 시뮬레이션을 수행할 수 있다. The fluid analysis
도 8은 본 발명의 일 실시예에 따른 유체 해석 시뮬레이션 장치에서 유체 해석 시뮬레이션 방법의 순서도이다. 도 8에 도시된 일 실시예에 따른 입자 기반의 유체 해석 시뮬레이션 방법은 도 2에 도시된 유체 해석 시뮬레이션 장치에서 시계열적으로 처리되는 단계들을 포함한다. 따라서, 이하 생략된 내용이라고 하더라도 도 8에 도시된 실시예에 따라 수행되는 입자 기반의 유체 해석 시뮬레이션 방법에도 적용된다. 8 is a flowchart of a fluid analysis simulation method in the fluid analysis simulation apparatus according to an embodiment of the present invention. The particle-based fluid analysis simulation method according to the embodiment shown in FIG. 8 includes steps that are time-series processed by the fluid analysis simulation device shown in FIG. 2 . Therefore, even if omitted below, it is also applied to the particle-based fluid analysis simulation method performed according to the embodiment shown in FIG. 8 .
단계 S810에서 유체 해석 시뮬레이션 장치(1)는 유체에 대한 복수의 유체 입자를 시뮬레이션 영역 내에 배치하고, 경계 조건을 계산하기 위한 복수의 가상 입자를 시뮬레이션 영역 외에 배치하여 시뮬레이션 영역을 모델링할 수 있다. In operation S810 , the fluid analysis and
단계 S820에서 유체 해석 시뮬레이션 장치(1)는 복수의 유체 입자의 물리량에 기초하여 시뮬레이션 영역의 경계에 위치한 복수의 가상 마커의 물리량을 산출할 수 있다. In operation S820 , the fluid analysis and
단계 S830에서 유체 해석 시뮬레이션 장치(1)는 복수의 가상 마커의 물리량에 기초하여 경계 조건을 계산할 수 있다. In operation S830, the fluid analysis and
단계 S840에서 유체 해석 시뮬레이션 장치(1)는 복수의 유체 입자의 물리량 및 경계 조건에 기초하여 유체 해석 시뮬레이션을 수행할 수 있다. In operation S840, the fluid
상술한 설명에서, 단계 S810 내지 S840은 본 발명의 구현예에 따라서, 추가적인 단계들로 더 분할되거나, 더 적은 단계들로 조합될 수 있다. 또한, 일부 단계는 필요에 따라 생략될 수도 있고, 단계 간의 순서가 전환될 수도 있다.In the above description, steps S810 to S840 may be further divided into additional steps or combined into fewer steps, according to an embodiment of the present invention. In addition, some steps may be omitted as necessary, and the order between the steps may be switched.
도 8을 통해 설명된 유체 해석 시뮬레이션 방법은 매체에 저장된 컴퓨터 프로그램의 형태로 구현되거나, 컴퓨터에 의해 실행되는 프로그램 모듈과 같은 컴퓨터에 의해 실행 가능한 명령어를 포함하는 기록 매체의 형태로도 구현될 수 있다. 컴퓨터 판독 가능 매체는 컴퓨터에 의해 액세스될 수 있는 임의의 가용 매체일 수 있고, 휘발성 및 비휘발성 매체, 분리형 및 비분리형 매체를 모두 포함한다. 또한, 컴퓨터 판독가능 매체는 컴퓨터 저장 매체를 포함할 수 있다. 컴퓨터 저장 매체는 컴퓨터 판독가능 명령어, 데이터 구조, 프로그램 모듈 또는 기타 데이터와 같은 정보의 저장을 위한 임의의 방법 또는 기술로 구현된 휘발성 및 비휘발성, 분리형 및 비분리형 매체를 모두 포함한다.The fluid analysis simulation method described with reference to FIG. 8 may be implemented in the form of a computer program stored in the medium, or may be implemented in the form of a recording medium including instructions executable by a computer, such as a program module executed by a computer. . Computer-readable media can be any available media that can be accessed by a computer and includes both volatile and nonvolatile media, removable and non-removable media. Also, computer-readable media may include computer storage media. Computer storage media includes both volatile and nonvolatile, removable and non-removable media implemented in any method or technology for storage of information such as computer readable instructions, data structures, program modules or other data.
전술한 본 발명의 설명은 예시를 위한 것이며, 본 발명이 속하는 기술분야의 통상의 지식을 가진 자는 본 발명의 기술적 사상이나 필수적인 특징을 변경하지 않고서 다른 구체적인 형태로 쉽게 변형이 가능하다는 것을 이해할 수 있을 것이다. 그러므로 이상에서 기술한 실시예들은 모든 면에서 예시적인 것이며 한정적이 아닌 것으로 이해해야만 한다. 예를 들어, 단일형으로 설명되어 있는 각 구성 요소는 분산되어 실시될 수도 있으며, 마찬가지로 분산된 것으로 설명되어 있는 구성 요소들도 결합된 형태로 실시될 수 있다. The description of the present invention described above is for illustration, and those of ordinary skill in the art to which the present invention pertains can understand that it can be easily modified into other specific forms without changing the technical spirit or essential features of the present invention. will be. Therefore, it should be understood that the embodiments described above are illustrative in all respects and not restrictive. For example, each component described as a single type may be implemented in a dispersed form, and likewise components described as distributed may also be implemented in a combined form.
본 발명의 범위는 상기 상세한 설명보다는 후술하는 특허청구범위에 의하여 나타내어지며, 특허청구범위의 의미 및 범위 그리고 그 균등 개념으로부터 도출되는 모든 변경 또는 변형된 형태가 본 발명의 범위에 포함되는 것으로 해석되어야 한다.The scope of the present invention is indicated by the following claims rather than the above detailed description, and all changes or modifications derived from the meaning and scope of the claims and their equivalents should be interpreted as being included in the scope of the present invention. do.
1: 유체 해석 시뮬레이션 장치
210: 모델링부
220: 물리량 산출부
230: 경계 조건 계산부
235: 탐색부
240: 유체 시뮬레이션 수행부1: Fluid analysis simulation device
210: modeling unit
220: physical quantity calculator
230: boundary condition calculator
235: search unit
240: fluid simulation performing unit
Claims (18)
유체에 대한 복수의 유체 입자를 시뮬레이션 영역 내에 배치하고, 경계 조건을 계산하기 위한 복수의 가상 입자를 상기 시뮬레이션 영역 외에 배치하여 상기 시뮬레이션 영역을 모델링하는 단계;
상기 복수의 유체 입자의 물리량에 기초하여 상기 시뮬레이션 영역의 경계에 위치한 복수의 가상 마커의 물리량을 산출하는 단계;
상기 복수의 가상 마커의 물리량에 기초하여 상기 경계 조건을 계산하는 단계; 및
상기 복수의 유체 입자의 물리량 및 상기 경계 조건에 기초하여 유체 해석 시뮬레이션을 수행하는 단계를 포함하고,
상기 복수의 가상 입자의 물리량은 각 가상 입자로부터 기설정된 범위 내의 복수의 마커의 물리량의 평균에 기초한 것인, 유체 해석 시뮬레이션 방법.
In the particle-based fluid analysis simulation method performed in the fluid analysis simulation device,
modeling the simulation region by arranging a plurality of fluid particles for a fluid in a simulation region and arranging a plurality of virtual particles for calculating a boundary condition outside the simulation region;
calculating physical quantities of a plurality of virtual markers located at a boundary of the simulation area based on the physical quantities of the plurality of fluid particles;
calculating the boundary condition based on the physical quantities of the plurality of virtual markers; and
Comprising the step of performing a fluid analysis simulation based on the physical quantity of the plurality of fluid particles and the boundary condition,
wherein the physical quantity of the plurality of virtual particles is based on an average of physical quantities of a plurality of markers within a predetermined range from each virtual particle.
상기 복수의 가상 마커의 물리량을 산출하는 단계는,
상기 복수의 가상 마커 및 상기 복수의 유체 입자 간의 위치 벡터의 차를 산출하는 단계; 및
상기 산출된 위치 벡터의 차에 기초하여 상기 복수의 유체 입자에 대한 상기 복수의 가상 마커 각각의 속도 및 압력을 산출하는 단계를 포함하는 것인, 유체 해석 시뮬레이션 방법.
The method of claim 1,
Calculating the physical quantities of the plurality of virtual markers comprises:
calculating a difference in position vectors between the plurality of virtual markers and the plurality of fluid particles; and
Comprising the step of calculating the velocity and pressure of each of the plurality of virtual markers for the plurality of fluid particles based on the difference between the calculated position vectors, fluid analysis simulation method.
상기 경계 조건을 계산하는 단계는,
제 1 가상 입자로부터 소정의 탐색 반경 내에 위치한 적어도 하나의 인접 가상 마커를 탐색하는 단계; 및
상기 탐색된 적어도 하나의 인접 가상 마커의 상기 경계에 대한 법선 벡터를 도출하는 단계를 포함하는 것인, 유체 해석 시뮬레이션 방법.
3. The method of claim 2,
Calculating the boundary condition comprises:
searching for at least one adjacent virtual marker located within a predetermined search radius from the first virtual particle; and
and deriving a normal vector to the boundary of the found at least one adjacent virtual marker.
상기 적어도 하나의 인접 가상 마커를 탐색하는 단계는
상기 경계로부터 상기 제 1 가상 입자까지의 거리에 기초하여 상기 탐색 반경을 결정하는 단계
를 포함하는 것인, 유체 해석 시뮬레이션 방법.
4. The method of claim 3,
The step of searching for the at least one adjacent virtual marker includes:
determining the search radius based on a distance from the boundary to the first virtual particle;
Which will include, a fluid analysis simulation method.
상기 경계에 대한 법선 벡터를 도출하는 단계는,
상기 도출된 적어도 하나의 인접 가상 마커의 상기 경계에 대한 법선 벡터의 평균을 산출하는 단계; 및
상기 산출된 평균에 기초하여 상기 제 1 가상 입자의 상기 경계에 대한 단위 법선 벡터를 도출하는 단계를 포함하는 것인, 유체 해석 시뮬레이션 방법.
4. The method of claim 3,
The step of deriving a normal vector for the boundary comprises:
calculating an average of normal vectors with respect to the boundary of the derived at least one adjacent virtual marker; and
and deriving a unit normal vector for the boundary of the first virtual particle based on the calculated average.
상기 경계 조건을 계산하는 단계는,
상기 도출된 적어도 하나의 인접 가상 마커의 속도 및 압력의 평균을 산출하는 단계; 및
상기 산출된 평균에 기초하여 상기 제 1 가상 입자의 속도 및 압력를 도출하는 단계를 포함하는 것인, 유체 해석 시뮬레이션 방법.
6. The method of claim 5,
Calculating the boundary condition comprises:
calculating an average of the velocity and pressure of the derived at least one adjacent virtual marker; and
and deriving the velocity and pressure of the first virtual particle based on the calculated average.
상기 제 1 가상 입자의 상기 경계에 대한 단위 법선 벡터 및 상기 제 1 가상 입자의 속도에 기초하여 제 1 경계 조건을 계산하는 단계; 및
상기 제 1 가상 입자의 압력 및 외력에 기초하여 제 2 경계 조건을 계산하는 단계를 포함하는 것인, 유체 해석 시뮬레이션 방법.
7. The method of claim 6,
calculating a first boundary condition based on a unit normal vector of the first virtual particle to the boundary and a velocity of the first virtual particle; and
and calculating a second boundary condition based on the pressure and external force of the first virtual particle.
상기 유체에 대한 복수의 입자를 포함하는 시뮬레이션 영역을 모델링하는 단계는,
지형 정보, 구조물 정보, 경계 조건 정보, 입자 물성 정보 및 중력 가속도 정보 중 적어도 하나를 입력받는 단계; 및
상기 입력받은 적어도 하나에 기초하여 상기 시뮬레이션 영역을 모델링하는 단계
를 포함하는 것인, 유체 해석 시뮬레이션 방법.
The method of claim 1,
Modeling a simulation region including a plurality of particles for the fluid comprises:
receiving at least one of terrain information, structure information, boundary condition information, particle property information, and gravitational acceleration information; and
modeling the simulation region based on the received at least one
Which will include, a fluid analysis simulation method.
상기 구조물 정보는 밀도, 반발 계수 및 마찰 계수 중 적어도 하나를 포함하는 것인, 유체 해석 시뮬레이션 방법.
9. The method of claim 8,
The structure information will include at least one of a density, a coefficient of restitution, and a coefficient of friction, fluid analysis simulation method.
상기 입자 물성 정보는 입자 반경, 밀도, 점성, 음속 및 초기 속도 중 적어도 하나를 포함하는 것인, 유체 해석 시뮬레이션 방법.
9. The method of claim 8,
The particle property information will include at least one of a particle radius, density, viscosity, speed of sound, and initial velocity, fluid analysis simulation method.
유체에 대한 복수의 유체 입자를 시뮬레이션 영역 내에 배치하고, 경계 조건을 계산하기 위한 복수의 가상 입자를 상기 시뮬레이션 영역 외에 배치하여 상기 시뮬레이션 영역을 모델링하는 모델링부;
상기 복수의 유체 입자의 물리량에 기초하여 상기 시뮬레이션 영역의 경계에 위치한 복수의 가상 마커의 물리량을 산출하는 물리량 산출부;
상기 복수의 가상 마커의 물리량에 기초하여 상기 경계 조건을 계산하는 경계 조건 계산부; 및
상기 복수의 유체 입자의 물리량 및 상기 경계 조건에 기초하여 유체 해석 시뮬레이션을 수행하는 유체 시뮬레이션 수행부를 포함하고,
상기 복수의 가상 입자의 물리량은 각 가상 입자로부터 기설정된 범위 내의 복수의 마커의 물리량의 평균에 기초한 것인, 유체 해석 시뮬레이션 장치.
In the particle-based fluid simulation fluid analysis simulation device,
a modeling unit that models the simulation area by arranging a plurality of fluid particles for a fluid in a simulation area and arranging a plurality of virtual particles for calculating a boundary condition outside the simulation area;
a physical quantity calculator configured to calculate physical quantities of a plurality of virtual markers located at a boundary of the simulation area based on the physical quantities of the plurality of fluid particles;
a boundary condition calculator configured to calculate the boundary condition based on physical quantities of the plurality of virtual markers; and
A fluid simulation performing unit for performing a fluid analysis simulation based on the physical quantities of the plurality of fluid particles and the boundary conditions,
wherein the physical quantity of the plurality of virtual particles is based on an average of physical quantities of a plurality of markers within a predetermined range from each virtual particle.
상기 물리량 산출부는 상기 복수의 가상 마커 및 상기 복수의 유체 입자 간의 위치 벡터의 차를 산출하고, 상기 산출된 위치 벡터의 차에 기초하여 상기 복수의 유체 입자에 대한 상기 복수의 가상 마커 각각의 속도 및 압력을 산출하는 것인, 유체 해석 시뮬레이션 장치.
12. The method of claim 11,
The physical quantity calculator calculates a difference between the position vectors between the plurality of virtual markers and the plurality of fluid particles, and based on the calculated difference between the position vectors, the velocity of each of the plurality of virtual markers with respect to the plurality of fluid particles and A fluid analysis simulation device that calculates a pressure.
제 1 가상 입자로부터 소정의 탐색 반경 내에 위치한 적어도 하나의 인접 가상 마커를 탐색하는 인접 가상 마커 탐색부를 더 포함하고,
상기 경계 조건 계산부는 상기 탐색된 적어도 하나의 인접 가상 마커의 상기 경계에 대한 법선 벡터를 도출하는 것인, 유체 해석 시뮬레이션 장치.
13. The method of claim 12,
Further comprising: a neighbor virtual marker search unit for searching for at least one adjacent virtual marker located within a predetermined search radius from the first virtual particle,
The boundary condition calculation unit is to derive a normal vector for the boundary of the at least one searched adjacent virtual marker, the fluid analysis simulation apparatus.
상기 인접 가상 마커 탐색부는 상기 경계로부터 상기 제 1 가상 입자까지의 거리에 기초하여 상기 탐색 반경을 결정하는 것인, 유체 해석 시뮬레이션 장치.
14. The method of claim 13,
The adjacent virtual marker search unit will determine the search radius based on the distance from the boundary to the first virtual particle, fluid analysis simulation apparatus.
상기 경계 조건 계산부는 상기 도출된 적어도 하나의 인접 가상 마커의 상기 경계에 대한 법선 벡터의 평균을 산출하고, 상기 산출된 평균에 기초하여 상기 제 1 가상 입자의 상기 경계에 대한 단위 법선 벡터를 도출하는 것인, 유체 해석 시뮬레이션 장치.
14. The method of claim 13,
The boundary condition calculation unit calculates an average of a normal vector with respect to the boundary of the at least one adjacent virtual marker, and derives a unit normal vector for the boundary of the first virtual particle based on the calculated average that, a fluid analysis simulation device.
상기 경계 조건 계산부는 상기 도출된 적어도 하나의 인접 가상 마커의 속도 및 압력의 평균을 산출하고, 상기 산출된 평균에 기초하여 상기 제 1 가상 입자의 속도 및 압력를 도출하는 것인, 유체 해석 시뮬레이션 장치.
16. The method of claim 15,
The boundary condition calculator calculates an average of the derived velocity and pressure of the at least one adjacent virtual marker, and derives the velocity and pressure of the first virtual particle based on the calculated average.
상기 제 1 가상 입자의 상기 경계에 대한 단위 법선 벡터 및 상기 제 1 가상 입자의 속도에 기초하여 제 1 경계 조건을 계산하고, 상기 제 1 가상 입자의 압력 및 외력에 기초하여 제 2 경계 조건을 계산하는 것인, 유체 해석 시뮬레이션 장치.
17. The method of claim 16,
A first boundary condition is calculated based on a unit normal vector of the first virtual particle to the boundary and a velocity of the first virtual particle, and a second boundary condition is calculated based on a pressure and an external force of the first virtual particle. which is, a fluid analysis simulation device.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
KR1020200091710A KR102436658B1 (en) | 2020-07-23 | 2020-07-23 | Fluid analysis simulation method and fluid simulation apparatus |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
KR1020200091710A KR102436658B1 (en) | 2020-07-23 | 2020-07-23 | Fluid analysis simulation method and fluid simulation apparatus |
Publications (2)
Publication Number | Publication Date |
---|---|
KR20220012662A true KR20220012662A (en) | 2022-02-04 |
KR102436658B1 KR102436658B1 (en) | 2022-08-26 |
Family
ID=80268506
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
KR1020200091710A KR102436658B1 (en) | 2020-07-23 | 2020-07-23 | Fluid analysis simulation method and fluid simulation apparatus |
Country Status (1)
Country | Link |
---|---|
KR (1) | KR102436658B1 (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR102612846B1 (en) * | 2023-05-02 | 2023-12-12 | 씨넘버스 주식회사 | Flow analysis simulation apparatus, method and computer program based on boundary layer grid generation |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR100963114B1 (en) * | 2008-12-23 | 2010-06-15 | 고려대학교 산학협력단 | Apparatus and method for simulating fluid |
JP2011179835A (en) * | 2010-02-26 | 2011-09-15 | Toray Ind Inc | Fluid analysis method and fluid analysis apparatus |
KR20120100419A (en) * | 2011-03-04 | 2012-09-12 | (주)에프엑스기어 | System and method for fluid simulation using interaction between grid and particle |
JP2019070597A (en) * | 2017-10-10 | 2019-05-09 | 富士通株式会社 | Fluid simulation program, fluid simulation method and fluid simulation device |
KR102181986B1 (en) * | 2019-12-13 | 2020-11-23 | 이에이트 주식회사 | Dummy particle based fluid analysis simulation method and fluid simulation apparatus |
-
2020
- 2020-07-23 KR KR1020200091710A patent/KR102436658B1/en active IP Right Grant
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR100963114B1 (en) * | 2008-12-23 | 2010-06-15 | 고려대학교 산학협력단 | Apparatus and method for simulating fluid |
JP2011179835A (en) * | 2010-02-26 | 2011-09-15 | Toray Ind Inc | Fluid analysis method and fluid analysis apparatus |
KR20120100419A (en) * | 2011-03-04 | 2012-09-12 | (주)에프엑스기어 | System and method for fluid simulation using interaction between grid and particle |
JP2019070597A (en) * | 2017-10-10 | 2019-05-09 | 富士通株式会社 | Fluid simulation program, fluid simulation method and fluid simulation device |
JP6897477B2 (en) * | 2017-10-10 | 2021-06-30 | 富士通株式会社 | Fluid simulation program, fluid simulation method and fluid simulation equipment |
KR102181986B1 (en) * | 2019-12-13 | 2020-11-23 | 이에이트 주식회사 | Dummy particle based fluid analysis simulation method and fluid simulation apparatus |
Non-Patent Citations (1)
Title |
---|
Mitsuteru ASAI 외 3명, "Slip and No-Slip Boundary Treatment for Particle Simulation Model with Incompatible Step-Shaped Boundaries by Using a Virtual Maker",2013.12., pp1-13. 1부.* * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR102612846B1 (en) * | 2023-05-02 | 2023-12-12 | 씨넘버스 주식회사 | Flow analysis simulation apparatus, method and computer program based on boundary layer grid generation |
Also Published As
Publication number | Publication date |
---|---|
KR102436658B1 (en) | 2022-08-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
KR102181989B1 (en) | Particle based fluid simulation method using a plurality of processors and fluid simulation apparatus | |
Zheng et al. | Incompressible SPH method based on Rankine source solution for violent water wave simulation | |
Korzani et al. | Parametric study on smoothed particle hydrodynamics for accurate determination of drag coefficient for a circular cylinder | |
KR102181986B1 (en) | Dummy particle based fluid analysis simulation method and fluid simulation apparatus | |
KR20090027142A (en) | Method for constructing surface of fluid-body simulation based on particle method, program for the same, and storage medium for storing program | |
Manik et al. | A generic algorithm for three-dimensional multiphase flows on unstructured meshes | |
Evje et al. | Numerical treatment of two-phase flow in capillary heterogeneous porous media by finite-volume approximations | |
Diehl et al. | Bond-based peridynamics: a quantitative study of mode i crack opening | |
KR102026154B1 (en) | The method for numerical simulation of shallow water waves in shallow flows | |
KR102436658B1 (en) | Fluid analysis simulation method and fluid simulation apparatus | |
KR102181988B1 (en) | Sph based on fluid simulation method and fluid analysis simulation apparatus | |
KR102200443B1 (en) | Apparatus, method and computer program for lattice boltzmann method based fluid analysis simulation | |
Rasool et al. | A strategy to interface isogeometric analysis with Lagrangian finite elements—Application to incompressible flow problems | |
Dolejší et al. | An anisotropic hp-mesh adaptation method for time-dependent problems based on interpolation error control | |
Angst et al. | Robust and efficient wave simulations on deforming meshes | |
JP6163897B2 (en) | Numerical calculation program, numerical calculation method, and information processing apparatus | |
KR102436664B1 (en) | Apparatus, method and computer program for weakly compressible smoothed particle hydrodynamics based fluid analysis simulation | |
KR102181979B1 (en) | Apparatus, method and computer program for sph based fluid analysis simulation | |
KR102181981B1 (en) | Apparatus, method and computer program for sph based fluid analysis simulation | |
Kostorz et al. | Semi‐analytical treatment of spherical solids in smoothed‐particle hydrodynamics fluid simulations | |
KR102200444B1 (en) | Apparatus, method and computer program for sph based fluid analysis simulation | |
JP7395456B2 (en) | Simulation device and program | |
KR102181985B1 (en) | Apparatus, method and computer program for sph based fluid analysis simulation | |
Atallah et al. | Weak boundary conditions for Lagrangian shock hydrodynamics: A high-order finite element implementation on curved boundaries | |
JP2019101939A (en) | Method for analyzing motion state of fibrous material in fluid and analysis device therefor |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
E701 | Decision to grant or registration of patent right | ||
GRNT | Written decision to grant |