Progress in Agricultural Engineering Sciences
Volume/Issue:
Volume 16: Issue 1
Comparing different CFD software with NACA 2412 airfoil
Authors:
Csaba Hetyei Óbuda University, Doctoral School on Safety and Security Sciences, Hungary

Search for other papers by Csaba Hetyei in
Current site
Google Scholar
PubMed Close
https://orcid.org/0000-0003-2915-4540
,
Ildikó Molnár Óbuda University, Donát Bánki Faculty of Mechanical and Safety Engineering, Hungary

Search for other papers by Ildikó Molnár in
Current site
Google Scholar
PubMed Close
, and
Ferenc Szlivka Óbuda University, Donát Bánki Faculty of Mechanical and Safety Engineering, Hungary

Search for other papers by Ferenc Szlivka in
Current site
Google Scholar
PubMed Close
https://orcid.org/0000-0002-3298-4142
Pages:
25–40
Online Publication Date:
22 Oct 2020
Publication Date:
12 Dec 2020
Article Category:
Research Article
DOI:
https://doi.org/10.1556/446.2020.00004
Keywords:
airfoil; CFD; finite volume method; simulation; validation
Open access

Abstract

The engineering application's design process starts with a concept, based on our theoretical knowledge and continues with a numerical simulation. In our paper, we review the finite volume method (FVM) which is used generally for heat and fluid dynamic simulations.

We compare three different computational fluid dynamics (CFD) software (based in the fine volume method) for validating a NACA airfoil, which can be used for example in the aerospace industry for an airplane's wing profile, and it can be used for example in the renewable industry for a wind turbine's blade or a water turbine's impeller profile. At the end of this paper, the result of our simulations will be compared with a validation case and the difference between the CFD software and the measured data will be presented.

Abstract

The engineering application's design process starts with a concept, based on our theoretical knowledge and continues with a numerical simulation. In our paper, we review the finite volume method (FVM) which is used generally for heat and fluid dynamic simulations.

We compare three different computational fluid dynamics (CFD) software (based in the fine volume method) for validating a NACA airfoil, which can be used for example in the aerospace industry for an airplane's wing profile, and it can be used for example in the renewable industry for a wind turbine's blade or a water turbine's impeller profile. At the end of this paper, the result of our simulations will be compared with a validation case and the difference between the CFD software and the measured data will be presented.

Introduction

In agricultural areas which are far from the electric grid, establishing a connection with the energy system is expensive and if the investment returns, it will be in a long time. For this region, an independent, off-grid renewable energy source is a good option, i.e., the wind and hydropower. The development process of this power generation equipment is the same as with other products. It starts with a preconcept which is based on our theoretical background or a previous similar product. After the preconcept, the design engineer creates a model which is tested with simulation products.

The simulation method can be different depending on the simulated physics, i.e., for stress simulation the engineer usually uses a finite element software, for fluid simulation they use a finite volume (FV) software. Besides, there are the meshless, finite difference, domain decomposition and other methods. With the simulation tools, the engineers can iteratively improve product performance and eliminate the possible errors.

For fluid dynamics collectively we called the simulation software computational fluid dynamics (CFD). In our article, we will use three different CFD software, which are based on the finite volume method (FVM). We will validate an airfoil (which is also called wing profile) with these tools. The airfoils are the principal parts of the modern lift force-based equipment, i.e., the previously indicated wind turbines' blade and water turbines' impeller. We will reflect the difference with this validation process within the same method with different software.

Briefly about airfoils

Nowadays the geometry of wind turbines blades is based on wing profiles used in the aerospace industry and the impellers use the lift force favorably. For this process, the turbines' blade is made of airfoils, which are known as wing profiles in the aerospace industry. If the working medium is water, the airfoil's name changes to hydrofoil referring to the liquid medium.

Perhaps one of the most well-known airfoil families is the National Advisory Committee for Aeronautics (NACAs) wing profiles, for which a large amount of measurement results are available. The NACA profiles have a four- or five-digit code, which represents the geometry. In the four-digit system,

  1. The first digit shows the maximum asymmetry (called chamber) between the two acting surfaces of an airfoil in the percentage of the cord line length.

  2. The second digit shows where the point defined with the first digit can be found on the cord line. Also, this point is where the cord line and the chamber mean line (green line on the next figure) have the greatest distance.

  3. The third and fourth show the percentage of the length and thickness (Abbott & von Doenhoff, 1959).

