Odour source identification in a complex flow environment using a particle trajectory model

am . . Zugriffe: 3227

P11Ip   This paper presents a methodology that generates complex flow fields based on the Navier Stokes equations for fine-mode trajectory modelling. Initial u- and v-component of wind is input to the model from an onsite weather station on a 2-minute resolution.

   The perturbated flow is then calculated in an iterative fashion until the numerical solution converges. The results are then be stored and the next initial u- and v-component of wind will be input from the anemometer. Each solution is stored and converted for use in a particle trajectory model.

P. D’Abreton1, R. Ormerod1, L. Marchant2

1 Envirosuite, Brisbane, Australia
2 Envirosuite, Santiago, Chile

   Competing interests: The author has declared that no competing interests exist.

   Academic editor: Carlos N Díaz.

   Content quality: This paper has been peer reviewed by at least two reviewers. See scientific committee here

   Citation: P. D’Abreton1, R. Ormerod1, L. Marchant2, 2019, Odour source identification in a complex flow environment using a particle trajectory model, OLORES19 Conference, Santiago, Chile, www.olores.org.

   Copyright: 2019 Olores.org. Open Content Creative Commons license. It is allowed to download, reuse, reprint, modify, distribute, and/or copy articles in olores.org website, as long as the original authors and source are cited. No permission is required from the authors or the publishers.

   ISBN: 978-84-09-22553-8

   Keyword: Navier Stokes, fast solver, real time, sources

Download presentation




   This paper presents a methodology that generates complex flow fields based on the Navier Stokes equations for fine-mode trajectory modelling. Initial u- and v-component of wind is input to the model from an onsite weather station on a 2-minute resolution. The perturbated flow is then calculated in an iterative fashion until the numerical solution converges. The results are then be stored and the next initial u- and v-component of wind will be input from the anemometer. Each solution is stored and converted for use in a particle trajectory model.

   Preliminary results indicate that the complex flow model approximates results from more sophisticated models. A real-world test case also demonstrates the usefulness of this methodology as a quick tool in identifying potential odour sources.


1. Introduction

   Back trajectories are a useful tool for facilities to identify odour sources for timely management actions to be implemented, thereby minimising potential offsite complaints. However, the most commonly used trajectory models such HySplit (Stein et al, 2015) cannot accurately depict particle motion in industrial situations where air flow is almost always subject to interference from structures. In addition, the spatial and temporal resolution of the data used in these models are often too coarse to resolve the flow in these situations. Computational Fluid Dynamics (CFD) models are often used to resolve flows in these situations. However, these models are computationally intensive and cannot be used on a near real-time basis.

   A finite difference fast solver for the steady state Navier Stokes equations is presented in this paper. The results of this model are compared to those obtained from the QUIC-URB fast response urban dispersion model (Pardyjak and Brown, 2002; 2003). Finally, a test case is examined for a waste water treatment facility.


2. Materials and methods

   The Navier-Stokes equation can be written in vector form as:


   where u is the fluid velocity, p is the fluid pressure, ρ is the fluid density and ϒ is dynamic viscosity. These equations describe how the velocity, pressure, temperature, and density of a moving fluid are related.

   The numerical solution of the Navier-Stokes equations involves discretisation of simple differential equations and applying it to the continuous (Navier-Stokes) equations. The finite difference method is used for the discretisation and follows the methodology of Griebel et al. (1998). Beginning at time t = n, with given initial values for u and v, time is increased by δt until a final time (t = n+Δn) is reached. The time stepping loop commences by first performing time discretisation of the terms ∂u/∂t and ∂v/∂t in equations 1 and 2, thereby obtaining the following form:




   The derivatives in Eqn. 2a and 2b are discretised in two dimensions using finite differencing:


   Substituting equations 2a and 2b into the continuity equation and rearranging, gives the Poisson equation for pressure


   The following boundary values are required for the discretised momentum equations:

  • Buildings and structures - No-slip conditions are set for obstacles within the flow field. This means that velocities equal zero on these boundaries.
  • Inflow – u and v velocities from the weather station are explicitly given as an inflow boundary condition.
  • Outflow – the normal derivatives of u and v are set to zero at the boundary, thereby setting velocity values along the boundary to equal to their neighbouring values inside the domain.

   The complex flow model described above requires the following initial inputs:

  • Fine temporal resolution u- and v- wind components;
  • 2-dimensional array with structures identified by their heights and non-structures (i.e. ground) identified as zero; and
  • Reynolds number for each model setup.

   Comparison of near-surface wind fields generated by the Envirosuite code with those generated by QUIC-Urban (Pardyjak and Brown, 2003) for initial southwesterly flow shows excellent agreement, with both models showing establishment of a twin wake eddy downwind of the obstacle (Figure 1).

