CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM

Time varying velocity causes unwanted pressure gradient

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   April 28, 2023, 02:25
Default Time varying velocity causes unwanted pressure gradient
  #1
New Member
 
Ben
Join Date: Feb 2023
Posts: 8
Rep Power: 3
Ben0000 is on a distinguished road
I'm trying to simulate the effect of wind on building, specifically the the pressure distribution on the surface of the building. I've had success with a constant wind velocity with respect to time, attached is a screenshot of the results, but when I try to make the wind velocity change with time I'm running into some issues. The problem is that I'm getting a pressure gradient across the entire domain (see other attachment). This gradient is then being applied onto the building surface, and essentially overwhelms any other pressure that occurs not as a result of the gradient. This gradient also changes significantly between different time steps.

At first I thought this was nonphysical behaviour, but looking at the Navier-Stokes Equations I can see that a non-zero du/dt will lead to some constant dp/dx, i.e. a pressure gradient. I'm using the Von Karman wind turbulence model to generate the time-varying wind velocity, and I've checked this against real-world data to confirm that the acceleration (du/dt) has a realistic magnitude.

My question is how would I go about modelling the time-varying wind around the building without the pressure gradient, or if the pressure gradient is necessary, without the pressure gradient applied onto the building? Currently, it seems that the entire domain changes velocity simultaneously, would it be possible to model some sort of gust generator that generates these varying wind speeds upwind of the building, and they propagate through the domain and around the building?

Below are my current boundary conditions for p and U (freeWall refers to the left, right and top boundaries of the domain).

p:
Code:
dimensions      [0 2 -2 0 0 0 0];

internalField   uniform 0;

boundaryField
{
    inlet
    {
        type            zeroGradient;
    }
    outlet
    {
        type            fixedValue;
        value            uniform 0;
    }
    freeWall
    {
        type            slip;
    }
    floor
    {
        type            zeroGradient;
    }
    buildingWall
    {
        type            zeroGradient;
    }
}
U:
Code:
dimensions      [0 1 -1 0 0 0 0];

internalField    uniform (45 0 0);

boundaryField
{
    inlet
    {
        type            timeVaryingMappedFixedValue;
        setAverage        off;
    }
    outlet
    {
        type            zeroGradient;
    }
    freeWall
    {
        type            slip;
    }
    floor
    {
        type            noSlip;
    }
    buildingWall
    {
        type            noSlip;
    }
}
Attached Images
File Type: jpg ConstantVelocityPressureDistribution.jpg (30.0 KB, 5 views)
File Type: jpg ChangingVelocityPressureGradient.jpg (31.9 KB, 5 views)
Ben0000 is offline   Reply With Quote

Old   May 17, 2023, 04:17
Default
  #2
Senior Member
 
Josh Williams
Join Date: Feb 2021
Location: Scotland
Posts: 112
Rep Power: 5
joshwilliams is on a distinguished road
A flow is driven by pressure gradients. For example, if you have incompressible flow in a pipe with pressure = 0 Pa at the outlet and 0 Pa at the inlet, the fluid within the pipe will be stationary. It is only a difference/gradient in pressure between the inlet and outlet that creates flow. So, there is no way to simulate your problem 'without a pressure gradient' if I understand you correctly.
joshwilliams 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
Postprocess: sampleDict works but creates no output folder shock77 OpenFOAM Post-Processing 14 November 15, 2021 08:27
Convergence problem of OF WUYing OpenFOAM Running, Solving & CFD 2 September 20, 2021 10:09
Fail to converge when solving with a fabricated solution zizhou FLUENT 0 March 22, 2021 06:33
AMI speed performance danny123 OpenFOAM 21 October 24, 2020 04:13
Floating point exception error lpz_michele OpenFOAM Running, Solving & CFD 53 October 19, 2015 02:50


All times are GMT -4. The time now is 10:02.