Fig. 1 presents a NACA 2412, which we chose for simulation.

Fig. 1.
Fig. 1.

NACA 2412 airfoil

Citation: Progress in Agricultural Engineering Sciences 16, 1; 10.1556/446.2020.00004

If we put the airfoil in the position seen in the next figures (Figs. 2 and 3) and the direction of the free flow is horizontal, the pressure on the upper part of the airfoil decreases, and the velocity increases. In the lower part, the pressure increases while wind velocity drops. The pressure difference between the sides generates the lift.

Fig. 2.
Fig. 2.

Pressure distribution near an airfoil (in a blue to red scale, the blue is a lowest and red is higher pressure) (For interpretation of the references to colour in this figure legend, the reader is referred to the web version of this article.)

Citation: Progress in Agricultural Engineering Sciences 16, 1; 10.1556/446.2020.00004

Fig. 3.
Fig. 3.

Velocity distribution near an airfoil (in a blue to red scale, the blue is a lowest and red is higher velocity) (For interpretation of the references to color in this figure legend, the reader is referred to the web version of this article.)

Citation: Progress in Agricultural Engineering Sciences 16, 1; 10.1556/446.2020.00004

The lift force can be calculated with the following equation:
FL=12×cL×ρ×A×v∞2

In the equation FL is the lift force, cL is the lift coefficient, ρ is the density, A is the area of the blade's cross-section, v∞2 is the velocity in free stream.

Shortly about numerical fluid dynamics

On the CFD software's market there are a large number of CFD software with a different mathematical background. We used three CFD software (CFX from ANSYS, OpenFOAM from OpenFOAM Foundation, and FloEFD from Mentor Graphics) based on FVMs.

The FVM based software divides the computational domain into FVs. Pressure, velocity, and temperature can be defined with the help of the conservation laws (mass, momentum, and energy). This calculation typically is based on the following transport equation:
∂∂t∫VUdV+∮AFdA−∫VSVdV−∮ASAdA=R

In the equation the ∂U/∂t is the time-dependent member (∂U/∂t=0 the steady-state), U is a conserved quantity (e.g., mass, impulse, energy), F is the same quantity's flux, Sv is the U quantity's volumetric source, SA is the U quantity's surface source, V is the arbitrary enclosed control volume, A is the surface of this control volume and R is the residual (error of the equation).

For the computational domain's discretized part (each mesh cells), the solver starts to compute with an initial value and iteratively recalculates the flow field. In this calculation, the software analyses the residual quantity (R), which should decrease and go towards 0 or a prescribed zero close value. After the computation, the solver calculates each physical value, in Reynold–Averaged Navier–Stokes method (RANS), the main quantities are the velocity, pressure, and temperature.

Turbulence

In the engineering application, the flow characteristics are usually turbulence phenomena. The turbulence is a transient, random, chaotic, diffusive, and spatial (3D) phenomenon. It can occur due to near the wall's boundary layer, in the mixing flow shear layer and the unstable density stratification. The development of a turbulence flow is shown in Fig. 4.

Fig. 4.
Fig. 4.

Laminar-turbulent transition (COMSOL Blog)

Citation: Progress in Agricultural Engineering Sciences 16, 1; 10.1556/446.2020.00004

In the CFD software, the turbulence can be modeled in several ways. The most computational resource-intensive method is not to use any turbulence model and wall function (which describe the flow in the near-wall region). In this case, the software calculates the viscous sublayer, the buffer layer and the turbulent region, where the mesh size should be smaller than the smallest eddy that influences the flow field and its properties.

There is a less computational resource method, which is to use a turbulent model. The most common RANS transport equation-based turbulence models are the two-equation models, for example, the k–ϵ and k–ω models. Each model has advantages and disadvantages, therefore Menter in 1993 created the SST (Shear Stress Transport) model, which includes the benefits of both models.

When we use a wall function to describe the flow field in the near-wall regions, the computational requirements are much less than in the previous case.

k–ϵ turbulence model

