Scribd
Upload
123 views10 pages

QQQQ

This document discusses using discrete element method (DEM) simulations in OpenFOAM to analyze the destruction of a solid snow block. It provides details on setting up DEM simulations in OpenFOAM, including defining particle properties and positions, simulation parameters, and post-processing results. It also references an example DEM case study of snow destruction and lists the basic steps of mesh generation, running the simulation, and visualizing results.

Uploaded by

Miguel Rumsfeld
Copyright
© Attribution Non-Commercial (BY-NC)
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
123 views10 pages

QQQQ

This document discusses using discrete element method (DEM) simulations in OpenFOAM to analyze the destruction of a solid snow block. It provides details on setting up DEM simulations in OpenFOAM, including defining particle properties and positions, simulation parameters, and post-processing results. It also references an example DEM case study of snow destruction and lists the basic steps of mesh generation, running the simulation, and visualizing results.

Uploaded by

Miguel Rumsfeld
Copyright
© Attribution Non-Commercial (BY-NC)
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 10

OpenFOAM Seminar for beginner 13th Interim report

Destruction of solid with DEM

Takahiro Shibata
[email protected]

1, Introduction
I want to analyse destruction of solid(snow block) with DEM(Discrete Element Method).
snow block

2, DEM(Discrete Element Method)


DEM is numerical methods for computing the motion of a large number of particles of micrometre-scale size and above.

http://www.openfoam.com/version2.0.0/lagrangian.php

3, CaseFile
Reffering to Mr. Kasugas PENGUINITIS

http://www.geocities.jp/penguinitis2002/study/OpenFOAM/dem/dem.html

square system controlDict fvSchemes fvSolution constant RASProperties g kinematicCloudPositions kinematicCloudProperties transportProperties turbulenceProperties polyMesh blockMeshDict boundary 0 Initial condtions & size of particles Parameter of DEM

kinematicCloudPositions
FoamFile { version format class object } // * * * * *

100 particles
2.0; ascii; vectorField; k inematicCloudPositions; * * * * * * * * * * * * * * * * * * * * //

diam eter 0.08

( (0.04595 0.07097 0.0051) (0.04709 0.0709 0.00497) (0.04792 0.07099 0.00497) (0.04894 0.07094 0.00499) (0.05009 0.07105 0.005) (0.05098 0.07106 0.00508) (0.05204 0.07102 0.00494) (0.05506 0.07998 0.0049) )

0.8 0.7

0.45

0.55

10 particles

10 particles

kinematicCloudProperties

FoamFile { version format class location object } // * * * * * subModels {

2.0; ascii; dictionary; "constant"; reactingCloud1Properties; * * * * * * * * * * * * * * * * * * * * * * * * * //

pairCollisionCoeffs { // Maximum possible particle diameter expected at any time maxInteractionDistance 0.0008; Size of particles writeReferredParticleCloud no; pairModel pairSpringSliderDashpot; pairSpringSliderDashpotCoeffs { useEquivalentSize no; alpha 0.12; b 1.5; mu 0.52; cohesionEnergyDensity 0; collisionResolutionSteps 12; }; wallModel

Param eter of DEM

manualInjectionCoeffs { sizeDistribution { type fixedValue; fixedValueDistribution { value 0.0008; } } }

wallLocalSpringSliderDashpot;

Size of particles

wallLocalSpringSliderDashpotCoeffs { useEquivalentSize no; collisionResolutionSteps 12; movingW all { youngsModulus 1e10; poissonsRatio 0.23; alpha 0.12; b 1.5; mu 0.43; } fixedW alls { youngsModulus 1e10; poissonsRatio 0.23; alpha 0.12; b 1.5; mu 0.43; } frontAndBack { youngsModulus 1e10; poissonsRatio 0.23; alpha 0.12; b 1.5; mu 0.1; } };

Param eter of DEM

standardW allInteractionCoeffs { type rebound; }

cloudFunctions {}

4, Calculation
Mesh generation
$ blockMesh

Running the code


$ icoUncoupledKinematicParcelFoam

Post-processing
$ foam ToVTK $ paraview

Open square/VTK/square_.vtk Open square/VTK/lagrangian/kinematicCloud /kinematicCloud_.vtk Glyph

Glyph Type Sphere

Radius 0.52

Set Scale Factor 0.002

Case1 Regular arrangement Case2 Irregular arrangement

You might also like