/ Journal File for GAMBIT 2.1.6
/ File opened for write Sat Jan 22 16:44:49 2005.
identifier name "Basicwing2" new nosaveprevious

/Parameters for Geometry
/Case 
$Alpha = 0

/Domain size (in chord lengths)

$Upstream=3
$Downstream=5
$Height=10
$Width=5
$Base=-0.0075

/Parameters for Meshing
/Domain

$ynodes = 50
$znodes = 50
$xnodes_upstream = 30
$xnodes_downstream = 50
$Stretch = 0.01

/Delta Wing
$xnodes_wing = 50
$znodes_wing = 50
$xnodes_bevel = 20
$wing_first = 0.5
$wing_last = 0.2

/Constants
$Chord=24.9849

/Mathematics for case
/Case
$AoA = -1 * $Alpha
/Side volume
$LEside = $Chord * $Width
$TEside = $LEside - 9.0938
/Front and rear volumes
$Front = -1 * $Chord * $Upstream
$Rear = $Chord * $Downstream
$High = $Chord * $Height * 0.5
$Low = -1 * $High
$High_front = $High - $Chord * SIN(-1 * $AoA)
$Low_front = $Low - $Chord * SIN(-1 * $AoA)
/Volumes above and below wing
$TE_length = SQRT( 0.1875 * 0.1875 + 0.5 * 0.5 )
$dx = $TE_length * COS (20.55605 - $AoA)
$x7 = $Chord-$dx
$x8 = $x7 - ($Chord - 1.4619 - 0.5) * COS(-1 * $AoA)

