scholarly journals Particulate Matter Dispersion Modeling in Agricultural Applications: Investigation of a Transient Open Source Solver

Agronomy ◽  
2021 ◽  
Vol 11 (11) ◽  
pp. 2246
Author(s):  
David Janke ◽  
Senthilathiban Swaminathan ◽  
Sabrina Hempel ◽  
Robert Kasper ◽  
Thomas Amon
Keyword(s):  
Particulate Matter ◽  
Open Source ◽  
Turbulence Models ◽  
Continuous Phase ◽  
Experimental Chamber ◽  
Dispersion Modeling ◽  
Reynolds Stresses ◽  
Agricultural Applications ◽  
Inflow Turbulence ◽  
Turbulence Generator

Agriculture is a major emitter of particulate matter (PM), which causes health problems and can act as a carrier of the pathogen material that spreads diseases. The aim of this study was to investigate an open-source solver that simulates the transport and dispersion of PM for typical agricultural applications. We investigated a coupled Eulerian–Lagrangian solver within the open source software package OpenFOAM. The continuous phase was solved using transient large eddy simulations, where four different subgrid-scale turbulence models and an inflow turbulence generator were tested. The discrete phase was simulated using two different Lagrangian solvers. For the validation case of a turbulent flow of a street canyon, the flowfield could be recaptured very well, with errors of around 5% for the non-equilibrium turbulence models (WALE and dynamicKeq) in the main regions. The inflow turbulence generator could create a stable and accurate boundary layer for the mean vertical velocity and vertical profile of the turbulent Reynolds stresses R11. The validation of the Lagrangian solver showed mixed results, with partly good agreements (simulation results within the measurement uncertainty), and partly high deviations of up to 80% for the concentration of particles. The higher deviations were attributed to an insufficient turbulence regime of the used validation case, which was an experimental chamber. For the simulation case of PM dispersion from manure application on a field, the solver could capture the influence of features such as size and density on the dispersion. The investigated solver is especially useful for further investigations into time-dependent processes in the near-source area of PM sources.

Implementation of an Efficient Synthetic Inflow Turbulence-Generator in the Open-Source Code OpenFOAM for 3D LES/DNS Applications

2021 ◽  
pp. 207-221
Author(s):  
Flavio Cesar Cunha Galeazzo ◽  
Feichi Zhang ◽  
Thorsten Zirwes ◽  
Peter Habisreuther ◽  
Henning Bockhorn ◽  
...  
Keyword(s):  
Open Source ◽  
Source Code ◽  
Open Source Code ◽  
Inflow Turbulence ◽  
Turbulence Generator

CFD Simulation of Particle Transport and Dispersion in Indoor Environment by Human Walking

2012 ◽  
Author(s):  
Iman Goldasteh ◽  
Goodarz Ahmadi ◽  
Andrea Ferro
Keyword(s):  
Particulate Matter ◽  
High Speed ◽  
Cfd Simulation ◽  
Three Dimensional ◽  
Study Data ◽  
Turbulent Dispersion ◽  
Experimental Chamber ◽  
Indoor Environments ◽  
Three Dimensional Model ◽  
Particle Resuspension

Particle resuspension is an important source of particulate matter in indoor environments that significantly affects the indoor air quality and could potentially have adverse effect on human health. Earlier efforts to investigate indoor particle resuspension hypothesized that high speed airflow generated at the floor level during the gate cycle is the main cause of particle resuspension. The resuspended particles are then assumed to be dispersed by the airflow in the room, which is impacted by both the ventilation and the occupant movement, leading to increased PM concentration. In this study, a three dimensional model of a room was developed using FLUENT™ CFD package. A RANS approach with the RNG k-ε turbulence model was used for simulating the airflow field in the room for different ventilation conditions. The trajectories of resuspended particulate matter were computed with a Lagrangian method by solving the equations of particle motion. The effect of turbulent dispersion was included with the use of the eddy lifetime model. The resuspension of particles due to gait cycle was estimated and included in the computational model. The dispersion and transport of particles resuspended from flooring as well as particle re-deposition on flooring and walls were simulated. Particle concentrations in the room generated by the resuspension process were evaluated and the results were compared with experimental chamber study data as well as simplified model predictions, and good agreement was found.

