Skip to content

Examples update for gas-solid fluidized and spouted-bed #1370

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
merged 5 commits into from
Nov 20, 2024
Merged

Examples update for gas-solid fluidized and spouted-bed #1370

Show file tree
Hide file tree
Changes from 3 commits
Commits
Show all changes
5 commits
Select commit Hold shift + click to select a range
9e48471
Updates gas-solid-fluidized-bed
OresteMarquis Nov 18, 2024
83c4583
update doc gas-solid-spouted-bed
OresteMarquis Nov 18, 2024
bdaae50
Prmindent on files changed
OresteMarquis Nov 18, 2024
c1feda2
Update following review comments
OresteMarquis Nov 19, 2024
f5e0882
Merge branch 'master' into examples_gas_solid
OresteMarquis Nov 19, 2024
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
71 changes: 38 additions & 33 deletions ...examples/unresolved-cfd-dem/gas-solid-fluidized-bed/gas-solid-fluidized-bed.rst
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -148,8 +148,8 @@ The ``insertion info`` subsection manages the insertion of particles. It allows

subsection insertion info
set insertion method = volume
set inserted number of particles at each time step = 2000
set insertion frequency = 2000
set inserted number of particles at each time step = 30000
set insertion frequency = 100000
set insertion box points coordinates = -0.018, -0.05, -0.018 : 0.018, 0.05, 0.018
set insertion distance threshold = 1.5
set insertion maximum offset = 0.2
Expand Down Expand Up @@ -224,11 +224,12 @@ The simulation is run for 1 s with a time step of 0.002 s. The time scheme chose

subsection simulation control
set method = bdf1
set output name = result_
set number mesh adapt = 0
set output name = cfd_dem
set output frequency = 10
set startup time scaling = 0.6
set time end = 1
set time step = 0.002
set time step = 0.001
set output path = ./output/
end

Expand Down Expand Up @@ -268,38 +269,42 @@ For the boundary conditions, we choose a slip boundary condition on the walls of

.. code-block:: text

subsection boundary conditions
set number = 5
subsection bc 0
set id = 0
set type = slip
end
subsection bc 1
set id = 1
set type = slip
subsection boundary conditions
set number = 6
subsection bc 0
set id = 0
set type = slip
end
subsection bc 1
set id = 1
set type = slip
end
subsection bc 2
set id = 2
set type = function
subsection u
set Function expression = 0
end
subsection bc 2
set id = 4
set type = slip
subsection v
set Function expression = 2
end
subsection bc 3
set id = 5
set type = slip
end
subsection bc 4
set id = 2
set type = function
subsection u
set Function expression = 0
end
subsection v
set Function expression = 2
end
subsection w
set Function expression = 0
end
subsection w
set Function expression = 0
end
end
subsection bc 3
set id = 3
set type = outlet
end
subsection bc 4
set id = 4
set type = slip
end
subsection bc 5
set id = 5
set type = slip
end
end

The additional sections for the CFD-DEM simulations are the void fraction subsection and the CFD-DEM subsection. These subsections are described in detail in the `CFD-DEM parameters <../../../parameters/unresolved-cfd-dem/unresolved-cfd-dem.html>`_ .

Expand All @@ -311,7 +316,7 @@ Since we are calculating the void fraction using the packed bed of the DEM simul
.. code-block:: text

subsection void fraction
set mode = dem
set mode = pcm
set read dem = true
set dem file name = dem
set l2 smoothing factor = 0.000005
Expand Down
15 changes: 5 additions & 10 deletions ...rce/examples/unresolved-cfd-dem/gas-solid-spouted-bed/gas-solid-spouted-bed.rst
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -130,13 +130,13 @@ The physical properties section of the particles allows us to specify the differ
set diameter = 0.0025
set number = 31050
set density particles = 2526
set young modulus particles = 10000000
set young modulus particles = 1000000
set poisson ratio particles = 0.25
set restitution coefficient particles = 0.97
set friction coefficient particles = 0.4
set rolling friction particles = 0.3
end
set young modulus wall = 10000000
set young modulus wall = 1000000
set poisson ratio wall = 0.25
set restitution coefficient wall = 0.33
set friction coefficient wall = 0.2
Expand All @@ -154,12 +154,7 @@ The ``insertion info`` subsection manages the insertion of particles. It allows
set insertion method = volume
set inserted number of particles at each time step = 31050
set insertion frequency = 2000
set insertion box minimum x = -0.075
set insertion box minimum y = 0.0
set insertion box minimum z = 0
set insertion box maximum x = 0.075
set insertion box maximum y = 0.3
set insertion box maximum z = 0.015
set insertion box points coordinates = -0.075, 0.0, 0 : 0.075, 0.3, 0.015
set insertion distance threshold = 1.05
set insertion maximum offset = 0.3
set insertion prn seed = 19
Expand Down Expand Up @@ -224,7 +219,7 @@ The simulation is run for 5 s with a time step of 0.0001 s. The time scheme chos
subsection simulation control
set method = bdf1
set number mesh adapt = 0
set output name = result_
set output name = cfd-dem
set output frequency = 50
set startup time scaling = 0.6
set time end = 5
Expand Down Expand Up @@ -346,7 +341,6 @@ We also enable grad-div stabilization in order to improve local mass conservatio
set pressure force = true
set saffman lift force = false
set drag model = rong
set post processing = true
set coupling frequency = 100
set implicit stabilization = false
set grad-div length scale = 0.005
Expand Down Expand Up @@ -390,6 +384,7 @@ Linear Solver
set ilu preconditioner absolute tolerance = 1e-12
set ilu preconditioner relative tolerance = 1
set verbosity = verbose
set max krylov vectors = 200
end
end

