CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Main CFD Forum

Eigenvalues & Eigenvectors in General Geometries for 3-D TVD Solver

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   October 6, 1999, 10:42
Default Eigenvalues & Eigenvectors in General Geometries for 3-D TVD Solver
  #1
Axel Rohde
Guest
 
Posts: n/a
As part of my dissertation I am developing a TVD solver for 3-D inviscid flow in finite volume format. I have been extensively searching TVD papers for the eigenvector-matrices of left and right eigenvectors in terms of an arbitrary unit normal direction (Nx, Ny, Nz), i.e. consistent with finite volume formulation. So far the only papers for 3-D flow that I have seen decompose the eigenvector matrix along the three Cartesian directions (1, 0, 0); (0, 1, 0); (0, 0, 1) and treat each matrix separately. Anyhow, I decided to derive the more general expression in terms of an arbitrary unit normal (Nx, Ny, Nz) myself with the help of Mathematica. This task was not that simple, because for repeated eigenvalues, the corresponding eigenvectors are not distinct. They form a subspace, and any vector in that subspace is also an eigenvector. I was wondering if anyone has ever seen the eigensystem (eigenvalues & eigenvectors) for 3-D inviscid flow (Euler equations) in general geometries, i.e. for any unit normal direction. If not, I am thinking about publishing my result.
  Reply With Quote

Old   October 6, 1999, 11:32
Default Re: Eigenvalues & Eigenvectors in General Geometries for 3-D TVD Solver
  #2
Frank Bramkamp
Guest
 
Posts: n/a
Sometimes, people solve riemann problems in the x,y,z direction (or i,j,k) and then multiply it by the normal vectors. E.g. F*n_x+G*n_y+H*N_z. So one solves 3 riemann problems. This is (of course) not required. A more general and better approach is (that is what I do) to use the rotational invariance of the euler equations, which reduces everything to a solely 1D problem.

It is simply,

F*n_x+G*n_y+H*N_z = T_R^{-1}*F(T_R*U)

where U is the vector of conserved variables and T_R a rotational matrix which rotates the equations for the momentum equations into a local coordinate system, alligned with each face (in fact for each gauss quadrature point). So one computes ONE Riemann problem and maps it back into the x,y,z coordinate system. The euler equations have this nice property.

In my opinion, that makes live much easier and one can use the standard eigenvalues and vectors. In this case e.g. u denotes the velocity, perpendicular to the face and v,w the appropriate tangential velocities.

The nice thing is: The difference between a 3D and 2D code gets in fact very small and can be treated via ifdefs in a nice way.

I already did this stuff for a 3D unstructured finite volume scheme, that works on arbitrary grids of any topology using 2D/3D reconstruction schemes, Vankatakrishnan limiter, Barth etc...works fine so far.

Hope this helps a bit.

Frank

  Reply With Quote

Old   October 6, 1999, 12:36
Default Re: Eigenvalues & Eigenvectors in General Geometries for 3-D TVD Solver
  #3
Patrick Godon
Guest
 
Posts: n/a
I think the group in Brown University has been working in curved geometry with spectral methods. Since spectral methods require the treatment of the boundary conditions on the characteristics of the flow, they probably solved for the eigenvector (characteristics of the flow) and eigenvalue in curved geometry using a Jacobian, etc..

You might want to try to contact them:

Jan Hesthaven: jansh@cfm.brown.edu

David Gottlieb: dig@cfm.brown.edu

David Gottlieb is probably pretty busy, so try Jan.

I have myself solved for the egienvalues and eigenvectors in spherical and cylindrical coordinates for the Normal vector in the radial direction. The matrices were similar to the one in cartesian coordinates (after linearization of the euler equations).

Patrick
  Reply With Quote

Old   October 6, 1999, 14:48
Default Re: Eigenvalues & Eigenvectors in General Geometries for 3-D TVD Solver
  #4
Bert Laney
Guest
 
Posts: n/a
You should consult papers on "truly multidimensional" numerical methods, such as those by Hirsch, Roe, and van Leer. Most ordinary "TVD" papers will not have what you want - it is very difficult to usefully exploit truly-multidimensional characteristics in a numerical method and thus few attempt it.

The last chapter of my book, "Computational Gasdynamics," includes some expressions for the multidimensional characteristics of the Euler equations that you might find helpful. Unfortunately, there was not room to include all of the derivations such as those involving the Jacobian matrices. I have quite lengthy formal notes with the full details, originally planned for inclusion in the book, which I may someday see fit to distribute.

Bascially, in multidimensions, instead of a Jacobian matrix (two-dimensional), you have a Jacobian tensor (three-dimensional, two deep for two-dimensions, and three deep for three-dimensions). Multiplying by unit vectors, such as coordinate unit axis vectors, brings the Jacobian tensor from three- to two-dimensions, where it is much easier to express and work with. On the other hand, the three-dimensional tensor is invariant to coordinate transformations, which is not true of the Jacobian matrices derived from the coordinate unit axis vectors.
  Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
general doubt on solver control arunraj CFX 5 July 31, 2011 07:51
why the solver reject it? Anyone with experience? bearcat CFX 6 April 28, 2008 14:08
Texts for eigenvalues, eigenvectors & pde's diaw (Des Aubery) Main CFD Forum 24 April 24, 2006 11:32
CFX 5.5 Roued CFX 1 October 2, 2001 16:49
Eigenvalues and Eigenvectors for 3-D TVD Finite Volume Solver Axel Rohde Main CFD Forum 3 April 2, 2000 17:42


All times are GMT -4. The time now is 18:43.