The k–ϵ turbulence model explains the flow field variability with two equations. The first equation is the k, which describes the kinetic energy, the second is the ϵ, for the turbulence dissipation energy.

The benefit of the model is the need for low computing resources, and it was one of the first turbulence models, so it was validated with a lot of experiments and its limits are well known.

The disadvantage of the model is that it shows poor performance for fully developed flows in noncircular pipes, for rotating flows, for flow with large strains, and for mixing layers (Jakubík & Feszty, 2016; Kristóf, 2019).

The most common k–ϵ models are the Standard, the Realisable, and the RNG versions. We used the Standard k–ϵ model for our simulation.

k–ω turbulence model

The k–ω model is a two-equation eddy-viscosity turbulence model like the k–ϵ. In this case, instead of ϵ, the solver computes ω, which is the specific rate of dissipation.

The model shows better performance near the walls, but its results highly depend on the turbulence intensity of the boundary conditions (BCs), which is a disadvantage compared to the k–ϵ (Jakubík & Feszty, 2016).

SST turbulence model

The SST turbulence model is a two-equation eddy-viscosity turbulence model, which used k–ω near the walls and k–ϵ for free stream. The model used the benefits of each model, therefore it describes the fluid field more accurately than the other two-equation models. Its known disadvantage is the overestimation of the turbulence intensity in the large strain regions (e.g., at the stagnation point which is shown in Fig. 5), but it is less like the overestimation of the k–ϵ (CFD Online).

Fig. 5.
Fig. 5.

Boundary-layer separation (Liberal Dictionary)

Citation: Progress in Agricultural Engineering Sciences 16, 1; 10.1556/446.2020.00004

Simulations properties

Computational domain

We used a 60 × 61 × 1 m domain with plain flow (2D) simplification for the simulation. Fig. 6 shows the whole domain with the airfoil.

Fig. 6.
Fig. 6.

Computational domain with the NACA airfoil

Citation: Progress in Agricultural Engineering Sciences 16, 1; 10.1556/446.2020.00004

Mesh for Mentor Graphics FloEFD

The FloEFD uses its own Cartesian mesh process, which generates structured hexahedron elements, whose edges are aligned with the global X, Y, and Z coordinate axes.

The mesher is using the cell mating technology for finer mesh size and it uses the cut-cell method for the elements near the boundary of the walls. The cell mating process finer the grid by splitting the cells into the half, in 3D it generates 8 cells from 1 cell (halved by the X, Y, and Z axis). The cut cell method divides the mesh elements near-wall in two pieces and generates polyhedral elements on the fluid and the solid sides, using Mentor Graphics' SmartCells technology to describe the flow field for this control volumes (Mentor Graphics, 2019). Fig. 7 shows these methods (the colors from the blue to red represent the cell mating's steps).

Fig. 7.
Fig. 7.

Cell mating method with three dividing steps near a wing profile

Citation: Progress in Agricultural Engineering Sciences 16, 1; 10.1556/446.2020.00004

The FloEFDs mesh grid is in each case (3D and 2D) a volumetric mesh, so in our study, the initial mesh (represented as blue in the previous figure) thickness has only one element.

The simulation may give different results depending on the mesh density, therefore for 0° angle of attack (AoA) we ran mesh dependency simulations, where we followed the convergence of the lift coefficient (shown in Fig. 8).

Fig. 8.
Fig. 8.

Mesh dependency analysis (cL versus mesh density)

Citation: Progress in Agricultural Engineering Sciences 16, 1; 10.1556/446.2020.00004

As a result of the analysis, we chose the 1,520,510 number element settings. The final mesh can be found in the next three figures (Figs. 9, 10 and 11).

Fig. 9.
Fig. 9.

Final finite volume (FV) mesh for the whole domain (FloEFD)

Citation: Progress in Agricultural Engineering Sciences 16, 1; 10.1556/446.2020.00004

Fig. 10.
Fig. 10.

Final FV mesh near the NACA airfoil (FloEFD)

Citation: Progress in Agricultural Engineering Sciences 16, 1; 10.1556/446.2020.00004

Fig. 11.
Fig. 11.

Final FV mesh near the airfoil's leading edge (FloEFD)

Citation: Progress in Agricultural Engineering Sciences 16, 1; 10.1556/446.2020.00004