Expand Down
50 changes: 29 additions & 21 deletions examples/unresolved-cfd-dem/gas-solid-fluidized-bed/gas-solid-fluidized-bed.prm
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -11,13 +11,14 @@ set dimension = 3
#---------------------------------------------------

subsection simulation control
set method = bdf1
set number mesh adapt = 0
set output name = cfd_dem
set output frequency = 10
set time end = 1.
set time step = 0.001
set output path = ./output/
set method = bdf1
set number mesh adapt = 0
set output name = cfd_dem
set output frequency = 10
set startup time scaling = 0.6
set time end = 1.
set time step = 0.001
set output path = ./output/
end

#---------------------------------------------------
Expand Down Expand Up @@ -62,17 +63,23 @@ subsection void fraction
set l2 smoothing factor = 0.000005
set l2 lower bound = 0
set l2 upper bound = 1
set bound void fraction = false
set bound void fraction = true
end

#---------------------------------------------------
# CFD-DEM
#---------------------------------------------------

subsection cfd-dem
set grad div = true
set drag model = difelice
set coupling frequency = 40
set grad div = true
set void fraction time derivative = true
set drag force = true
set buoyancy force = true
set shear force = false
set pressure force = false
set drag model = difelice
set coupling frequency = 100
set vans model = modelB
end

#---------------------------------------------------
Expand Down Expand Up @@ -100,14 +107,6 @@ subsection boundary conditions
set type = slip
end
subsection bc 2
set id = 4
set type = slip
end
subsection bc 3
set id = 5
set type = slip
end
subsection bc 4
set id = 2
set type = function
subsection u
Expand All @@ -120,10 +119,18 @@ subsection boundary conditions
set Function expression = 0
end
end
subsection bc 5
subsection bc 3
set id = 3
set type = outlet
end
subsection bc 4
set id = 4
set type = slip
end
subsection bc 5
set id = 5
set type = slip
end
end

#---------------------------------------------------
Expand All @@ -149,7 +156,7 @@ end
subsection model parameters
subsection contact detection
set contact detection method = dynamic
set neighborhood threshold = 1.8
set neighborhood threshold = 1.3
end
set particle particle contact force method = hertz_mindlin_limit_overlap
set particle wall contact force method = nonlinear
Expand Down Expand Up @@ -232,5 +239,6 @@ subsection linear solver
set ilu preconditioner absolute tolerance = 1e-14
set ilu preconditioner relative tolerance = 1.00
set verbosity = verbose
set max krylov vectors = 200
end
end
12 changes: 6 additions & 6 deletions examples/unresolved-cfd-dem/gas-solid-spouted-bed/gas-solid-spouted-bed.prm
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -13,11 +13,11 @@ set dimension = 3
subsection simulation control
set method = bdf1
set number mesh adapt = 0
set output name = result
set output name = cfd-dem
set output frequency = 50
set startup time scaling = 0.6
set time end = 5
set time step = 0.001
set time step = 0.0001
set subdivision = 1
set log precision = 10
set output path = ./output/
Expand Down Expand Up @@ -144,10 +144,6 @@ subsection boundary conditions
set type = slip
end
subsection bc 1
set id = 2
set type = outlet
end
subsection bc 2
set id = 1
set type = function
subsection u
Expand All @@ -160,6 +156,10 @@ subsection boundary conditions
set Function expression = 0
end
end
subsection bc 2
set id = 2
set type = outlet
end
subsection bc 3
set id = 3
set type = function
Expand Down
Loading