P16I10 P16I11
 Fig.1.: Streamlines of flow out the south west around a single cube using QUIC-URB (after Pardyjak and Brown 2003) (Left) and flow vectors out the south west around a single cube using the complex flow model.

   The trajectory model is based on the kinematic trajectory model presented in D’Abreton (1996) and D’Abreton and Tyson (1996). The model is Lagrangian, with atmospheric motion being described in terms of individual air parcels moving with air streams. The model uses the explicit method of integration defined by


Where x(t+dt) is the new three-dimensional particle position at t+dt, x(t) is the old position and V(t) is the parcel velocity vector (V=[u,v,w]) which is the sum of the mean and turbulent components:


The mean component (V) is obtained from the gridded three-dimensional derived wind fields, interpolated to produce u-, v- and w-components at the precise location of the particle for each time step. The trajectory model estimates the turbulent component (V’) from the standard deviations of the fluctuating velocities (σu, σv and σw). These were derived from the definition of TKE (where TKE was produced as an output from the complex flow model). Solving for σu and σv gives:


   The turbulence partition between the vertical and horizontal components is defined through the turbulence anisotropy ratio. A factor of 0.18 is assumed from the literature.

   The time for the model to converge (reach a solution) depends on the size of the domain and the spatial resolution of the mesh. Typically, a domain of 200 m x 200m (at 4 m resolution) takes approximately 5 minutes to converge.


3. Results and discussion

   A wastewater treatment facility has been selected to test whether this methodology can be used to identify odour sources in an area with complex flow. Figure 2 shows the facility with potential odour sources identifies as well as three e-nose odour monitors, which were developed and deployed by Odotech (Canada). Each e-nose incorporates 16 metal oxide sensors, some of which are in duplicate, and which respond to the presence of different groups of odorous compounds. The net result is that the e-noses respond to a broad range of odour types. Those deployed in this wastewater treatment facility had been in place for over 3 years and were not configured to provide odour fingerprints, which can help identify odour types and sources, unlike more recent units.

   Building heights were entered into the model and 2-minute wind data (between 0.9 m/s and 4.5 m/s) recorded by means of 2D sonic anemometer at the on-site meteorological station, was used to determine inflow conditions for the complex flow model. Note that turbulence measurements were not used in the model, thereby allowing the model to generate spatially varying turbulence fields.

   Figure 3 shows backward trajectories for four separate times. Trajectories for 10 April 22:42 (upper left) show that the high odour concentrations measured at EN03 is most likely due to aeration while the slightly lower measured odour at EN01 is most likely due to the secondary clarifier. Trajectories for 10 April 23:34 (upper right) show that sludge storage is responsible for the odour measured at EN01. By contrast, odour at EN03 may be due to emissions from the secondary clarifier, sludge thickening and/or sludge storage. Trajectories for 11 April 15:38 (lower left) show that the elevated odour concentration measured at EN03 may not due to emissions from the WWTP but rather a metal casting facility to the west. The odour measured at EN01 may be due to emissions from the secondary clarifier and/or aeration. The odour concentration measured at EN03 on 6 May 17:32 (lower right) may be partly due to offsite emissions to the west and emissions from the secondary clarifier. By contrast, odour measured at EN01 may be due to emissions from sludge thickening or an unidentified offsite source. Measured odour at EN02 is not due to WWTP emissions, but rather from manufacturing facilities to the west.

   More recent e-nose models available to the authors are configured to readily provide sensor-specific data and to analyse signals to identify a broad odour type – the specificity will depend on the nature of the odour, its intensity and the extent to which odours of different types are mixed in plumes.

Fig. 2.: Location of potential odour sources, weather station and H2S monitors. Buildings considered in the model are also indicated.


P16I16 P16I17
P16I18 P16I19
Fig. 3.: Backward trajectories for 10 April 22:42 (upper left), 10 April 23:34 (upper right), 11 April 15:38 (lower left) and 6 May 17:32 (lower right). Measured odour concentrations at each monitor is indicated on the figures. 


4. Conclusions

   A finite difference fast solver to estimate flow in a complex regime was presented in this paper. This model allows near real-time solution of complex flow from a single onsite weather station, with execution times within a few minutes. In contrast to existing trajectory models, the complex flow model can be used the drive a particle trajectory model to determine near-field transport in disturbed flow.

   A case study shows that in principle, odour sources can be identified from a combination of trajectories under complex flow conditions. In particular, offsite odour sources may be identified, thereby averting the need for mitigatory actions and even potential legal issues.


