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

Divergence free polynomials

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

Reply
 
LinkBack Thread Tools Display Modes
Old   January 22, 2011, 07:40
Default Divergence free polynomials
  #1
Member
 
kostas
Join Date: Jun 2010
Posts: 34
Rep Power: 7
pankos is on a distinguished road
I want to use divergence-free basis in finite element framework for discretizing the Maxwell equations due to divergence free magnetic field.
All paper which i read, include vector free divergence basis for approximation field. Scalar function for approximation of Hx,Hy distinct, how i will carry out?
pankos is offline   Reply With Quote

Old   January 22, 2011, 11:25
Default
  #2
Member
 
Jonas T. Holdeman, Jr.
Join Date: Mar 2009
Location: Knoxville, Tennessee
Posts: 90
Rep Power: 8
Jonas Holdeman is on a distinguished road
Quote:
Originally Posted by pankos View Post
I want to use divergence-free basis in finite element framework for discretizing the Maxwell equations due to divergence free magnetic field.
All paper which i read, include vector free divergence basis for approximation field. Scalar function for approximation of Hx,Hy distinct, how i will carry out?
You might want to look at two recent publications:

[1] J. T. Holdeman. "A Hermite finite element method for incompressible fluid flow." Int. J. Numer. Meth. Fluids, 64:376-408 (2010)
[2] J. T. Holdeman. "Computation of incompressible thermal flows using Hermite elements." Comput. Meth. Appl. Mech. Engrg., 199:3297-3304 (2010)

The finite elements described are strictly divergence-free and application is straight-forward. The trick is that the polynomials used are necessarily Hermite. The papers describe applications in 2D, but a paper for 3D is under development and should be available in preprint form within a year. Irrotational fields are even simpler (mentioned very briefly in [2]) as is extension to 3D.

Last edited by Jonas Holdeman; January 22, 2011 at 11:27. Reason: spelling
Jonas Holdeman is offline   Reply With Quote

Old   January 22, 2011, 12:25
Default
  #3
Member
 
kostas
Join Date: Jun 2010
Posts: 34
Rep Power: 7
pankos is on a distinguished road
Thanks. I will read carefully these papers for 2d case. Finally, i must construct the 3d version of maxwell solver and thus i will need solenoidal basis for 3d elements. Initially, i read the paper "Locally divergence-free discontinuous galerkin methods for the Maxwell equations." by Cockburn, Shu and Li and i saw that very easy can be constructed div-free vector bases (i have 2d and 3d vector basis for rect and hexa). What is your opinion for these basis? Have you worked with these?
pankos is offline   Reply With Quote

Old   January 22, 2011, 16:54
Default
  #4
Member
 
Jonas T. Holdeman, Jr.
Join Date: Mar 2009
Location: Knoxville, Tennessee
Posts: 90
Rep Power: 8
Jonas Holdeman is on a distinguished road
The velocity elements I work with and described in the papers are the curl of a stream function or vector potential element and hence are strictly divergence-free locally. In addition they have continuous normal components at element boundaries and hence are divergence-free there (tangential continuity is not a requirement). Discontinuous elements are obviously not divergence-free at the element interfaces except perhaps on average. Both yield satisfactory results, but implementation is simpler with the Hermite elements.

The 3D element I referred to is again the curl of a vector potential element and has necessary normal continuity at element interfaces, but some VP degrees-of-freedom that do not contribute to the velocity (curl) have been eliminated to produce a simpler velocity element. That VP form is probably not suitable for vector potential formulations in electromagnetics, but that can be remedied.
Jonas Holdeman is offline   Reply With Quote

Old   January 22, 2011, 19:25
Default
  #5
Member
 
kostas
Join Date: Jun 2010
Posts: 34
Rep Power: 7
pankos is on a distinguished road
The sense of vector div-free basis is to be zero the divergence of each of these. Therefore, independenly of DOFs, the div of approximation is set zero.
Consequently, in case of vector basis, i must use only one set of DOFs for approximation of solution vector field. If i use a set of DOFs for approximation of each solution vector componets(scalar), distinctly, i can not carry out the desirable property. For example:

Hx=Σ ci*Pi
Hy=Σ gi*Fi
Hx,Hy are componets of vec(H), ci,gi two sets of DOFs (i=0..n) and Pi,Fi scalar basis.

div(vec(H))=d(Hx)/dx + d(Hy)/dy =Σ ci *d(Pi)/dx + Σ gi *d(Fi)/dy =0 if all Pi=Fi=0 or if d(Pi)/dx=-d(Fi)/dy with special values of DOFs ci,gi.