Mesh for ANSYS CFX and OpenFOAM

For the CFX and OpenFOAM, we were able to use the same mesh. The mesh can be converted from the ANSYS products to the Cfx4ToFoam, fluentMeshToFoam and fluent3DMeshToFoam commands to OpenFOAM, and with the FoamMeshToFluent we convert the grid from OpenFOAM to ANSYS Fluent, which can be used for CFX.

For this grid, we used an unstructured, deformable surface followed mesh, which contained hexahedron and prism elements. Each software uses a volumetric mesh for 3D and 2D simulations, so the final mesh was a one element thick mesh.

Like previously, we started the validation process with a mesh dependency analysis for the AoA = 0° case, whose results are shown in Fig. 12.

Fig. 12.
Fig. 12.

Mesh dependency analysis for CFXs and OpenFOAMs mesh

Citation: Progress in Agricultural Engineering Sciences 16, 1; 10.1556/446.2020.00004

The lift coefficient started converging very quickly, and in the last five cases, the difference was below 1% compared to the previous simulation. We did not see this fast convergence in the case of the drag coefficient, therefore we continue to finer the cells near the airfoil's region. For the grid which contain 152,578 elements the drag coefficient change was decreased to 1.7% (compared with the previous grid), and its value does not change significantly for the last grid, which contains 377,644 cells. As a result of these simulations, we chose the 152,578 elements mesh, which can be seen in the next three figures (Figs. 13, 14 and 15).

Fig. 13.
Fig. 13.

Final finite volume (FV) mesh for the whole domain (OpenFOAM, CFX)

Citation: Progress in Agricultural Engineering Sciences 16, 1; 10.1556/446.2020.00004

Fig. 14.
Fig. 14.

Final FV mesh near the NACA airfoil (OpenFOAM, CFX)

Citation: Progress in Agricultural Engineering Sciences 16, 1; 10.1556/446.2020.00004

Fig. 15.
Fig. 15.

Final FV mesh near the airfoil's leading edge (OpenFOAM, CFX)

Citation: Progress in Agricultural Engineering Sciences 16, 1; 10.1556/446.2020.00004

Boundary conditions (BCs)

The air entered the flow field from the left side. The BCs of the other sides were open, where there is atmospheric pressure. The BC of the wing's surface was no-slip wall and the other two surfaces were symmetric planes.

In the three different software, the open condition was prescribed in two different ways, in FloEFD, we defined the absolute environment pressure (101,325 Pa), for the other software we defined a relative pressure (0 Pa). In the FloEFD and CFX, the symmetry was made with the symmetric condition. For the OpenFOAM, in 2D we have to use the “empty” BC instead of symmetry, which is in this CFD software a 3D BC.

The inlet velocity was 42.89 m/s. We changed the inlet vector between 0 and 10°. The Reynold number was RE = 2.85×106.

In FloEFD and CFX, we used the “Air” fluid from the software material database. In the OpenFOAM, there is no material database, so we defined the material with the following properties: ρ = 1.2041 kg/m3, v = 1.511083 × 10−5 m2/s.

Stopping criteria

The software solves Eq. (2) for each mesh elements and monitors the U quantity. Based on the continuity laws, the entering flow should leave the flow fields and no residual fluid should remain in the system. The FVM based CFD software iteratively monitor this residual U value and the calculations run again and again until the residual value decreases near the required 0 value.

For each CFD software, we use the default quantities (velocity [u, v, and w], pressure, turbulence [k and ϵ or ω] and the temperature) and the lift and drag coefficients and forces. For the lift coefficient, we used the following equation (which is the (1) equation rearranged).
cL=FL12ρAv∞2

In Eq. (3) cL is the lift coefficient, FL is the lift force, ρ is the density, A is the area of the wing's cross-section (which is, in this case, the cord length multiplied by the wings' thickness (Gulyás,2016)), v∞ is the air velocity in the free-steam. For the AoA = 6°, the changing of the lift coefficient is shown in the next figure (Fig. 16) during the first iteration to the last iteration.

Fig. 16.
Fig. 16.

cL's changing versus the iterations (AoA = 6°)

Citation: Progress in Agricultural Engineering Sciences 16, 1; 10.1556/446.2020.00004

