Boundary conditions in computational fluid dynamics: Difference between revisions

Revision as of 10:14, 8 August 2016

Latest revision as of 14:03, 19 July 2025

(12 intermediate revisions by 7 users not shown)

Line 1:

{{Main|Boundary conditions in fluid dynamics}}

[[File:Fig 1 Formation of grid in cfd.JPG|thumb|Fig 1 Formation of grid in cfd]]Almost every [[computational fluid dynamics]] problem is defined under the limits of initial and boundary conditions. For implementation of boundary conditions when we construct a staggered grid we add an extra node across the physical boundary in order to get,

[[File:Fig 1 Formation of grid in cfd.JPG|thumb|Fig 1 Formation of grid in cfd]]Almost every [[computational fluid dynamics]] problem is defined under the limits of initial and boundary conditions. When constructing a staggered grid, it is common to implement boundary conditions by adding an extra node across the physical boundary. The nodes just outside the inlet of the system are used to assign the inlet conditions and the physical boundaries can coincide with the [[scalar control]] volume boundaries. This makes it possible to introduce the boundary conditions and achieve discrete equations for nodes near the boundaries with small modifications.

'''MostThe most common boundary conditions used in [[computational fluid dynamics]] are'''

*The nodes just outside the inlet of the system are used to assign the inlet conditions

*The physical boundaries can coincide with the scalar control volume boundaries.

This allow us to introduce the boundary conditions and achieve discretion equations for nodes near boundary with small modifications.

'''Most common boundary conditions used in [[computational fluid dynamics]] are'''

*Intake conditions

Line 12 ⟶ 8:

*Physical boundary conditions

*Cyclic conditions

*pressurePressure conditions

*exitExit conditions

==Intake boundary conditions==

We are consideringConsider the case of an inlet perpendicular to the x- direction -[[File:Fig.2 u-velocity cell at intake boundary.JPG|left|thumb|Fig.2 u-velocity cell at intake boundary]]

[[File:Fig.4 pressure correction cell at intake boundary.JPG|thumb|Fig.4 pressure correction cell at intake boundary]]

[[File:Fig.3 v-velocity cell at intake boundary.JPG|center|thumb|Fig.3 v-velocity cell at intake boundary]]

[[File:Fig. 5 scalar cell at intake boundary.JPG|left|thumb|Fig. 5 scalar cell at intake boundary]]

*For the first u, v, φ-cell all links to neighboring nodes are active, so there is no need of any modifications to discretion equations.

*At one of the inlet node absolute pressure is fixed and made pressure correction to zero at that node.

*Generally [[computational fluid dynamics]] codes estimate k and ε with approximate formulate based on turbulent intensity between 1 and 6% and length scale

{|

|+

|-

| [[File:Fig.2 u-velocity cell at intake boundary.JPG|thumb|Fig.2 u-velocity cell at intake boundary|alt=|none]] || [[File:Fig.3 v-velocity cell at intake boundary.JPG|thumb|Fig.3 v-velocity cell at intake boundary|alt=|none]] || [[File:Fig.4 pressure correction cell at intake boundary.JPG|thumb|Fig.4 pressure correction cell at intake boundary|alt=|none]] || [[File:Fig. 5 scalar cell at intake boundary.JPG|thumb|Fig. 5 scalar cell at intake boundary|alt=|none]]

|}

==Symmetry boundary condition==

Line 36:

In this type of situations values of properties just adjacent to the solution ___domain are taken as values at the nearest node just inside the ___domain.

==physicalPhysical boundary conditions==

Consider situation solid wall parallel to the x-direction:

[[File:Fig.6 u-velocity cell at a physical boundary.JPG|left|thumb|Fig.6 u-velocity cell at a physical boundary]]

[[File:Fig.8 v-cell at physical boundary j=NJ.JPG|thumb|Fig.8 v-cell at physical boundary j=NJ]]

[[File:Fig.7 v-cell at physical boundary j=3.JPG|center|thumb|Fig.7 v-cell at physical boundary j=3]]

