# Classification of boundary conditions

### From Thermal-FluidsPedia

Yuwen Zhang (Talk | contribs) |
|||

Line 1: | Line 1: | ||

+ | For any problem to be well defined, there are boundary/initial conditions that must be applied. There are three basic boundary conditions for second order PDEs. These boundary conditions are the Dirichlet [<math>\Phi = f\left( {\eta ,\zeta } \right)</math>], the Neumann [<math>\partial \Phi /\partial \eta = f\left( \zeta \right)</math>, or <math>\partial \Phi /\partial \zeta = f\left( \eta \right)</math>], and the mixed [<math>a\,(\partial \Phi /\partial \eta ) + b\Phi = f(\zeta ) | ||

+ | </math>, or <math>a(\partial \Phi /\partial \zeta ) + b\Phi = f(\eta ) </math>] type. | ||

+ | |||

+ | The number of boundary conditions that must be applied for each independent variable is equal to the highest order of that variable. For example, if the PDE is second order in both <math>\eta</math> and <math>\zeta</math>, then two boundary conditions are needed for each independent variable for a total of four boundary conditions. If <math>\zeta</math> is time and first order, then an initial Dirichlet-type boundary condition is needed, and if it is second order, then both a Dirichlet and Neumann boundary condition are needed. | ||

+ | |||

+ | There are implications to experimental measurements and numerical analysis based on the classification of the governing PDE. For example, experimental measurements in an incompressible flow field with moderate to low Reynolds numbers are very difficult, because the governing equations are highly elliptic. The elliptic nature means that disturbances downstream greatly affect the upstream flow field. Therefore, any measurements that require a device in the flow field may change the nature of that flow field, and are inherently inaccurate. Numerical simulations are very reliable for these cases when the flow is in the laminar regime. In the laminar flow regime, the full Navier-Stokes equations can be directly solved as an elliptic problem with no approximations. | ||

+ | |||

+ | In the compressible flow regime with high Reynolds numbers, the characteristic of the flow is parabolic or hyperbolic, depending on the Mach number. In these cases, disturbances produced downstream do not affect the upstream flow field; therefore measurement devices in the flow will give an accurate depiction of the flow field without the device. In numerical simulations, it is computationally efficient to solve a parabolic flow field, because disturbances propagate in one direction. However, a flow field that is truly parabolic usually has a high Reynolds number and therefore is turbulent, which means it is three-dimensional and transient. Turbulence modeling either involves a very fine mesh, which is very computationally expensive, or an averaging technique, which requires a closure problem. The classification of a particular problem will help reduce the time it takes to get reliable results. | ||

+ | |||

The conservation equations introduced above can be applied within each phase and up to an interface. However, they are not valid across the interface, where sharp changes in various properties occur. Appropriate boundary conditions at the interface must be specified in order to solve the governing equations for heat, mass, and momentum transfer in the two adjoining phases. The interface conditions will serve as boundary conditions for the transport equations | The conservation equations introduced above can be applied within each phase and up to an interface. However, they are not valid across the interface, where sharp changes in various properties occur. Appropriate boundary conditions at the interface must be specified in order to solve the governing equations for heat, mass, and momentum transfer in the two adjoining phases. The interface conditions will serve as boundary conditions for the transport equations | ||

## Revision as of 09:31, 28 June 2010

For any problem to be well defined, there are boundary/initial conditions that must be applied. There are three basic boundary conditions for second order PDEs. These boundary conditions are the Dirichlet [], the Neumann [, or ], and the mixed [, or ] type.

The number of boundary conditions that must be applied for each independent variable is equal to the highest order of that variable. For example, if the PDE is second order in both η and ζ, then two boundary conditions are needed for each independent variable for a total of four boundary conditions. If ζ is time and first order, then an initial Dirichlet-type boundary condition is needed, and if it is second order, then both a Dirichlet and Neumann boundary condition are needed.

There are implications to experimental measurements and numerical analysis based on the classification of the governing PDE. For example, experimental measurements in an incompressible flow field with moderate to low Reynolds numbers are very difficult, because the governing equations are highly elliptic. The elliptic nature means that disturbances downstream greatly affect the upstream flow field. Therefore, any measurements that require a device in the flow field may change the nature of that flow field, and are inherently inaccurate. Numerical simulations are very reliable for these cases when the flow is in the laminar regime. In the laminar flow regime, the full Navier-Stokes equations can be directly solved as an elliptic problem with no approximations.

In the compressible flow regime with high Reynolds numbers, the characteristic of the flow is parabolic or hyperbolic, depending on the Mach number. In these cases, disturbances produced downstream do not affect the upstream flow field; therefore measurement devices in the flow will give an accurate depiction of the flow field without the device. In numerical simulations, it is computationally efficient to solve a parabolic flow field, because disturbances propagate in one direction. However, a flow field that is truly parabolic usually has a high Reynolds number and therefore is turbulent, which means it is three-dimensional and transient. Turbulence modeling either involves a very fine mesh, which is very computationally expensive, or an averaging technique, which requires a closure problem. The classification of a particular problem will help reduce the time it takes to get reliable results.

The conservation equations introduced above can be applied within each phase and up to an interface. However, they are not valid across the interface, where sharp changes in various properties occur. Appropriate boundary conditions at the interface must be specified in order to solve the governing equations for heat, mass, and momentum transfer in the two adjoining phases. The interface conditions will serve as boundary conditions for the transport equations

in the adjacent phases. Jump conditions at the interface can be obtained by applying the basic laws (conservation of mass, momentum, energy, and the second law of thermodynamics) at the interface. It is the objective of this subsection to specify mass, momentum, and energy balance at a non-flat liquid-vapor interface (see Fig. 1), as well as species balance in solid-liquid-vapor interfaces. For solid-liquid or solid-vapor interfaces, these jump conditions can be significantly simplified.

## Contents |