1
0
Fork 0
fgdata/Aircraft/c172p/dev/fdm/gerris/buoyancy3D.gfs
Stuart Buchanan 60d7073794 Initial checking of c172p-detailed.
From https://github.com/Juanvvc/c172p-detailed.git

Commit df6b2c5f4f75d000c5f5ef9c6cc4446784348f20
2015-06-16 22:20:54 +01:00

82 lines
2.4 KiB
Text

1 0 GfsSimulation GfsBox GfsGEdge {} {
# One object immersed in a fluid.
# Coordinate system
# X-axis goes towards the back.
# Y-axis goes to the left looking back.
# Z-axis goes up.
# Origo is located in the center of the domain 0.5 units behind the
# the initial boundary face.
# Simulation scaling parameters
# Reference length (GfsBox side) L = 10 m
# Reference speed U = 0 m/s
# Reference density Rho = 1000 kg/m^3
# Reference gravity G = 9.81 m/s^2
#
# length scale 1/L = 0.10
# time scale L/U = inf
# Force and moment scaling.
# F_r = F_s * U^3 * Rho * G
# M_r = M_s * U^4 * Rho * G
# We don't care about Reynolds number or viscosity for this scenario.
# We are interested in the static situation for buoyancy computation.
Time { iend = 1 }
Global {
#define VAR(T,min,max) (min + CLAMP(T,0,1)*(max - min))
#define RATIO (1.2/1000.)
}
ApproxProjectionParams { tolerance = 1e-12 }
ProjectionParams { tolerance = 1e-12 }
# Keep the domain simple.
Refine 3
# Fluid density.
Init {} { rho = 1.0 }
# Gravity.
Source W -rho
# Viscosity is uninteresting for buoyancy.
#GfsSourceViscosity 1e-3
# FIXME: Make sure moments and rotations are around HYDRORP.
# Is it really ok to set the position and orientation of the solid
# by moving its origin?
Solid floats.gts { scale=0.10 }
RefineSolid 9
# Output pressure force and moment on the object
OutputSolidForce { istart = 1 istep = 1 } f
# Refine the water surface to four levels
RefineSurface { return 4; } (1e-4 - z)
VariableTracerVOF T
# For high-density ratios we cannot use the volume fraction field
# directly to define the density. We need a smoother version.
VariableFiltered T1 T 1
# Initialise the water surface at z = 1e-4 (i.e. pretty much in the middle)
InitFraction T (1e-4 - z)
# air/water density ratio
PhysicalParams { alpha = 1./VAR(T1,RATIO,1.) }
OutputTime { istep = 1 } stderr
OutputBalance { istep = 1 } stderr
OutputProjectionStats { istep = 1 } stderr
OutputTiming { istep = 10 } stderr
# Output and compress the saved simulation files
OutputSimulation { start = 0 step = 1.0 } sim-%g.gfs
EventScript { start = 0 step = 1 } { gzip -f -q sim-*.gfs }
}
GfsBox {}