deforming mesh and freestream preservation

syoo · April 12, 2011, 03:22

Hello,

I am currently trying to code an 2D Euler-eq solver using the Beam-Warming scheme for a given deforming mesh. So far, I can't preserve freestream using my solver (embarrassing ...) for deforming grids. The freestream is conserved for a non-deforming grid.

My current understanding is that the mesh deformation (mesh velocity) is accounted for in the solution through the Jacobian which transforms the grid from the cartesian coordinate to generalized body-fitted coordinate. I also satisfy the Geometric Conservation Law by calculating the grid metrics written the conservative form, suggested by Thomas and Lombard.

Should I be including the grid velocity in the flux terms? I have went to many literatures that use Beam-Warming scheme and none of them mention this.

Any help, suggestion, or comment will be very much appreciated.

thank you.

-Paul

agd · April 12, 2011, 15:11

You need to include the grid velocity terms in the fluxes, since the flux across the face depends on the fluid velocity relative to the face. Additionally, Thomas and Lombard describe a procedure for ensuring geometric conservation when the grid is moving that involves the proper construction of the time derivative of the Jacobian.

syoo · April 14, 2011, 01:05

thank you. I'll try it.

-paul

syoo · April 18, 2011, 12:12

I would like some clarification.

I have seen many papers using ALE formulation which includes the grid velocities in the fluxes to accommodate for non-rigidly deforming grids. But I also see papers (ex. Visbal, "On the Use of Higher-ORder Finite-Difference Schemes on Curviilinear and Deforming Meshes") that does not mention including the grid velocities in the flux. In Visbal paper, he uses the same time-implicit method that I use which the Beam-Warming scheme.

so is including the grid velocity necessary? or does Visbal imply that grid velocities are already included in the velocity terms?

thank you.

-Paul

agd · April 18, 2011, 13:42

Sorry - I should have been clearer in my original post. The grid velocities are included in the flux terms, but they appear in the form of the grid time metrics (the time derivatives of the grid transformation terms). The Jacobian relationship in Thomas-Lombard connects the time derivative of the Jacobian to the spatial variation of the time metric terms. All of this is necessary for GCL.

April 12, 2011, 03:22	deforming mesh and freestream preservation	#1
syoo New Member Paul Join Date: Apr 2011 Posts: 12 Rep Power: 15	Hello, I am currently trying to code an 2D Euler-eq solver using the Beam-Warming scheme for a given deforming mesh. So far, I can't preserve freestream using my solver (embarrassing ...) for deforming grids. The freestream is conserved for a non-deforming grid. My current understanding is that the mesh deformation (mesh velocity) is accounted for in the solution through the Jacobian which transforms the grid from the cartesian coordinate to generalized body-fitted coordinate. I also satisfy the Geometric Conservation Law by calculating the grid metrics written the conservative form, suggested by Thomas and Lombard. Should I be including the grid velocity in the flux terms? I have went to many literatures that use Beam-Warming scheme and none of them mention this. Any help, suggestion, or comment will be very much appreciated. thank you. -Paul

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
mesh generation for supersonic flow followed by freestream	padmanathan	FLUENT	6	September 14, 2010 08:53

April 12, 2011, 15:11		#2
agd Senior Member Join Date: Jul 2009 Posts: 358 Rep Power: 19	You need to include the grid velocity terms in the fluxes, since the flux across the face depends on the fluid velocity relative to the face. Additionally, Thomas and Lombard describe a procedure for ensuring geometric conservation when the grid is moving that involves the proper construction of the time derivative of the Jacobian.

April 14, 2011, 01:05		#3
syoo New Member Paul Join Date: Apr 2011 Posts: 12 Rep Power: 15	thank you. I'll try it. -paul

April 18, 2011, 12:12		#4
syoo New Member Paul Join Date: Apr 2011 Posts: 12 Rep Power: 15	I would like some clarification. I have seen many papers using ALE formulation which includes the grid velocities in the fluxes to accommodate for non-rigidly deforming grids. But I also see papers (ex. Visbal, "On the Use of Higher-ORder Finite-Difference Schemes on Curviilinear and Deforming Meshes") that does not mention including the grid velocities in the flux. In Visbal paper, he uses the same time-implicit method that I use which the Beam-Warming scheme. so is including the grid velocity necessary? or does Visbal imply that grid velocities are already included in the velocity terms? thank you. -Paul

April 18, 2011, 13:42		#5
agd Senior Member Join Date: Jul 2009 Posts: 358 Rep Power: 19	Sorry - I should have been clearer in my original post. The grid velocities are included in the flux terms, but they appear in the form of the grid time metrics (the time derivatives of the grid transformation terms). The Jacobian relationship in Thomas-Lombard connects the time derivative of the Jacobian to the spatial variation of the time metric terms. All of this is necessary for GCL.