Having compared the simulations, we found that the change of the cLs value decreased below 10−5 between 400th and 500th iterations, below 10−6 between 550th and 600th iterations and below 10−7 between 700th and 750th iterations. Therefore, we set the maximum iteration to 1,100.

Comparing steady-state with transient simulation

Validating the NACA 2412 airfoil is a time-dependent process, so it should be modeled with a time-dependent simulation. We ran a steady-state and a transient simulation for 0 and 10° AoA. The results of the simulations are shown in the next figures (Figs. 17, 18, 19 and 20) and in Table 1.

Fig. 17.
Fig. 17.

Mach number distribution near the wing for 0° (left) and 10° (right) AoA (steady-state)

Citation: Progress in Agricultural Engineering Sciences 16, 1; 10.1556/446.2020.00004

Fig. 18.
Fig. 18.

Mach number distribution near the wing for 0° (left) and 10° (right) AoA (transient; t = 1 s)

Citation: Progress in Agricultural Engineering Sciences 16, 1; 10.1556/446.2020.00004

Fig. 19.
Fig. 19.

Pressure distribution near the wing for 0° (left) and 10° (right) AoA (steady-state)

Citation: Progress in Agricultural Engineering Sciences 16, 1; 10.1556/446.2020.00004

Fig. 20.
Fig. 20.

Pressure distribution near the wing for 0° (left) and 10° (right) AoA (transient; t = 1 s)

Citation: Progress in Agricultural Engineering Sciences 16, 1; 10.1556/446.2020.00004

Table 1.

Steady-state simulation versus transient simulation results

0°0°Difference (%)10°10°Difference (%)
Stead-stateTransient (t = 1 s)Steady-stateTransient (t = 1 s)
Fx [N]0.0083970.0085792.17−2.317583−2.3123760.22
Fy [N]2.6166332.5991460.6713.93539113.9160180.14
FL [N]2.6159232.5984370.6714.12729314.1073660.14
cL0.2321700.2306180.671.2538361.2520680.14

In Table 1, Fx is the force on the wing's surface (X component), Fy is the force on the wing's surface (Y component), FL is the lift force and cL is the lift coefficient.

Comparing the flow field and the results in Table 1 (especially the lift coefficient), the differences between the steady-state and transient results are negligible. The consequence of running the simulation in steady-state is that we are unable to examine the turbulent phenomenon's time dependency, only the average value of the eddies can be examined.

Simulations' results

We ran our simulations between 0° and 10° with 2° steps. The next figures (Figs. 21–24) show the AoA = 0° and the AoA = 8° the cases' velocity field.

Fig. 21.
Fig. 21.

Velocity filed near the airfoil for AoA = 0° with k–ϵ turbulence model (CFX, OpenFOAM, FloEFD)

Citation: Progress in Agricultural Engineering Sciences 16, 1; 10.1556/446.2020.00004

Fig. 22.
Fig. 22.

Velocity field near the airfoil for AoA = 0° with SST turbulence model (CFX, OpenFOAM)

Citation: Progress in Agricultural Engineering Sciences 16, 1; 10.1556/446.2020.00004

Fig. 23.
Fig. 23.

Velocity field near the airfoil for AoA = 8° with k–ϵ turbulence model (CFX, OpenFOAM, FloEFD)

Citation: Progress in Agricultural Engineering Sciences 16, 1; 10.1556/446.2020.00004

Fig. 24.
Fig. 24.

Velocity field near the airfoil for AoA = 8° with SST turbulence model (CFX, OpenFOAM)

Citation: Progress in Agricultural Engineering Sciences 16, 1; 10.1556/446.2020.00004

The lift coefficients can be found in Fig. 25.

Fig. 25.
Fig. 25.

Lift coefficients versus AoA

Citation: Progress in Agricultural Engineering Sciences 16, 1; 10.1556/446.2020.00004

Examining the diagram, the lift coefficients are approximately the same at 0–6°, above the CFXs, and OpenFOAMs results, remain the same, while the Mentor Graphics FloEFDs cL values were a bit higher than the others.

Comparison of results

We would like to compare our results with the measured data. We chose three different measurements and we made one averaged value about them. These values are shown in Fig. 26.

