I am using a cylinder region to add external force in a domain. The script for the domain is a s follows
region a cylinder z -0.343 -0.004 0.6 -0.4025 -0.2375 units box
fix push all addforce 0.0 0.0 -0.1 region a
Is there any error in the script. Please help me out.
Thanks in advance.
| Attachment | Size |
|---|---|
 Please rename it to Frustum-1.stl after downloading | 45.67 KB |
ckloss | Thu, 10/30/2014 - 14:44
Hi b.saprativ,
Hi b.saprativ,
can you post a full script so one can run it?
Cheers
Christoph
b.saprativ | Fri, 10/31/2014 - 11:08
input script for the flowfield
#Problem Definition/Name
######################
# Physics/Generic Details
######################
units si
atom_style sphere
atom_modify map array sort 0 0
newton off
communicate single vel yes
processors 2 2 2
neighbor 0.001 bin
neigh_modify delay 0
boundary f f f
region reg block -0.11 0.11 -0.11 0.11 -0.0525 1.1025 units box
create_box 2 reg
region aa cylinder z 0.1 0 0 0.1 0 1 units box
fix push all addforce 0.0 0.0 -0.2 region aa
fix 1 all nve/sphere
fix gravi all gravity 9.81 vector 0.0 0.0 -1.0
pair_style gran model hertz tangential history cohesion sjkr rolling_friction cdt
pair_coeff * *
######################
# Material Data
######################
fix m1 all property/global youngsModulus peratomtype 5000000.0 70000000000.0
fix m2 all property/global poissonsRatio peratomtype 0.45 0.45
#Material properties required for new pair styles
fix m3 all property/global coefficientRestitution peratomtypepair 2 0.5 0.5 0.5 0.5
fix m4 all property/global coefficientFriction peratomtypepair 2 0.5 0.5 0.5 0.5
fix m5 all property/global coefficientRollingFriction peratomtypepair 2 0.0 0.0 0.0 0.0
fix m7 all property/global cohesionEnergyDensity peratomtypepair 2 0.0 0.0 0.0 0.0
######################
# Geometry Details
######################
fix cad1 all mesh/surface file Frustum-1.stl type 2 scale 1.0 curvature 1e-05
fix granwall all wall/gran model hertz tangential history mesh n_meshes 1 meshes cad1
####################
# Liquid bridging
####################
######################
# Particle Details
######################
#distributions for injection
fix pts1_1 all particletemplate/sphere 1 atom_type 1 density constant &
2000.0 radius constant 0.001
fix pdd1 all particledistribution/discrete 1. 1 pts1_1 1.0
######################
# Factory Details
######################
#region and insertion
region fact1 cylinder z 0.0 0.0 0.08 0.95 1.05 units box
######################
# Recipe Details
######################
timestep 2e-05
fix inj1 all insert/rate/region seed 1234 distributiontemplate pdd1 &
maxattempt 10 nparticles 10000.0 massrate 0.837758040957 all_in yes &
insert_every 1000 overlapcheck yes &
vel constant 0.0 0.0 -0.1 omega constant 0.0 0.0 0.0 region fact1 ntry_mc 10000
######################
# Applied flow field
######################
######################
# Output Controls
######################
restart 1000 Restart/test_flow.*.res
#output settings, include total thermal energy
thermo_style custom step atoms ke cpu
thermo 1000
thermo_modify lost ignore norm no
compute_modify thermo_temp dynamic yes
dump dmp all custom 1000 Post/dump.test_flow_ &
id type x y z ix iy iz vx vy vz fx fy fz omegax omegay omegaz radius
dump dmp_geom1 all mesh/vtk 1000 Post/test_flow_cad1_*.vtk id cad1
dump dmp_part all atom/vtk 1000 Post/test_flow_*.vtu
#######################
# Events
#######################
run 500000
######################
# Solve
######################
write_restart test_flow.res
b.saprativ | Fri, 10/31/2014 - 12:33
Additional Info
Hi,
The input script requires a Frustum-1.stl file which I have attached as a txt file in the main post. Please rename it to Frustum-1.stl after downloading