Development of High Order LES Solver for Heat Transfer Applications Based on the Open Source OpenFOAM Framework

2015 ◽  
Author(s):  
Luca Mangani ◽  
David Roos Launchbury ◽  
Ernesto Casartelli ◽  
Giulio Romanelli
Keyword(s):  
Heat Transfer ◽  
Open Source ◽  
Gas Turbines ◽  
Incompressible Flows ◽  
Stokes Equations ◽  
Turbulence Models ◽  
Local Heat Transfer ◽  
Local Heat ◽  
High Order ◽  
Computational Time

The computation of heat transfer phenomena in gas turbines plays a key role in the continuous quest to increase performance and life of both component and machine. In order to assess different cooling approaches computational fluid dynamics (CFD) is a fundamental tool. Until now the task has often been carried out with RANS simulations, mainly due to the relatively short computational time. The clear drawback of this approach is in terms of accuracy, especially in those situations where averaged turbulence-structures are not able to capture the flow physics, thus under or overestimating the local heat transfer. The present work shows the development of a new explicit high-order incompressible solver for time-dependent flows based on the open source C++ Toolbox OpenFOAM framework. As such, the solver is enabled to compute the spatially filtered Navier-Stokes equations applied in large eddy simulations for incompressible flows. An overview of the development methods is provided, presenting numerical and algorithmic details. The solver is verified using the method of manufactured solutions, and a series of numerical experiments is performed to show third-order accuracy in time and low temporal error levels. Typical cooling devices in turbomachinery applications are then investigated, such as the flow over a turbulator geometry involving heated walls and a film cooling application. The performance of various sub-grid-scale models are tested, such as static Smagorinsky, dynamic Lagrangian, dynamic one-equation turbulence models, dynamic Smagorinsky, WALE and sigma-model. Good results were obtained in all cases with variations among the individual models.

Comparative Assessment of Turbulence Models for Prediction of Flow-Induced Corrosion Damages

2011 ◽  
Author(s):  
Kaushik Das ◽  
Debashis Basu ◽  
Todd Mintz
Keyword(s):  
Corrosion Rate ◽  
Cross Sections ◽  
Volume Fraction ◽  
Solid Phase ◽  
Turbulence Models ◽  
Reynolds Stress Model ◽  
Continuous Phase ◽  
Comparative Assessment ◽  
Solid Particles ◽  
Protective Film

The present study makes a comparative assessment of different turbulence models in simulating the flow-assisted corrosion (FAC) process for pipes with noncircular cross sections and bends, features regularly encountered in heat exchangers and other pipeline networks. The case study investigates material damage due to corrosion caused by dissolved oxygen (O2) in a stainless steel pipe carrying an aqueous solution. A discrete solid phase is also present in the solution, but the transport of the solid particles is not explicitly modeled. It is assumed that the volume fraction of the solid phase is low, so it does not affect the continuous phase. Traditional two-equation models are compared, such as isotropic eddy viscosity, standard k-ε and k-ω models, shear stress transport (SST) k-ω models, and the anisotropic Reynolds Stress Model (RSM). Computed axial and radial velocities, and turbulent kinetic energy profiles predicted by the turbulence models are compared with available experimental data. Results show that all the turbulence models provide comparable results, though the RSM model provided better predictions in certain locations. The convective and diffusive motion of dissolved O2 is calculated by solving the species transport equations. The study assumes that solid particle impingement on the pipe wall will completely remove the protective film formed by corrosion products. It is also assumed that the rate of corrosion is controlled by diffusion of O2 through the mass transfer boundary layer. Based on these assumptions, corrosion rate is calculated at the internal pipe walls. Results indicate that the predicted O2 corrosion rate along the walls varies for different turbulence models but show the same general trend and pattern.

On the skin friction due to turbulence in ducts of various shapes