Fig. 26.
Fig. 26.

Lift coefficients depending on the angle of attack (Abbott & von Doenhoff, 1959; Matsson et al., 2016; Seetharam et al., 1977)

Citation: Progress in Agricultural Engineering Sciences 16, 1; 10.1556/446.2020.00004

The three cL series were measured with different free stream velocities, but they could be compared without any correction because the lift coefficient was independent of Reynold's number (Abbott & von Doenhoff, 1959). We can conclude that the measured value is similar to Abbot's and Seetharam's measurements. There was only a significant difference over 10°. Matsson's results below 10° (except for AoA = 8°) are similar, but the values are less than in the two other cases.

In Fig. 27, we can see the comparison of the simulations' results and the average of the measured lift coefficients.

Fig. 27.
Fig. 27.

Simulated results and the average of the measured lift coefficients

Citation: Progress in Agricultural Engineering Sciences 16, 1; 10.1556/446.2020.00004

In Fig. 28, we can see the difference between the simulations' results and the measured average lift coefficients.

Fig. 28.
Fig. 28.

Difference between the simulated results and the average of the measured values

Citation: Progress in Agricultural Engineering Sciences 16, 1; 10.1556/446.2020.00004

In the two previous figures, the lift coefficient deviated with a couple of percentage in the 0–6° AoA range for the CFX and the OpenFOAM, while for the FloEFD the difference in this range was between 2.47 and 7.44%. In the 8 and 10° for each software and each turbulence model, the deviation increased up to 11.17%.

In the case of FloEFD, the lift value was higher, for OpenFOAM with SST, the lift value was always lower in each case than the averaged measured values. In the case of other software and turbulence models, we could not find any regularities.

Conclusions

In this article, we reviewed how to calculate lift coefficients for wing profiles. Then, we briefly described the FVM and our chosen software's meshers and then we presented the results of our simulations.

Our simulation was about to run validation for a NACA 2412 profile in the range of 0–10°. Comparing the results in steady-state (stationery) and in transient (time-dependent, unsteady state) for AoA = 0° and AoA = 10°, the difference was negligible, therefore we ran our simulation in steady-state.

With our simulation, we established that with the chosen method and software in the case of the lower AoA, the lifting coefficient can be simulated with lower differences, in the case of the larger AoAs (in this case 8° and above) transient simulation is necessary for accurate results.

