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

pressure neumann boundary condition for invicid flow around moving boundaries

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   April 17, 2022, 21:51
Default pressure neumann boundary condition for invicid flow around moving boundaries
  #1
Member
 
Anders Aamodt Resell
Join Date: Dec 2021
Location: Oslo, Norway
Posts: 66
Rep Power: 4
ander is on a distinguished road
Hi

I'm struggling to understand how a a correct neumann bc for the pressure is supposed to be derived when the Euler equations are used, and moving boundaries are permitted. The impermeability condition reads \mathbf{u}\cdot\mathbf{n} = \mathbf{u}_\Gamma\cdot\mathbf{n}, where \mathbf{u} is the flow velocity, \mathbf{u}_\Gamma is the velocity of the solid boundary, and \mathbf{n} is the normal vector pointing out of the solid domain. The momentum equation projected in the normal direction gives \frac{\partial p}{\partial n} = -\rho\frac{D\mathbf{u}}{Dt}\cdot\mathbf{n}, so the challenge is to find \frac{D\mathbf{u}}{Dt}\cdot\mathbf{n}. I've seen a paper with the name of "A high order moving boundary treatment for compressible inviscid flows" take the material derivative of the impermeability condition stated above, which by application of the product rule gives \frac{D\mathbf{u}}{Dt}\cdot\mathbf{n} = \frac{D(\mathbf{u}_\Gamma\cdot\mathbf{n})}{Dt} - \mathbf{u}\cdot\frac{D\mathbf{n}}{Dt}. A source of confusion in the mentioned paper was that the term \frac{D(\mathbf{u}_\Gamma\cdot\mathbf{n})}{Dt} was written as \frac{d(\mathbf{u}_\Gamma\cdot\mathbf{n})}{dt}.

I guess my primary issue is to find out how the material derivative of the normal velocity \frac{D\mathbf{n}}{Dt} and the term \frac{D(\mathbf{u}_\Gamma\cdot\mathbf{n})}{Dt} (=?)\; \frac{d(\mathbf{u}_\Gamma\cdot\mathbf{n})}{dt} are computed. I'm confused on whether for instance the material normal vector should be "advected with the flow" or if it's described relative to the solid body. In the case of a curved, static boundary
\frac{D\mathbf{n}}{Dt} would yield a nonzero term for a nonzero tangential velocity, while \frac{d\mathbf{n}}{dt} would yield zero anyway, if I'm not mistaken.

If someone could explain or possibly point me to a nice source on the topic, it would be greatly appreciated.
ander is offline   Reply With Quote

Reply


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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Question about different kinds of Boundaries and Boundary Conditions granzer Main CFD Forum 17 April 12, 2022 18:27
JOB ABORT invoked by rank 17 ltl CONVERGE 3 January 25, 2021 16:05
Match Pressure Inlet/Outlet Boundary Condition Mass Flow Rate MSchneid Fluent UDF and Scheme Programming 3 February 23, 2019 07:00
Problem with SIMPLEC-like finite volume channel flow boundary conditions ghobold Main CFD Forum 3 June 15, 2015 12:14
An error has occurred in cfx5solve: volo87 CFX 5 June 14, 2013 18:44


All times are GMT -4. The time now is 11:39.