2018 ◽  
Vol 838 ◽  
pp. 369-378 ◽  
Author(s):  
P. R. Spalart ◽  
A. Garbaruk ◽  
A. Stabnikov
Keyword(s):  
Reynolds Number ◽  
Skin Friction ◽  
Numerical Solutions ◽  
Mathematical Reasoning ◽  
Turbulence Models ◽  
Turbulence Theory ◽  
Reynolds Stresses ◽  
Cross Sectional ◽  
Law Of The Wall ◽  
Secondary Motion

We consider fully developed turbulence in straight ducts of non-circular cross-sectional shape, for instance a square. A global friction velocity $\overline{u}_{\unicode[STIX]{x1D70F}}$ is defined from the streamwise pressure gradient $|\text{d}p/\text{d}x|$ and a single characteristic length $h$, half the hydraulic diameter (shapes with disparate length scales, due to high aspect ratio, are excluded). We reason that as the Reynolds number $Re$ reaches high values, outside the viscous region the streamwise velocity differences and the secondary motion scale with $\overline{u}_{\unicode[STIX]{x1D70F}}$ and the Reynolds stresses with $\overline{u}_{\unicode[STIX]{x1D70F}}^{2}$. This extends the classical defect-law argument, associated with Townsend and many others, and is successful in channel and pipe flows. We then posit matched asymptotic expansions with overlap of the law of the wall and the behaviour we assumed in the core region. The wall may be smooth, or have a Nikuradse roughness $k_{S}$ (such that it is fully rough, with $k_{S}^{+}\gg 1$). The consequences include the familiar logarithmic behaviour of the velocity profile, but also the surprising prediction that the skin friction tends to uniformity all around the duct, except near possible corners, asymptotically as $Re\rightarrow \infty$ or $k_{S}/h\rightarrow 0$. This is confirmed by numerical solutions for a square and two ellipses, using a conventional turbulence model, albeit the trend with Reynolds number is slow. The magnitude of the secondary motion also scales as expected, and the skin-friction coefficient follows the logarithm of the appropriate Reynolds number. This is a validation of the mathematical reasoning, but is by no means independent physical evidence, because the turbulence models embody the same assumptions as the theory. The uniformity of the skin friction appears to be a new and falsifiable deduction from turbulence theory, and a candidate for high-Reynolds-number experiments.

Improved Prediction of Turbomachinery Flows Using Near-Wall Reynolds-Stress Model

2001 ◽  
Vol 124 (1) ◽  
pp. 86-99 ◽  
Author(s):  
G. A. Gerolymos ◽  
J. Neubauer ◽  
V. C. Sharma ◽  
I. Vallet
Keyword(s):  
Experimental Data ◽  
Turbulent Flows ◽  
Reynolds Stress ◽  
Turbulence Models ◽  
Reynolds Stress Model ◽  
Stress Model ◽  
Reynolds Stresses ◽  
Mean Flow ◽  
Flow Equations ◽  
Near Wall

In this paper an assessment of the improvement in the prediction of complex turbomachinery flows using a new near-wall Reynolds-stress model is attempted. The turbulence closure used is a near-wall low-turbulence-Reynolds-number Reynolds-stress model, that is independent of the distance-from-the-wall and of the normal-to-the-wall direction. The model takes into account the Coriolis redistribution effect on the Reynolds-stresses. The five mean flow equations and the seven turbulence model equations are solved using an implicit coupled OΔx3 upwind-biased solver. Results are compared with experimental data for three turbomachinery configurations: the NTUA high subsonic annular cascade, the NASA_37 rotor, and the RWTH 1 1/2 stage turbine. A detailed analysis of the flowfield is given. It is seen that the new model that takes into account the Reynolds-stress anisotropy substantially improves the agreement with experimental data, particularily for flows with large separation, while being only 30 percent more expensive than the k−ε model (thanks to an efficient implicit implementation). It is believed that further work on advanced turbulence models will substantially enhance the predictive capability of complex turbulent flows in turbomachinery.