References

  • Abbott, I.H. and von Doenhoff, A.E. (1959). Theory of wing sections: including a summary of airfoil data, Dover Publications Inc., New York, ISBN: .

    • Search Google Scholar
    • Export Citation

  • CFD Online: SST k-omega model, https://www.cfd-online.com/Wiki/SST_k-omega_model(09 May 2019).

  • COMSOL Blog: Which Turbulence Model Should I Choose for My CFD Application? https://www.comsol.com/blogs/which-turbulence-model-should-choose-cfd-application/(22 August 2018).

    • Search Google Scholar
    • Export Citation

  • Gulyás, A. (2016). H10 – Mérési segédlet (Egyszerűsített szárnyprofil geometriák áramlási paramétereinek vizsgálata) [In English: Gulyás, A. (2016) H10 -Measurement guide (Examination of flow parameters of simplified wing profile geometries)], BME, Áramlástan Tanszék, http://www.ara.bme.hu/neptun/BMEGEATMG01/2016-2017-II/labor/MSc_H10.pdf. p. 5 (26 August 2018).

    • Search Google Scholar
    • Export Citation

  • Jakubík, T. and Feszty, D. (2016). Computational fluid dynamics – 6. Turbulence modelling in CFD, Széchenyi University, Audi Department of Vehicle Engineering, https://www.feszty.com/uploads/3/6/7/8/3678219/alkalmazott%C3%81raml%C3%A1stan_week7_eng.pdf(25 June 2019).

    • Search Google Scholar
    • Export Citation

  • Kristóf, G. (2019). Áramlások numerikus modellezése [In English: Kristóf, G. (2019) Numerical modeling of flows], Akadémia Kiadó. ISBN: , https://doi.org/10.1556/9789634544128.

    • Crossref
    • Search Google Scholar
    • Export Citation

  • Liberal Dictionary: Boundary-layer separation, http://www.liberaldictionary.com/wp-content/uploads/2019/01/adverse-pressure-gradient-0472.jpg(22 June 2019).

    • Search Google Scholar
    • Export Citation

  • Matsson, J.E., Voth, J.A., McCain, C.A., and McGraw, C. (2016). Aerodynamic performance of the NACA 2412 airfoil at Low Reynolds Number, 2016 A SEE Annual Conference & Exposition, New Orleans.

    • Crossref
    • Search Google Scholar
    • Export Citation

  • Mentor Graphics (2019). FloEFD 2019.1 Technical Reference, pp. 70, 89.

  • Seetharam, H.C., Rodgers, E.J., and Wentz, W.H. Jr. (1977). Experimental studies of flow separation of the NACA 2412 airfoil at low speeds, NASA-CR-197497 19950002355, https://ntrs.nasa.gov/search.jsp?R=19950002355(25 August 2018).

    • Search Google Scholar
    • Export Citation

  • Abbott, I.H. and von Doenhoff, A.E. (1959). Theory of wing sections: including a summary of airfoil data, Dover Publications Inc., New York, ISBN: .

    • Search Google Scholar
    • Export Citation

  • CFD Online: SST k-omega model, https://www.cfd-online.com/Wiki/SST_k-omega_model(09 May 2019).

  • COMSOL Blog: Which Turbulence Model Should I Choose for My CFD Application? https://www.comsol.com/blogs/which-turbulence-model-should-choose-cfd-application/(22 August 2018).

    • Search Google Scholar
    • Export Citation

  • Gulyás, A. (2016). H10 – Mérési segédlet (Egyszerűsített szárnyprofil geometriák áramlási paramétereinek vizsgálata) [In English: Gulyás, A. (2016) H10 -Measurement guide (Examination of flow parameters of simplified wing profile geometries)], BME, Áramlástan Tanszék, http://www.ara.bme.hu/neptun/BMEGEATMG01/2016-2017-II/labor/MSc_H10.pdf. p. 5 (26 August 2018).

    • Search Google Scholar
    • Export Citation

  • Jakubík, T. and Feszty, D. (2016). Computational fluid dynamics – 6. Turbulence modelling in CFD, Széchenyi University, Audi Department of Vehicle Engineering, https://www.feszty.com/uploads/3/6/7/8/3678219/alkalmazott%C3%81raml%C3%A1stan_week7_eng.pdf(25 June 2019).

    • Search Google Scholar
    • Export Citation

  • Kristóf, G. (2019). Áramlások numerikus modellezése [In English: Kristóf, G. (2019) Numerical modeling of flows], Akadémia Kiadó. ISBN: , https://doi.org/10.1556/9789634544128.

    • Crossref
    • Search Google Scholar
    • Export Citation

  • Liberal Dictionary: Boundary-layer separation, http://www.liberaldictionary.com/wp-content/uploads/2019/01/adverse-pressure-gradient-0472.jpg(22 June 2019).

    • Search Google Scholar
    • Export Citation

  • Matsson, J.E., Voth, J.A., McCain, C.A., and McGraw, C. (2016). Aerodynamic performance of the NACA 2412 airfoil at Low Reynolds Number, 2016 A SEE Annual Conference & Exposition, New Orleans.

    • Crossref
    • Search Google Scholar
    • Export Citation

  • Mentor Graphics (2019). FloEFD 2019.1 Technical Reference, pp. 70, 89.

  • Seetharam, H.C., Rodgers, E.J., and Wentz, W.H. Jr. (1977). Experimental studies of flow separation of the NACA 2412 airfoil at low speeds, NASA-CR-197497 19950002355, https://ntrs.nasa.gov/search.jsp?R=19950002355(25 August 2018).

    • Search Google Scholar
    • Export Citation
  • Collapse
  • Expand
Cover Progress in Agricultural Engineering Sciences
Progress in Agricultural Engineering Sciences
Print ISSN:
1786-335X
Online ISSN:
1787-0321

Senior editors

Editor(s)-in-Chief: Felföldi, József

Chair of the Editorial Board Szendrő, Péter

