CFDEM DEM DEM(MPI) LIGGGHTS CFD CFD

5) 5) 5)

11) 10) β D n = βd (1) D n β D 10) 10) β = 0.2 0.5 β β β = 0.2 0.5 β = 0.2 β = 0.5

35 30 25 ( ) 20 15 10 5 0 0 0.1 0.2 0.3 0.4 0.5 0.6 β (-) β β 1) Zhu, H.P., Zhou, Z.Y., Yang, R.Y., Yu, A.B. Discrete particle simulation of particulate systems: A review of major applications and findings, Chemical Engineering Science, 63, (2008), 5728. 2) Goniva, C., Kloss, C., Deen, N.G., Kuipers, J.A.M., Pirker, S.Influence of rolling friction modelling on single spout fluidized bed simulations, Particuology, 10, 5, (2012). 3) Weller, H. G., Tabor, G., Jasak, H., Fureby, C.A tensorial approach to computational continuum mechanics using object-oriented techniques, Computers in Physics, 12(6), (1998), 620. 4) Kloss, C., Goniva, C., Hager, A., Amberger, S., Stefan Pirker, S.Models, algorithms and validation for opensource DEM and CFD-DEM, Progress in Computational Fluid Dynamics, An Int. J. 12, (2012), 140. 5) Feng, Y. Q., Yu, A. B.Assessment of model formulations in the discrete particle Simulation of gas-solid flow, Industrial and Engineering Chemistry Research, 43, (2004), 8378. 6) Cleary, P.W.Industrial particle flow modeling using discrete element method, Eng. Comput., 26, 6 (2009), 698. 7) Lillie, C., Wriggers, P.Three-dimensional modeling of discrete particles by superellipsoids, Proc. Appl. Math. Mech., 6 (2006), 101. 8) Zhao, D., Nezami, E. G., Hashash, Y. M. A. Three-dimensional discrete element simulation for granular materials., Eng. Comput. Int. J. Comput. Aided Eng. Softw., 23, 7 (2006), 749. 9) Bierwisch, C., Kübler, R., Kleer, G., Moseler, M.Modelling of contact regimes in wire sawing with dissipative particle dynamics, Phil. Trans. R. Soc. A, 369 (2011), 2422. 10) 10 (2011), 43. 11) 59A (2013), 208.

( ( ( (Discrete) v d x F( x 3 ) 1 2 x x d v d u d x d v d 1 12,2017

pbesti gbest v i1 wvi r1 c1( pbesti xi ) r c ( gbest x ) 2 2 x x v i1 i i1 x i x i1 xi xi i r 0.5 r 0.5 E A exp T Ti 1 T i 1 x1x 2 f ( x) min 6.931 x3x4 12 x 60 i 1,2,3,4 i x1, x2, x3, x4 16,19,43,49 19,16,43,49 16,19,49,43 19,16,43,49 1) Sandgren, E., Nonliear Integer and Discrete Programming in Mechanical Design Optimization, ASME/J. of Mechanical Design, Vol.112, (1990), pp.223-229. 2 7 12,2017

FluentANSYS HyperStudy Altair 2 Fluent Fluent HyperStudy Fluent HyperStudy 1 Gambit Fluent HyperMesh HyperStudy 1. 2 HyperMesh 2 3 1000 3 500mm 8 12

2. HyperStudy 9 3 1000mm 500mm 125mm 4 1000mm ARSM Adaptive Response Surface Method ARSM 4 1 ARSM 5 ARSM 998.6mm 12.188Pa 5. ARSM 500mm 125mm 3. Altair HyperStudy 4. 12

12

0.687 C 24(1 0.15Re ) / Re Re 1000 D 0.44 Re 1000

SSB 3) Q=FV [Nm/s] =FV [Watt] FLUX MARC F V 3 4 Q 3 R 4000 GLUE F V Q 12 12

Q Q =FV Watt H p H=RH0.01[W/( 2 K)] H=1.0[W/( 2 K)] 8 RH H 4) 1), 405 2016 p124-128, NS-SSB 2),Vol.6 2015,p28-35, NS-SSB 3),2015,p465-466,(SSB) 4),2009,p17 13 12

u J ( Wdy T ds) x 1 J 2-1 J A u q W1 i ) x j 1 ij x1 i ( da 12

2-2 MARC 2-3 2-4 J 2-5 J K E Bb (3 =1.9) (3) 12

1) 20 2m 2m 16 12

1000 900 800 mm/sec 700 600 500 400 300 200 100 0 0 0.5 1 1.5 2 2.5 3 sec 1.6m 6.6m 2.3m 0.3m 0.9m 1) Abaqus 6.14 17 12

MATLAB Simulink - etc etc x =u y=x+c Qin T Qout Text Qin [W] 5000 Tgoal [] 60 Text [] 20 rho [kg/l] 1 Cp [J/kgK] 4217 V [L] 5 A [m 2 ] 2 T0 [] 20 s Alpha(s) Qout dq Q T [W/m 2 K] =-0.0875s 2 + 7.875s + 23.75 [W] =AAlpha(s)(T-Text) =Qin-Qout Tgoal [] = dqdt+t0rhocpv [] =Q/(rhoCpV) 18 12

-2 Simulink - Heat Model Heat Model Level2 MATLAB S-Function API *1 Simulink MATLAB Simulink MATLAB MATLAB MATLAB *1 Application Program Interface: Level2 MATLAB S-Function Simulink Level2 MATLAB S-Function % function HeatModel( block ) block.regblockmethod('initializeconditions', @InitConditions); % block.regblockmethod('outputs', @Output); % block.regblockmethod('derivatives', @Derivative); % end % function InitConditions(block) % Q Q=T0*rho*Cp*V; block.contstates.data=q; end % function Derivative (block) % dq Alpha=@(s) -0.0875*s^2 + 7.875*s + 23.75; Qin=block.InputPort(2).Data; s=block.inputport(1).data; Qout=A*Alpha(s)*(T-Text); dq=qin-qout; block.derivatives.data=dq; end % function Output (block) % Q T Q=block.ContStates.Data T=Q/(rho*Cp*V); block.outputport(1).data =T; end Level2 MATLAB S-Function - 180 300 [] Level2 MATLAB S-Function 19 12

Drag&Drop 1 2

2 0M1000 1000 7.86 10 9.175 10 0.722 1 V=0.0m/s 2 V=1.0m/s 5.0 10 A. 1 2 5.0 10 0.727 B. 1 2 TEL: 093-588-7233 E-mailnpd.cae-sales@eng.nssmc.com

22 12