[[File:Fig.9 scalar cell at a physical boundary.JPG|left|thumb|Fig.9 scalar cell at a physical boundary]]

'''Assumptions made and relations considered'''-

Line 49 ⟶ 45:

*No slip condition: u = v = 0.

*In this we are applying the “wall functions” instead of the mesh points.

{|

|+

|-

| [[File:Fig.6 u-velocity cell at a physical boundary.JPG|left|thumb|Fig.6 u-velocity cell at a physical boundary]] || [[File:Fig.7 v-cell at physical boundary j=3.JPG|center|thumb|Fig.7 v-cell at physical boundary j=3]] || [[File:Fig.8 v-cell at physical boundary j=NJ.JPG|thumb|Fig.8 v-cell at physical boundary j=NJ]] ||

[[File:Fig.49 pressure correctionscalar cell at intakea physical boundary.JPG|left|thumb|Fig.49 pressure correctionscalar cell at intakea physical boundary]]

|}

'''[[Turbulent flow]]''':

Line 54 ⟶ 57:

<math> y^+ > 11.63\,</math>.

in the log-law region of a turbulent [[boundary layer]].

'''[[Laminar flow]]''' :

Line 67 ⟶ 70:

*No chemical reactions at the wall

==cyclicCyclic boundary condition==

*We take flux of flow leaving the outlet cycle boundary equal to the flux entering the inlet cycle boundary

*Values of each variable at the nodes at upstream and downstream of the inlet of the inlet plane are equal to values at the nodes at upstream and downstream of the outlet plane.

==Pressure boundary condition==

[[File:Fig.10 p’-cell at an intake boundary.JPG|left|thumb|Fig.10 p’-cell at an intake boundary]]

{|

[[File:Fig. 11 p’-cell at an exit boundary.JPG|thumb|Fig. 11 p’-cell at an exit boundary]]

|+

|-

| [[File:Fig.10 p’-cell at an intake boundary.JPG|left|thumb|Fig.10 p’-cell at an intake boundary]] || [[File:Fig. 11 p’-cell at an exit boundary.JPG|thumb|Fig. 11 p’-cell at an exit boundary]]

|}

These conditions are used when we don’t know the exact details of flow distribution but boundary values of pressure are known

For example: external flows around objects, internal flows with multiple outlets, [[buoyancy]]-driven flows, [[free surface]] flows, etc.

*The pressure corrections are taken zero at the nodes.

==exitExit boundary conditions==

Considering the case of an outlet perpendicular to the x-direction -

[[File:Fig.12 A control volume at an exit boundary.JPG|left|thumb|Fig.12 A control volume at an exit boundary]]

{|

[[File:Fig. 14 pressure correction cell at an exit boundary.JPG|thumb|Fig. 14 pressure correction cell at an exit boundary]]

|+

[[File:Fig. 13 v-control volume at an exit boundary.JPG|center|thumb|Fig. 13 v-control volume at an exit boundary]]

|-

[[File:Fig.15 scalar cell at an exit boundary.JPG|left|thumb|Fig.15 scalar cell at an exit boundary]]

| [[File:Fig.12 A control volume at an exit boundary.JPG|left|thumb|Fig.12 A control volume at an exit boundary]] || [[File:Fig. 13 v-control volume at an exit boundary.JPG|center|thumb|Fig. 13 v-control volume at an exit boundary]] || [[File:Fig. 14 pressure correction cell at an exit boundary.JPG|thumb|Fig. 14 pressure correction cell at an exit boundary]] || [[File:Fig.15 scalar cell at an exit boundary.JPG|left|thumb|Fig.15 scalar cell at an exit boundary]]

|}

In fully developed flow no changes occurs in flow direction, gradient of all variables except pressure are zero in flow direction

Line 99 ⟶ 108:

==References==

<references/>

*An introduction to computational fluid dynamics by Versteeg, PEARSON.

[[Category:Computational fluid dynamics]]

[[Category:Boundary conditions|computational fluid dynamics in]]