KR102640561B1 - Apparatus for simulating aerosol and method thereof - Google Patents
Apparatus for simulating aerosol and method thereof Download PDFInfo
- Publication number
- KR102640561B1 KR102640561B1 KR1020200142924A KR20200142924A KR102640561B1 KR 102640561 B1 KR102640561 B1 KR 102640561B1 KR 1020200142924 A KR1020200142924 A KR 1020200142924A KR 20200142924 A KR20200142924 A KR 20200142924A KR 102640561 B1 KR102640561 B1 KR 102640561B1
- Authority
- KR
- South Korea
- Prior art keywords
- aerosol
- simulation
- equation
- generating article
- calculation
- Prior art date
Links
- 239000000443 aerosol Substances 0.000 title claims abstract description 157
- 238000000034 method Methods 0.000 title claims abstract description 38
- 238000004088 simulation Methods 0.000 claims abstract description 207
- 238000012546 transfer Methods 0.000 claims abstract description 63
- 239000012530 fluid Substances 0.000 claims abstract description 25
- 238000004364 calculation method Methods 0.000 claims description 120
- 238000009826 distribution Methods 0.000 claims description 42
- 238000010438 heat treatment Methods 0.000 claims description 24
- 238000012800 visualization Methods 0.000 claims description 19
- 238000004590 computer program Methods 0.000 claims description 10
- 230000003252 repetitive effect Effects 0.000 claims description 3
- 230000004044 response Effects 0.000 claims description 3
- 238000012827 research and development Methods 0.000 abstract description 12
- 230000006870 function Effects 0.000 abstract description 9
- 239000000126 substance Substances 0.000 description 32
- 239000000463 material Substances 0.000 description 31
- 239000000758 substrate Substances 0.000 description 20
- 238000006243 chemical reaction Methods 0.000 description 18
- 230000008859 change Effects 0.000 description 17
- 230000000875 corresponding effect Effects 0.000 description 17
- 235000019504 cigarettes Nutrition 0.000 description 13
- 230000000391 smoking effect Effects 0.000 description 12
- 238000004891 communication Methods 0.000 description 10
- 238000001816 cooling Methods 0.000 description 10
- 241000208125 Nicotiana Species 0.000 description 8
- 235000002637 Nicotiana tabacum Nutrition 0.000 description 8
- 238000004519 manufacturing process Methods 0.000 description 8
- 238000013461 design Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 7
- 239000007788 liquid Substances 0.000 description 7
- 238000003860 storage Methods 0.000 description 7
- SNICXCGAKADSCV-JTQLQIEISA-N (-)-Nicotine Chemical compound CN1CCC[C@H]1C1=CC=CN=C1 SNICXCGAKADSCV-JTQLQIEISA-N 0.000 description 6
- 238000000889 atomisation Methods 0.000 description 6
- 238000001704 evaporation Methods 0.000 description 6
- 230000008020 evaporation Effects 0.000 description 6
- 229960002715 nicotine Drugs 0.000 description 6
- SNICXCGAKADSCV-UHFFFAOYSA-N nicotine Natural products CN1CCCC1C1=CC=CN=C1 SNICXCGAKADSCV-UHFFFAOYSA-N 0.000 description 6
- 238000009792 diffusion process Methods 0.000 description 5
- 230000004907 flux Effects 0.000 description 5
- 210000000214 mouth Anatomy 0.000 description 5
- 230000035699 permeability Effects 0.000 description 5
- 238000011144 upstream manufacturing Methods 0.000 description 5
- 238000004134 energy conservation Methods 0.000 description 4
- 239000000203 mixture Substances 0.000 description 4
- 239000012071 phase Substances 0.000 description 4
- 230000000704 physical effect Effects 0.000 description 4
- DNIAPMSPPWPWGF-UHFFFAOYSA-N Propylene glycol Chemical compound CC(O)CO DNIAPMSPPWPWGF-UHFFFAOYSA-N 0.000 description 3
- 229920002301 cellulose acetate Polymers 0.000 description 3
- 230000007423 decrease Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 239000007791 liquid phase Substances 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 239000007787 solid Substances 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- PEDCQBHIVMGVHV-UHFFFAOYSA-N Glycerine Chemical compound OCC(O)CO PEDCQBHIVMGVHV-UHFFFAOYSA-N 0.000 description 2
- 239000004909 Moisturizer Substances 0.000 description 2
- 230000004913 activation Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 239000000470 constituent Substances 0.000 description 2
- 230000002596 correlated effect Effects 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 239000003906 humectant Substances 0.000 description 2
- 210000004072 lung Anatomy 0.000 description 2
- 230000001333 moisturizer Effects 0.000 description 2
- 239000011148 porous material Substances 0.000 description 2
- 238000000926 separation method Methods 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 240000003483 Leersia hexandra Species 0.000 description 1
- 239000003795 chemical substances by application Substances 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 239000003571 electronic cigarette Substances 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 239000000796 flavoring agent Substances 0.000 description 1
- 235000019634 flavors Nutrition 0.000 description 1
- 229940055390 glycerin / propylene glycol Drugs 0.000 description 1
- 235000011187 glycerol Nutrition 0.000 description 1
- 230000006698 induction Effects 0.000 description 1
- 238000009413 insulation Methods 0.000 description 1
- 210000003928 nasal cavity Anatomy 0.000 description 1
- 210000001331 nose Anatomy 0.000 description 1
- 230000001151 other effect Effects 0.000 description 1
- 239000002994 raw material Substances 0.000 description 1
- 241000894007 species Species 0.000 description 1
- 238000000003 thermogravimetry coupled to Fourier transform infrared spectroscopy Methods 0.000 description 1
- 238000009834 vaporization Methods 0.000 description 1
- 230000008016 vaporization Effects 0.000 description 1
- 238000007794 visualization technique Methods 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A24—TOBACCO; CIGARS; CIGARETTES; SIMULATED SMOKING DEVICES; SMOKERS' REQUISITES
- A24F—SMOKERS' REQUISITES; MATCH BOXES; SIMULATED SMOKING DEVICES
- A24F40/00—Electrically operated smoking devices; Component parts thereof; Manufacture thereof; Maintenance or testing thereof; Charging means specially adapted therefor
- A24F40/50—Control or monitoring
-
- A—HUMAN NECESSITIES
- A24—TOBACCO; CIGARS; CIGARETTES; SIMULATED SMOKING DEVICES; SMOKERS' REQUISITES
- A24F—SMOKERS' REQUISITES; MATCH BOXES; SIMULATED SMOKING DEVICES
- A24F40/00—Electrically operated smoking devices; Component parts thereof; Manufacture thereof; Maintenance or testing thereof; Charging means specially adapted therefor
- A24F40/60—Devices with integrated user interfaces
-
- A—HUMAN NECESSITIES
- A24—TOBACCO; CIGARS; CIGARETTES; SIMULATED SMOKING DEVICES; SMOKERS' REQUISITES
- A24F—SMOKERS' REQUISITES; MATCH BOXES; SIMULATED SMOKING DEVICES
- A24F40/00—Electrically operated smoking devices; Component parts thereof; Manufacture thereof; Maintenance or testing thereof; Charging means specially adapted therefor
- A24F40/65—Devices with integrated communication means, e.g. wireless communication means
-
- 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]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Optimization (AREA)
- Computer Hardware Design (AREA)
- Computing Systems (AREA)
- Mathematical Analysis (AREA)
- Algebra (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Fluid Mechanics (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Human Computer Interaction (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
에어로졸 시뮬레이션 장치 및 그 방법이 제공된다. 본 개시의 몇몇 실시예들에 따른 에어로졸 시뮬레이션 장치는, 타깃 에어로졸 발생 물품의 구조데이터를 포함하는 시뮬레이션 기초데이터를 입력받는 데이터 입력부, 시뮬레이션 기초데이터에 기반하여 CFD(Computational Fluid Dynamics) 모델을 구성하는 모델 구성부 및 구성된 CFD 모델을 이용하여 타깃 에어로졸 발생 물품에서 나타나는 열전달 현상과 에어로졸의 유동 현상을 시뮬레이션하는 시뮬레이션부를 포함할 수 있다. 이러한 에어로졸 시뮬레이션 장치는 CFD 모델을 통해 물품 내부의 열전달 현상과 에어로졸 유동 현상에 대한 시뮬레이션 기능을 제공함으로써, 에어로졸 발생 물품의 연구 개발에 투입되는 인적 비용과 시간 비용을 절감시킬 수 있다.An aerosol simulation device and method are provided. An aerosol simulation device according to some embodiments of the present disclosure includes a data input unit that receives simulation basic data including structural data of a target aerosol-generating article, and a model that constructs a CFD (Computational Fluid Dynamics) model based on the simulation basic data. It may include a simulation unit that simulates the heat transfer phenomenon and aerosol flow phenomenon that appear in the target aerosol-generating article using the configuration unit and the configured CFD model. This aerosol simulation device can reduce human and time costs invested in research and development of aerosol-generating products by providing simulation functions for heat transfer phenomena and aerosol flow phenomena inside the product through a CFD model.
Description
본 개시는 에어로졸 시뮬레이션 장치 및 그 방법에 관한 것이다. 보다 자세하게는, 에어로졸 발생 물품 내에서 발생되는 에어로졸 현상을 시뮬레이션할 수 있는 장치 및 그 장치에서 수행되는 방법에 관한 것이다.The present disclosure relates to an aerosol simulation device and method. More specifically, it relates to a device capable of simulating aerosol phenomena occurring within an aerosol-generating article and a method performed in the device.
근래에 일반적인 궐련의 단점들을 극복하는 대체 방법에 관한 수요가 증가하고 있다. 예를 들어, 궐련을 전기적으로 가열함으로써 에어로졸을 발생시키는 장치(e.g. 전자담배 디바이스)와 그 장치에 적용되는 새로운 구조의 궐련 스틱에 관한 수요가 증가하고 있다. 이에 따라, 궐련형 에어로졸 발생 장치와 궐련 스틱에 대한 연구 개발이 활발하게 이루어지고 있다.In recent years, there has been an increasing demand for alternative methods to overcome the disadvantages of regular cigarettes. For example, demand is increasing for devices that generate aerosols by electrically heating cigarettes (e.g. electronic cigarette devices) and cigarette sticks with new structures applied to the devices. Accordingly, research and development on cigarette-type aerosol generating devices and cigarette sticks is being actively conducted.
통상적으로, 궐련 스틱에 대한 연구 개발은 요구사항 설정, 설계, 제조 및 검증 순으로 진행된다. 가령, 궐련 스틱의 무화량 개선을 위한 연구 개발의 경우에는, 목표 무화량 설정, 궐련 스틱의 구조 설계, 제조 및 흡연 실험을 통한 무화량 검증 순서로 연구 개발이 진행될 수 있다. 그런데, 이러한 연구 개발 프로세스는 때때로 반복적인 제조와 실제 실험을 요구하여, 연구 개발에 투입되는 인적 비용과 시간 비용이 상당히 증가시킬 수 있다. 가령, 설계 상의 결함이 마지막 검증 단계에서 비로소 발견됨으로써, 다시 궐련 스틱을 설계하고 제조하는 일이 빈번하게 발생될 수 있다.Typically, research and development for cigarette sticks proceeds in the following order: requirement setting, design, manufacturing and verification. For example, in the case of research and development to improve the atomization amount of cigarette sticks, research and development may be conducted in the following order: setting a target atomization amount, designing the structure of the cigarette stick, manufacturing, and verifying the atomization amount through smoking experiments. However, this research and development process sometimes requires repetitive manufacturing and actual experiments, which can significantly increase the human and time costs invested in research and development. For example, defects in the design are discovered only in the final verification stage, which may frequently result in the design and manufacture of cigarette sticks again.
따라서, 연구 개발에 투입되는 인적 비용 및 시간 비용을 절감하고, 연구 개발의 효율성을 강화하기 위해서는, 궐련 스틱 내의 다양한 에어로졸 현상을 시뮬레이션할 수 있는 방법이 절실히 요구된다.Therefore, in order to reduce human and time costs invested in research and development and enhance the efficiency of research and development, a method that can simulate various aerosol phenomena within cigarette sticks is urgently needed.
본 개시의 몇몇 실시예들을 통해 해결하고자 하는 기술적 과제는, 에어로졸 발생 물품에서 나타나는 열전달 현상과 에어로졸의 유동 현상을 시뮬레이션할 수 있는 장치 및 그 장치에서 수행되는 방법을 제공하는 것이다.The technical problem to be solved through several embodiments of the present disclosure is to provide a device that can simulate the heat transfer phenomenon and aerosol flow phenomenon that appear in an aerosol-generating article and a method performed in the device.
본 개시의 기술적 과제들은 이상에서 언급한 기술적 과제들로 제한되지 않으며, 언급되지 않은 또 다른 기술적 과제들은 아래의 기재로부터 본 개시의 기술분야에서의 통상의 기술자에게 명확하게 이해될 수 있을 것이다.The technical problems of the present disclosure are not limited to the technical problems mentioned above, and other technical problems not mentioned can be clearly understood by those skilled in the art from the description below.
상기 기술적 과제를 해결하기 위한, 본 개시의 몇몇 실시예들에 따른 에어로졸 시뮬레이션 장치는, 타깃 에어로졸 발생 물품의 구조데이터를 포함하는 시뮬레이션 기초데이터를 입력받는 데이터 입력부, 상기 시뮬레이션 기초데이터에 기반하여 CFD(Computational Fluid Dynamics) 모델을 구성하는 모델 구성부 및 상기 구성된 CFD 모델을 이용하여 상기 타깃 에어로졸 발생 물품에서 나타나는 열전달 현상과 에어로졸의 유동 현상을 시뮬레이션하는 시뮬레이션부를 포함할 수 있다.In order to solve the above technical problem, an aerosol simulation device according to some embodiments of the present disclosure includes a data input unit that receives simulation basic data including structural data of a target aerosol-generating article, and a CFD (CFD) based on the simulation basic data. It may include a model component that configures a Computational Fluid Dynamics (Computational Fluid Dynamics) model and a simulation unit that simulates the heat transfer phenomenon and aerosol flow phenomenon that appear in the target aerosol-generating article using the configured CFD model.
몇몇 실시예들에서, 상기 모델 구성부는, 상기 구조데이터에 기반하여 계산격자를 생성하는 격자 생성부, 상기 CFD 모델의 지배방정식을 설정하는 지배방정식 설정부 및 상기 CFD 모델의 초기조건 및 경계조건을 설정하는 조건 설정부를 포함할 수 있다.In some embodiments, the model constructor includes a grid generator that generates a calculation grid based on the structural data, a governing equation setting unit that sets the governing equation of the CFD model, and initial conditions and boundary conditions of the CFD model. It may include a condition setting unit to set.
몇몇 실시예들에서, 상기 조건 설정부는 상기 구조데이터에 기반하여 상기 에어로졸의 유동과 관련된 벽면의 경계조건을 설정할 수 있다.In some embodiments, the condition setting unit may set boundary conditions of the wall related to the flow of the aerosol based on the structural data.
몇몇 실시예들에서, 상기 시뮬레이션 기초데이터는 퍼프데이터를 더 포함하고, 상기 조건 설정부는 상기 퍼프데이터에 기반하여 상기 에어로졸에 대한 유동경계조건을 설정하되, 상기 유동경계조건은 상기 타깃 에어로졸 발생 물품으로부터 상기 에어로졸이 토출되는 하류 말단 영역에 설정되는 경계조건일 수 있다.In some embodiments, the simulation basic data further includes puff data, and the condition setting unit sets a flow boundary condition for the aerosol based on the puff data, wherein the flow boundary condition is set from the target aerosol-generating article. It may be a boundary condition set at the downstream end area where the aerosol is discharged.
몇몇 실시예들에서, 상기 시뮬레이션 기초데이터는 상기 타깃 에어로졸 발생 물품을 가열하는 히터의 온도 프로파일을 더 포함하고, 상기 조건 설정부는 상기 온도 프로파일에 기반하여 열 경계조건을 설정할 수 있다.In some embodiments, the simulation basic data further includes a temperature profile of a heater that heats the target aerosol-generating article, and the condition setting unit may set a thermal boundary condition based on the temperature profile.
몇몇 실시예들에서, 상기 시뮬레이션 기초데이터는 상기 히터의 가열 구조에 관한 데이터를 더 포함하고, 상기 조건 설정부는 상기 히터의 가열 구조에 기반하여 상기 히터로부터 타깃 에어로졸 발생 물품으로 열이 전달되는 경계영역을 결정하고, 상기 결정된 경계영역에 상기 열 경계조건을 설정할 수 있다.In some embodiments, the simulation basic data further includes data regarding the heating structure of the heater, and the condition setting unit is a boundary area where heat is transferred from the heater to the target aerosol-generating article based on the heating structure of the heater. can be determined, and the thermal boundary condition can be set in the determined boundary area.
몇몇 실시예들에서, 상기 지배방정식 설정부는 제1 지배방정식 및 제2 지배방정식을 상기 CFD 모델의 지배방정식으로 설정하고, 상기 시뮬레이션부는 상기 설정된 제1 지배방정식 및 제2 지배방정식을 연계하여 계산함으로써 상기 에어로졸의 유동 현상과 상기 열전달 현상을 시뮬레이션하되, 상기 제1 지배방정식은 상기 에어로졸의 유동에 관한 방정식이고, 상기 제2 지배방정식은 상기 타깃 에어로졸 발생 물품 내의 열전달에 관한 방정식일 수 있다.In some embodiments, the governing equation setting unit sets the first governing equation and the second governing equation as the governing equation of the CFD model, and the simulation unit calculates the set first governing equation and the second governing equation by linking them. Simulating the flow phenomenon of the aerosol and the heat transfer phenomenon, the first governing equation may be an equation related to the flow of the aerosol, and the second governing equation may be an equation related to heat transfer within the target aerosol-generating article.
몇몇 실시예들에서, 상기 시뮬레이션부는, 상기 생성된 계산격자 내에서 상기 설정된 초기조건 및 상기 설정된 경계조건에 따라 상기 설정된 지배방정식을 계산함으로써 상기 타깃 에어로졸 발생 물품 내의 온도 분포를 도출하고, 상기 도출 결과에 기반하여 온도가 기준치 이하인 계산격자 영역을 검출할 수 있다.In some embodiments, the simulation unit derives a temperature distribution within the target aerosol-generating article by calculating the set governing equation according to the set initial conditions and the set boundary conditions within the generated calculation grid, and the derived result. Based on this, areas of the calculation grid where the temperature is below the standard value can be detected.
몇몇 실시예들에서, 상기 시뮬레이션의 결과를 시각화하는 시각화부를 더 포함할 수 있다.In some embodiments, a visualization unit may be further included to visualize the results of the simulation.
상술한 기술적 과제를 해결하기 위한 본 개시의 몇몇 실시예들에 따른 에어로졸 시뮬레이션 방법은, 컴퓨팅 장치에서 수행되는 에어로졸 시뮬레이션 방법에 있어서, 타깃 에어로졸 발생 물품의 구조 데이터를 포함하는 시뮬레이션 기초데이터를 입력받는 단계, 상기 시뮬레이션 기초데이터에 기반하여 CFD(Computational Fluid Dynamics) 모델을 구성하는 단계 및 상기 구성된 CFD 모델을 이용하여 상기 타깃 에어로졸 발생 물품에서 나타나는 열전달 현상과 에어로졸의 유동 현상을 시뮬레이션하는 단계를 포함할 수 있다.The aerosol simulation method according to some embodiments of the present disclosure for solving the above-mentioned technical problems includes the step of receiving basic simulation data including structural data of a target aerosol-generating article in an aerosol simulation method performed on a computing device. , configuring a CFD (Computational Fluid Dynamics) model based on the simulation basic data, and simulating the heat transfer phenomenon and aerosol flow phenomenon that appear in the target aerosol-generating article using the constructed CFD model. .
상술한 기술적 과제를 해결하기 위한 본 개시의 몇몇 실시예들에 따른 컴퓨터 프로그램은, 컴퓨팅 장치와 결합되어, 타깃 에어로졸 발생 물품의 구조 데이터를 포함하는 시뮬레이션 기초데이터를 입력받는 단계, 상기 시뮬레이션 기초데이터에 기반하여 CFD(Computational Fluid Dynamics) 모델을 구성하는 단계 및 상기 구성된 CFD 모델을 이용하여 상기 타깃 에어로졸 발생 물품에서 나타나는 열전달 현상과 에어로졸의 유동 현상을 시뮬레이션하는 단계를 실행시키기 위하여 컴퓨터 판독가능한 기록매체에 저장될 수 있다.A computer program according to some embodiments of the present disclosure for solving the above-mentioned technical problem includes the steps of combining with a computing device and receiving simulation basic data including structural data of a target aerosol-generating article, and adding the simulation basic data to the simulation basic data. Storing in a computer-readable recording medium to execute the steps of constructing a CFD (Computational Fluid Dynamics) model based on the CFD model and simulating the heat transfer phenomenon and aerosol flow phenomenon that appear in the target aerosol-generating article using the constructed CFD model. It can be.
상술한 본 개시의 몇몇 실시예들에 따르면, CFD 모델을 통해 타깃 에어로졸 발생 물품에서 발생되는 에어로졸 현상이 시뮬레이션될 수 있다. 이에 따라, 타깃 에어로졸 발생 물품에 대한 제조 및 실험 없이도 설계 결과가 용이하게 검증될 수 있으며, 에어로졸 발생 물품의 연구 개발에 투입되는 인적 비용 및 시간 비용이 크게 절감될 수 있다.According to some embodiments of the present disclosure described above, an aerosol phenomenon generated from a target aerosol-generating article can be simulated through a CFD model. Accordingly, design results can be easily verified without manufacturing and testing target aerosol-generating products, and human and time costs invested in research and development of aerosol-generating products can be greatly reduced.
또한, 에어로졸 유동에 관한 지배방정식과 열전달에 관한 지배방정식을 연계하여 풀이함으로써, 타깃 에어로졸 발생 물품에서 나타나는 열전달 현상과 에어로졸 유동 현상이 정확하게 시뮬레이션될 수 있다.In addition, by linking and solving the governing equations for aerosol flow and the governing equations for heat transfer, the heat transfer phenomenon and aerosol flow phenomenon that appear in the target aerosol-generating article can be accurately simulated.
또한, 에어로졸의 유동, 열전달, 화학반응 및 물질이송에 관한 복수의 지배방정식을 연계하여 풀이함으로써, 에어로졸 발생 물품 내의 다양한 에어로졸 현상이 정확하게 시뮬레이션될 수 있다.Additionally, by linking and solving a plurality of governing equations for aerosol flow, heat transfer, chemical reaction, and material transfer, various aerosol phenomena within aerosol-generating articles can be accurately simulated.
또한, 시뮬레이션 변수값의 변동 정도에 기초하여 계산격자의 조밀도가 동적으로 조절될 수 있다. 가령, 시뮬레이션 변수값의 변동 정도가 큰 계산격자 영역의 조밀도는 증가하고, 변동 정도가 작은 계산격자 영역의 조밀도는 감소될 수 있다, 이에 따라, 중요하지 않은 계산격자 영역에 투입되는 컴퓨팅 비용이 절감될 수 있고, 중요한 계산격자 영역에 대해서는 시뮬레이션의 정확도가 향상될 수 있다.Additionally, the density of the computational grid can be dynamically adjusted based on the degree of variation in simulation variable values. For example, the density of calculation grid areas where the degree of variation in simulation variable values is large may increase, and the density of calculation grid areas where the degree of variation is small may be reduced. Accordingly, the computing cost invested in unimportant calculation grid areas may be increased. This can be reduced, and the accuracy of simulation can be improved for important computational grid areas.
본 개시의 기술적 사상에 따른 효과들은 이상에서 언급한 효과들로 제한되지 않으며, 언급되지 않은 또 다른 효과들은 아래의 기재로부터 통상의 기술자에게 명확하게 이해될 수 있을 것이다.The effects according to the technical idea of the present disclosure are not limited to the effects mentioned above, and other effects not mentioned can be clearly understood by those skilled in the art from the description below.
도 1은 본 개시의 몇몇 실시예들에 따른 에어로졸 시뮬레이션 장치를 나타내는 예시적인 구성도이다.
도 2는 본 개시의 몇몇 실시예들에 따른 격자 생성부의 동작을 설명하기 위한 예시적인 도면이다.
도 3은 본 개시의 몇몇 실시예들에 따른 지배방정식 설정부의 동작과 지배방정식 간의 연계 관계를 설명하기 위한 예시적인 도면이다.
도 4 및 도 5는 본 개시의 몇몇 실시예들에 따른 조건 설정부의 동작을 설명하기 위한 예시적인 도면이다.
도 6 및 도 7은 본 개시의 몇몇 실시예들에 따른 시뮬레이션부의 동작을 설명하기 위한 예시적인 도면이다.
도 8 내지 도 10은 본 개시의 몇몇 실시예들에 따른 시각화부에 의해 시각화된 시뮬레이션 결과를 예시한다.
도 11은 본 개시의 몇몇 실시예들에 따른 에어로졸 시뮬레이션 방법을 나타내는 예시적인 흐름도이다.
도 12는 본 개시의 몇몇 실시예들에 따른 시뮬레이션 장치를 구현할 수 있는 예시적인 컴퓨팅 장치를 도시한다.1 is an exemplary configuration diagram showing an aerosol simulation device according to some embodiments of the present disclosure.
FIG. 2 is an exemplary diagram for explaining the operation of a grid generator according to some embodiments of the present disclosure.
FIG. 3 is an exemplary diagram illustrating the linkage between the operation of the governing equation setting unit and the governing equations according to some embodiments of the present disclosure.
4 and 5 are exemplary diagrams for explaining the operation of a condition setting unit according to some embodiments of the present disclosure.
FIGS. 6 and 7 are exemplary diagrams for explaining the operation of a simulation unit according to some embodiments of the present disclosure.
8 to 10 illustrate simulation results visualized by a visualization unit according to some embodiments of the present disclosure.
11 is an example flowchart illustrating an aerosol simulation method according to some embodiments of the present disclosure.
12 illustrates an example computing device that can implement a simulation device according to some embodiments of the present disclosure.
이하, 첨부된 도면을 참조하여 본 개시의 바람직한 실시예들을 상세히 설명한다. 본 개시의 이점 및 특징, 그리고 그것들을 달성하는 방법은 첨부되는 도면과 함께 상세하게 후술되어 있는 실시예들을 참조하면 명확해질 것이다. 그러나 본 개시의 기술적 사상은 이하의 실시예들에 한정되는 것이 아니라 서로 다른 다양한 형태로 구현될 수 있으며, 단지 이하의 실시예들은 본 개시의 기술적 사상을 완전하도록 하고, 본 개시가 속하는 기술분야에서 통상의 지식을 가진 자에게 본 개시의 범주를 완전하게 알려주기 위해 제공되는 것이며, 본 개시의 기술적 사상은 청구항의 범주에 의해 정의될 뿐이다.Hereinafter, preferred embodiments of the present disclosure will be described in detail with reference to the attached drawings. The advantages and features of the present disclosure and methods for achieving them will become clear by referring to the embodiments described in detail below along with the accompanying drawings. However, the technical idea of the present disclosure is not limited to the following embodiments and may be implemented in various different forms. The following examples are merely intended to complete the technical idea of the present disclosure and to be used in the technical field to which the present disclosure belongs. It is provided to fully inform those skilled in the art of the scope of the present disclosure, and the technical idea of the present disclosure is only defined by the scope of the claims.
각 도면의 구성요소들에 참조부호를 부가함에 있어서, 동일한 구성요소들에 대해서는 비록 다른 도면상에 표시되더라도 가능한 한 동일한 부호를 가지도록 하고 있음에 유의해야 한다. 또한, 본 개시를 설명함에 있어, 관련된 공지 구성 또는 기능에 대한 구체적인 설명이 본 개시의 요지를 흐릴 수 있다고 판단되는 경우에는 그 상세한 설명은 생략한다.When adding reference numerals to components in each drawing, it should be noted that identical components are given the same reference numerals as much as possible even if they are shown in different drawings. Additionally, in describing the present disclosure, if it is determined that a detailed description of a related known configuration or function may obscure the gist of the present disclosure, the detailed description will be omitted.
다른 정의가 없다면, 본 명세서에서 사용되는 모든 용어(기술 및 과학적 용어를 포함)는 본 개시가 속하는 기술분야에서 통상의 지식을 가진 자에게 공통적으로 이해될 수 있는 의미로 사용될 수 있다. 또 일반적으로 사용되는 사전에 정의되어 있는 용어들은 명백하게 특별히 정의되어 있지 않는 한 이상적으로 또는 과도하게 해석되지 않는다. 본 명세서에서 사용된 용어는 실시예들을 설명하기 위한 것이며 본 개시를 제한하고자 하는 것은 아니다. 본 명세서에서, 단수형은 문구에서 특별히 언급하지 않는 한 복수형도 포함한다.Unless otherwise defined, all terms (including technical and scientific terms) used in this specification may be used with meanings that can be commonly understood by those skilled in the art to which this disclosure pertains. Additionally, terms defined in commonly used dictionaries are not interpreted ideally or excessively unless clearly specifically defined. The terminology used herein is for the purpose of describing embodiments and is not intended to limit the disclosure. As used herein, singular forms also include plural forms, unless specifically stated otherwise in the context.
또한, 본 개시의 구성요소를 설명하는 데 있어서, 제1, 제2, A, B, (a), (b) 등의 용어를 사용할 수 있다. 이러한 용어는 그 구성요소를 다른 구성요소와 구별하기 위한 것일 뿐, 그 용어에 의해 해당 구성요소의 본질이나 차례 또는 순서 등이 한정되지 않는다. 어떤 구성요소가 다른 구성요소에 "연결", "결합" 또는 "접속"된다고 기재된 경우, 그 구성요소는 그 다른 구성요소에 직접적으로 연결되거나 또는 접속될 수 있지만, 각 구성요소 사이에 또 다른 구성요소가 "연결", "결합" 또는 "접속"될 수도 있다고 이해되어야 할 것이다.Additionally, in describing the components of the present disclosure, terms such as first, second, A, B, (a), and (b) may be used. These terms are only used to distinguish the component from other components, and the nature, order, or order of the component is not limited by the term. When a component is described as being “connected,” “coupled,” or “connected” to another component, that component may be directly connected or connected to that other component, but there is another component between each component. It will be understood that elements may be “connected,” “combined,” or “connected.”
본 개시에서 사용되는 "포함한다(comprises)" 및/또는 "포함하는(comprising)"은 언급된 구성요소, 단계, 동작 및/또는 소자는 하나 이상의 다른 구성요소, 단계, 동작 및/또는 소자의 존재 또는 추가를 배제하지 않는다.As used in this disclosure, “comprises” and/or “comprising” refers to a referenced component, step, operation and/or element that includes one or more other components, steps, operations and/or elements. Does not exclude presence or addition.
본 개시의 다양한 실시예들에 대한 설명에 앞서, 실시예들에서 사용되는 몇몇 용어들에 대하여 명확하게 하기로 한다.Before describing various embodiments of the present disclosure, some terms used in the embodiments will be clarified.
이하의 실시예들에서, "에어로졸 형성 기재"는 에어로졸(aerosol)을 형성할 수 있는 물질을 의미할 수 있다. 에어로졸은 휘발성 화합물을 포함할 수 있다. 에어로졸 형성 기재는 고체 또는 액상일 수 있다. 예를 들면, 고체의 에어로졸 형성 기재는 판상엽 담배, 각초, 재구성 담배 등 담배 원료를 기초로 하는 고체 물질을 포함할 수 있으며, 액상의 에어로졸 형성 기재는 니코틴, 담배 추출물, 보습제 및/또는 다양한 향미제를 기초로 하는 액상 조성물을 포함할 수 있다. 그러나, 본 개시의 범위가 상기 열거된 예시에 한정되는 것은 아니다.In the following examples, “aerosol-forming substrate” may mean a material capable of forming an aerosol. Aerosols may contain volatile compounds. The aerosol-forming substrate may be solid or liquid. For example, the solid aerosol-forming substrate may include solid substances based on tobacco raw materials such as leaf tobacco, cut tobacco, reconstituted tobacco, etc., and the liquid aerosol-forming substrate may include nicotine, tobacco extract, humectants, and/or various flavors. It may include a liquid composition based on an agent. However, the scope of the present disclosure is not limited to the examples listed above.
이하의 실시예들에서, "에어로졸 발생 물품"은 에어로졸을 발생시킬 수 있는 물품(article)을 의미할 수 있다. 에어로졸 발생 물품은 에어로졸 형성 기재를 포함할 수 있다. 에어로졸 발생 물품의 대표적인 예로는 궐련 스틱을 들 수 있을 것이나, 본 개시의 범위가 이러한 예시에 한정되는 것은 아니다.In the following examples, “aerosol-generating article” may refer to an article capable of generating an aerosol. Aerosol-generating articles can include an aerosol-forming substrate. Representative examples of aerosol-generating articles include cigarette sticks, but the scope of the present disclosure is not limited to these examples.
이하의 실시예들에서, "에어로졸 발생 장치"는 사용자의 입을 통해 사용자의 폐로 직접적으로 흡입 가능한 에어로졸을 발생시키기 위해 에어로졸 형성 기재 또는 에어로졸 발생 물품을 이용하여 에어로졸을 발생시키는 장치를 의미할 수 있다. 에어로졸 발생 장치는 예를 들어 궐련 스틱을 통해 에어로졸을 발생시키는 장치를 포함할 수 있으나, 이에 한정되는 것은 아니다.In the following embodiments, “aerosol-generating device” may refer to a device that generates an aerosol using an aerosol-forming substrate or aerosol-generating article to generate an aerosol that can be inhaled directly through the user's mouth and into the user's lungs. The aerosol generating device may include, but is not limited to, a device that generates an aerosol through, for example, a cigarette stick.
이하의 실시예들에서, "상류"(upstream) 또는 "상류 방향"은 흡연자의 구부로부터 멀어지는 방향을 의미하고, "하류"(downstream) 또는 "하류 방향"은 흡연자의 구부로부터 가까워지는 방향을 의미할 수 있다. 상류 및 하류라는 용어는 흡연 물품을 구성하는 요소들의 상대적 위치를 설명하기 위해 이용될 수 있다. 예를 들어, 도 2에 예시된 에어로졸 발생 물품(200)에서, 마우스피스 세그먼트(240)는 냉각 세그먼트(230)의 하류 또는 하류 방향에 위치하고, 냉각 세그먼트(230)는 마우스피스 세그먼트(240)의 상류 또는 상류 방향에 위치한다.In the following embodiments, “upstream” or “upstream direction” means a direction away from the smoker's mouth, and “downstream” or “downstream direction” means a direction closer to the smoker's mouth. can do. The terms upstream and downstream may be used to describe the relative positions of the elements that make up the smoking article. For example, in the aerosol-generating
이하의 실시예들에서, "퍼프(puff)"는 사용자의 흡입(inhalation)을 의미하며, 흡입이란 사용자의 입이나 코를 통해 사용자의 구강 내, 비강 내 또는 폐로 끌어 당기는 상황을 의미할 수 있다.In the following embodiments, “puff” refers to the user's inhalation, and inhalation may refer to the situation of pulling the product into the user's oral cavity, nasal cavity, or lungs through the user's mouth or nose. .
이하, 본 개시의 몇몇 실시예들에 대하여 첨부된 도면에 따라 상세하게 설명한다.Hereinafter, several embodiments of the present disclosure will be described in detail with reference to the attached drawings.
도 1은 본 개시의 몇몇 실시예들에 따른 에어로졸 시뮬레이션 장치(100)를 나타내는 예시적인 구성도이다.1 is an exemplary configuration diagram showing an
에어로졸 시뮬레이션 장치(100)는 타깃(target) 에어로졸 발생 물품에서 발생되는 에어로졸 현상을 시뮬레이션할 수 있는 컴퓨팅 장치일 수 있다. 여기서, 타깃 에어로졸 발생 물품은 시뮬레이션 대상이 되는 에어로졸 발생 물품을 의미할 수 있다. 또한, 에어로졸 발생 물품은 에어로졸 발생 장치를 통해 에어로졸을 발생시키는 물품으로서, 예를 들어 궐련 스틱이 될 수 있을 것이나, 이에 한정되는 것은 아니다. 또한, 에어로졸 발생 장치는 예를 들어 삽입된 궐련 스틱을 가열하여 에어로졸을 발생시키는 장치일 수 있을 것이나, 이에 한정되는 것은 아니다. 이하에서는, 설명의 편의상, 에어로졸 시뮬레이션 장치(100)를 "시뮬레이션 장치"(100)로 약칭하도록 한다.The
시뮬레이션 장치(100)는 하나 이상의 컴퓨팅 장치로 구현될 수 있다. 다시 말해, 시뮬레이션 장치(100)는 단일 컴퓨팅 장치로 구현될 수도 있고, 복수의 컴퓨팅 장치로 구현될 수도 있다. 예컨대, 시뮬레이션 장치(100)의 제1 기능은 제1 컴퓨팅 장치에서 구현되고, 제2 기능은 제2 컴퓨팅 장치에서 구현될 수 있다. 또는, 시뮬레이션 장치(100)의 특정 기능이 복수의 컴퓨팅 장치에서 구현될 수도 있다.
상기 컴퓨팅 장치는 예를 들어 노트북, 데스크톱(desktop), 랩탑(laptop), 등이 될 수 있으나, 이에 국한되는 것은 아니며 컴퓨팅 기능이 구비된 모든 종류의 장치를 포함할 수 있다. 컴퓨팅 장치의 일 예시에 관하여서는 도 12를 참조하도록 한다.The computing device may be, for example, a laptop, desktop, laptop, etc., but is not limited thereto and may include all types of devices equipped with a computing function. Refer to FIG. 12 for an example of a computing device.
도 1에 도시된 바와 같이, 실시예에 따른 시뮬레이션 장치(100)는 데이터 입력부(110), 모델 구성부(120), 시뮬레이션부(130), 시각화부(140) 및 DB(150)를 포함할 수 있다. 단, 이는 본 개시의 목적을 달성하기 위한 바람직한 실시예일뿐이며, 필요에 따라 일부 구성요소가 추가되거나 생략될 수 있음은 물론이다. 또한, 도 1에 도시된 시뮬레이션 장치(100)의 각각의 구성요소들(e.g. 110)은 기능적으로 구분되는 기능요소들을 나타낸 것으로서, 복수의 구성요소(e.g. 110)가 실제 물리적 환경에서는 서로 통합되는 형태로 구현될 수도 있음에 유의한다. 물론, 실제 물리적 환경에서 상기 각각의 구성요소들(e.g. 110)은 복수의 기능요소로 세분화되는 형태로 구현될 수도 있다. 이하, 시뮬레이션 장치(100)의 각 구성요소들(e.g. 110)에 대하여 설명하도록 한다.As shown in FIG. 1, the
데이터 입력부(110)는 타깃 에어로졸 발생 물품에 관한 시뮬레이션 기초데이터를 입력받을 수 있다. 시뮬레이션 기초데이터는 예를 들어 에어로졸 발생 물품의 종류, 에어로졸 발생 물품의 구조데이터, 에어로졸 발생 물품을 구성하는 재료(e.g. 에어로졸 형성 기재, 래퍼 등)에 관한 데이터(이하, "재료관련데이터"), 에어로졸 발생 물품을 가열하는 히터에 관한 데이터(이하, "히터관련데이터"), 퍼프데이터, 시뮬레이션 설정과 관련된 데이터(이하, "설정관련데이터") 등을 포함할 수 있을 것이나, 이에 한정되는 것은 아니며, 시뮬레이션과 관련된 다양한 데이터를 더 포함할 수 있다. 데이터 입력부(110)를 통해 입력된 시뮬레이션 기초데이터는 DB(150)에 저장될 수 있다.The
상기 구조데이터는 에어로졸 발생 물품의 구조를 나타내거나 정의할 수 있는 다양한 데이터를 포함할 수 있다. 예를 들어, 구조데이터는 에어로졸 발생 물품에 대한 모델링데이터(e.g. 2D 모델링 데이터, 3D 모델링 데이터), 에어로졸 발생 물품의 체제데이터(e.g. 길이, 직경, 내부 구조 등), 에어로졸 발생 물품에 형성된 천공 구조(e.g. 천공 크기, 개수, 위치, 배열 형태 등)에 관한 데이터(이하, "천공관련데이터") 등을 포함할 수 있다. 등을 포함할 수 있다. 그러나, 이에 한정되는 것은 아니다.The structural data may include various data that can indicate or define the structure of the aerosol-generating article. For example, structural data includes modeling data (e.g. 2D modeling data, 3D modeling data) for the aerosol-generating article, regime data (e.g. length, diameter, internal structure, etc.) of the aerosol-generating article, and perforation structure formed in the aerosol-generating article (e.g. It may include data on perforation size, number, location, arrangement type, etc.) (hereinafter referred to as “perforation-related data”). It may include etc. However, it is not limited to this.
또한, 상기 재료관련데이터는 예를 들어 에어로졸 형성 기재(e.g. 니코틴, 글리세린/프로필렌글리콜 등의 보습제, 각초/판상엽 등의 담배 물질 등)의 종류, 함량 및 물성(e.g. 점도, 다공도, 투수율, 밀도 등), 필터 물질(e.g. 셀룰로오스 아세테이트 섬유)의 함량과 물성(e.g. 다공도, 투수율 등), 래퍼의 사양(e.g. 두께, 길이 등)과 물성(e.g. 기공도, 열전도율 등) 등을 포함할 수 있다. 그러나, 이에 한정되는 것은 아니다.In addition, the material-related data includes, for example, the type, content, and physical properties (e.g. viscosity, porosity, water permeability, density) of the aerosol-forming base (e.g. nicotine, humectants such as glycerin/propylene glycol, tobacco substances such as cut grass/lamellae, etc.) etc.), the content and physical properties (e.g. porosity, water permeability, etc.) of the filter material (e.g. cellulose acetate fiber), and the specifications of the wrapper (e.g. thickness, length, etc.) and physical properties (e.g. porosity, thermal conductivity, etc.). However, it is not limited to this.
또한, 상기 히터관련데이터는 히터의 온도 프로파일(e.g. 시간에 따른 가열 온도), 히터의 가열 구조에 관한 데이터(e.g. 가열 방식, 가열 위치, 형상, 길이, 두께 등), 히터의 유형(e.g. 내부 가열식, 외부 가열식, 전기 저항성 히터, 유도 가열식 히터) 등을 포함할 수 있다. 그러나, 이에 한정되는 것은 아니다. 히터의 가열 방식은 내부 가열식과 외부 가열식으로 구분될 수 있을 것이나, 이에 한정되는 것은 아니다.In addition, the heater-related data includes the temperature profile of the heater (e.g. heating temperature over time), data on the heating structure of the heater (e.g. heating method, heating position, shape, length, thickness, etc.), and the type of heater (e.g. internal heating type). , external heating type, electrical resistance heater, induction heating type heater), etc. However, it is not limited to this. The heating method of the heater may be divided into internal heating type and external heating type, but is not limited thereto.
또한, 상기 퍼프데이터는 사용자의 퍼프를 나타내거나 정의할 수 있는 다양한 데이터를 포함할 수 있다. 예를 들어, 퍼프데이터는 퍼프 횟수, 퍼프 강도, 퍼프 간격 등에 관한 데이터를 포함할 수 있다. 그러나, 이에 한정되는 것은 아니다. 참고로, 퍼프데이터의 형태는 다양할 수 있는데, 예를 들어 시간에 따른 부피유량의 형태로 정의될 수 있다(도 5의 360 참조).Additionally, the puff data may include various data that can represent or define the user's puff. For example, puff data may include data on the number of puffs, puff intensity, puff interval, etc. However, it is not limited to this. For reference, the form of puff data may vary. For example, it may be defined in the form of volumetric flow rate over time (see 360 in FIG. 5).
또한, 상기 시뮬레이션 설정관련데이터는 예를 들어 시뮬레이션 수행 시간(또는 시간 스텝의 개수), 시뮬레이션 정확도, 계산격자(calculation grid or computation grid)의 형태(e.g. 2D, 3D, 정렬형, 비정렬형), 계산격자의 조밀도, 계산셀(cell)의 개수, 시뮬레이션 목표요소, 지배방정식의 종류, 초기조건, 경계조건, 관심 영역 등을 포함할 수 있다. 그러나, 이에 한정되는 것은 아니다.In addition, the simulation setting-related data includes, for example, simulation execution time (or number of time steps), simulation accuracy, type of calculation grid or computation grid (e.g. 2D, 3D, aligned, unaligned), It may include the density of the calculation grid, the number of calculation cells, simulation target elements, type of governing equation, initial conditions, boundary conditions, area of interest, etc. However, it is not limited to this.
상기 시뮬레이션 목표요소는 시뮬레이션을 통해 도출하고자 요소를 의미하는 것으로서, 예를 들어 온도 분포, 농도 분포(e.g. 니코틴, 보습제 등의 농도 분포), 속도 분포(e.g. 에어로졸의 속도 분포), 에어로졸 이행량(즉, 무화량), 흡인저항 등을 포함할 수 있으나, 이에 한정되는 것은 아니다. 시뮬레이션을 통해 목표요소의 값을 도출하는 방법에 관하여서는 후술하도록 한다.The simulation target elements refer to elements to be derived through simulation, such as temperature distribution, concentration distribution (e.g. concentration distribution of nicotine, moisturizer, etc.), velocity distribution (e.g. aerosol velocity distribution), and aerosol transfer amount (i.e. , atomization amount), suction resistance, etc., but is not limited thereto. The method of deriving the value of the target element through simulation will be described later.
한편, 데이터 입력부(110)가 시뮬레이션 기초데이터를 입력받는 방식은 어떠한 방식이 되더라도 무방하다. 예를 들어, 데이터 입력부(110)는 웹 인터페이스(e.g. 웹 GUI)를 통해 시뮬레이션 기초데이터를 입력받을 수 있다. 또는, 데이터 입력부(110)는 로컬 저장소에 저장된 파일로부터 시뮬레이션 기초데이터를 입력받을 수도 있다. 물론, 데이터 입력부(110)는 웹 인터페이스를 통해 파일의 형태로 시뮬레이션 기초데이터를 입력받을 수도 있다.Meanwhile, the
몇몇 실시예들에서는, 데이터 입력부(110)가 시뮬레이션 가능한 에어로졸 발생 물품의 리스트를 제공하고, 사용자로부터 타깃 에어로졸 발생 물품을 선택받을 수도 있다. 그러면, 데이터 입력부(110)는 DB(150)를 조회하여 타깃 에어로졸 발생 물품에 대한 시뮬레이션 기초데이터를 획득할 수 있다.In some embodiments, the
다음으로, 모델 구성부(120)는 입력된 시뮬레이션 기초데이터에 기반하여 타깃 에어로졸 발생 물품의 에어로졸 현상을 시뮬레이션하기 위한 CFD(Computational Fluid Dynamics) 모델을 구성할 수 있다. CFD 모델은 유체 현상을 해석하기 위한 모델을 의미하는데, 당해 기술 분야의 종사자라면 CFD 모델에 대하여 자명하게 이해할 수 있을 것인 바, 이에 대한 자세한 설명은 생략하도록 한다.Next, the
도시된 바와 같이, 모델 구성부(120)는 격자 생성부(121), 지배방정식 설정부(122) 및 조건 설정부(123)를 포함하도록 구성될 수 있다. 이하, 모델 구성부(120)의 구성요소에 대하여 설명하도록 한다.As shown, the
격자 생성부(121)는 시뮬레이션 기초데이터에 기반하여 계산격자를 생성할 수 있다. 계산격자는 격자 형태로 이루어진 계산영역을 의미하는데, 계산격자를 구성하는 각 계산셀이 시뮬레이션 계산의 최소 단위가 될 수 있다. 따라서, 계산셀의 개수가 증가할수록(즉, 계산격자가 조밀할수록) 시뮬레이션 정확도가 증가되고, 시뮬레이션에 소요되는 컴퓨팅 비용 또한 증가될 수 있다. 반대로, 계산셀의 개수가 감소할수록(즉, 계산격자가 조밀하지 않을수록) 시뮬레이션 정확도는 떨어지고 컴퓨팅 비용도 감소될 수 있다. 계산격자는 정렬 격자 또는 비정렬 격자의 형태를 가질 수 있으며, 2차원의 격자일 수 있고 3차원의 격자일 수도 있다. 당해 기술 분야의 종사자라면, CFD 모델의 계산격자에 대하여 자명하게 이해할 수 있을 것인 바, 이에 대한 더 이상의 자세한 설명은 생략하도록 한다.The
구체적으로, 격자 생성부(121)는 타깃 에어로졸 발생 물품의 구조데이터에 기반하여 타깃 에어로졸 발생 물품의 물리적 구조에 대응되는 계산격자를 생성할 수 있다. 격자 생성부(121)는 계산격자에 관한 시뮬레이션 설정관련데이터(e.g. 계산격자의 형태, 계산셀의 개수, 조밀도 등)에 더 기반하여 계산격자를 생성할 수도 있다. 보다 이해의 편의를 제공하기 위해, 도 2에 도시된 예를 참조하여 부연 설명하도록 한다.Specifically, the
도 2는 복수의 세그먼트(210 내지 240)로 구성된 타깃 에어로졸 발생 물품(200)에 대한 계산격자(300)를 생성하는 것을 예시하고 있으며, 도 2에 도시된 화살표는 에어로졸(A)의 이동 경로를 나타내고 있다.FIG. 2 illustrates creating a
도 2에 도시된 바와 같이, 타깃 에어로졸 발생 물품(200)이 전단 필터 세그먼트(210), 에어로졸 형성 기재 세그먼트(220), 냉각 세그먼트(230) 및 마우스피스 세그먼트(240) 및 래퍼(250)로 구성되었다고 가정하자. 또한, 전단 필터 세그먼트(210)는 중공 채널을 포함하는 셀룰로오스 아세테이트 필터이고, 에어로졸 형성 기재 세그먼트(220)는 에어로졸 형성 기재를 포함하며, 냉각 세그먼트(230)는 에어로졸을 냉각시키기 위한 중공 채널을 포함하는 지관 필터이고, 마우스피스 세그먼트(240)는 셀룰로오스 아세테이트 필터라고 가정하자. 본 개시의 논지를 흐리지 않기 위해, 각 세그먼트(210 내지 240)에 관한 자세한 설명은 생략하도록 한다. 또한, 타깃 에어로졸 발생 물품(200)은 에어로졸 발생 장치(미도시)의 외부 가열식 히터(260)에 의해 가열되어 에어로졸을 발생시킨다고 가정하자.As shown in Figure 2, the target aerosol-generating
위와 같은 경우, 격자 생성부(121)는 타깃 에어로졸 발생 물품(200)의 구조데이터(e.g. 2D or 3D 모델링 데이터, 체제 데이터)에 기반하여 타깃 에어로졸 발생 물품(200)의 물리적 구조에 대응되는 계산격자(300)를 생성할 수 있다. 가령, 격자 생성부(121)에 의해 생성된 계산격자(300)는 타깃 에어로졸 발생 물품(200)을 구성하는 각 세그먼트(210 내지 240)에 대응되는 계산격자 영역(310 내지 340)을 포함할 수 있다.In the above case, the
다시 도 1을 참조하여 설명한다.Description will be made again with reference to FIG. 1 .
몇몇 실시예들에서는, 격자 생성부(121)가 시뮬레이션 장치(100)의 컴퓨팅 성능 및/또는 가용 컴퓨팅 자원에 기반하여 계산셀의 개수 및/또는 계산격자의 조밀도를 조절할 수 있다. 예를 들어, 가용 컴퓨팅 자원(또는 컴퓨팅 성능)이 기준치 이상인 경우, 격자 생성부(121)는 3차원의 계산격자를 생성하거나 조밀도가 높은 2차원의 계산격자를 생성하는 등 계산셀의 개수를 증가시킬 수 있다. 반대의 경우, 격자 생성부(121)는 2차원의 계산격자를 생성하거나 조밀도가 낮은 3차원의 계산격자를 생성하는 등 계산셀의 개수를 감소시킬 수 있다. 이러한 경우, 시뮬레이션 장치(100)의 컴퓨팅 성능 및/또는 가용 컴퓨팅 자원을 고려하여 시뮬레이션이 효율적으로 수행될 수 있다.In some embodiments, the
또한, 몇몇 실시예들에서는, 격자 생성부(121)가 관심 영역에 대응되는 계산격자 영역의 조밀도(또는 계산셀 개수)를 증가시킬 수 있다. 예를 들어, 타깃 에어로졸 발생 물품(e.g. 200)의 마우스피스 세그먼트(e.g. 240)가 관심 영역으로 지정된 경우, 격자 생성부(121)는 마우스피스 세그먼트(e.g. 240)에 대응되는 계산격자 영역(e.g. 340)의 조밀도를 증가시킬 수 있다. 이러한 경우, 관심 영역에 대한 시뮬레이션 정확도가 향상될 수 있다.Additionally, in some embodiments, the
또한, 몇몇 실시예들에서는, 격자 생성부(121)가 에어로졸의 이동(유동) 경로와 연관성이 높은 제1 계산격자 영역(e.g. 320, 330, 340)의 조밀도(또는 계산셀의 개수)를 증가시키고, 연관성이 낮은 제2 계산격자 영역(e.g. 310)의 조밀도를 감소시킬 수 있다. 여기서, 에어로졸의 이동 경로와 연관성이 높은 영역은 예를 들어 에어로졸의 생성에 영향을 미치는 영역, 에어로졸의 이동 경로와 그 인접 영역을 포함할 수 있을 것이나, 이에 한정되는 것은 아니며, 다양하게 정의될 수 있다. 가령, 격자 생성부(121)는 제1 계산격자 영역에 포함된 계산셀의 개수가 제2 계산격자 영역보다 많도록 계산격자를 생성할 수 있다. 이러한 경우, 에어로졸의 이동 경로와 그 인접 영역에 대한 시뮬레이션 정확도가 향상될 수 있다.Additionally, in some embodiments, the
또한, 몇몇 실시예들에서는, 격자 생성부(121)가 시뮬레이션 수행 중에 동적으로 계산격자의 조밀도(또는 계산셀의 개수)를 조절할 수도 있다. 가령, 격자 생성부(121)는 특정 계산셀을 분할하거나 복수의 계산셀을 병합함으로써 계산셀의 개수를 조절할 수 있는데, 본 실시예와 관련하여서는 추후 도 6을 참조하여 부연 설명하도록 한다.Additionally, in some embodiments, the
다음으로, 지배방정식 설정부(122)는 CFD 모델의 지배방정식(governing equation)을 설정할 수 있다. 상기 지배방정식은 유체 현상을 정의하거나 유체 현상에 영향을 미치는 다양한 변수(e.g. 온도, 압력 등)와의 관계를 기술하고 있는 방정식을 의미하는데, 당해 기술 분야의 종사자라면 지배방정식의 개념에 대하여 이미 숙지하고 있을 것인 바 이에 대한 설명은 생략하도록 한다. 이하에서는, 보다 이해의 편의를 제공하기 위해, 실시예에 따른 CFD 모델에 적용될 수 있는 지배방정식의 종류와 연계 관계에 대하여 도 3을 참조하여 설명하도록 한다.Next, the governing
도 3에 예시된 바와 같이, 상기 지배방정식은 에어로졸의 유동 현상을 시뮬레이션하기 위한 유동 방정식(410)과 열전달 현상을 시뮬레이션하기 위한 열전달 방정식(420)을 포함하고, 그 외에 화학반응 방정식(430) 및 물질이송 방정식(440)을 더 포함할 수 있다. 그러나, 이에 한정되는 것은 아니다. 화학반응 방정식(430) 및 물질이송 방정식(440)은 에어로졸의 유동 현상과 열전달 현상을 보다 정확하게 시뮬레이션하기 위해 또는 에어로졸의 유동 현상과 열전달 현상 외에 물리 화학적 현상들(e.g. 에어로졸의 발생, 물질 이송 등)까지 복합적으로 시뮬레이션하기 위해 이용될 수 있다. 이하, 도 3에 예시된 지배방정식(410 내지 440)에 관하여 설명하도록 한다.As illustrated in FIG. 3, the governing equations include a
유동 방정식(410)은 에어로졸의 유동 현상을 기술하는 방정식으로, 예를 들어 에어로졸 속도(u)와 압력(p) 간의 관계를 나타내거나, 에어로졸 속도(u)와 압력(p) 중 적어도 하나를 변수로 갖는 수식일 수 있다. 보다 구체적인 예로서, 유동 방정식(410)은 질량 보존 방정식(또는 연속 방정식)과 나비에-스토크스 방정식(Navier-Stokes'equation)과 같은 운동량 방정식(momentum equation) 등을 포함할 수 있다. 그러나, 이에 한정되는 것은 아니다.The
질량 보존 방정식과 운동량 방정식의 일 예시는 하기의 수학식 1 및 수학식 2에 기재되어 있다. 하기의 수학식 1은 질량 보존의 법칙으로부터 유도된 방정식이고, 하기의 수학식 2는 운동량 보존의 법칙으로부터 유도된 방정식(또는 운동량 보존 방정식)으로서 유체의 속도, 압력 및 점도 등의 관계를 나타내는 수식이다.An example of the mass conservation equation and momentum equation is described in
[수학식 1][Equation 1]
, ,
[수학식 2][Equation 2]
상기 수학식 1 및 2에서, i는 에어로졸을 구성하는 개별 유체 물질(e.g. 니코틴, 글리세린, 프로필렌글리콜, 타르, 수분 등과 같은 에어로졸 형성 기재가 증기화된 것)을 지칭하고, eja는 개별 유체 물질의 혼합체인 에어로졸을 지칭한다. 여기서, 에어로졸은 에어로졸과 공기의 혼합체를 포함할 수 있다. 그리고, ρi는 물질(i)의 유체 밀도(fluid density)를 의미하고, ni는 물질(i)의 질량 플럭스(mass flux)를 의미하며, ρeja는 에어로졸(eja)의 밀도를 의미하며, u는 유체 속도(fluid velocity)를 의미한다. 그리고, p는 압력을 의미하고, μeja는 에어로졸(eja)의 점도(viscosity)를 의미하며, I는 관련 항이 텐서(tensor) 형태임(또는 텐서 형태로 변환됨)을 의미한다. 이하의 설명에서는, 본 개시의 명료함을 위해, 앞서 설명한 수학식과 중복되는 파라미터(parameter; 변수, 상수 등)에 대한 설명은 생략하도록 하고, 동일한 기호의 파라미터가 다른 의미로 사용되는 경우에는 관련 수식에서 새롭게 설명하도록 한다. 또한, 본 개시에서 언급된 수학식들의 파라미터 중 일부는 시뮬레이션 과정에서 계산될 수 있고, 다른 일부에는 실험적으로 얻어진 값이 대입될 수 있다.In
참고로, 상기 수학식 2는 좌변의 2개 항과 우변의 2개 항으로 구성되는데, 좌변의 첫번째 항은 시간에 따른 속도 변화항이고, 좌변의 두번째 항은 대류항이며, 우변의 첫번째 항은 압력 변화항이고, 우변의 두번째 항은 점도에 관한 유체 저항항이다. 또한, 에어로졸을 구성하는 물질(i)이 복수개로 설정되면, 상기 수학식 1 또한 복수개가 될 수 있으며(마찬가지로, 아래의 수학식들도 복수개가 될 수 있음), 시뮬레이션 수행 시 상기 수학식 1 및 2는 서로 연계되어 계산될 수 있다.For reference,
한편, 도 2에 예시된 타깃 에어로졸 발생 물품(200)의 구조에서, 냉각 세그먼트(230)에서의 에어로졸 유동은 에어로졸 형성 기재 세그먼트(220) 또는 마우스피스 세그먼트(240)와 달라질 수 밖에 없다. 냉각 세그먼트(230)에서의 에어로졸 유동은 장애물이 없는 자유 흐름 영역에서의 기체 이동과 다를 바 없으나, 에어로졸 형성 기재 세그먼트(220) 또는 마우스피스 세그먼트(240)에는 에어로졸의 유동을 방해하는 다공성의 물질(e.g. 필터 물질, 담배 물질)이 존재하기 때문이다. 따라서, 에어로졸 형성 기재 세그먼트(220) 또는 마우스피스 세그먼트(240)와 같이 다공성 물질이 포함된 영역에서의 에어로졸 유동을 보다 정확하게 모사하기 위해, 유동 방정식(410)에 다르시의 법칙(Darcy's law)이 더 적용될 수 있다. 다르시의 법칙은 다공성 매체를 통과하는 유체의 흐름을 관찰하여 얻어진 경험식으로부터 유도된 법칙으로, 당해 기술 분야의 종사자라면 다르시의 법칙에 대해 이미 숙지하고 있을 것인 바, 이에 대한 자세한 설명은 생략하도록 한다. 다르시의 법칙이 적용되는 경우, 상기 수학식 1 및 2는 하기의 수학식 3 및 4의 형태로 변형될 수 있다.On the other hand, in the structure of the target aerosol-generating
[수학식 3][Equation 3]
, ,
[수학식 4][Equation 4]
상기 수학식 3에서, Qi는 물질(i)에 관한 생성항으로, 증발에 의한 기체상 물질(i)의 생성(양수 형태) 또는 액상 물질(i)의 소진(음수 형태)을 의미할 수 있으며, 질량 플럭스의 단위를 가질 수 있다. 그리고, uPM은 다공성 매체(porous medium; PM) 내에서의 유체 속도를 의미하고, ki는 물질(i)에 대한 다공성 매체의 투수율(permeability)을 의미하며, Di는 물질(i)의 확산계수를 의미하고, ci는 물질(i)의 농도를 의미한다.In Equation 3, Q i is a production term related to substance (i), which can mean the creation of gaseous substance (i) by evaporation (positive number form) or the exhaustion (negative form) of liquid substance (i). and can have units of mass flux. And, u PM refers to the fluid velocity within the porous medium (PM), k i refers to the permeability of the porous medium to material (i), and D i refers to the permeability of material (i). It refers to the diffusion coefficient, and c i refers to the concentration of substance (i).
상기 수학식 4는 운동량 보존 방정식인 상기 수학식 2의 속도 관련항(e.g. 속도 변화항, 대류항 등)에 다공성 매체의 다공도를 반영한 수식으로 이해될 수 있으며, k는 다공성 매체를 지칭하고, φk는 다공성 매체(k)의 다공도를 의미하며, Keja는 에어로졸(eja)에 대한 다공성 매체의 투수율을 의미한다.Equation 4 can be understood as an equation that reflects the porosity of the porous medium in the velocity-related terms (eg velocity change term, convection term, etc.) of
다음으로, 열전달 방정식(420)은 에어로졸 발생 물품 내부의 열 전달을 기술하는 방정식으로, 예를 들어 에어로졸의 유동 관련 변수(e.g. 속도 u, 압력 p)와 온도(T)와의 관계를 나타내거나, 에어로졸 유동 관련 변수(e.g. 속도 u, 압력 p)와 온도(T) 중 적어도 하나를 변수로 갖는 수식일 수 있다. 예를 들어, 열전달 방정식(420)은 에너지 보존 방정식을 포함할 수 있다. 그러나, 이에 한정되는 것은 아니다.Next, the
에너지 보존 방정식의 일 예시는 하기의 수학식 5 및 6에 기재되어 있다. 하기의 수학식 5는 자유 흐름 영역에 적용될 수 있는 에너지 보존 방정식을 나타내고, 하기의 수학식 6은 다공성 매체 영역에 적용될 수 있는 에너지 보존 방정식을 나타낸다.An example of an energy conservation equation is shown in Equations 5 and 6 below. Equation 5 below represents the energy conservation equation applicable to the free flow region, and Equation 6 below represents the energy conservation equation applicable to the porous media region.
[수학식 5][Equation 5]
[수학식 6][Equation 6]
상기 수학식 5 및 6에서, i는 기체(증기)상의 개별 유체 물질(즉, 기체 상의 에어로졸 형성 기재)을 지칭하고, j는 액체상의 개별 유체 물질(즉, 액체상의 에어로졸 형성 기재)을 지칭하며, tot는 복수의 개별 물질(i, j)을 포함하는 전체 물질을 지칭한다. 그리고, Cp,eja는 에어로졸(eja)의 열용량(Cp)을 의미하고, T는 온도를 의미하며, hi는 물질(i)의 열전달 계수를 의미한다. 그리고, a는 공기(air)를 지칭하고, na는 공기(a)의 질량 플럭스를 의미하며, ha는 공기(a)의 열전달 계수를 의미하고, keja는 에어로졸(eja)의 열전도율을 의미한다. 그리고, eff는 다공성 매체 내에서 기체상과 액상이 혼재하는 상태(e.g. 액상의 니코틴과 기체상의 니코틴이 혼재하는 상태, 기체 에어로졸과 액화된 에어로졸이 혼재하는 상태)를 지칭하고, Cp,eff는 혼재 상태(eff)에서 전체 물질(tot)의 열용량을 의미하며, nj는 물질(j)의 질량 플럭스를 의미하고, hj는 물질(j)의 열전달 계수를 의미하며, keff,tot는 혼재 상태(eff)에서 전체 물질(tot)의 열전도율을 의미한다. 그리고, Hevap,i는 물질(i)의 기화 열용량을 의미하고, q는 온도에 관한 생성항으로 에어로졸 발생 물품 내의 별도의 발열 물질이 존재하지 않는 경우에는 0으로 설정될 수 있다.In equations 5 and 6 above, i refers to an individual fluid substance in the gas (vapor) phase (i.e., an aerosol-forming substrate in the gas phase), and j refers to an individual fluid substance in the liquid phase (i.e., an aerosol-forming substrate in the liquid phase), , tot refers to the entire substance including a plurality of individual substances (i, j). And, C p,eja means the heat capacity (C p ) of the aerosol (eja), T means the temperature, and h i means the heat transfer coefficient of the substance (i). And, a refers to air, n a refers to the mass flux of air (a), h a refers to the heat transfer coefficient of air (a), and k eja refers to the thermal conductivity of aerosol (eja). it means. And, eff refers to a state in which gas and liquid phases coexist in a porous medium (eg, a state in which liquid nicotine and gaseous nicotine coexist, gas aerosol and liquefied aerosol coexist), and C p,eff is It refers to the heat capacity of the entire material (tot) in the mixed state (eff), n j refers to the mass flux of material (j), h j refers to the heat transfer coefficient of material (j), and k eff,tot is It refers to the thermal conductivity of the entire material (tot) in the mixed state (eff). In addition, H evap,i refers to the heat capacity of vaporization of substance (i), and q is a generation term related to temperature and can be set to 0 when no separate heating material is present in the aerosol-generating article.
다음으로, 화학반응 방정식(430)은 에어로졸 생성에 관한 화학반응을 기술하는 방정식으로, 예를 들어 에어로졸과 연관된 화학반응속도(K)와 물질 농도(c) 및 온도(T)와의 관계를 나타내거나, 상기 화학반응속도(K)와 물질 농도(c) 및 온도(T) 중 적어도 하나를 변수로 갖는 수식일 수 있다. 예를 들어, 화학반응 방정식(430)은 아레니우스 방정식(Arrhenius equation), 증발량 계산식, 이상기체 방정식을 포함할 수 있다. 그러나, 이에 한정되는 것은 아니다. 아레니우스 방정식, 증발량 계산식, 이상기체 방정식의 일 예시는 각각 하기의 수학식 7, 8 및 9에 기재되어 있다.Next, the
[수학식 7][Equation 7]
, ,
[수학식 8][Equation 8]
[수학식 9][Equation 9]
상기 수학식 7 내지 9에서, Kf는 반응속도상수를 의미하고, Af는 빈도인자를 의미하며, Ea는 활성화에너지를 의미하고, R은 기체상수를 의미하며, αi는 물질(i)의 반응속도 변수(kinetic parameter)를 의미하고, E0i는 물질(i)의 초기 활성화에너지를 의미한다. 그리고, EVi와 mevap,i는 물질(i)의 증발량(e.g. Kg/(m3*s) or mol/(m3*s))을 의미하고, Ki는 물질(i)의 증발속도상수를 의미하며, Csat,i는 물질(i)의 포화농도를 의미하고, Ci는 물질(i)의 농도를 의미한다. 그리고, ai는 물질(i)에 대한 포화농도(Csat,i)를 실제 환경(또는 시뮬레이션 목적)에 맞게 조절(제어)하는 인자로서, 적절한 값으로 설정될 수 있다. 포화농도(Csat,i)는 수학식 9의 포화농도(csat) 계산식에 기초하여 계산될 수 있으며, 수학식 7 내지 9는 연계되어 계산될 수 있다.In Equations 7 to 9, K f means the reaction rate constant, A f means the frequency factor, E a means the activation energy, R means the gas constant, and α i means the substance (i ), and E 0i refers to the initial activation energy of substance (i). And, EV i and m evap,i mean the evaporation amount of substance (i) (eg Kg/(m 3 *s) or mol/(m 3 *s)), and K i is the evaporation rate of substance (i) It means a constant, C sat,i means the saturation concentration of substance (i), and C i means the concentration of substance (i). In addition, a i is a factor that adjusts (controls) the saturation concentration (C sat,i ) for substance (i) to suit the actual environment (or simulation purpose), and can be set to an appropriate value. The saturation concentration (C sat,i ) can be calculated based on the saturation concentration (c sat ) calculation formula in Equation 9, and Equations 7 to 9 can be calculated in conjunction.
참고로, 빈도인자(Af)와 반응속도 변수(αi)의 값은 실험적으로 얻을 수 있는데, 예를 들어 TGA-DSC, TG-FTIR, GC-mass, LC-mass 등의 기법을 통해 얻을 수 있다. 또한, 수학식 7의 반응속도상수(Kf)의 값은 수학식 8의 증발속도상수(Ki)에 적용(대입)될 수 있고, 수학식 9의 포화농도(csat) 값은 수학식 8의 포화농도(Csat)에 적용(대입)될 수 있다.For reference, the values of the frequency factor (A f ) and reaction rate variable (α i ) can be obtained experimentally, for example, through techniques such as TGA-DSC, TG-FTIR, GC-mass, and LC-mass. You can. In addition, the value of the reaction rate constant (K f ) in Equation 7 can be applied (substituted) for the evaporation rate constant (K i ) in Equation 8, and the saturation concentration (c sat ) value in Equation 9 is Equation It can be applied (substituted) for a saturation concentration (C sat ) of 8.
다음으로, 물질이송 방정식(440)은 에어로졸 발생 물품 내에서의 물질 이송을 기술하는 방정식으로, 예를 들어 에어로졸 속도(u)와 이송량과의 관계를 나타내거나, 에어로졸 속도(u)와 이송량 중 적어도 하나를 변수로 갖는 수식일 수 있다. 예를 들어, 물질이송 방정식(440)은 종이송 방정식(species transport equation)을 포함할 수 있다. 그러나, 이에 한정되는 것은 아니다.Next, the
종이송 방정식의 일 예시는 하기의 수학식 10 및 11에 기재되어 있다. 하기의 수학식 10은 자유 흐름 영역에 적용될 수 있는 이송 방정식을 나타내고, 하기의 수학식 11은 다공성 매체 영역에 적용될 수 있는 이송 방정식을 나타낸다.An example of the paper song equation is shown in
[수학식 10][Equation 10]
[수학식 11][Equation 11]
상기 수학식 10 및 11에서, ξi는 물질(i)의 이송량으로 질량 플럭스 또는 농도의 단위를 가질 수 있고, Di는 물질(i)의 확산계수를 의미하며, Deff,i는 다공성 매체 내에서 액체상과 기체상이 혼재된 상태(eff)에 있는 물질(i)의 확산계수를 의미한다. 기타 다른 파라미터의 의미는 앞선 수학식을 참조하도록 한다.In
참고로, 수학식 11의 속도항(uPM 관련항)은 수학식 10의 속도항(u 관련항)에 다르시의 법칙이 적용된 것으로 이해될 수 있고, 수학식 11의 확산항(Deff,i 관련항)은 액체상과 기체상의 혼합으로 이루어진 유체의 특성을 표현하는 식으로 변형된 것으로 이해될 수 있다.For reference, the velocity term (u PM- related term) in Equation 11 can be understood as Darcy's law applied to the velocity term (u-related term) in
지금까지 설명한 지배방정식들은 도 3에 도시된 바와 같이 서로 연계되어 계산될 수 있다. 예를 들어, 상기 수학식 1 내지 11이 CFD 모델의 지배방정식으로 설정된 경우, 수학식 1 내지 4와 같은 유동 방정식(e.g. 410)을 통해 계산된 속도(u) 값 또는 압력(p) 값이 수학식 10 또는 11과 같은 물질이송 방정식(e.g. 440)에 적용(대입)될 수 있다. 다른 예로서, 수학식 5 또는 6과 같은 열전달 방정식(e.g. 420)을 통해 계산된 온도(T) 값이 수학식 7 내지 9와 같은 화학반응 방정식(e.g. 430)에 적용(대입)될 수 있다. 또 다른 예로서, 수학식 10 또는 11와 같은 물질이송 방정식(e.g. 440)을 통해 계산된 물질 이송량(ξ) 값이 수학식 7 내지 9와 같은 화학반응 방정식(e.g. 430)의 농도(C) 값에 영향을 미칠 수 있다. 또 다른 예로서, 수학식 5 또는 6과 같은 열전달 방정식(e.g. 420)을 통해 계산된 온도(T) 값이 수학식 1 내지 4와 같은 유동 방정식(e.g. 410)의 밀도(ρ), 점도(μ), 확산계수(D) 등에 영향을 미칠 수 있다.The governing equations described so far can be calculated in conjunction with each other as shown in FIG. 3. For example, when
이와 같이, 특정 지배방정식을 통해 계산된 변수의 값은 다른 지배방정식에 적용(대입)되거나 영향을 미침으로써, 복수의 지배방정식이 서로 연계되어 계산될 수 있으며, 이에 따라 타깃 에어로졸 발생 물품에서 나타나는 에어로졸 유동 현상과 열전달 현상뿐만 아니라 다른 복합적인 물리 화학 현상(e.g. 에어로졸의 발생, 물질 이송 등)까지 정확하게 시뮬레이션될 수 있다. 가령, 에어로졸의 발생과 유동 현상은 화학반응, 열전달, 유체의 유동 등이 복잡하게 얽혀있는 현상인데, 복수의 지배방정식을 연계하여 계산함으로써 이러한 현상들이 정확하게 시뮬레이션될 수 있다. 참고로, 특정 지배방정식의 변수 값은 동일 시간 스텝에서 다른 지배방정식의 변수에 적용(대입)될 수도 있고, 다음 시간 스텝에서 적용(대입)될 수도 있다.In this way, the values of variables calculated through a specific governing equation can be applied (substituted) or influenced by other governing equations, so that multiple governing equations can be calculated in conjunction with each other, thereby allowing the aerosol that appears in the target aerosol-generating article to be calculated. Not only flow and heat transfer phenomena, but also other complex physical and chemical phenomena (e.g. aerosol generation, material transport, etc.) can be accurately simulated. For example, the generation and flow of aerosols is a phenomenon in which chemical reactions, heat transfer, and fluid flow are complexly intertwined, and these phenomena can be accurately simulated by calculating multiple governing equations in conjunction. For reference, the variable values of a specific governing equation may be applied (substituted) to variables of another governing equation at the same time step, or may be applied (substituted) at the next time step.
한편, 지배방정식 설정부(122)는 사용자로부터 직접 입력받은 수식을 CFD 모델의 지배방정식으로 설정할 수 있고, 시뮬레이션 목표요소에 기초하여 자동으로 지배방정식을 설정할 수도 있다. 가령, 시뮬레이션 목표요소가 무화량인 경우, 지배방정식 설정부(122)는 유동 방정식(410), 열전달 방정식(420), 화학반응 방정식(430) 및 물질이송 방정식(440)을 CFD 모델의 지배방정식으로 설정할 수 있다. 다른 예로서, 시뮬레이션 목표요소가 흡인저항인 경우, 지배방정식 설정부(122)는 유동 방정식(410)만을 CFD 모델의 지배방정식으로 설정할 수 있다. 또 다른 예로서, 시뮬레이션 목표 요소가 온도 분포인 경우, 지배방정식 설정부(122)는 유동 방정식(410)과 열전달 방정식(420)을 CFD 모델의 지배방정식으로 설정할 수 있다.Meanwhile, the governing
다시 도 1을 참조하여 시뮬레이션 장치(100)의 다른 구성요소에 대한 설명을 이어가도록 한다.The description of other components of the
조건 설정부(123)는 CFD 모델의 초기조건 및/또는 경계조건을 설정할 수 있다. 초기조건 및/또는 경계조건은 시뮬레이션을 정확하게 수행하기 위해 설정되는 조건들로서, 수치 해석 기법을 통해 지배방정식을 풀이하는 데 이용될 수 있다.The
상기 초기조건은 다양한 시뮬레이션 변수의 초기값을 설정하는 조건을 포함할 수 있고, 조건 설정부(123)는 시뮬레이션 기초데이터 및/또는 입력된 초기값(DB 150에 미리 저장된 초기값)에 따라 초기조건을 설정할 수 있다. 예를 들어, 타깃 에어로졸 발생 물품의 시뮬레이션 변수(e.g. 지배방정식의 변수들)가 속도(에어로졸의 속도), 압력, 온도 및 농도(e.g. 구성 물질의 농도)를 포함하는 경우, 도 4에 도시된 바와 같이, 조건 설정부(123)는 계산격자(300)를 구성하는 각 계산셀(e.g. 321, 341)에 대해 속도, 압력, 온도 및 농도에 대한 초기조건(e.g. 322, 342)을 설정할 수 있다.The initial conditions may include conditions for setting initial values of various simulation variables, and the
보다 구체적인 예로서, 조건 설정부(123)는 재료관련데이터(e.g. 구성 물질의 함량 데이터)에 기반하여 물질 농도에 대한 초기조건을 설정할 수 있다. 이를테면, 조건 설정부(123)는 물질(e.g. 담배 물질, 보습제 등)의 함량(e.g. 충진량)을 상기 물질이 포함된 부분(e.g. 에어로졸 형성 기재 세그먼트 260)의 부피로 나누어 농도값을 구하고, 구해진 농도값을 초기조건으로 설정할 수 있다. 또한, 조건 설정부(123)는 [속도 = 0], [압력 = 1기압], [온도 = 상온]과 같은 형태로 속도, 압력 및 온도에 대한 초기조건을 설정할 수 있다. 흡연 시작 전에는 히터에서 열이 발생되지 않고 에어로졸 발생 물품 내에 에어로졸의 유동도 없을 것이기 때문이다.As a more specific example, the
상기 경계조건은 에어로졸의 유동, 히터의 열, 물질의 전달 등에 관한 다양한 경계조건을 포함할 수 있고, 조건 설정부(123)는 시뮬레이션 기초데이터 또는 입력된 경계조건에 따라 CFD 모델의 경계조건을 설정할 수 있다.The boundary conditions may include various boundary conditions related to aerosol flow, heater heat, material transfer, etc., and the
예를 들어, 조건 설정부(123)는 퍼프데이터에 기반하여 에어로졸 토출(흡입)에 관한 유동 경계조건을 설정할 수 있다. 또한, 조건 설정부(123)는 상기 구조 데이터에 기반하여 에어로졸의 유동을 제한하거나 차단하기 위한 벽면 경계조건(e.g. non-slip)을 설정할 수 있다. 가령, 조건 설정부(123)는 타깃 에어로졸 발생 물품의 하류 말단이 아닌 다른 방향으로 에어로졸이 이동하는 것을 제한하거나 차단하기 위한 벽면 경계조건을 설정할 수 있다. 또한, 조건 설정부(123)는 히터관련데이터에 기반하여 열 경계조건을 설정할 수 있으며, 타깃 에어로졸 발생 물품(200)의 재료관련데이터에 기반하여 유동 또는 물질 전달에 관한 경계 조건(e.g. 벽면 경계 조건)을 설정할 수 있다. 이외에도, 조건 설정부(123)는 재료관련데이터(e.g. 래퍼의 종류, 물성 등)에 기초하여 단열 경계 조건을 설정할 수도 있다. 보다 이해의 편의를 제공하기 위해, 도 5를 참조하여 경계조건을 설정하는 예에 관하여 부연 설명하도록 한다.For example, the
도 5는 퍼프데이터(360)와 히터의 온도 프로파일(350)에 기반하여 유동 경계조건과 열 경계조건을 설정하는 예를 도시하고 있다. 참고로, 도 5는 퍼프데이터(360)가 시간(t)에 따른 부피유량(Qv; ml/s)의 변화로 정의된 것을 예로써 도시하고 있으며, 부피유량의 변화가 특정 시점(t1) 이후에 발생되는 것은 예열이 완료된 이후(t1 시점)에 퍼프가 발생되는 것을 모사한 것으로 이해될 수 있다.Figure 5 shows an example of setting flow boundary conditions and thermal boundary conditions based on
도 5에 도시된 바와 같이, 조건 설정부(123)는 히터의 온도 프로파일(350)과 가열 구조에 관한 데이터에 기반하여 열 경계조건을 설정할 수 있다. 가령, 조건 설정부(123)는 히터의 가열 구조에 기반하여 히터로부터 타깃 에어로졸 발생 물품(e.g. 200)으로 열이 전달되는 경계영역(370, 380)을 결정하고, 결정된 경계영역(370, 380)에 온도 프로파일(350)에 따른 경계조건을 설정할 수 있다. 이러한 경우, 히터의 가열 구조 및 온도 프로파일에 의한 가열이 CFD 모델에 정확하게 반영될 수 있다.As shown in FIG. 5, the
참고로, 도 6은 타깃 에어로졸 발생 물품(e.g. 200)을 가열하는 히터(e.g. 260)가 외부 가열식이기 때문에(도 2 참조), 경계영역(370, 380)이 타깃 에어로졸 발생 물품(e.g. 200)의 외곽부에 형성되는 것을 예로써 도시하고 있을 뿐이고, 히터의 가열 구조가 변경되면 열전달과 관련된 경계영역(e.g. 370, 380)은 다른 위치에 형성될 수 있다. 가령, 대상 히터가 내부 가열식인 경우, 상기 경계영역(e.g. 370, 380)은 타깃 에어로졸 발생 물품의 내부에 형성될 수 있다.For reference, in FIG. 6, since the heater (e.g. 260) that heats the target aerosol-generating article (e.g. 200) is an external heating type (see FIG. 2), the
또한, 조건 설정부(123)는 퍼프데이터(360)에 기반하여 유동 경계조건을 설정할 수 있다. 가령, 조건 설정부(1230)는 퍼프에 의해 에어로졸이 토출되는 하류 말단 영역(390)에 유동 경계조건을 설정할 수 있다. 이러한 경우, 퍼프에 의한 에어로졸의 흡입이 CFD 모델에 정확하게 반영될 수 있다.Additionally, the
다시 도 1을 참조하여 설명한다.Description will be made again with reference to FIG. 1 .
다른 예로서, 조건 설정부(123)는 천공관련데이터에 기반하여 천공에 관한 유동경계조건을 설정할 수 있다. 이때, 천공관련데이터는 타깃 에어로졸 발생 물품에 형성되어 있는 천공의 위치, 개수, 배열 형태, 천공을 통해 유입되는 외기의 부피유량, 외기의 온도 등에 관한 데이터를 포함할 수 있고, 상기 유동경계조건은 상기 천공관련데이터에 기반하여 천공을 통해 외기가 유입되는 경계영역에 설정될 수 있다. 본 예시에서, 천공에 관한 유동경계조건이 설정되면, 천공을 통해 유입되는 외기의 유량과 온도를 고려하여 에어로졸 현상(e.g. 에어로졸의 유동 현상, 열전달 현상)이 시뮬레이션될 수 있으며, 시뮬레이션 결과는 천공 설계에 효과적으로 활용될 수 있다.As another example, the
다음으로, 시뮬레이션부(130)는 구성된 CFD 모델을 이용하여 타깃 에어로졸 발생 물품에 대한 에어로졸 현상을 시뮬레이션할 수 있다. 구체적으로, 시뮬레이션부(130)는 수치 해석 기법(e.g. FDM, FEM, FVM 등)을 통해 지배방정식을 대수방정식으로 변환하고, 초기조건 및 경계조건에 따라 대수방정식을 풀이함으로써 에어로졸 현상과 연관된 시뮬레이션 변수의 값을 계산할 수 있다. 가령, 도 6에 도시된 바와 같이, 시뮬레이션부(130)는 계산격자(300)에 포함된 각 계산셀(e.g. 321, 341)에 대하여 복수의 시간 스텝에 걸쳐 시뮬레이션 변수(e.g. 속도, 압력, 온도, 농도 등)의 값을 반복적으로 계산할 수 있다. 이러한 계산 과정에 따라, 흡연에 의해 타깃 에어로졸 발생 물품 내부에서 발생되는 다양한 에어로졸 현상이 정확하게 시뮬레이션될 수 있다.Next, the
일 예로서, 시뮬레이션부(130)는 유동 방정식(e.g. 410)을 계산함으로써, 타깃 에어로졸 발생 물품 내의 에어로졸 유동 현상을 시뮬레이션할 수 있다. 그리고 시뮬레이션 결과로, 타깃 에어로졸 발생 물품 내의 압력 분포와 속도 분포(e.g. 에어로졸의 속도 분포)가 도출될 수 있다. 이때, 시뮬레이션부(130)는 타깃 에어로졸 발생 물품(e.g. 200)의 자유 흐름 영역 세그먼트(e.g. 냉각 세그먼트 230)에 대응되는 계산격자 영역(e.g. 330)에 대해서는 자유 흐름 영역에 대한 유동방정식(e.g. 수학식 1 및 2)에 따라 계산을 수행하고, 다공성 매체 영역 세그먼트(e.g. 에어로졸 형성 기재 세그먼트 220, 마우스피스 세그먼트 240)에 대응되는 계산격자 영역(e.g. 320, 340)에 대해서는 다공성 매체 영역에 대한 유동방정식(e.g. 수학식 3 및 4)에 따라 계산을 수행할 수 있다. 본 예시에서, 시뮬레이션부(130)는 흡연 완료 후의 압력 분포, 특정 흡연 시점의 압력 분포, 퍼프에 따른 압력 분포 변화, 흡연 동안의 압력 분포 변화, 특정 위치의 시간에 따른 압력 변화 등을 도출할 수 있다. 또한, 시뮬레이션부(130)는 퍼프에 따른 에어로졸 이동 속도 변화, 이동 경로 등을 도출할 수도 있다. 또한, 시뮬레이션부(130)는 타깃 에어로졸 발생 물품의 양단부의 압력 차이(즉, 양단부에 대응되는 계산격자 영역 또는 계산셀의 압력값 차이)에 기초하여 타깃 에어로졸 발생 물품 전체에 대한 흡인저항을 도출할 수 있다. 유사한 방식으로, 시뮬레이션부(130)는 타깃 에어로졸 발생 물품의 일부분에 대한 흡인저항을 도출할 수도 있다. 가령, 시뮬레이션부(130)는 타깃 에어로졸 발생 물품(e.g. 200)을 구성하는 마우스피스 세그먼트(e.g. 240)의 양단부에 대응되는 계산격자 영역(또는 계산셀)의 압력값 차이에 기초하여 마우스피스 세그먼트(240)에 대한 흡인저항을 도출할 수 있다.As an example, the
다른 예로서, 시뮬레이션부(130)는 유동 방정식(e.g. 410)과 열전달 방정식(e.g. 420)을 연계하여 계산함으로써, 타깃 에어로졸 발생 물품 내의 에어로졸 유동 현상과 열전달 현상을 시뮬레이션할 수 있다. 그리고 시뮬레이션 결과로, 타깃 에어로졸 발생 물품 내의 압력 분포, 속도 분포 및 온도 분포가 도출될 수 있다. 이때, 시뮬레이션부(130)는 타깃 에어로졸 발생 물품(e.g. 200)의 자유 흐름 영역 세그먼트(e.g. 냉각 세그먼트 230)에 대응되는 계산격자 영역(e.g. 330)에 대해서는 자유 흐름 영역에 대한 지배방정식(e.g. 수학식 1, 2, 5)에 따라 계산을 수행하고, 다공성 매체 영역 세그먼트(e.g. 에어로졸 형성 기재 세그먼트 220, 마우스피스 세그먼트 240)에 대응되는 계산격자 영역(e.g. 320, 340)에 대해서는 다공성 매체 영역에 대한 지배방정식(e.g. 수학식 3, 4, 6)에 따라 계산을 수행할 수 있다. 본 예시에서, 시뮬레이션부(130)는 흡연 완료 후의 온도 분포, 특정 흡연 시점의 온도 분포, 퍼프에 따른 온도 분포 변화, 흡연 동안의 온도 분포 변화, 특정 위치의 시간에 따른 온도 변화 등을 도출할 수 있다.As another example, the
또 다른 예로서, 시뮬레이션부(130)는 유동 방정식(e.g. 410), 열전달 방정식(e.g. 420), 화학반응 방정식(e.g. 430) 및 물질이송 방정식(e.g. 440)을 연계하여 계산함으로써, 에어로졸의 발생, 유동, 열전달 및 물질 이송을 포함하는 다양한 에어로졸 현상을 시뮬레이션할 수 있다. 그리고 시뮬레이션 결과로, 타깃 에어로졸 발생 물품 내의 압력 분포, 속도 분포, 온도 분포, 농도 분포, 화학반응속도, 외부로 토출되는 에어로졸 이송량 등이 도출될 수 있다. 이때, 시뮬레이션부(130)는 타깃 에어로졸 발생 물품(e.g. 200)의 자유 흐름 영역 세그먼트(e.g. 냉각 세그먼트 230)에 대응되는 계산격자 영역(e.g. 330)에 대해서는 자유 흐름 영역에 대한 지배방정식(e.g. 수학식 1, 2, 5, 7 내지 10)에 따라 계산을 수행하고, 다공성 매체 영역 세그먼트(e.g. 에어로졸 형성 기재 세그먼트 220, 마우스피스 세그먼트 240)에 대응되는 계산격자 영역(e.g. 320, 340)에 대해서는 다공성 매체 영역에 대한 지배방정식(e.g. 수학식 3, 4, 6, 7 내지 9, 11)에 따라 계산을 수행할 수 있다. 본 예시에서, 시뮬레이션부(130)는 타깃 에어로졸 발생 물품으로부터 에어로졸이 토출되는 하류 말단 영역에 대응되는 계산격자 영역의 에어로졸 이송량을 종합하여(즉, 상기 계산격자 영역에 포함된 계산셀의 에어로졸 이송량을 종합하여) 에어로졸 이행량(무화량)을 도출할 수 있다. 이를테면, 계산격자가 2D 격자인 경우 시뮬레이션부(130)는 면적분(surface integral)을 통해 에어로졸 이송량을 종합할 수 있고, 3D 격자인 경우 부피적분(volume integral)을 통해 에어로졸 이송량을 종합할 수 있다. 또한, 시뮬레이션부(130)는 퍼프에 따른 농도 분포 변화, 흡연 동안의 농도 분포 변화, 특정 위치의 시간에 따른 농도 변화, 물질 별 농도 변화, 퍼프에 따른 에어로졸 이행량 변화 등을 도출할 수 있다.As another example, the
몇몇 실시예들에서는, 도 6에 도시된 바와 같이, 시뮬레이션부(130)가 복수의 시간 스텝에 걸쳐 시뮬레이션 변수를 계산하는 도중에, 계산격자(300)의 적어도 일부 영역의 조밀도가 격자 생성부(121)에 의해 동적으로 조절될 수 있다. 가령, 격자 생성부(121)는 시뮬레이션 변수(e.g. 온도, 압력 등)의 값이 기준치 이상으로 변동되는 계산셀을 탐지하고, 탐지된 계산셀을 복수개로 분할할 수 있다(즉, 조밀도를 증가시킬 수 있다). 또는, 격자 생성부(121)는 시뮬레이션 변수(e.g. 온도, 압력 등)의 값이 기준치 이하로 변동되는 복수개의 계산셀을 탐지하고, 탐지된 복수개의 계산셀을 병합할 수도 있다(즉, 조밀도를 감소시킬 수 있다). 시뮬레이션 변수의 값이 크게 변동하는 계산격자 영역은 에어로졸 현상과 밀접하게 연관된 중요한 영역일 가능성이 높다. 따라서, 본 실시예에 따르면, 중요 계산격자 영역에 대한 시뮬레이션 정확도가 향상되고, 중요하지 않은 계산격자 영역에 투입되는 컴퓨팅 비용은 절감되어, 보다 효율적으로 시뮬레이션이 수행될 수 있다.In some embodiments, as shown in FIG. 6, while the
또한, 몇몇 실시예들에서는, 시뮬레이션부(130)가 시뮬레이션 결과를 토대로 타깃 에어로졸 발생 물품 내의 데드존(dead zone)을 검출할 수 있다. 여기서, 데드존은 흡연 시에 타깃 에어로졸 발생 물품 내에서 잘 활용되지 않는 영역을 의미할 수 있다. 구체적으로, 시뮬레이션부(130)는 전체 계산격자 영역에서 화학반응속도, 온도 또는 물질의 농도 감소 정도가 기준치 이하인 영역(또는 계산셀들)을 데드존으로 검출할 수 있다. 이렇게 검출된 데드존은 에어로졸 발생 물품, 히터 등의 구조를 재설계하는데 효과적으로 활용될 수 있다. 보다 이해의 편의를 제공하기 위해, 도 7에 도시된 예를 참조하여 부연 설명하도록 한다.Additionally, in some embodiments, the
도 7은 에어로졸 형성 기재부(520)와 필터부(510)를 포함하는 타깃 에어로졸 발생 물품(500)의 온도 분포에 관한 시뮬레이션 결과(600)를 예시하고 있다. 특히, 도 7은 계산격자(600)의 온도값을 히트맵 형태로 시각화한 것을 예시하고 있다.FIG. 7 illustrates
도 7을 참조하면, 시뮬레이션부(130)는 계산격자(600) 내에서 온도가 기준치 이하인 영역(610, 620)을 데드존으로 검출할 수 있다. 이때, 상기 기준치는 계산격자 영역 별 또는 계산셀 별로 다른 값으로 설정될 수 있다. 가령, 히터로부터 가까이 위치한 계산격자 영역의 기준치는 높은 값으로 설정되고, 멀리 위치한 계산격자 영역의 기준치는 낮은 값으로 설정될 수 있다. 또는, 에어로졸 형성 기재부(520)에 대응되는 계산격자 영역의 기준치는 높은 값으로 설정되고, 필터부(510)에 대응되는 계산격자 영역의 기준치는 낮은 값으로 설정될 수 있다.Referring to FIG. 7 , the
다시 도 1을 참조하여 시뮬레이션 장치(100)의 다른 구성요소에 대한 설명을 이어가도록 한다.The description of other components of the
시각화부(140)는 타깃 에어로졸 발생 물품에 대한 시뮬레이션 결과를 시각화할 수 있다. 예를 들어, 시각화부(140)는 타깃 에어로졸 발생 물품을 나타내는 2D 또는 3D 객체 상에 시뮬레이션 결과(e.g. 온도 분포, 속도 분포, 압력 분포, 농도 분포 등)를 히트맵 형태로 시각화할 수 있다(도 7 내지 도 10 참조). 또한, 시각화부(140)는 에어로졸의 이동 경로 및/또는 이동 방향을 상기 객체 상에 표시하거나, 데드존을 강조하여 표시할 수도 있다.The
보다 구체적인 예로서, 도 8에 도시된 바와 같이, 시각화부(140)는 타깃 에어로졸 발생 물품을 나타내는 2D 객체 상에 물품 내부의 온도 분포가 히트맵 형태로 도시된 시각화 정보(700)를 생성할 수 있다. 물론, 시각화부(140)는 다른 시뮬레이션 결과(e.g. 속도 분포, 농도 분포, 압력 분포 등)에 대해서도 도 8 등에 도시된 바와 같은 시각화 정보를 생성할 수 있다.As a more specific example, as shown in FIG. 8, the
다른 예로서, 도 9에 도시된 바와 같이, 시각화부(140)는 타깃 에어로졸 발생 물품을 나타내는 3D 객체 상에 물품 내부의 온도 분포가 히트맵 형태로 도시된 시각화 정보(800)를 생성할 수도 있다.As another example, as shown in FIG. 9, the
또 다른 예로서, 도 10에 도시된 바와 같이, 시각화부(140)는 타깃 에어로졸 발생 물품의 특정 위치(e.g. 포인트 A, B, C)의 시간에 따른 온도 변화가 그래프 형태로 도시된 시각화 정보(900)를 생성할 수도 있다.As another example, as shown in FIG. 10, the
또한, 도시되어 있지는 않으나, 시각화부(140)는 에어로졸의 이동 방향을 화살표의 형태로 시각화할 수 있다. 이때, 시각화부(140)는 화살표의 굵기(또는 개수)를 통해 에어로졸의 이동량을 표시할 수 있고, 화살표의 색상을 통해 이동 속도를 표시할 수도 있다.Additionally, although not shown, the
다음으로, DB(150)에는 에어로졸 시뮬레이션과 관련된 각종 데이터가 저장될 수 있다. 가령, DB(150)에는 에어로졸 발생 물품의 유형에 따른 시뮬레이션 기초데이터가 저장되어 있을 수 있다. 또는, DB(150)에는 에어로졸 발생 물품의 유형(또는 시뮬레이션 기초데이터)에 따른 시뮬레이션 결과가 저장되어 있을 수 있다.Next, various data related to aerosol simulation may be stored in the
몇몇 실시예들에서, 시뮬레이션부(130)는 시뮬레이션 수행 전에 입력된 시뮬레이션 기초데이터에 대응되는 시뮬레이션 결과가 DB(150)에 저장되어 있는지 여부를 판단할 수 있다. 그리고, 저장되어 있다는 판단에 응답하여, 시뮬레이션부(130)는 시뮬레이션을 수행하지 않고 기 저장된 시뮬레이션 결과를 제공할 수 있다. 이러한 경우, 시뮬레이션에 수행되는 컴퓨팅 비용이 크게 절감되는 효과가 달성될 수 있다.In some embodiments, the
한편, 본 개시의 몇몇 실시예들에 따르면, 시뮬레이션 장치(100)가 사용자에게 에어로졸 시뮬레이션 서비스를 제공할 수도 있다. 예컨대, 시뮬레이션 장치(100)는 웹 인터페이스 등을 통해 사용자로부터 시뮬레이션 요청을 수신하고, 상기 요청에 응답하여 에어로졸 시뮬레이션 서비스를 제공할 수 있다. 이때, 시뮬레이션 장치(100)는 사용자에게 서비스 이용 요금을 요청할 수도 있다. 서비스 이용 요금은 시뮬레이션에 투입된 컴퓨팅 자원에 기반하여 산정될 수 있을 것이나, 다른 방식으로 산정되더라도 무방하다. 몇몇 예들에서는, 시뮬레이션 요청 시 사용자가 시뮬레이션 정확도(또는 서비스 이용 요금)을 선택할 수도 있다. 그러면, 시뮬레이션 장치(100)가 선택된 정확도(또는 서비스 이용 요금)에 따른 계산격자를 생성하여 시뮬레이션을 수행할 수 있다. 가령, 사용자가 높은 정확도를 선택한 경우, 시뮬레이션 장치(100)는 계산셀의 개수가 많고 조밀한 계산격자를 생성할 수 있다.Meanwhile, according to some embodiments of the present disclosure, the
도 1의 각 구성요소(e.g. 110)는 소프트웨어(Software) 또는, FPGA(Field Programmable Gate Array)나 ASIC(Application-Specific Integrated Circuit)과 같은 하드웨어(Hardware)를 의미할 수 있다. 그렇지만, 상기 구성요소들은 소프트웨어 또는 하드웨어에 한정되는 의미는 아니며, 어드레싱(Addressing)할 수 있는 저장 매체에 있도록 구성될 수도 있고, 하나 또는 그 이상의 프로세서들을 실행시키도록 구성될 수도 있다. 상기 구성요소들 안에서 제공되는 기능은 더 세분화된 구성요소에 의하여 구현될 수 있으며, 복수의 구성요소들을 합하여 특정한 기능을 수행하는 하나의 구성요소로 구현될 수도 있다.Each component (e.g. 110) in FIG. 1 may refer to software or hardware such as FPGA (Field Programmable Gate Array) or ASIC (Application-Specific Integrated Circuit). However, the above components are not limited to software or hardware, and may be configured to reside in an addressable storage medium, and may be configured to execute one or more processors. The functions provided within the above components may be implemented by more detailed components, or may be implemented as a single component that performs a specific function by combining multiple components.
지금까지, 도 1 내지 도 10을 참조하여 본 개시의 몇몇 실시예들에 따른 시뮬레이션 장치(100)의 구성 및 동작에 대하여 설명하였다. 상술한 바에 따르면, CFD 모델을 통해 타깃 에어로졸 발생 물품에서 발생되는 에어로졸 현상이 정확하게 시뮬레이션될 수 있다. 이에 따라, 타깃 에어로졸 발생 물품에 대한 제조 및 실험 없이도 설계 결과가 용이하게 검증될 수 있으며, 에어로졸 발생 물품의 연구 개발에 투입되는 인적 비용 및 시간 비용이 크게 절감될 수 있다.So far, the configuration and operation of the
이하에서는, 도 11 이하의 도면을 참조하여 본 개시의 몇몇 실시예들에 따른 에어로졸 시뮬레이션 방법에 대하여 상세하게 설명하도록 한다.Hereinafter, an aerosol simulation method according to some embodiments of the present disclosure will be described in detail with reference to the drawings of FIG. 11 and below.
상기 방법의 각 단계는 컴퓨팅 장치에 의해 수행될 수 있다. 다시 말하면, 상기 방법의 각 단계는 컴퓨팅 장치의 프로세서에 의해 실행되는 하나 이상의 인스트럭션들로 구현될 수 있다. 상기 방법에 포함되는 모든 단계는 하나의 물리적인 컴퓨팅 장치에 의하여 실행될 수도 있을 것이나, 복수의 물리적인 컴퓨팅 장치에 의해 분산되어 실행될 수도 있다. 예를 들면, 상기 방법의 제1 단계들은 제1 컴퓨팅 장치에 의하여 수행되고, 상기 방법의 제2 단계들은 제2 컴퓨팅 장치에 의하여 수행될 수도 있다. 이하에서는, 이해의 편의를 제공하기 위해, 상기 방법의 각 단계가 도 1에 예시된 시뮬레이션 장치(100)에 의해 수행되는 것을 가정하여 설명을 이어가도록 한다. 따라서, 아래의 설명에서 각 동작의 주어가 생략된 경우, 상기 예시된 장치(100)에 의하여 수행되는 것으로 이해될 수 있다.Each step of the method may be performed by a computing device. In other words, each step of the method may be implemented as one or more instructions executed by a processor of a computing device. All steps included in the method may be executed by one physical computing device, or may be distributed and executed by a plurality of physical computing devices. For example, the first steps of the method may be performed by a first computing device and the second steps of the method may be performed by a second computing device. Hereinafter, for convenience of understanding, the description will be continued assuming that each step of the method is performed by the
도 11은 본 개시의 몇몇 실시예들에 따른 에어로졸 시뮬레이션 방법을 나타내는 예시적인 흐름도이다. 단, 이는 본 개시의 목적을 달성하기 위한 바람직한 실시예일뿐이며, 필요에 따라 일부 단계가 추가되거나 삭제될 수 있음은 물론이다.11 is an example flowchart illustrating an aerosol simulation method according to some embodiments of the present disclosure. However, this is only a preferred embodiment for achieving the purpose of the present disclosure, and of course, some steps may be added or deleted as needed.
도 11에 도시된 바와 같이, 상기 에어로졸 시뮬레이션 방법은 타깃 에어로졸 발생 물품에 관한 시뮬레이션 기초데이터를 입력받는 단계 S100에서 시작될 수 있다. 본 단계에 관하여서는 데이터 입력부(110)에 대한 설명을 더 참조하도록 한다.As shown in FIG. 11, the aerosol simulation method may begin at step S100 of receiving basic simulation data regarding the target aerosol-generating article. For further information on this step, please refer to the description of the
단계 S200에서, 입력된 시뮬레이션 기초데이터를 기초로 CFD 모델이 구성될 수 있다. 가령, 시뮬레이션 장치(100)는 입력된 시뮬레이션 기초데이터를 기초로 계산격자를 생성하고, 지배방정식을 설정하며, 초기조건과 경계조건을 설정함으로써 CFD 모델을 구성할 수 있다. 본 단계에 관하여서는 모델 구성부(120)에 대한 설명을 더 참조하도록 한다.In step S200, a CFD model may be constructed based on the input basic simulation data. For example, the
단계 S300에서, 구성된 CFD 모델을 이용하여 타깃 에어로졸 발생 물품에 대한 에어로졸 현상이 시뮬레이션될 수 있다. 가령, 시뮬레이션 장치(100)는 설정된 초기조건과 경계조건에 따라 지배방정식을 계산하되, 복수의 시간 스텝에 따라 반복적으로 지배방정식을 계산함으로써 에어로졸 현상을 시뮬레이션할 수 있다. 이때, 지배방정식에 대한 계산은 계산격자를 구성한 계산셀 단위로 수행될 수 있다. 보다 구체적인 예로서, 시뮬레이션 장치(100)는 복수의 시간 스텝에 따라 유동 방정식(e.g. 410)과 열전달 방정식(e.g. 420)을 연계하여 반복 계산함으로써 타깃 에어로졸 발생 물품에서 나타나는 열전달 현상과 에어로졸 유동 현상을 시뮬레이션할 수 있다. 본 단계에 관하여서는 시뮬레이션부(130)에 대한 설명을 더 참조하도록 한다.In step S300, aerosol phenomena for the target aerosol-generating article may be simulated using the constructed CFD model. For example, the
단계 S400에서, 시뮬레이션 결과가 시각화될 수 있다. 가령, 시뮬레이션 장치(100)는 타깃 에어로졸 발생 물품을 가리키는 2차원 또는 3차원 객체 상에 시뮬레이션 결과를 시각화할 수 있다. 이를테면, 시뮬레이션 장치(100)는 온도 분포, 압력 분포, 속도 분포, 농도 분포 등을 히트맵 등의 형태로 시각화하거나, 상기 객체 상에 에어로졸의 이동 경로를 표시하거나, 데드존을 강조하여 표시할 수 있다. 다만, 시각화 방식이 이에 한정되는 것은 아니다. 본 단계에 관하여서는 시각화부(140)에 대한 설명을 더 참조하도록 한다.In step S400, the simulation results can be visualized. For example, the
참고로, 도 11예 예시된 단계 S100 내지 단계 S400은 각각 데이터 입력부(110), 모델 구성부(120), 시뮬레이션부(130) 및 시각화부(140)에 의해 수행될 수 있다.For reference, steps S100 to S400 illustrated in FIG. 11 may be performed by the
지금까지, 도 11을 참조하여 본 개시의 몇몇 실시예들에 따른 에어로졸 시뮬레이션 방법에 대하여 설명하였다. 상술한 바에 따르면, CFD 모델을 통해 타깃 에어로졸 발생 물품에서 발생되는 에어로졸 현상이 정확하게 시뮬레이션될 수 있다. 이에 따라, 타깃 에어로졸 발생 물품에 대한 제조 및 실험 없이도 설계 결과가 용이하게 검증될 수 있으며, 에어로졸 발생 물품의 연구 개발에 투입되는 인적 비용 및 시간 비용이 크게 절감될 수 있다.So far, the aerosol simulation method according to some embodiments of the present disclosure has been described with reference to FIG. 11. According to the above, the aerosol phenomenon generated from the target aerosol-generating article can be accurately simulated through the CFD model. Accordingly, design results can be easily verified without manufacturing and testing target aerosol-generating products, and human and time costs invested in research and development of aerosol-generating products can be greatly reduced.
이하에서는, 도 1에 예시된 시뮬레이션 장치(100)를 구현할 수 있는 예시적인 컴퓨팅 장치(1000)에 대하여 설명하도록 한다.Hereinafter, an
도 12는 본 개시의 컴퓨팅 장치(1000)를 나타내는 예시적인 하드웨어 구성도이다.FIG. 12 is an exemplary hardware configuration diagram showing the
도 12에 도시된 바와 같이, 컴퓨팅 장치(1000)는 하나 이상의 프로세서(1100), 버스(1300), 통신 인터페이스(1400), 프로세서(1100)에 의하여 수행되는 컴퓨터 프로그램을 로드(load)하는 메모리(1200)와, 컴퓨터 프로그램(1600)를 저장하는 스토리지(1500)를 포함할 수 있다. 다만, 도 12에는 본 개시의 실시예와 관련 있는 구성요소들만이 도시되어 있다. 따라서, 본 개시가 속한 기술분야의 통상의 기술자라면 도 12에 도시된 구성요소들 외에 다른 범용적인 구성요소들이 더 포함될 수 있음을 알 수 있다. 즉, 컴퓨팅 장치(1000)에는, 도 12에 도시된 구성요소 이외에도 다양한 구성요소가 더 포함될 수 있다.As shown in FIG. 12, the
프로세서(1100)는 컴퓨팅 장치(1000)의 각 구성의 전반적인 동작을 제어한다. 프로세서(1100)는 CPU(Central Processing Unit), MPU(Micro Processor Unit), MCU(Micro Controller Unit), GPU(Graphic Processing Unit) 또는 본 개시의 기술 분야에 잘 알려진 임의의 형태의 프로세서 중 적어도 하나를 포함하여 구성될 수 있다. 또한, 프로세서(1100)는 본 개시의 실시예들에 따른 방법/동작을 실행하기 위한 적어도 하나의 애플리케이션 또는 프로그램에 대한 연산을 수행할 수 있다. 컴퓨팅 장치(1000)는 하나 이상의 프로세서를 구비할 수 있다.The
메모리(1200)는 각종 데이터, 명령 및/또는 정보를 저장한다. 메모리(1200)는 본 개시의 실시예들에 따른 방법/동작을 실행하기 위하여 스토리지(1500)로부터 하나 이상의 프로그램(1600)을 로드할 수 있다. 가령, 메모리(1200)에 컴퓨터 프로그램(1600)이 로드되면, 도 1에 도시된 바와 같은 모듈이 메모리(1200) 상에 구현될 수 있다. 메모리(1200)는 RAM과 같은 휘발성 메모리로 구현될 수 있을 것이나, 이에 한정되는 것은 아니다.The
버스(1300)는 컴퓨팅 장치(1000)의 구성요소 간 통신 기능을 제공한다. 버스(1300)는 주소 버스(Address Bus), 데이터 버스(Data Bus) 및 제어 버스(Control Bus) 등 다양한 형태의 버스로 구현될 수 있다.
통신 인터페이스(1400)는 컴퓨팅 장치(1000)의 유무선 인터넷 통신을 지원한다. 또한, 통신 인터페이스(1400)는 인터넷 통신 외의 다양한 통신 방식을 지원할 수도 있다. 이를 위해, 통신 인터페이스(1400)는 본 개시의 기술 분야에 잘 알려진 통신 모듈을 포함하여 구성될 수 있다. 몇몇 실시예들에서, 통신 인터페이스(1400)는 생략될 수도 있다.The communication interface 1400 supports wired and wireless Internet communication of the
스토리지(1500)는 하나 이상의 프로그램(1600)을 비임시적으로 저장할 수 있다. 스토리지(1500)는 ROM(Read Only Memory), EPROM(Erasable Programmable ROM), EEPROM(Electrically Erasable Programmable ROM), 플래시 메모리 등과 같은 비휘발성 메모리, 하드 디스크, 착탈형 디스크, 또는 본 개시가 속하는 기술 분야에서 잘 알려진 임의의 형태의 컴퓨터로 읽을 수 있는 기록 매체를 포함하여 구성될 수 있다.
컴퓨터 프로그램(1600)은 메모리(1200)에 로드될 때 프로세서(1100)로 하여금 본 개시의 다양한 실시예에 따른 방법/동작을 수행하도록 하는 하나 이상의 인스트럭션들(instructions)을 포함할 수 있다. 즉, 프로세서(1100)는 상기 하나 이상의 인스트럭션들을 실행함으로써, 본 개시의 실시예들에 따른 방법/동작을 수행할 수 있다.The
예를 들어, 컴퓨터 프로그램(1600)은 타깃 에어로졸 발생 물품의 구조데이터를 포함하는 시뮬레이션 기초데이터를 입력받는 동작, 입력된 시뮬레이션 기초데이터에 기반하여 CFD 모델을 구성하는 동작 및 구성된 CFD 모델을 이용하여 타깃 에어로졸 발생 물품에서 나타나는 열전달 현상과 에어로졸의 유동 현상을 시뮬레이션하는 동작을 수행하도록 하는 인스트럭션들을 포함할 수 있다. 이러한 경우, 컴퓨팅 장치(1000)를 통해 본 개시의 몇몇 실시예들에 따른 시뮬레이션 장치(100)가 구현될 수 있다.For example, the
지금까지 도 1 내지 도 12를 참조하여 설명된 본 개시의 기술적 사상은 컴퓨터가 읽을 수 있는 매체 상에 컴퓨터가 읽을 수 있는 코드로 구현될 수 있다. 상기 컴퓨터로 읽을 수 있는 기록 매체는, 예를 들어 이동형 기록 매체(CD, DVD, 블루레이 디스크, USB 저장 장치, 이동식 하드 디스크)이거나, 고정식 기록 매체(ROM, RAM, 컴퓨터 구비 형 하드 디스크)일 수 있다. 상기 컴퓨터로 읽을 수 있는 기록 매체에 기록된 상기 컴퓨터 프로그램은 인터넷 등의 네트워크를 통하여 다른 컴퓨팅 장치에 전송되어 상기 다른 컴퓨팅 장치에 설치될 수 있고, 이로써 상기 다른 컴퓨팅 장치에서 사용될 수 있다.The technical idea of the present disclosure described so far with reference to FIGS. 1 to 12 may be implemented as computer-readable code on a computer-readable medium. The computer-readable recording medium may be, for example, a removable recording medium (CD, DVD, Blu-ray disk, USB storage device, removable hard disk) or a fixed recording medium (ROM, RAM, computer-equipped hard disk). You can. The computer program recorded on the computer-readable recording medium can be transmitted to another computing device through a network such as the Internet, installed on the other computing device, and thus used on the other computing device.
이상에서, 본 개시의 실시예를 구성하는 모든 구성요소들이 하나로 결합되거나 결합되어 동작하는 것으로 설명되었다고 해서, 본 개시의 기술적 사상이 반드시 이러한 실시예에 한정되는 것은 아니다. 즉, 본 개시의 목적 범위 안에서라면, 그 모든 구성요소들이 하나 이상으로 선택적으로 결합하여 동작할 수도 있다.In the above, even though all the components constituting the embodiments of the present disclosure have been described as being combined or operated in combination, the technical idea of the present disclosure is not necessarily limited to these embodiments. That is, within the scope of the purpose of the present disclosure, all of the components may be operated by selectively combining one or more of them.
도면에서 동작들이 특정한 순서로 도시되어 있지만, 반드시 동작들이 도시된 특정한 순서로 또는 순차적 순서로 실행되어야만 하거나 또는 모든 도시 된 동작들이 실행되어야만 원하는 결과를 얻을 수 있는 것으로 이해되어서는 안 된다. 특정 상황에서는, 멀티태스킹 및 병렬 처리가 유리할 수도 있다. 더욱이, 위에 설명한 실시예들에서 다양한 구성들의 분리는 그러한 분리가 반드시 필요한 것으로 이해되어서는 안 되고, 설명된 프로그램 컴포넌트들 및 시스템들은 일반적으로 단일 소프트웨어 제품으로 함께 통합되거나 다수의 소프트웨어 제품으로 패키지 될 수 있음을 이해하여야 한다.Although operations are shown in the drawings in a specific order, it should not be understood that the operations must be performed in the specific order shown or sequential order or that all illustrated operations must be performed to obtain the desired results. In certain situations, multitasking and parallel processing may be advantageous. Moreover, the separation of the various components in the embodiments described above should not be construed as necessarily requiring such separation, and the program components and systems described may generally be integrated together into a single software product or packaged into multiple software products. You must understand that it exists.
이상 첨부된 도면을 참조하여 본 개시의 실시예들을 설명하였지만, 본 개시가 속하는 기술분야에서 통상의 지식을 가진 자는 그 기술적 사상이나 필수적인 특징을 변경하지 않고서 본 개시가 다른 구체적인 형태로도 실시될 수 있다는 것을 이해할 수 있다. 그러므로 이상에서 기술한 실시예들은 모든 면에서 예시적인 것이며 한정적인 것이 아닌 것으로 이해해야만 한다. 본 개시의 보호 범위는 아래의 청구범위에 의하여 해석되어야 하며, 그와 동등한 범위 내에 있는 모든 기술 사상은 본 개시에 의해 정의되는 기술적 사상의 권리범위에 포함되는 것으로 해석되어야 할 것이다.Although embodiments of the present disclosure have been described above with reference to the attached drawings, those skilled in the art will understand that the present disclosure can be implemented in other specific forms without changing the technical idea or essential features. I can understand that there is. Therefore, the embodiments described above should be understood in all respects as illustrative and not restrictive. The scope of protection of this disclosure should be interpreted in accordance with the claims below, and all technical ideas within the equivalent scope should be interpreted as being included in the scope of rights of the technical ideas defined by this disclosure.
Claims (13)
상기 시뮬레이션 기초데이터에 기반하여 CFD(Computational Fluid Dynamics) 모델을 구성하는 모델 구성부; 및
상기 구성된 CFD 모델을 이용하여 상기 타깃 에어로졸 발생 물품에서 나타나는 열전달 현상과 에어로졸의 유동 현상을 시뮬레이션하는 시뮬레이션부를 포함하되,
상기 모델 구성부는,
상기 구조데이터에 기반하여 계산격자를 생성하는 격자 생성부;
상기 CFD 모델의 지배방정식을 설정하는 지배방정식 설정부; 및
상기 CFD 모델의 초기조건 및 경계조건을 설정하는 조건 설정부를 포함하고,
상기 격자 생성부는 관심 영역에 대응하여 계산 격자 영역의 조밀도를 조절하고,
상기 지배방정식 설정부는 제1 지배방정식 및 제2 지배방정식을 상기 CFD 모델의 지배방정식으로 설정하고,
상기 시뮬레이션부는 상기 설정된 제1 지배방정식 및 제2 지배방정식을 연계하여 계산함으로써 상기 에어로졸의 유동 현상과 상기 열전달 현상을 시뮬레이션하되,
상기 제1 지배방정식은 상기 에어로졸의 유동에 관한 방정식이고,
상기 제2 지배방정식은 상기 타깃 에어로졸 발생 물품 내의 열전달에 관한 방정식인,
에어로졸 시뮬레이션 장치.A data input unit that receives simulation basic data including structural data of the target aerosol-generating article;
A model configuration unit that configures a CFD (Computational Fluid Dynamics) model based on the simulation basic data; and
It includes a simulation unit that simulates the heat transfer phenomenon and aerosol flow phenomenon that appear in the target aerosol-generating article using the configured CFD model,
The model component part is,
a grid generator that generates a calculation grid based on the structural data;
A governing equation setting unit that sets the governing equation of the CFD model; and
It includes a condition setting unit that sets initial conditions and boundary conditions of the CFD model,
The grid generator adjusts the density of the calculation grid area in response to the area of interest,
The governing equation setting unit sets the first governing equation and the second governing equation as the governing equation of the CFD model,
The simulation unit simulates the flow phenomenon of the aerosol and the heat transfer phenomenon by calculating in conjunction with the set first and second governing equations,
The first governing equation is an equation for the flow of the aerosol,
The second governing equation is an equation related to heat transfer within the target aerosol-generating article,
Aerosol simulation device.
상기 시뮬레이션부는 상기 설정된 지배방정식, 상기 설정된 초기조건 및 상기 설정된 경계조건에 따라 상기 에어로졸의 유동 현상과 연관된 시뮬레이션 변수의 값을 계산을 계산하되, 상기 생성된 계산격자에 포함된 계산셀에 대하여 복수의 시간 스텝에 걸쳐 반복적인 계산을 수행하고,
상기 격자 생성부는 상기 반복적인 계산 도중에 상기 시뮬레이션 변수의 값이 기준치 이상으로 변동되는 계산셀을 탐지하고, 상기 탐지된 계산셀을 복수개로 분할하는,
에어로졸 시뮬레이션 장치.According to claim 1,
The simulation unit calculates the values of simulation variables related to the flow phenomenon of the aerosol according to the set governing equations, the set initial conditions, and the set boundary conditions, and calculates a plurality of calculation cells included in the generated calculation grid. Perform iterative calculations over time steps,
The grid generator detects a calculation cell in which the value of the simulation variable changes more than a standard value during the repetitive calculation, and divides the detected calculation cell into a plurality of cells.
Aerosol simulation device.
상기 시뮬레이션부는 상기 설정된 지배방정식, 상기 설정된 초기조건 및 상기 설정된 경계조건에 따라 상기 에어로졸의 유동 현상과 연관된 시뮬레이션 변수의 값을 계산하되, 상기 생성된 계산격자에 포함된 계산셀에 대하여 복수의 시간 스텝에 걸쳐 반복적인 계산을 수행하고,
상기 격자 생성부는 상기 반복적인 계산 도중에 상기 시뮬레이션 변수의 값이 기준치 이하로 변동되는 복수개의 계산셀을 탐지하고, 상기 탐지된 복수개의 계산셀을 병합하는,
에어로졸 시뮬레이션 장치.According to claim 1,
The simulation unit calculates the values of simulation variables associated with the flow phenomenon of the aerosol according to the set governing equations, the set initial conditions, and the set boundary conditions, and performs a plurality of time steps for the calculation cells included in the generated calculation grid. Perform iterative calculations over,
The grid generator detects a plurality of calculation cells in which the value of the simulation variable changes below a standard value during the repetitive calculation, and merges the detected plurality of calculation cells.
Aerosol simulation device.
상기 조건 설정부는 상기 구조데이터에 기반하여 상기 에어로졸의 유동과 관련된 벽면의 경계조건을 설정하는,
에어로졸 시뮬레이션 장치.According to claim 1,
The condition setting unit sets boundary conditions of the wall related to the flow of the aerosol based on the structural data,
Aerosol simulation device.
상기 시뮬레이션 기초데이터는 퍼프데이터를 더 포함하고,
상기 조건 설정부는 상기 퍼프데이터에 기반하여 상기 에어로졸에 대한 유동경계조건을 설정하되,
상기 유동경계조건은 상기 타깃 에어로졸 발생 물품으로부터 상기 에어로졸이 토출되는 하류 말단 영역에 설정되는 경계조건인 것인,
에어로졸 시뮬레이션 장치.According to claim 1,
The simulation basic data further includes puff data,
The condition setting unit sets flow boundary conditions for the aerosol based on the puff data,
The flow boundary condition is a boundary condition set at the downstream end area where the aerosol is discharged from the target aerosol-generating article,
Aerosol simulation device.
상기 시뮬레이션 기초데이터는 상기 타깃 에어로졸 발생 물품을 가열하는 히터의 온도 프로파일을 더 포함하고,
상기 조건 설정부는 상기 온도 프로파일에 기반하여 열 경계조건을 설정하는,
에어로졸 시뮬레이션 장치.According to claim 1,
The simulation basic data further includes a temperature profile of a heater that heats the target aerosol-generating article,
The condition setting unit sets thermal boundary conditions based on the temperature profile,
Aerosol simulation device.
상기 시뮬레이션 기초데이터는 상기 히터의 가열 구조에 관한 데이터를 더 포함하고,
상기 조건 설정부는 상기 히터의 가열 구조에 기반하여 상기 히터로부터 타깃 에어로졸 발생 물품으로 열이 전달되는 경계영역을 결정하고, 상기 결정된 경계영역에 상기 열 경계조건을 설정하는,
에어로졸 시뮬레이션 장치.According to clause 7,
The simulation basic data further includes data regarding the heating structure of the heater,
The condition setting unit determines a boundary area where heat is transferred from the heater to the target aerosol-generating article based on the heating structure of the heater, and sets the thermal boundary condition in the determined boundary area,
Aerosol simulation device.
상기 시뮬레이션부는,
상기 생성된 계산격자 내에서 상기 설정된 초기조건 및 상기 설정된 경계조건에 따라 상기 설정된 지배방정식을 계산함으로써 상기 타깃 에어로졸 발생 물품 내의 온도 분포를 도출하고,
상기 도출 결과에 기반하여 온도가 기준치 이하인 계산격자 영역을 검출하는,
에어로졸 시뮬레이션 장치.According to claim 1,
The simulation unit,
Derive the temperature distribution within the target aerosol-generating article by calculating the set governing equation according to the set initial conditions and the set boundary conditions within the generated calculation grid,
Detecting a calculation grid area where the temperature is below the standard value based on the derived results,
Aerosol simulation device.
상기 시뮬레이션의 결과를 시각화하는 시각화부를 더 포함하는,
에어로졸 시뮬레이션 장치.According to claim 1,
Further comprising a visualization unit that visualizes the results of the simulation,
Aerosol simulation device.
타깃 에어로졸 발생 물품의 구조 데이터를 포함하는 시뮬레이션 기초데이터를 입력받는 단계;
상기 시뮬레이션 기초데이터에 기반하여 CFD(Computational Fluid Dynamics) 모델을 구성하는 단계; 및
상기 구성된 CFD 모델을 이용하여 상기 타깃 에어로졸 발생 물품에서 나타나는 열전달 현상과 에어로졸의 유동 현상을 시뮬레이션하는 단계를 포함하되,
상기 CFD 모델을 구성하는 단계는,
상기 구조 데이터에 기반하여 계산 격자를 생성하는 단계,
상기 CFD 모델의 지배방정식을 설정하는 단계, 및
상기 시뮬레이션 기초데이터에 기반하여 초기조건 및 경계조건을 설정하는 단계를 포함하고,
상기 계산 격자를 생성하는 단계는 관심 영역에 대응하여 계산 격자 영역의 조밀도를 조절하되,
상기 CFD 모델의 지배방정식을 설정하는 단계는 제1 지배방정식 및 제2 지배방정식을 상기 CFD 모델의 지배방정식으로 설정하는 단계를 포함하고,
상기 시뮬레이션하는 단계는 상기 설정된 제1 지배방정식 및 상기 설정된 제2 지배방정식을 연계하여 계산함으로써 상기 에어로졸의 유동 현상과 상기 열전달 현상을 시뮬레이션하는 단계를 포함하되,
상기 제1 지배방정식은 상기 에어로졸의 유동에 관한 방정식이고,
상기 제2 지배방정식은 상기 타깃 에어로졸 발생 물품 내의 열전달에 관한 방정식인,
에어로졸 시뮬레이션 방법.In an aerosol simulation method performed on a computing device,
Receiving basic simulation data including structural data of the target aerosol-generating article;
Constructing a CFD (Computational Fluid Dynamics) model based on the simulation basic data; and
Including the step of simulating heat transfer phenomena and aerosol flow phenomena appearing in the target aerosol-generating article using the configured CFD model,
The step of constructing the CFD model is,
generating a computational grid based on the structural data;
Setting the governing equations of the CFD model, and
It includes setting initial conditions and boundary conditions based on the simulation basic data,
In the step of generating the calculation grid, the density of the calculation grid area is adjusted according to the area of interest,
The step of setting the governing equation of the CFD model includes setting the first governing equation and the second governing equation as the governing equation of the CFD model,
The simulating step includes simulating the flow phenomenon of the aerosol and the heat transfer phenomenon by calculating in conjunction with the set first governing equation and the set second governing equation,
The first governing equation is an equation for the flow of the aerosol,
The second governing equation is an equation related to heat transfer within the target aerosol-generating article,
Aerosol simulation method.
상기 컴퓨터 프로그램은,
타깃 에어로졸 발생 물품의 구조 데이터를 포함하는 시뮬레이션 기초데이터를 입력받는 단계;
상기 시뮬레이션 기초데이터에 기반하여 CFD(Computational Fluid Dynamics) 모델을 구성하는 단계; 및
상기 구성된 CFD 모델을 이용하여 상기 타깃 에어로졸 발생 물품에서 나타나는 열전달 현상과 에어로졸의 유동 현상을 시뮬레이션하는 단계를 실행하고,
상기 CFD 모델을 구성하는 단계는,
상기 구조 데이터에 기반하여 계산 격자를 생성하는 단계,
상기 CFD 모델의 지배방정식을 설정하는 단계, 및
상기 시뮬레이션 기초데이터에 기반하여 초기조건 및 경계조건을 설정하는 단계를 포함하고,
상기 계산 격자를 생성하는 단계는 관심 영역에 대응하여 계산 격자 영역의 조밀도를 조절하되,
상기 CFD 모델의 지배방정식을 설정하는 단계는 제1 지배방정식 및 제2 지배방정식을 상기 CFD 모델의 지배방정식으로 설정하는 단계를 포함하고,
상기 시뮬레이션하는 단계는 상기 설정된 제1 지배방정식 및 상기 설정된 제2 지배방정식을 연계하여 계산함으로써 상기 에어로졸의 유동 현상과 상기 열전달 현상을 시뮬레이션하는 단계를 포함하되,
상기 제1 지배방정식은 상기 에어로졸의 유동에 관한 방정식이고,
상기 제2 지배방정식은 상기 타깃 에어로졸 발생 물품 내의 열전달에 관한 방정식인,
컴퓨터 프로그램.In a computer program combined with a computing device and stored on a computer-readable recording medium,
The computer program is,
Receiving basic simulation data including structural data of the target aerosol-generating article;
Constructing a CFD (Computational Fluid Dynamics) model based on the simulation basic data; and
Executing the step of simulating the heat transfer phenomenon and aerosol flow phenomenon that appear in the target aerosol-generating article using the configured CFD model,
The step of constructing the CFD model is,
generating a computational grid based on the structural data;
Setting the governing equations of the CFD model, and
It includes setting initial conditions and boundary conditions based on the simulation basic data,
In the step of generating the calculation grid, the density of the calculation grid area is adjusted according to the area of interest,
The step of setting the governing equation of the CFD model includes setting the first governing equation and the second governing equation as the governing equation of the CFD model,
The simulating step includes simulating the flow phenomenon of the aerosol and the heat transfer phenomenon by calculating in conjunction with the set first governing equation and the set second governing equation,
The first governing equation is an equation for the flow of the aerosol,
The second governing equation is an equation related to heat transfer within the target aerosol-generating article,
computer program.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
KR1020200142924A KR102640561B1 (en) | 2020-10-30 | 2020-10-30 | Apparatus for simulating aerosol and method thereof |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
KR1020200142924A KR102640561B1 (en) | 2020-10-30 | 2020-10-30 | Apparatus for simulating aerosol and method thereof |
Publications (2)
Publication Number | Publication Date |
---|---|
KR20220057847A KR20220057847A (en) | 2022-05-09 |
KR102640561B1 true KR102640561B1 (en) | 2024-02-23 |
Family
ID=81582287
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
KR1020200142924A KR102640561B1 (en) | 2020-10-30 | 2020-10-30 | Apparatus for simulating aerosol and method thereof |
Country Status (1)
Country | Link |
---|---|
KR (1) | KR102640561B1 (en) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2008108147A (en) * | 2006-10-26 | 2008-05-08 | Riam-Compact Co Ltd | Fluid analysis support system and program |
KR101363964B1 (en) * | 2010-04-22 | 2014-02-18 | 윈챵 슈우 | Combined multifunctional electronic simulated cigarette |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20090000644A (en) * | 2007-03-14 | 2009-01-08 | 쎄딕(주) | System and method for analyzing flow field and flow noise automatically |
GB201805268D0 (en) * | 2018-03-29 | 2018-05-16 | Nicoventures Trading Ltd | Apaaratus for generating aerosol from an aerosolisable medium, an article of aerosolisable medium and a method of operating an aerosol generating apparatus |
-
2020
- 2020-10-30 KR KR1020200142924A patent/KR102640561B1/en active IP Right Grant
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2008108147A (en) * | 2006-10-26 | 2008-05-08 | Riam-Compact Co Ltd | Fluid analysis support system and program |
KR101363964B1 (en) * | 2010-04-22 | 2014-02-18 | 윈챵 슈우 | Combined multifunctional electronic simulated cigarette |
Also Published As
Publication number | Publication date |
---|---|
KR20220057847A (en) | 2022-05-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Zhang et al. | Size-change and deposition of conventional and composite cigarette smoke particles during inhalation in a subject-specific airway model | |
Zhang et al. | Vapor deposition during cigarette smoke inhalation in a subject-specific human airway model | |
Walters et al. | Computational fluid dynamics simulations of particle deposition in large-scale, multigenerational lung models | |
Robinson et al. | A framework to investigate the impact of topography and product characteristics on electronic cigarette emissions | |
Asgari et al. | Multispecies aerosol evolution and deposition in a human respiratory tract cast model | |
Carr et al. | A variable-stepsize Jacobian-free exponential integrator for simulating transport in heterogeneous porous media: Application to wood drying | |
KR102640561B1 (en) | Apparatus for simulating aerosol and method thereof | |
KR102640559B1 (en) | Apparatus for simulating aerosol and method thereof | |
KR102640560B1 (en) | Apparatus for simulating aerosol and method thereof | |
Bao et al. | Thermal and moisture analysis for tobacco leaf flue-curing with heat pump technology | |
CN110121275A (en) | Electrically operated aerosol generates system | |
CN105138736A (en) | Numerical simulation method of lignite microwave drying furnace flow field | |
Talih et al. | Carbonyl emissions and heating temperatures across 75 nominally identical electronic nicotine delivery system products: do manufacturing variations drive pulmonary toxicant exposure? | |
Panigrahi et al. | A 3D transient CFD model to predict heat and moisture transfer in on-farm stored grain silo through parallel computing using compiler directives: Impact of discretization methods on solution efficacy | |
Asgharian et al. | Simulation modeling of air and droplet temperatures in the human respiratory tract for inhaled tobacco products | |
Qi et al. | Mechanisms of trace water vapor desublimation over airfoil in transonic cryogenic wind tunnels | |
Forest et al. | Considerations on dosimetry for in vitro assessment of e-cigarette toxicity | |
Tang et al. | Research on cigarette during smoking based on reverse engineering and numerical simulation | |
Kulkarni et al. | High-temperature effects on the mucus layers in a realistic human upper airway model | |
Sun et al. | An improved CFD model of gas flow and particle interception in a fiber material | |
JPH10132807A (en) | Evaluation method for smoke of cigarette | |
Kim et al. | Numerical analysis of the coupling between heat transfer and pyrolysis in heat-not-burn tobacco using computational fluid dynamics | |
CN112504902A (en) | Electronic cigarette expansion experiment table and building method and testing method thereof | |
Zhan et al. | Design and development of CYO with adjustable smoke volume | |
CN110793669B (en) | Method for rapidly measuring temperature of filter stick in suction process |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
E902 | Notification of reason for refusal | ||
E902 | Notification of reason for refusal | ||
E90F | Notification of reason for final refusal | ||
E701 | Decision to grant or registration of patent right | ||
GRNT | Written decision to grant |