This shows you the differences between two versions of the page.
Both sides previous revision Previous revision Next revision | Previous revision | ||
docs:math:symmetry [2009/05/04 05:37] predrag |
docs:math:symmetry [2014/12/04 11:53] (current) gibson [Symmetry of flows in channel geometries] |
||
---|---|---|---|
Line 3: | Line 3: | ||
The symmetry group of 3D fields in channel geometries is generated by | The symmetry group of 3D fields in channel geometries is generated by | ||
- | <latex> $ \begin{align*} | + | \begin{eqnarray*} |
- | [u,v,w](x,y,z) &\rightarrow [-u, v, w](-x,y,z) \\ | + | [ u,v,w ](x,y,z) & \rightarrow [-u, v, w](-x,y,z) \\ |
- | [u,v,w](x,y,z) &\rightarrow [ u, -v, w](x,-y,z) \\ | + | [ u,v,w ](x,y,z) & \rightarrow [ u, -v, w](x,-y,z) \\ |
- | [u,v,w](x,y,z) &\rightarrow [ u, v, -w](x, y,-z) \\ | + | [ u,v,w ](x,y,z) & \rightarrow [ u, v, -w](x, y,-z) \\ |
- | [u,v,w](x,y,z) &\rightarrow [-u,-v,-w](x,y,z) \\ | + | [ u,v,w ](x,y,z) & \rightarrow [-u,-v,-w](x,y,z) \\ |
- | [u,v,w](x,y,z) &\rightarrow [ u, v, w](x+\ell_x, y, z+\ell_z) | + | [ u,v,w ](x,y,z) & \rightarrow [ u, v, w](x+\ell_x, y, z+\ell_z) |
- | \end{align*} $ </latex> | + | \end{eqnarray*} |
- | "Channel geometry" means a domain that is periodic or infinite in x and z | + | "Channel geometry" means a domain that is periodic or infinite in //x// and //z// |
- | and bounded in y, with %%-Ly/2 <= y <= Ly/2%% and Dirichlet or Neumann | + | and bounded in //y//, with $-L_y/2 \leq y \leq L_y/2$ and Dirichlet or Neumann |
- | boundary conditions at the bounds in y. The symmetry groups of velocity | + | boundary conditions at the bounds in //y//. The symmetry groups of velocity |
fields for specific flows, with constraints such as incompressibility | fields for specific flows, with constraints such as incompressibility | ||
and specific boundary conditions, are subgroups of the group | and specific boundary conditions, are subgroups of the group | ||
generated by the above symmetries. | generated by the above symmetries. | ||
- | |||
====== Symmetry of plane Couette flow======= | ====== Symmetry of plane Couette flow======= | ||
+ | |||
+ | For the full description of 67 isotropy subgroups of plane Couette, see J. Halcrow, J. F. Gibson, and P. Cvitanović, | ||
+ | //Equilibrium and traveling-wave solutions of plane Couette flow//, [[http://arxiv.org/abs/0808.3375|arXiv:0808.3375]], J. Fluid Mech. (to appear, 2009), and [[http://chaosbook.org/projects/Halcrow/thesis.pdf|J. Halcrow, "Charting the state space of plane Couette flow: Equilibria, relative equilibria, and heteroclinic connections"]] (Georgia Tech Ph.D. thesis, Aug 2008). Here are some highlights. | ||
===== Invariance ===== | ===== Invariance ===== | ||
Line 24: | Line 26: | ||
Plane Couette flow is invariant under the following symmetries | Plane Couette flow is invariant under the following symmetries | ||
- | <latex> | + | \begin{eqnarray*} |
- | $ \begin{align*} | + | |
\sigma_x \, [u,v,w](x,y,z) &= [-u,-v,w](-x,-y,z) \\ | \sigma_x \, [u,v,w](x,y,z) &= [-u,-v,w](-x,-y,z) \\ | ||
\sigma_z \, [u,v,w](x,y,z) &= [u, v,-w](x,y,-z) \\ | \sigma_z \, [u,v,w](x,y,z) &= [u, v,-w](x,y,-z) \\ | ||
\tau(\ell_x, \ell_z) [u,v,w](x,y,z) &= [u, v,-w](x+\ell_x, y, z+\ell_z) \\ | \tau(\ell_x, \ell_z) [u,v,w](x,y,z) &= [u, v,-w](x+\ell_x, y, z+\ell_z) \\ | ||
- | \end{align*} | + | \end{eqnarray*} |
- | $ | + | |
- | </latex> | + | |
That is, if f^t(u) is the time-t map of plane Couette flow, then | That is, if f^t(u) is the time-t map of plane Couette flow, then | ||
- | <latex> | + | \begin{eqnarray*} |
f^t(s u) = s f^t(u) | f^t(s u) = s f^t(u) | ||
- | </latex> | + | \end{eqnarray*} |
for any s in group G generated by <latex> \{\sigma_x, \sigma_z, \tau(\ell_x, \ell_z)\}. </latex> | for any s in group G generated by <latex> \{\sigma_x, \sigma_z, \tau(\ell_x, \ell_z)\}. </latex> | ||
Line 43: | Line 42: | ||
Let u(t) be a solution of Navier Stokes with initial condition u(0), | Let u(t) be a solution of Navier Stokes with initial condition u(0), | ||
- | <latex> | + | \begin{eqnarray*} |
u(t) = f^t(u(0)) | u(t) = f^t(u(0)) | ||
- | </latex> | + | \end{eqnarray*} |
then | then | ||
- | <latex> | + | \begin{eqnarray*} |
s u(t) = s f^t(u(0)) = f^t(s u(0)) | s u(t) = s f^t(u(0)) = f^t(s u(0)) | ||
- | </latex> | + | \end{eqnarray*} |
is a solution of Navier-Stokes with initial condition s u(0). | is a solution of Navier-Stokes with initial condition s u(0). | ||
Line 60: | Line 59: | ||
Suppose //u(0)// is invariant under a symmetry //s// in //G//, i.e. | Suppose //u(0)// is invariant under a symmetry //s// in //G//, i.e. | ||
- | <latex> | + | \begin{eqnarray*} |
s u(0) = u(0) | s u(0) = u(0) | ||
- | </latex> | + | \end{eqnarray*} |
Then //u(t)// satisfies that symmetry for all //t//, since | Then //u(t)// satisfies that symmetry for all //t//, since | ||
- | <latex> | + | \begin{eqnarray*} |
s u(t) = s f^t(u(0)) = f^t(s u(0)) = f^t(u(0)) = u(t) | s u(t) = s f^t(u(0)) = f^t(s u(0)) = f^t(u(0)) = u(t) | ||
- | </latex> | + | \end{eqnarray*} |
The set of all symmetries //s// in //G// satisfied by u forms a subgroup //H ⊂ G//, | The set of all symmetries //s// in //G// satisfied by u forms a subgroup //H ⊂ G//, | ||
Line 77: | Line 76: | ||
===== Isotropy groups of known solutions ===== | ===== Isotropy groups of known solutions ===== | ||
- | Most known equilibria and periodic orbits of plane Couette flow have the | + | The isotropy group most known equilibria and periodic orbits of plane Couette flow is |
- | same isotropy group | + | |
- | <latex> | + | \begin{eqnarray*} |
S = \{1, s_1, s_2, s_3 \} | S = \{1, s_1, s_2, s_3 \} | ||
- | </latex> | + | \end{eqnarray*} |
where | where | ||
- | <latex> | + | \begin{eqnarray*} |
- | $ \begin{align*} | + | |
s_1 \, [u, v, w](x,y,z) &= [u, v, -w](x+L_x/2, y, -z) \\ | s_1 \, [u, v, w](x,y,z) &= [u, v, -w](x+L_x/2, y, -z) \\ | ||
s_2 \, [u, v, w](x,y,z) &= [-u, -v, w](-x+L_x/2,-y,z+L_z/2) \\ | s_2 \, [u, v, w](x,y,z) &= [-u, -v, w](-x+L_x/2,-y,z+L_z/2) \\ | ||
s_3 \, [u, v, w](x,y,z) &= [-u,-v,-w](-x, -y, -z+L_z/2) \\ | s_3 \, [u, v, w](x,y,z) &= [-u,-v,-w](-x, -y, -z+L_z/2) \\ | ||
- | \end{align*} | + | \end{eqnarray*} |
- | $ | + | It is helpful to express these symmetries in terms of //σ<sub>x</sub>, σ<sub>z</sub>,// and translations. Let |
- | </latex> | + | |
- | It's helpful to express break these symmetries into these σ<sub>x</sub>, σ<sub>z</sub>, | + | \begin{eqnarray*} |
- | and translations. Let | + | |
- | + | ||
- | <latex> $ \begin{align*} | + | |
\tau_x &= \tau(L_x/2, 0) \\ | \tau_x &= \tau(L_x/2, 0) \\ | ||
\tau_z &= \tau(0, L_z/2) \\ | \tau_z &= \tau(0, L_z/2) \\ | ||
\tau_{xz} &= \tau_x \tau_z | \tau_{xz} &= \tau_x \tau_z | ||
- | \end{align*} $ </latex> | + | \end{eqnarray*} |
then | then | ||
- | <latex> $ \begin{align*} | + | \begin{eqnarray*} |
S = \{1, \, \tau_x \sigma_z, \, \tau_{xz} \sigma_x, \, \tau_z \sigma_{xz} \} | S = \{1, \, \tau_x \sigma_z, \, \tau_{xz} \sigma_x, \, \tau_z \sigma_{xz} \} | ||
- | \end{align*} $ </latex> | + | \end{eqnarray*} |
===== Fun facts ===== | ===== Fun facts ===== | ||
1. If u has isotropy group S, then | 1. If u has isotropy group S, then | ||
- | <latex> | + | \begin{eqnarray*} |
\tau_x u, \, \tau_z u, \, \text{ and } \, \tau_{xz} u | \tau_x u, \, \tau_z u, \, \text{ and } \, \tau_{xz} u | ||
- | </latex> | + | \end{eqnarray*} |
also have isotropy group S. Thus for each equilibrium or periodic orbit with isotropy group S, there are four half-box shifted partners. | also have isotropy group S. Thus for each equilibrium or periodic orbit with isotropy group S, there are four half-box shifted partners. | ||
Line 131: | Line 123: | ||
relative periodic orbits only of the form | relative periodic orbits only of the form | ||
- | <latex> | + | \begin{eqnarray*} |
\tau f^t(u) - u = 0 \text{ for } \tau \in T = \{1, \, \tau_x, \, \tau_z, \, \tau_{xz} \} | \tau f^t(u) - u = 0 \text{ for } \tau \in T = \{1, \, \tau_x, \, \tau_z, \, \tau_{xz} \} | ||
- | </latex> | + | \end{eqnarray*} |
- | ===== Isotropy groups and invariant solutions ===== | + | |
- | Our recent submission to JFM has a more complete treatment of the | + | ===== Isotropy groups and invariant solutions ===== |
- | isotropy groups of plane Couette flow and the types of solutions they admit. | + | |
- | For simplicity we have so far restricted most of our attention to the | + | So far we have restricted most of our attention to the |
- | solutions with S isotropy. We have a few solutions with other isotropies. | + | solutions with //S// isotropy. We have a few solutions with other isotropies. |
- | One of the main simplifications of the restriction to S is that reduces | + | One of the main simplifications of the restriction to //S// is that reduces |
the number of free parameters in the search for good initial guesses for | the number of free parameters in the search for good initial guesses for | ||
invariant solutions. E.g. we don't have to provide a guess for the wave | invariant solutions. E.g. we don't have to provide a guess for the wave | ||
speed of traveling waves, and for periodic orbits, there are only four | speed of traveling waves, and for periodic orbits, there are only four | ||
- | choices for the the symmetry σ in | + | choices for the symmetry //σ// in |
- | <latex> | + | \begin{eqnarray*} |
\sigma f^t u - u = 0 | \sigma f^t u - u = 0 | ||
- | </latex> | + | \end{eqnarray*} |
namely, <latex> \sigma = 1, \tau_x, \, \tau_z, \, \text{ or } \tau_{xz} </latex>, rather | namely, <latex> \sigma = 1, \tau_x, \, \tau_z, \, \text{ or } \tau_{xz} </latex>, rather | ||
than the continuum <latex> \tau(\ell_x, \ell_z) </latex>. | than the continuum <latex> \tau(\ell_x, \ell_z) </latex>. | ||
Line 158: | Line 147: | ||
close recurrence within a trajectory u(t) by | close recurrence within a trajectory u(t) by | ||
- | <latex> $ \begin{align*} | + | \begin{eqnarray*} |
r(t,T) &= min_{\tau} \| \sigma f^T u(t) - u(t) \| \\ | r(t,T) &= min_{\tau} \| \sigma f^T u(t) - u(t) \| \\ | ||
&= min_{\tau} \| \sigma u(t+T) - u(t) \| | &= min_{\tau} \| \sigma u(t+T) - u(t) \| | ||
- | \end{align*} $ </latex> | + | \end{eqnarray*} |
- | + | for $ \tau \in \{1, \, \tau_x, \, \tau_z, \tau_{xz}\}$. | |
- | for <latex> \tau \in \{1, \, \tau_x, \, \tau_z, \tau_{xz}\}</latex>. | + | |
We can compute r(t,T) from a time series of u(t) and look for places where r(t,T) << 1 for | We can compute r(t,T) from a time series of u(t) and look for places where r(t,T) << 1 for | ||
stretches of t and constant T. Those will be good guesses for periodic orbits. | stretches of t and constant T. Those will be good guesses for periodic orbits. | ||
- | |||
- |