Otherwise if i use one set of DOFs for approximation of vector field (not each scalar componets) i have:

div(vec(H))=Σ ci*vec(Pi)=ci*Σ div(vec(Pi))=0 for all cases.

The above problem i have when i take as scalar basis for Hx,Hy the componets of vector div-basis which i read in paper of Shu.
pankos is offline   Reply With Quote

Old   January 22, 2011, 21:30
Default
  #6
Member
 
Jonas T. Holdeman, Jr.
Join Date: Mar 2009
Location: Knoxville, Tennessee
Posts: 90
Rep Power: 8
Jonas Holdeman is on a distinguished road
Perhaps the simplest solenoidal vector element on a rectangle has 12 DOFs. At each corner node the DOFs are a stream function and two components of the velocity. For magnetics these would be the z-component of the vector potential and Hx and Hy. Each vector shape function for the 12 DOFs is solenoidal. Thus for any choice of coefficients the resulting field on the element will be solenoidal. So you only have to find values for the coefficients which fit your problem.

There are also triangular solenoidal elements.

Perhaps we could continue this discussion offline.

Last edited by Jonas Holdeman; January 23, 2011 at 21:24. Reason: Remove email address, available on papers
Jonas Holdeman is offline   Reply With Quote

Old   January 23, 2011, 09:05
Default
  #7
Member
 
kostas
Join Date: Jun 2010
Posts: 34
Rep Power: 7
pankos is on a distinguished road
It is necessary to put the potential component? Using only components
(Hx,Hy) as DOFs, can not i carry out the div zero? Having the stream
function as addition DOFs you raise the total number of DOFs.
Also, the expansion (approximation) of each DOF as you say what scalar
basis will be used? the component of vector basis?
pankos is offline   Reply With Quote

Old   January 23, 2011, 11:15
Default
  #8
Member
 
Jonas T. Holdeman, Jr.
Join Date: Mar 2009
Location: Knoxville, Tennessee
Posts: 90
Rep Power: 8
Jonas Holdeman is on a distinguished road
Is it necessary to put the potential component?

The answer is yes. Failure to recognize this is, I think, why this type of element was not recognized 25 years ago. The existance of a divergence-free field v necessarily implies the existence of a potential a such that v=curl a.

For 2D fluids this is no problem. The Navier-Stokes equation can be orthogonally decomposed into divergence-free and curl-free parts. The divergence-free part is a governing equation for the velocity independent of pressure. The curl-free part is related to the pressure Poisson equation.

So how do you control flow through a duct, as an example of pressure-driven flow, when there is no pressure in the equation? The answer is by boundary conditions on the stream function, the net flow is equal to the difference in stream function across the duct. Is this additional work? No, because the stream function replaces the pressure.

What about electromagnetics? If you are constructing divergence-free fields, you still must introduce an additional (potential?) function to project out the non-solenoidal part introduced by using uncorrelated scalar functions for Hx & Hy. So there is no additional work, and you don't have to solve a Poisson equation at each step.

I don't have enough experience with 3D solenoidal elements to understand what economies might be involved. But, you are not constrained by elements that satisfy the LBB condition. Things are complicated by the fact that the dimension of the kernel (null space) of the curl operator is relatively large.

Gresho & Sani point out in their 1998 book (p.468) that there is roughly one momentum conservation and one mass conservation equation at each point in the mesh. They argue that roughly the constraint ratio (mass/momentum) should be roughly one for "good" elements. So with solenoidal Hermite elements, the extra 3 VP DOFs and 3 velocity DOFs result in a constraint ratio of one. There is, in fact, no penalty for introducing the VP DOFs.

Again, in electromagnetics, you don't have a pressure but you must introduce an extraneous potential function to satisfy the divergence-free constraint.
Jonas Holdeman is offline   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
Quarter Burner mesh with periosic condition SamCanuck FLUENT 2 August 31, 2011 11:34
Divergence free initial field maka OpenFOAM Pre-Processing 2 October 26, 2005 08:52
Divergence Free on the Magnetic Field Valdemir Main CFD Forum 1 March 19, 2003 19:33
CFX 4.4 New free surface option Viatcheslav Anissimov CFX 0 April 3, 2002 06:27
help:spectral methods & divergence free functionsn D. Puigjaner Main CFD Forum 1 August 28, 2000 10:06


All times are GMT -4. The time now is 16:04.