Comments (9)
from amr-wind.
@StormMata have you tried adjusting Actuator.UniformCtDisk.diameters_to_sample
?
This variable tells it how far upstream to sample for computing the free-stream velocity. 1 Diameter is very close. From your plot it looks like you want to set this to at least 3 or 4.
Another parameter you can adjust is the number of points in r. This really just determines how much blending there is radially since the uniform ct disk has no span-wise variations.
With 5 points in R you are going to have forcing applied over an extra 24m beyond the disk (linearly decaying)
from amr-wind.
Hi @StormMata,
Also, what you're showing with AMR-Wind is generally the correct behavior in front of the turbine. There's an induction zone ahead of the rotor which causes the flow to decelerate before it hits the turbine. There's a formula from Medici et al (2011)
which gives the centerline velocity in terms of the induction factor UniformCtDisk
model, with minor deviations at high Ct:
On the downstream side there isn't a perfect analytical formula which describes the wake velocity, but it does seem like the wake deficit is a bit stronger than it should be. However, that is sensitive to the numerical setup of the problem. I should note that your mesh is relatively coarse for this problem:
amr.n_cell = 64 32 32 # Grid cells at coarsest AMRlevel
amr.max_level = 0 # Max AMR level in hierarchy
geometry.prob_lo = -1200.0 -1200.0 -1200.0
geometry.prob_hi = 2400.0 1200.0 1200.0
geometry.is_periodic = 0 1 0 # Periodicity x y z (0/1)
So that yields cell sizes of 56.25 x 75 x 75, which for a 240m rotor is pretty coarse.
Lawrence
from amr-wind.
Hi @psakievich and @lawrenceccheung,
Thanks to you both for your help. We've increased the mesh resolution (now 200 x 128 x 128) and the sampling distance (now 4.0D) and attached the updated input file for reference.
We still find with this setup that the disk velocity is about 6.4% lower than expected (both from Medici and standard momentum theory, which give the same result). Moreover, we see that the wake velocity is around 35% off of the expected value from
Input file
from amr-wind.
hi @StormMata,
Regarding the induction values, the uniform Ct disk model should get a pretty good value for the rotor disk velocities, but it won't be perfect, especially for higher Ct values. Here's the comparison that we've observed:
where the BEM theory curve is given by
Note that the rotor disks velocities depend on the numerical parameters used in AMR-Wind. One possible rule of thumb is to keep the num_points_r
parameter approximately the same as rotor diameter/mesh spacing. If you use too many or too few points then there are effects that can be observed at the edge of the rotor disk.
Another thing to note is that things also heavily depend on the epsilon parameter. I would set epsilon to be smaller than the mesh size, maybe in the range
Lawrence
from amr-wind.
@StormMata another parameter you can toy with is the actual spreading function being used. We offer two additional varieties beyond the default you are using.
A description of what the spreading functions are is here: https://www.osti.gov/servlets/purl/1870580
Slide 12 left figure is the default you are currently using
Slide 12 right figure is the Actuator.UniformCtDisk.spreading = LinearBasis
Slide 11 is the same method and is Actuator.UniformCtDisk.spreading = UniformGaussian
For the uniform_ct
disk the first two options should give the same results. These spreading functions and improving them is a line of research that has been put on the back burner for a bit. Happy to discuss further if you find any notable differences in behavior.
from amr-wind.
This issue is stale because it has been open 30 days with no activity.
from amr-wind.
@StormMata were you able to try some of the proposed fixes? If so, can you close the issue if you feel like it has been resolved?
from amr-wind.
This issue is stale because it has been open 30 days with no activity.
from amr-wind.
Related Issues (20)
- Discussion Question: Forested Areas and Terrain - Drag Based Models HOT 2
- Segfaults on initial MAC projection HOT 24
- Underestimated power predictions from OpenFAST/ALM turbine models in unsteady flows HOT 48
- Averaging (Post Processing AMR Wind) HOT 3
- Enable changing velocity in turbine simulations from precurors HOT 9
- Cannot find field: tke HOT 2
- Geostrophic Forcing flaws HOT 5
- FLLC Reintroduced bug HOT 1
- Feature request: ABL statistics on level >0 HOT 1
- Segfault of AMR-Wind + OpenFAST with oneapi HOT 2
- MMC forcing fails with GPU HOT 3
- Add time-vary & spatially varying body force for MMC simulations HOT 3
- GPU Unit test Failues: DiagnosticsTest.Max_Vel_MultiLevel and Max_MACvel_MultiLevel HOT 1
- Segmentation fault when running more than one turbine HOT 11
- Encountered NetCDF error; aborting HOT 6
- MAC velocity not getting filled correctly at the low boundary HOT 18
- Instability in some simulation HOT 10
- storing postprocessing files in another directory
- I encountered the following problems in the process of learning to use TurbineFASTLine. Is this due to the version of TurbineFastline HOT 4
- ow_w2a test is non deterministic on CPU
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from amr-wind.