5. References

   D’Abreton P.C. 1996. Lagrangian kinematic and isentropic trajectory models for aerosol and trace gas transport studies in southern Africa. S. Afr. Jour. Sci. 92, 157-160.

   D’Abreton P.C. and Tyson P.D. 1996. Three-dimensional kinematic modelling of water vapour transport over southern Africa. Water SA, 22, 297-306

   Stein, A.F., Draxler, R.R, Rolph, G.D., Stunder, B.J.B., Cohen, M.D., and Ngan, F. 2015. NOAA's HYSPLIT atmospheric transport and dispersion modeling system. Bull. Amer. Meteor. Soc., 96, 2059-2077.

   Griebel, M., Dornseifer, T. and Neunhoeffer, T. 1998. Numerical Simulation in Fluid Dynamics: A Practical Introduction, Society for Industrial and Applied Mathematics. Philadelphia.

   Pardyjak, E.R. and M.J. Brown, 2002. Fast Response Modeling of two Building Street Canyon. 4th AMS Urb. Env. Symp., Norfolk, VA. Paper J1.4

   Pardyjak, E.R. and Brown, M. (2003) QUIC-URB v1.1 Theory and User’s Guide, Los Alamos National Laboratory, LA-UR- 07-3181



Kommentare (0)

There are no comments posted here yet

Einen Kommentar verfassen

  1. Posting comment as a guest.
Anhänge (0 / 3)
Share Your Location

Please note that this site uses cookies in order to work properly.

See more about our cookie policy Learn more

I understand

Please read the following to learn more about our cookies policy:


What are cookies?

   A cookie is a text file stored in a user’s web browser on any device they use to access a website that holds information regarding the user’s visit, such as preferences. When the user returns, the browser provides the cookie with the stored information to the site.

What cookies are used for?

   Cookies are used for adjusting a website’s content to fit a user’s preferences and optimize the website. They store useful information that improve the user’s experience of a website. They are most commonly used for:

  •     Recognizing the type of device a user is browsing with and any preferences applied to optimize the website specifically for the device.
  •     Creating statistics that help website owners to understand how their users interact with their website, which allows them to improve their structure and content.

What types of cookies are used?

   There are two types of cookies: persistent cookies and session cookies. Persistent cookies remain on your hard drive for a period of time specified in the cookie’s file parameters or until removed manually. When you return to a website and it requires you to login again despite previously storing your login information, it is usually because the persistent cookie expired; this helps to increase security while maintaining accessibility.

   Session cookies, on the other hand, are used temporarily and expire once the website or browser is closed. They are used to track user activity on a website during a single visit. When a website requires that you verify your age or location once every visit before allowing you to view content and without requiring additional personal details, that is a session cookie at work.

Do cookies include personal data?

   If there is a need for the collection of personal information, such as for creating accounts, then cookies may store personal information. However, it is required by data protection law that users are informed of the collection of personal data. This data will also be encrypted to render it inaccessible for unauthorized users.

Managing cookies

   By default, browsers are configured to accept cookies. However, these settings may be changed to block cookies entirely, or to inform the user each time they are used. Detailed information about cookies and the options associated with them are available in each browsers’ settings.

Which cookies does collect olores.org?

   Olores.org collect cookies for 2 purposes:

  • Register statistical data.
  • Set language preferences.

   In addition we use third party cookies through Statcounter to collect different data.

StatCounter Analytics Cookies

   StatCounter is a web analytics service. We use StatCounter to track activity on our website. These stats help us to understand how people are interacting with our website and to improve the design and functionality of our site so that we can offer a better online experience to our visitors. If you visit olores.org, a StatCounter analytics cookie (called "is_unique") may  be placed in your browser.  This cookie is used only to determine whether you are a first-time or returning visitor and to estimate unique visits to the site. No personal information is stored in the cookie.

Refuse Statcounter cookies.

You may set your browser to refuse/accept StatCounter analytics cookies by clicking here.


    • Your decision to refuse/accept StatCounter analytics cookies applies to all websites which use the StatCounter service (including the StatCounter site itself).
    • If you refuse all StatCounter analytics cookie, a refusal cookie (called "refusal_cookie") will be set to remember this preference and any existing StatCounter analytics cookies in your browser will be destroyed.
    • If you delete/remove/destroy the refusal cookie, you must revisit this page in order to re-set your preference.
    • The refusal cookie is set only for your current browser and machine. If you use multiple browsers/machines, you must set a refusal cookie in each case.
    • You can also change your cookie settings directly in your browser. Learn more about cookies and how to manage them here: http://www.allaboutcookies.org/cookies/index.html
    • Or you can learn about how to adjust cookie settings for specific browsers here: