docs:utils:findorbit
====== Differences ====== This shows you the differences between two versions of the page.
| Both sides previous revision Previous revision Next revision | Previous revision | ||
|
docs:utils:findorbit [2009/02/11 08:39] gibson |
docs:utils:findorbit [2010/02/02 07:55] (current) |
||
|---|---|---|---|
| Line 18: | Line 18: | ||
| -Ng --Ngmres <int> default == 40 max number of GMRES iterations per restart | -Ng --Ngmres <int> default == 40 max number of GMRES iterations per restart | ||
| -d --delta <real> default == 0.01 initial radius of trust region | -d --delta <real> default == 0.01 initial radius of trust region | ||
| - | -sx --x-sign change u,x sign | + | -sigma --sigma <string> file containing symmetry of relative eqb/orb |
| - | -sy --y-sign change v,y sign | + | -symms --symmetries <string> constrain u(t) to invariant subspace generated by symmetries in listed file, arg is filename |
| - | -sz --z-sign change w,z sign | + | |
| - | -ax --axshift <real> default == 0 translate by ax*Lx | + | |
| - | -az --azshift <real> default == 0 translate by az*Lz | + | |
| -R --Reynolds <real> default == 400 Reynolds number | -R --Reynolds <real> default == 400 Reynolds number | ||
| -vdt --variableDt adjust dt to keep CFLmin<=CFL<CFLmax | -vdt --variableDt adjust dt to keep CFLmin<=CFL<CFLmax | ||
| Line 39: | Line 36: | ||
| * -eqb leaves the integration time fixed | * -eqb leaves the integration time fixed | ||
| * -orb adds the integration time T to the search space | * -orb adds the integration time T to the search space | ||
| - | * -rel modifies -orb / -eqb, makes the search for a relative periodic orbit or relative equilibrium (traveling wave), by adding the finite-time phase shift variables ''ax'' and ''az'' to the search space | + | * -rel modifies -orb / -eqb, makes the search for a relative periodic orbit or relative equilibrium (traveling wave), by adding finite-time phase shift variables ''ax'' and ''az'' to the search space |
| - | * -ax and -az specify a fixed or initial guess finite-time phase shifts for the solution (depending on whether -rel is specified) | + | * -sigma <filename> sets the initial guess (if -rel) or fixed value (if not -rel) of the symmetry \sigma |
| + | * -symms <filename> constrains the search to an invariant subspace. The file is a list of generators for the isotropy group. | ||
| In all cases we search for a solution of | In all cases we search for a solution of | ||
| <latex> | <latex> | ||
| - | \tau_{a_x, a_z} f^T(u) - u = 0 | + | \sigma f^T(u) - u = 0 |
| </latex> | </latex> | ||
| Line 51: | Line 49: | ||
| <latex> | <latex> | ||
| - | \tau_{a_x, a_z} [u,v,w](x,y,z) \rightarrow [u,v,w](x+a_x/L_x, \, y, \, z+a_z/L_z) | + | \sigma [u,v,w](x,y,z) \rightarrow s [s_x u, \, s_y v, \, s_z w](s_x x + a_x/L_x, \, s_y y, \, s_z z + a_z/L_z) |
| </latex> | </latex> | ||
| - | Further symmetry operations are available but not yet documented. | + | (will continue later today.... /2009-06-08 jfg) |
| ====== Usage examples ====== | ====== Usage examples ====== | ||
docs/utils/findorbit.1234370396.txt.gz · Last modified: 2009/02/11 08:39 by gibson
Page Tools
Except where otherwise noted, content on this wiki is licensed under the following license: CC Attribution-Share Alike 3.0 Unported
