popiol Posted March 10, 2012 Share Posted March 10, 2012 I'm looking for a numerical stability and error estimation of a finite element approximation of Navier-Stokes equations (with combustion). I define variables and operators on a domain that has both space and time axes (Ω = Ωs x [0,tmax]), so the transport equation looks generally like this div ( u [ v; 1 ] - D [ su; 0 ] ) = Q, where u is either density (of one of the species), velocity, temperature or pressure; v is velocity; D is diffusion coefficient; s is gradient over the space domain Ωs; and Q is either reaction rate, pressure gradient plus buoyancy force (-sp + fb), energy release or 0. The approximation scheme is ui = Σj<i ( D ( dij-1 - di.-1) - [ vj; 1 ] . ( xj - xi ) / dij ) wij uj where ui is a value in the i-th mesh node; dij = ||xj - xi||; di.-1= Σj<i wij / dij; Σj<i wij = 1; and xi is a mesh node in Ω. Mesh nodes are sorted by time, so t(xi) > t(xj) => i > j. It is a bit difficult to define stability in this case, but the following condition seems reasonable lim i -> ∞ ( ui - Σj<i ujwij ) = 0, which implies lim i -> ∞D ( dij-1 - di.-1) - [ vj; 1 ] . ( xj - xi ) / dij = 1. It should also be true that for each j > 0 Σi>j ( D ( dij-1 - di.-1) - [ vj; 1 ] . ( xj - xi ) / dij ) wij = 1 Any idea how to verify those conditions? Link to comment Share on other sites More sharing options...
Recommended Posts
Create an account or sign in to comment
You need to be a member in order to leave a comment
Create an account
Sign up for a new account in our community. It's easy!
Register a new accountSign in
Already have an account? Sign in here.
Sign In Now