Windtunnel/Boundary Conditions: Difference between revisions

From Arbeitsgruppe Kuiper
Jump to navigation Jump to search
(moved from WT main page)
(+keywords)
Line 1: Line 1:
[[Windtunnel online|up]]
[[Windtunnel online|up]]


= New Keywords =


= Boundary condition for the walls =
== in pluto.ini ==
 
* <code>ObjectType_int</code>: 0=nothing, 1=cylinder/sphere, 2=square/cube of size <code>ObjectDiameter_cm</code>
* <code>WindPressure_Pa</code>, <code>WindPressure_mbar</code> (mutually exclusive): pressure at entry (left)
* <code>PressureGradient_mbar_per_cm</code> positive (even though pressure drops from left to right)
* Wind velocity is ramped up to <code>WindVelocity_m_per_s</code> over time <code>InjectionTime_s</code>
* <code>WindTemperature_C</code>
* <code>KinematicViscosity_m2_per_s</code>, <code>DynamicViscosity_Pa_s</code>,
* <code>Wall_BoundaryCondition</code> for tangential walls, value from {0,1,2,3}
 
(ignored if given values <0)
 
== in init.c ==
 
* macro <code>SOLID</code> to enable an object
* array <code>solid</code> containing value 1 or 0 for object matter being present/absent
 
= Walls =
[[File:WT- Boundary Condition Skizze.jpg|thumb]]
[[File:WT- Boundary Condition Skizze.jpg|thumb]]
== tangential ==
{| class="wikitable" id="ganz_oben"
{| class="wikitable" id="ganz_oben"
|+ id="Ü-id" |  
|+ id="Ü-id" |  
Line 20: Line 41:
|}
|}


= Boundary condition for the object =
== entry/exit ==
 
=== prescribed velocity at entry (left) ===
 
=== prescribe pressure drop (left to right) ===
 
* current simulation result: v-profile slightly asymmetric


= Object =


= Open Questions =
= Open Questions =

Revision as of 09:03, 18 September 2024

up

New Keywords

in pluto.ini

  • ObjectType_int: 0=nothing, 1=cylinder/sphere, 2=square/cube of size ObjectDiameter_cm
  • WindPressure_Pa, WindPressure_mbar (mutually exclusive): pressure at entry (left)
  • PressureGradient_mbar_per_cm positive (even though pressure drops from left to right)
  • Wind velocity is ramped up to WindVelocity_m_per_s over time InjectionTime_s
  • WindTemperature_C
  • KinematicViscosity_m2_per_s, DynamicViscosity_Pa_s,
  • Wall_BoundaryCondition for tangential walls, value from {0,1,2,3}

(ignored if given values <0)

in init.c

  • macro SOLID to enable an object
  • array solid containing value 1 or 0 for object matter being present/absent

Walls

WT- Boundary Condition Skizze.jpg

tangential

Boundary condition for the wall
vx1 vx2
no-shear zero-gradient reflective
no-slip reflective reflective
no-wall zero-gradient zero-gradient
one-way wall zero-gradient zero-gradient & no-inflow

entry/exit

prescribed velocity at entry (left)

prescribe pressure drop (left to right)

  • current simulation result: v-profile slightly asymmetric

Object

Open Questions

Reflective cells

Setting the vector \(\vec v_g\) in a ghost cell as \(\vec v_g=-\vec v\), with \(\vec v\) being the value in the adjacent real cell, yields \(\vec 0\) as interpolation right at the boundary. This works for plane walls. What to do in cases where a ghost cell has more than one real cell as nearest neighbor? This calls for discussion!