This is an old revision of the document!
====== findsoln ====== find an equilibrium, traveling wave, periodic orbit, or relative periodic orbit of plane Couette or channel flow ====== main options ====== <code> -eqb --equilibrium search for equilibrium or relative equilibrium (trav wave) -orb --periodicorbit search for periodic orbit or relative periodic orbit -xrel --xrelative search over x phase shift for relative orbit or eqb -zrel --zrelative search over z phase shift for relative orbit or eqb -T --maptime <real> default == 20 initial guess for orbit period or time of eqb/reqb map f^T(u) -R --Reynolds <real> default == 400 Reynolds number -sigma --sigma <string> file containing symmetry of relative solution (default == identity) -symms --symmetries <string> file containing generators of isotropy group for symmetry-constrained search ... -o --outdir <string> default == ./ output directory -log --logfile <string> default == stdout output log (filename or "stdout") <flowfield> (trailing arg 1) initial guess for solution </code> ====== solution type ====== ''findsoln'' finds solutions of //σ f^T(u) - u = 0// where //f^T// is the time-//T// forward map of the Navier-Stokes equations plus boundary conditions and σ is a symmetry of the flow. This equation has four kinds of solutions * equilibrium: //σ = 1///, all //T// * relative equilibrium: //σ = // a translation that scales linearly in //T//, all //T// (also called a traveling wave) * periodic orbit: //σ = 1// and some finite value of //T// (and integer multiples of it) * relative periodic orbit: some element //σ// of the appropriate symmetry group, and some finite value of //T// (and integer multiples of it) You specify which kind of solution to find with the <code> -eqb, -orb, -xrel, -zrel</code> options. The symmetry σ for the relative equilibrium and relative periodic orbit solutions involves unknown phase shifts (translations) in //x//, //z//, or both //x,z//, depending on the symmetry of the initial guess //u//, and you specify which to treat as unknowns by using the ''-xrel'' and ''-zrel'' options. Here are some examples of searches of various types ^command^search type^ | ''findsoln -eqb uguess.h5'' | equilibrium solution | | ''findsoln -eqb -xrel uguess.h5'' | traveling wave in //x// | | ''findsoln -eqb -xrel -zrel uguess.h5'' | traveling wave in //x,z// | | ''findsoln -eqb -zrel -sigma sigma.asc -T 20 uguess.h5'' | traveling wave in //z// with an initial guess for wavespeed | | ''findsoln -orb -T 64.3 uguess.asc'' | periodic orbit with initial guess for period //T// | | ''findsoln -orb -T 34.8 -sigma sigma.asc uguess.h5'' | "pre-periodic orbit" | | ''findsoln -orb -xrel -T 49.1 -sigma sigma.asc uguess.h5 | relative periodic orbit with phase shift in //x// | ([[docs:math:newton_krylov_hookstep|mathematical details of algorithm]]