Editorial Board

  • Beke, János (Szent István University, Faculty of Mechanical Engineerin, Gödöllő – Hungary)
  • Fenyvesi, László (Szent István University, Faculty of Mechanical Engineering, Gödöllő – Hungary)
  • Szendrő, Péter (Szent István University, Faculty of Mechanical Engineering, Gödöllő – Hungary)
  • Felföldi, József (Szent István University, Faculty of Food Science, Budapest – Hungary)

 

Advisory Board

  • De Baerdemaeker, Josse (KU Leuven, Faculty of Bioscience Engineering, Leuven - Belgium)
  • Funk, David B. (United States Department of Agriculture | USDA • Grain Inspection, Packers and Stockyards Administration (GIPSA), Kansas City – USA
  • Geyer, Martin (Leibniz Institute for Agricultural Engineering and Bioeconomy (ATB), Department of Horticultural Engineering, Potsdam - Germany)
  • Janik, József (Szent István University, Faculty of Mechanical Engineering, Gödöllő – Hungary)
  • Kutzbach, Heinz D. (Institut für Agrartechnik, Fg. Grundlagen der Agrartechnik, Universität Hohenheim – Germany)
  • Mizrach, Amos (Institute of Agricultural Engineering. ARO, the Volcani Center, Bet Dagan – Israel)
  • Neményi, Miklós (Széchenyi University, Department of Biosystems and Food Engineering, Győr – Hungary)
  • Schulze-Lammers, Peter (University of Bonn, Institute of Agricultural Engineering (ILT), Bonn – Germany)
  • Sitkei, György (University of Sopron, Institute of Wood Engineering, Sopron – Hungary)
  • Sun, Da-Wen (University College Dublin, School of Biosystems and Food Engineering, Agriculture and Food Science, Dublin – Ireland)
  • Tóth, László (Szent István University, Faculty of Mechanical Engineering, Gödöllő – Hungary)

Prof. Felföldi, József
Institute: MATE - Hungarian University of Agriculture and Life Sciences, Institute of Food Science and Technology, Department of Measurements and Process Control
Address: 1118 Budapest Somlói út 14-16
E-mail: felfoldi.jozsef@uni-mate.hu

Indexing and Abstracting Services:

  • CABI
  • ERIH PLUS
  • SCOPUS

2023  
Scopus  
CiteScore 1.8
CiteScore rank Q2 (General Agricultural and Biological Sciences)
SNIP 0.497
Scimago  
SJR index 0.258
SJR Q rank Q3

Progress in Agricultural Engineering Sciences
Publication Model Hybrid
Submission Fee none
Printed Color Illustrations 40 EUR (or 10 000 HUF) + VAT / piece
Article Processing Charge 900 EUR/article (only for OA publications)
Regional discounts on country of the funding agency World Bank Lower-middle-income economies: 50%
World Bank Low-income economies: 100%
Further Discounts Editorial Board / Advisory Board members: 50%
Corresponding authors, affiliated to an EISZ member institution subscribing to the journal package of Akadémiai Kiadó: 100%
Subscription fee 2025 Online subsscription: 173 EUR / 190 USD
Print + online subscription: 200 EUR / 220 USD
Subscription Information Online subscribers are entitled access to all back issues published by Akadémiai Kiadó for each title for the duration of the subscription, as well as Online First content for the subscribed content.
Purchase per Title Individual articles can be purchased at the prices indicated.

Progress in Agricultural Engineering Sciences
Language English
Size B5
Year of
Foundation		 2004
Volumes
per Year		 1
Issues
per Year		 1
Founder Magyar Tudományos Akadémia  
Founder's
Address		 H-1051 Budapest, Hungary, Széchenyi István tér 9.
Publisher Akadémiai Kiadó
Publisher's
Address		 H-1117 Budapest, Hungary 1516 Budapest, PO Box 245.
Responsible
Publisher		 Chief Executive Officer, Akadémiai Kiadó
ISSN 1786-335X (Print)
ISSN 1787-0321 (Online)

Monthly Content Usage

Abstract Views Full Text Views PDF Downloads
Apr 2024 0 103 80
May 2024 0 130 48
Jun 2024 0 158 47
Jul 2024 0 53 34
Aug 2024 0 126 53
Sep 2024 0 102 34
Oct 2024 0 99 29