On Application of Nonlinear k-ε Models for Internal Combustion Engine Flows

2002 ◽  
Vol 124 (3) ◽  
pp. 668-677 ◽  
Author(s):  
G. M. Bianchi ◽  
G. Cantore ◽  
P. Parmeggiani ◽  
V. Michelassi
Keyword(s):  
Internal Combustion Engine ◽  
Eddy Viscosity ◽  
Combustion Engine ◽  
Turbulence Models ◽  
Internal Combustion ◽  
Moment Closure ◽  
Reynolds Stresses ◽  
Mean Velocity ◽  
Numerical Predictions ◽  
Viscosity Models

The linear k-ε model, in its different formulations, still remains the most widely used turbulence model for the solutions of internal combustion engine (ICE) flows thanks to the use of only two scale-determining transport variables and the simple constitutive relation. This paper discusses the application of nonlinear k-ε turbulence models for internal combustion engine flows. Motivations to nonlinear eddy viscosity models use arise from the consideration that such models combine the simplicity of linear eddy-viscosity models with the predictive properties of second moment closure. In this research the nonlinear k-ε models developed by Speziale in quadratic expansion, and Craft et al. in cubic expansion, have been applied to a practical tumble flow. Comparisons between calculated and measured mean velocity components and turbulence intensity were performed for simple flow structure case. The effects of quadratic and cubic formulations on numerical predictions were investigated too, with particular emphasis on anisotropy and influence of streamline curvature on Reynolds stresses.

A Review of Wind Turbine Polar Data and its Effect on Fatigue Loads Simulation Accuracy

2015 ◽  
Author(s):  
Matthew Lennie ◽  
Georgios Pechlivanoglou ◽  
David Marten ◽  
Christian Navid Nayeri ◽  
Oliver Paschereit
Keyword(s):  
Wind Turbine ◽  
Turbine Blades ◽  
Turbulence Models ◽  
Leading Edge ◽  
Sensitivity Study ◽  
Wind Turbine Blades ◽  
Simulation Code ◽  
Simulation Accuracy ◽  
Inflow Turbulence ◽  
Lift And Drag

To certify a Wind Turbine the standard processes set out by the GL guidelines and the IEC61400 demand a large number of simulations in order to justify the safe operation of the machine in all reasonably probable scenarios. The result of this rather demanding process is that the simulations rely on lower fidelity methods such as the Blade Element Momentum (BEM) method. The BEM method relies on a number of simplified inputs including the coefficient of lift and drag polar data (usually referred to as polars). These polars are usually either measured experimentally, generated using tools such as XFoil or, in some cases obtained using 2D CFD. It is typical to then modify these polars in order to make them suitable for aeroelastic simulations. Some of these modifications include 360° angle of attack extrapolation methods and polar modifications to account for 3D effects. Many of these modifications can be perceived to be a black art due to the manual selection of coefficients. The polars can misrepresent reality for many reasons, for example, inflow turbulence can affect measurements obtained in wind tunnels. Furthermore, on real wind turbine blades leading edge erosion can reduce performance. Simulated polars can even vary significantly due to the choice of turbulence models. Stack these effects on top of the uncertainties caused by yaw error, pitch error and dynamic stall and one can clearly see an operating environment hostile to accurate simulations. Colloquial evidence suggests that experienced designers would account for all of these sources of errors methodically, however, this is not reflected by the certification process. A review of experimental data and literature was performed to identify some of the inaccuracies in wind turbine polars. Significant variations were found between a range of 2D polar techniques and wind tunnel measurements. A sensitivity study was conducted using the aeroelastic simulation code FAST (National Renewable Energy Laboratory) with lift and drag polars sourced using different methods. The results were post-processed to give comparisons the rotor blade fatigue damage; variations in accumulated damages reached levels of 164%. This variation is not disastrous but is certainly enough to motivate a new approach for certifying the aerodynamic performance of wind turbines. Such an approach would simply see the source of polar data and all post-processing steps documented and included in the checks performed by certification bodies.

Sign in / Sign up

Export Citation Format

Share Document