/********Volume generation***************

/Wing Vertices at 0deg Alpha
vertex create "p1" coordinates 0 0 0
vertex create "p2" coordinates 1.4619 0.1875 0
vertex create "p3" coordinates 24.4849 0.1875 0
vertex create "p4" coordinates 24.9849 0 0
vertex create "p5" coordinates 24.9849 0 9.0938
vertex create "p6" coordinates 24.4849 0.1875 8.37970
/Wing edges
edge create "e12" straight "p1" "p2"
edge create "e23" straight "p2" "p3"
edge create "e34" straight "p3" "p4"
edge create "e14" straight "p1" "p4"
edge create "e45" straight "p4" "p5"
edge create "e36" straight "p3" "p6"
edge create "e56" straight "p5" "p6"
edge create "e26" straight "p2" "p6"
edge create "e15" straight "p1" "p5"
/Wing Faces
face create "Topface" wireframe "e23" "e36" "e26" real
face create "Botface" wireframe "e15" "e14" "e45" real
face create "LEface" wireframe "e12" "e15" "e56" "e26" real
face create "TEface" wireframe "e34" "e36" "e56" "e45" real
face create "Root" wireframe "e34" "e36" "e56" "e45" real
/Rotate via edge45 to $AoA angle
face move "Topface" "Botface" "LEface" "TEface" "Root" dangle $AoA vector 0 0 \
  9.0938 origin 24.9849 0 0
/Create Wingbase for H-Type topology
vertex cmove "p1" "p4" "p5" multiple 1 offset 0 $Base 0
/Renaming to fit my lables
vertex modify "vertex.7" label "p1L"
vertex modify "vertex.9" label "p5L"
vertex modify "vertex.8" label "p4L"
edge create "e15L" straight "p1L" "p5L"
edge create "e45L" straight "p4L" "p5L"
edge create "e41L" straight "p1L" "p4L"
edge create "e4L4" straight "p4L" "p4"
edge create "e5L5" straight "p5L" "p5"
edge create "e1L1" straight "p1L" "p1"
face create "LEbase" wireframe "e15" "e1L1" "e15L" "e5L5"
face create "TEbase" wireframe "e4L4" "e45" "e5L5" "e45L" real
face create "Botbase" wireframe "e41L" "e45L" "e15L" real

/Starting on thin plates
/Side of wing first, since that defines so much of the domain

face create translate "e1L1" vector 0 0 $LEside
edge modify "edge.17" label "e1side"
edge modify "edge.16" label "e1Lside"
edge modify "edge.18" label "e1L1side"
vertex modify "vertex.11" label "p1side"
vertex modify "vertex.10" label "p1Lside"
face modify "face.9" label "Wingsidefront"

face create translate "e5L5" vector 0 0 $TEside
edge modify "edge.20" label "e5side"
edge modify "edge.19" label "e5Lside"
edge modify "edge.21" label "e5L5side"
vertex modify "vertex.13" label "p5side"
vertex modify "vertex.12" label "p5Lside"
face modify "face.10" label "Wingsideback"
edge create "e15side" straight "p1side" "p5side"
edge create "e1L5Lside" straight "p1Lside" "p5Lside"
face create "BotbaseSide" wireframe "e15side" "e1L1side" "e1L5Lside" \
  "e5L5side" real
face create "WingsideTop" wireframe "e5side" "e15" "e1side" "e15side" real
face create "WingsideBot" wireframe "e1L5Lside" "e1Lside" "e15L" "e5Lside" \
  real
volume create "Wingside" stitch "Wingsideback" "BotbaseSide" "Wingsidefront" \
  "WingsideTop" "WingsideBot" "LEbase" real

/Sweeping volumes to create upstream and downstream volumes

volume create "Upstream_thin" translate "Wingsidefront" vector $Front 0 0
volume create "Downstream__thin_side" translate "Wingsideback" vector $Rear 0 0
volume create "Downstream_thin" translate "TEbase" vector $Rear 0 0

/Sweeping Rear Volumes

volume create "Downstream_top" translate "face.27" "face.22" vector 0 $High 0
volume create "Downstream_bot" translate "face.25" "face.20" vector 0 $Low 0 

/Sweeping Front Volume
volume create "Upstream_top" translate "face.17" vector 0 $High_front 0
volume create "Upstream_bot" translate "face.15" vector 0 $Low_front 0 

/ Creating vertices for Volume directly above wing

vertex create "p7" coordinates $x7 $High 0
vertex create "p8" coordinates $x8 $High 0
vertex create "p9" coordinates $x7 $High 8.37970
/ Creating edges and faces for Volume about wing
edge create "e78" straight "p7" "p8"
edge create "e89" straight "p8" "p9"
edge create "e79" straight "p7" "p9"
edge create "e12top" straight "vertex.41" "p8"
edge create "e15top" straight "vertex.41" "vertex.26"
edge create "e56top" straight "vertex.26" "p9"
edge create "e34top" straight "p7" "vertex.24"
edge create "e37" straight "p3" "p7"
edge create "e69" straight "p6" "p9"
edge create "e18" straight "p2" "p8"

/ Creating horizontal top faces to build volume above wing
face create "TopFace_top" wireframe "e78" "e89" "e79" real
face create "LEface_top" wireframe "e89" "e12top" "e15top" "e56top" real
face create "TEface_top" wireframe "e79" "e56top" "edge.48" "e34top" real
/ Creating vertical side faces to build volume above wing
face create wireframe "e37" "e34top" "edge.45" "e34" real
face create wireframe "e37" "e79" "e69" "e36" real
face create wireframe "edge.50" "e56top" "e69" "e56" real
face create wireframe "e15top" "edge.78" "e15" "edge.50" real
face create wireframe "e89" "e18" "e26" "e69" real
face create wireframe "edge.78" "e12top" "e18" "e12" real
face create wireframe "e78" "e37" "e23" "e18" real
/Building Volumes
volume create "Wingtop_volume" stitch "TopFace_top" "face.67" "face.62" \
  "face.65" "Topface" real
volume create "WingTE_volume" stitch "face.29" "TEface_top" "face.62" \
  "face.61" "face.63" "TEface" real
volume create "WingLE_volume" stitch "face.65" "LEface_top" "face.64" \
  "face.66" "face.63" "LEface" real

/ Creating edges,faces and volumes for volume underneath the wing (nostep)

edge create straight "vertex.45" "vertex.34"
edge create straight "vertex.45" "vertex.32"
face create "Botface_bot" wireframe "edge.103" "edge.104" "edge.64" real
face create "Botface_bot" wireframe "e15L" "edge.87" "edge.103" "edge.67" \
  real
face create "Botface_bot" wireframe "edge.104" "edge.87" "e41L" "edge.62" \
  real
volume create "Wingbot_volume" stitch "Botbase" "face.40" "face.69" \
  "Botface_bot" "face.70" real

/ Creating volumes off to the side of the wing

edge create straight "vertex.43" "vertex.30"
edge create straight "vertex.47" "vertex.38"
face create wireframe "edge.81" "edge.105" "edge.56" "e15top" real
face create wireframe "edge.58" "edge.105" "edge.83" "e15side" real
face create wireframe "e1L5Lside" "edge.92" "edge.106" "edge.75" real
face create wireframe "edge.106" "edge.89" "edge.103" "edge.72" real
volume create "Wing_side_top_volume" stitch "face.71" "face.50" "face.34" \
  "face.72" "face.64" "WingsideTop" real
volume create "Wing_side_bot_volume" stitch "WingsideBot" "face.73" "face.56" \
  "face.45" "face.69" "face.74" real

/Cleaning up face used for construction
face delete "Botface" lowertopology
face delete "Root" lowertopology

/********The Volume is now complete***************

/********Starting edge meshing********************

/All the vertical edges first

edge mesh "edge.45" "edge.46" "edge.50" "edge.51" "edge.58" "edge.59" \
  "edge.82" "e37" "e69" "e18" "edge.77" "edge.78" "edge.83" "edge.62" \
  "edge.63" "edge.67" "edge.68" "edge.75" "edge.76" "edge.86" "edge.87" \
  "edge.91" "edge.92" firstlength ratio1 $Stretch intervals $ynodes

/Thinplate verticals (single node only)

edge mesh "edge.41" "edge.33" "e5L5" "e4L4" "edge.38" "e5L5side" "e1L1side" \
  "e1L1" "edge.30" "edge.25" successive ratio1 1 intervals 1

/Now for the horizontal edges

edge mesh "edge.57" "edge.81" "e1Lside" "e1side" \
  "e5Lside" "e5side" "edge.28" "edge.35" "edge.72" \
  "edge.88" firstlength ratio1 $Stretch  intervals $znodes

/Fixing horizontal ones that were backwards

edge delete "edge.80" "edge.56" "edge.27" "edge.89" "edge.73" "edge.36" \
  keepsettings onlymesh
edge modify "edge.80" "edge.56" "edge.27" "edge.89" "edge.73" "edge.36" \
  backward
edge mesh "edge.80" "edge.56" "edge.27" "edge.89" "edge.73" "edge.36" \
  firstlength ratio1 $Stretch intervals $znodes

/Downstream Volume (important for wake)

edge mesh "edge.47" "edge.52" "edge.60" "edge.32" "edge.34" "edge.40" \
  "edge.42" "edge.37" "edge.39" "edge.61" "edge.66" "edge.74" firstlength \
  ratio1 $Stretch intervals $xnodes_downstream
/Fixing horizontal ones that were backwards
edge delete "edge.52" "edge.61" "edge.60" keepsettings onlymesh
edge modify "edge.52" forward
edge modify "edge.61" backward
edge modify "edge.60" backward
edge mesh "edge.52" "edge.61" "edge.60" firstlength ratio1 $Stretch \
intervals $xnodes_downstream
/Horizontal edges at rear of thin wake volume
edge mesh "edge.43" "edge.44" successive ratio1 1 intervals $znodes_wing
/Horizontal edges in domain (tied to TE nodes)

edge mesh "edge.48" "e79" "edge.49" "edge.65" "edge.64" successive ratio1 1 \
  intervals $znodes_wing

/Upstream Volume

edge mesh "edge.24" "edge.26" "edge.29" "edge.31" "edge.79" "edge.84" \
  "edge.85" "edge.90" firstlength ratio1 $Stretch intervals $xnodes_upstream
/Fixing horizontal ones that were backwards
edge delete "edge.79" "edge.90" keepsettings onlymesh
edge modify "edge.79" "edge.90" backward
edge mesh "edge.79" "edge.90" firstlength ratio1 $Stretch intervals $xnodes_upstream

/Wing Edges (all double sided except TE)

edge mesh "e23" "e26" "e15" "e15L" "e41L" "e15side" "e1L5Lside" "e78" "e89" \
  "e15top" "edge.103" "edge.104" "edge.105" "edge.106" lastlength ratio1 \
$wing_first ratio2 $wing_last intervals $xnodes_wing
edge delete "e78" keepsettings onlymesh
edge modify "e78" backward
edge mesh "e78" lastlength ratio1 $wing_first ratio2 $wing_last \
intervals $xnodes_wing
/TE of the wing
edge mesh "e45" "e36" "e45L" successive ratio1 1 intervals $znodes_wing

/Beveled edges (no stretching, same number of nodes)

/LE 
edge mesh "e12top" "e12" successive ratio1 1 intervals $xnodes_bevel
/TE
edge mesh "e34top" "e56top" "e34" "e56" successive ratio1 1 intervals \
$xnodes_bevel


/********Starting Volume Meshing*********************

/Volume diretly above and below wing (Cooper/Hex)

volume mesh "Wingtop_volume" cooper
volume mesh "Wingbot_volume" cooper

/LE and TE Volumes above wing (Hex)

volume mesh "WingTE_volume" "WingLE_volume" map

/Thinplate volumes (Hex)

volume mesh "Downstream__thin_side" "Downstream_thin" "Wingside" \
  "Upstream_thin" map

/Rest of Domain (Hex)

volume mesh "Upstream_top" "Downstream_top" "volume.6" "Upstream_bot" \
  "Wing_side_bot_volume" "volume.8" "Wing_side_top_volume" \
  "Downstream_bot" map

/********Volume Meshing Complete*********************

/********Starting Setting B.C.***********************

solver select "FLUENT 5/6"

/Symmetry
physics create "Plane of Symmetry" btype "SYMMETRY" face "face.28" "face.38" \
  "face.24" "face.61" "face.67" "face.66" "face.48" "face.14" "face.53" \
  "face.70"

/Inlet

physics create "Inlet" btype "PRESSURE_FAR_FIELD" face "face.49" "face.16" \
  "face.55"

/Outlet

physics create "Outlet" btype "PRESSURE_FAR_FIELD" face "face.35" "face.46" \
  "face.30" "face.41" "face.26" "face.21"

/Pressure Far Field

physics create "Freestream" btype "PRESSURE_FAR_FIELD" face "face.36" \
  "face.37" "face.31" "face.71" "face.72" "LEface_top" "TopFace_top" \
  "TEface_top" "face.51" "face.52" "face.44" "face.47" "face.73" "face.74" \
  "face.54" "face.57" "face.39" "Botface_bot" "face.23" "BotbaseSide" \
  "face.18"

/Wing

physics create "Wing" btype "WALL" face "Topface" "LEface" "TEface" "LEbase" \
  "TEbase" "Botbase"

/********Running Face and Volume Checks***********************

volume check "Wingside" "Upstream_thin" "Downstream__thin_side" \
  "Downstream_thin" "Downstream_top" "volume.6" "Downstream_bot" "volume.8" \
  "Upstream_top" "Upstream_bot" "Wingtop_volume" "WingTE_volume" \
  "WingLE_volume" "Wingbot_volume" "Wing_side_top_volume" \
  "Wing_side_bot_volume" topology geometry
face check "Topface" "LEface" "TEface" "LEbase" "TEbase" "Botbase" \
  "Wingsidefront" "Wingsideback" "BotbaseSide" "WingsideTop" "WingsideBot" \
  "face.14" "face.15" "face.16" "face.17" "face.18" "face.19" "face.20" \
  "face.21" "face.22" "face.23" "face.24" "face.25" "face.26" "face.27" \
  "face.28" "face.29" "face.30" "face.31" "face.32" "face.34" "face.35" \
  "face.36" "face.37" "face.38" "face.39" "face.40" "face.41" "face.42" \
  "face.44" "face.45" "face.46" "face.47" "face.48" "face.49" "face.50" \
  "face.51" "face.52" "face.53" "face.54" "face.55" "face.56" "face.57" \
  "TopFace_top" "LEface_top" "TEface_top" "face.61" "face.62" "face.63" \
  "face.64" "face.65" "face.66" "face.67" "Botface_bot" "face.69" "face.70" \
  "face.71" "face.72" "face.73" "face.74" topology geometry