Second international workshop on theory and applications of process visualization

The representation of business process models has been a continuing research topic for many years now. However, many process model representations have not developed beyond minimally interactive 2D icon-based representations of directed graphs and networks, with little or no annotation for information overlays. In addition, very few of these representations have undergone a thorough analysis or design process with reference to psychological theories on data and process visualization. This dearth of visualization research, we believe, has led to problems with BPM uptake in some organizations, as the representations can be difficult for stakeholders to understand, and thus remains an open research question for the BPM community. In addition, business analysts and process modeling experts themselves need visual representations that are able to assist with key BPM life cycle tasks in the process of generating optimal solutions. With the rise of desktop computers and commodity mobile devices capable of supporting rich interactive 3D environments, we believe that much of the research performed in computer human interaction, virtual reality, games and interactive entertainment have much potential in areas of BPM; to engage, provide insight, and to promote collaboration amongst analysts and stakeholders alike. We believe this is a timely topic, with research emerging in a number of places around the globe, relevant to this workshop. This is the second TAProViz workshop being run at BPM. The intention this year is to consolidate on the results of last year's successful workshop by further developing this important topic, identifying the key research topics of interest to the BPM visualization community.

Detecting anomalous events at railway level crossings

Collisions between pedestrians and vehicles continue to be a major problem throughout the world. Pedestrians trying to cross roads and railway tracks without any caution are often highly susceptible to collisions with vehicles and trains. Continuous financial, human and other losses have prompted transport related organizations to come up with various solutions addressing this issue. However, the quest for new and significant improvements in this area is still ongoing. This work addresses this issue by building a general framework using computer vision techniques to automatically monitor pedestrian movements in such high-risk areas to enable better analysis of activity, and the creation of future alerting strategies. As a result of rapid development in the electronics and semi-conductor industry there is extensive deployment of CCTV cameras in public places to capture video footage. This footage can then be used to analyse crowd activities in those particular places. This work seeks to identify the abnormal behaviour of individuals in video footage. In this work we propose using a Semi-2D Hidden Markov Model (HMM), Full-2D HMM and Spatial HMM to model the normal activities of people. The outliers of the model (i.e. those observations with insufficient likelihood) are identified as abnormal activities. Location features, flow features and optical flow textures are used as the features for the model. The proposed approaches are evaluated using the publicly available UCSD datasets, and we demonstrate improved performance using a Semi-2D Hidden Markov Model compared to other state of the art methods. Further we illustrate how our proposed methods can be applied to detect anomalous events at rail level crossings.

Porous graphene and nanomeshes

This chapter briefly introduces the concepts and modeling of gas/isotope separation by two dimensional carbon frameworks, i.e. porous graphene and carbon nanomeshes, on the basis of reviewing recent literatures. The small size of evenly distributed pores on these carbon frameworks make them ideal not only for the separation of small gas molecules but also for isotope separation by utilizing the different zero point energies induced by confinement of the pores. The related simulations were treated by transition state theory, an affordable yet precise method that could be adopted in combination with different levels of theory. Such method could be employed to evaluate the performance, as well as to aid the design, of other 2D carbon frameworks toward the goal of gas/isotope separation in the future.

Computational analysis of laminated glass panels under blast loads: A comparison of two dimensional and three dimensional modelling approaches

This paper illustrates the use of finite element (FE) technique to investigate the behaviour of laminated glass (LG) panels under blast loads. Two and three dimensional (2D and 3D) modelling approaches available in LS-DYNA FE code to model LG panels are presented. Results from the FE analysis for mid-span deflection and principal stresses compared well with those from large deflection plate theory. The FE models are further validated using the results from a free field blast test on a LG panel. It is evident that both 2D and 3D LG models predict the experimental results with reasonable accuracy. The 3D LG models give slightly more accurate results but require considerably more computational time compared to the 2D LG models.

Separable and non-separable discrete wavelet transform based texture features and image classification of breast thermograms

Highly sensitive infrared cameras can produce high-resolution diagnostic images of the temperature and vascular changes of breasts. Wavelet transform based features are suitable in extracting the texture difference information of these images due to their scale-space decomposition. The objective of this study is to investigate the potential of extracted features in differentiating between breast lesions by comparing the two corresponding pectoral regions of two breast thermograms. The pectoral regions of breastsare important because near 50% of all breast cancer is located in this region. In this study, the pectoral region of the left breast is selected. Then the corresponding pectoral region of the right breast is identified. Texture features based on the first and the second sets of statistics are extracted from wavelet decomposed images of the pectoral regions of two breast thermograms. Principal component analysis is used to reduce dimension and an Adaboost classifier to evaluate classification performance. A number of different wavelet features are compared and it is shown that complex non-separable 2D discrete wavelet transform features perform better than their real separable counterparts.

The function and mechanisms of action of ghrelin and obestatin in ovarian cancer

In this study, we have demonstrated that the preproghrelin derived hormones, ghrelin and obestatin, may play a role in ovarian cancer. Ghrelin and obestatin stimulated an increase in cell migration in ovarian cancer cell lines and may play a role in cancer progression. Ovarian cancer is the leading cause of death among gynaecological cancers and is the sixth most common cause of cancer-related deaths in women in developed countries. As ovarian cancer is difficult to diagnose at a low tumour grade, two thirds of ovarian cancers are not diagnosed until the late stages of cancer development resulting in a poor prognosis for the patient. As a result, current treatment methods are limited and not ideal. There is an urgent need for improved diagnostic markers, as well better therapeutic approaches and adjunctive therapies for this disease. Ghrelin has a number of important physiological effects, including roles in appetite regulation and the stimulation of growth hormone release. It is also involved in regulating the immune, cardiovascular and reproductive systems and regulates sleep, memory and anxiety, and energy metabolism. Over the last decade, the ghrelin axis, (which includes the hormones ghrelin and obestatin and their receptors), has been implicated in the pathogenesis of many human diseases and it may t may also play an important role in the development of cancer. Ghrelin is a 28 amino acid peptide hormone that exists in two forms. Acyl ghrelin (usually referred to as ghrelin), has a unique n-octanoic acid post-translational modification (which is catalysed by ghrelin O-acyltransferase, GOAT), and desacyl ghrelin, which is a non-octanoylated form. Octanoylated ghrelin acts through the growth hormone secretagogue receptor type 1a (GHSR1a). GHSR1b, an alternatively spliced isoform of GHSR, is C-terminally truncated and does not bind ghrelin. Ghrelin has been implicated in the pathophysiology of a number of diseases Obestatin is a 23 amino acid, C-terminally amidated peptide which is derived from preproghrelin. Although GPR39 was originally thought to be the obestatin receptor this has been disproven, and its receptor remains unknown. Obestatin may have as diverse range of roles as ghrelin. Obestatin improves memory, inhibits thirst and anxiety, increases pancreatic juice secretion and has cardioprotective effects. Obestatin also has been shown to regulate cell proliferation, differentiation and apoptosis in some cell types. Prior to this study, little was known regarding the functions and mechanisms of action ghrelin and obestatin in ovarian cancer. In this study it was demonstrated that the full length ghrelin, GHSR1b and GOAT mRNA transcripts were expressed in all of the ovarian-derived cell lines examined (SKOV3, OV-MZ-6 and hOSE 17.1), however, these cell lines did not express GHSR1a. Ovarian cancer tissue of varying stages and normal ovarian tissue expressed the coding region for ghrelin, obestatin, and GOAT, but not GHSR1a, or GHSR1b. No correlations between cancer grade and the level of expression of these transcripts were observed. This study demonstrated for the first time that both ghrelin and obestatin increase cell migration in ovarian cancer cell lines. Treatment with ghrelin (for 72 hours) significantly increased cell migration in the SKOV3 and OV-MZ-6 ovarian cancer cell lines. Ghrelin (100 nM) stimulated cell migration in the SKOV3 (2.64 +/- 1.08 fold, p <0.05) and OV-MZ-6 (1.65 +/- 0.31 fold, p <0.05) ovarian cancer cell lines, but not in the representative normal cell line hOSE 17.1. This increase in migration was not accompanied by an increase in cell invasion through Matrigel. In contrast to other cancer types, ghrelin had no effect on proliferation. Ghrelin treatment (10nM) significantly decreased attachment of the SKOV3 ovarian cancer cell line to collagen IV (24.7 +/- 10.0 %, p <0.05), however, there were no changes in attachment to the other extracellular matrix molecules (ECM) tested (fibronectin, vitronectin and collagen I), and there were no changes in attachment to any of the ECM molecules in the OV-MZ-6 or hOSE 17.1 cell lines. It is, therefore, unclear if ghrelin plays a role in cell attachment in ovarian cancer. As ghrelin has previously been demonstrated to signal through the ERK1/2 pathway in cancer, we investigated ERK1/2 signalling in ovarian cancer cell lines. In the SKOV3 ovarian cancer cell line, a reduction in ERK1/2 phosphorylation (0.58 fold +/- 0.23, p <0.05) in response to 100 nM ghrelin treatment was observed, while no significant change in ERK1/2 signalling was seen in the OV-MZ-6 cell line with treatment. This suggests that this pathway is unlikely to be involved in mediating the increased migration seen in the ovarian cancer cell lines with ghrelin treatment. In this study ovarian cancer tissue of varying stages and normal ovarian tissue expressed the coding region for obestatin, however, no correlation between cancer grade and level of obestatin transcript expression was observed. In the ovarian-derived cell lines studied (SKOV3, OV-MZ-6 and hOSE 17.1) it was demonstrated that the full length preproghrelin mRNA transcripts were expressed in all cell lines, suggesting they have the ability to produce mature obestatin. This is the first study to demonstrate that obestatin stimulates cell migration and cell invasion. Obestatin induced a significant increase in migration in the SKOV3 ovarian cancer cell line with 10 nM (2.80 +/- 0.52 fold, p <0.05) and 100 nM treatments (3.12 +/- 0.68 fold, p <0.05) and in the OV-MZ-6 cancer cell line with 10 nM (2.04 +/- 0.10 fold, p <0.01) and 100 nM treatments (2.00 +/- 0.37 fold, p <0.05). Obestatin treatment did no affect cell migration in the hOSE 17.1normal ovarian epithelial cell line. Obestatin treatment (100 nM) also stimulated a significant increase in cell invasion in the OV-MZ-6 ovarian cancer cell line (1.45 fold +/- 0.13, p <0.05) and in the hOSE17.1 normal ovarian cell line cells (1.40 fold +/- 0.04 and 1.55 fold +/- 0.05 respectively, p <0.01) with 10 nM and 100 nM treatments. Obestatin treatment did not stimulate cell invasion in the SKOV3 ovarian cancer cell line. This lack of obestatin-stimulated invasion in the SKOV3 cell line may be a cell line specific result. In this study, obestatin did not stimulate cell proliferation in the ovarian cell lines and it has previously been shown to have no effect on cell proliferation in the BON-1 pancreatic neuroendocrine and GC rat somatotroph tumour cell lines. In contrast, obestatin has been shown to affect cell proliferation in gastric and thyroid cancer cell lines, and in some normal cell lines. Obestatin also had no effect on attachment of any of the cell lines to any of the ECM components tested (fibronectin, vitronectin, collagen I and collagen IV). The mechanism of action of obestatin was investigated further using a two dimensional-difference in gel electrophoresis (2D-DIGE) proteomic approach. After treatment with obestating (0, 10 and 100 nM), SKOV3 ovarian cancer and hOSE 17.1 normal ovarian cell lines were collected and 2D-DIGE analysis and mass spectrometry were performed to identify proteins that were differentially expressed in response to treatment. Twenty-six differentially expressed proteins were identified and analysed using Ingenuity Pathway Analysis (IPA). This linked 16 of these proteins in a network. The analysis suggested that the ERK1/2 MAPK pathway was a major mediator of obestatin action. ERK1/2 has previously been shown to be associated with obestatin-stimulated cell proliferation and with the anti-apoptotic effects of obestatin. Activation of the ERK1/2 signalling pathway by obestatin was, therefore, investigated in the SKOV3 and OV-MZ-6 ovarian cancer cell lines using anti-active antibodies and Western immunoblots. Obestatin treatment significantly decreased ERK1/2 phosphorylation at higher obestatin concentrations in both the SKOV3 (100 nM and 1000 nM) and OV-MZ-6 (1000 nM) cell lines compared to the untreated controls. Currently, very little is known about obestatin signalling in cancer. This thesis has demonstrated for the first time that the ghrelin axis may play a role in ovarian cancer migration. Ghrelin and obestatin increased cell migration in ovarian cancer cell lines, indicating that they may be a useful target for therapies that reduce ovarian cancer progression. Further studies investigating the role of the ghrelin axis using in vivo ovarian cancer metastasis models are warranted.

Visualising experimental flow fields through a stormwater gross pollutant trap

Abstract An experimental dataset representing a typical flow field in a stormwater gross pollutant trap (GPT) was visualised. A technique was developed to apply the image-based flow visualisation (IBFV) algorithm to the raw dataset. Particle image velocimetry (PIV) software was previously used to capture the flow field data by tracking neutrally buoyant particles with a high speed camera. The dataset consisted of scattered 2D point velocity vectors and the IBFV visualisation facilitates flow feature characterisation within the GPT. The flow features played a pivotal role in understanding stormwater pollutant capture and retention behaviour within the GPT. It was found that the IBFV animations revealed otherwise unnoticed flow features and experimental artefacts. For example, a circular tracer marker in the IBFV program visually highlighted streamlines to investigate the possible flow paths of pollutants entering the GPT. The investigated flow paths were compared with the behaviour of pollutants monitored during experiments.

Real - time prediction of rainfall induced instability of residual soil slopes associated with deep cracks

The early warning based on real-time prediction of rain-induced instability of natural residual slopes helps to minimise human casualties due to such slope failures. Slope instability prediction is complicated, as it is influenced by many factors, including soil properties, soil behaviour, slope geometry, and the location and size of deep cracks in the slope. These deep cracks can facilitate rainwater infiltration into the deep soil layers and reduce the unsaturated shear strength of residual soil. Subsequently, it can form a slip surface, triggering a landslide even in partially saturated soil slopes. Although past research has shown the effects of surface-cracks on soil stability, research examining the influence of deep-cracks on soil stability is very limited. This study aimed to develop methodologies for predicting the real-time rain-induced instability of natural residual soil slopes with deep cracks. The results can be used to warn against potential rain-induced slope failures. The literature review conducted on rain induced slope instability of unsaturated residual soil associated with soil crack, reveals that only limited studies have been done in the following areas related to this topic: - Methods for detecting deep cracks in residual soil slopes. - Practical application of unsaturated soil theory in slope stability analysis. - Mechanistic methods for real-time prediction of rain induced residual soil slope instability in critical slopes with deep cracks. Two natural residual soil slopes at Jombok Village, Ngantang City, Indonesia, which are located near a residential area, were investigated to obtain the parameters required for the stability analysis of the slope. A survey first identified all related field geometrical information including slope, roads, rivers, buildings, and boundaries of the slope. Second, the electrical resistivity tomography (ERT) method was used on the slope to identify the location and geometrical characteristics of deep cracks. The two ERT array models employed in this research are: Dipole-dipole and Azimuthal. Next, bore-hole tests were conducted at different locations in the slope to identify soil layers and to collect undisturbed soil samples for laboratory measurement of the soil parameters required for the stability analysis. At the same bore hole locations, Standard Penetration Test (SPT) was undertaken. Undisturbed soil samples taken from the bore-holes were tested in a laboratory to determine the variation of the following soil properties with the depth: - Classification and physical properties such as grain size distribution, atterberg limits, water content, dry density and specific gravity. - Saturated and unsaturated shear strength properties using direct shear apparatus. - Soil water characteristic curves (SWCC) using filter paper method. - Saturated hydraulic conductivity. The following three methods were used to detect and simulate the location and orientation of cracks in the investigated slope: (1) The electrical resistivity distribution of sub-soil obtained from ERT. (2) The profile of classification and physical properties of the soil, based on laboratory testing of soil samples collected from bore-holes and visual observations of the cracks on the slope surface. (3) The results of stress distribution obtained from 2D dynamic analysis of the slope using QUAKE/W software, together with the laboratory measured soil parameters and earthquake records of the area. It was assumed that the deep crack in the slope under investigation was generated by earthquakes. A good agreement was obtained when comparing the location and the orientation of the cracks detected by Method-1 and Method-2. However, the simulated cracks in Method-3 were not in good agreement with the output of Method-1 and Method-2. This may have been due to the material properties used and the assumptions made, for the analysis. From Method-1 and Method-2, it can be concluded that the ERT method can be used to detect the location and orientation of a crack in a soil slope, when the ERT is conducted in very dry or very wet soil conditions. In this study, the cracks detected by the ERT were used for stability analysis of the slope. The stability of the slope was determined using the factor of safety (FOS) of a critical slip surface obtained by SLOPE/W using the limit equilibrium method. Pore-water pressure values for the stability analysis were obtained by coupling the transient seepage analysis of the slope using finite element based software, called SEEP/W. A parametric study conducted on the stability of an investigated slope revealed that the existence of deep cracks and their location in the soil slope are critical for its stability. The following two steps are proposed to predict the rain-induced instability of a residual soil slope with cracks. (a) Step-1: The transient stability analysis of the slope is conducted from the date of the investigation (initial conditions are based on the investigation) to the preferred date (current date), using measured rainfall data. Then, the stability analyses are continued for the next 12 months using the predicted annual rainfall that will be based on the previous five years rainfall data for the area. (b) Step-2: The stability of the slope is calculated in real-time using real-time measured rainfall. In this calculation, rainfall is predicted for the next hour or 24 hours and the stability of the slope is calculated one hour or 24 hours in advance using real time rainfall data. If Step-1 analysis shows critical stability for the forthcoming year, it is recommended that Step-2 be used for more accurate warning against the future failure of the slope. In this research, the results of the application of the Step-1 on an investigated slope (Slope-1) showed that its stability was not approaching a critical value for year 2012 (until 31st December 2012) and therefore, the application of Step-2 was not necessary for the year 2012. A case study (Slope-2) was used to verify the applicability of the complete proposed predictive method. A landslide event at Slope-2 occurred on 31st October 2010. The transient seepage and stability analyses of the slope using data obtained from field tests such as Bore-hole, SPT, ERT and Laboratory tests, were conducted on 12th June 2010 following the Step-1 and found that the slope in critical condition on that current date. It was then showing that the application of the Step-2 could have predicted this failure by giving sufficient warning time.

Three-dimensional modeling of steel portal frame buildings

The realistic strength and deflection behavior of industrial and commercial steel portal frame buildings are understood only if the effects of rigidity of end frames and profiled steel claddings are included. The conventional designs ignore these effects and are very much based on idealized two-dimensional (2D) frame behavior. Full-scale tests of a 1212 m steel portal frame building under a range of design load cases indicated that the observed deflections and bending moments in the portal frame were considerably different from those obtained from a 2D analysis of frames ignoring these effects. Three-dimensional (3D) analyses of the same building, including the effects of end frames and cladding, were carried out, and the results agreed well with full-scale test results. Results clearly indicated the need for such an analysis and for testing to study the true behavior of steel portal frame buildings. It is expected that such a 3D analysis will lead to lighter steel frames as the maximum moments and deflections are reduced.

The effect of organ motion of radiotherapy plans

Organ motion as a result of respiration is an important field of research for medical physics. Knowledge of magnitude and direction of this motion is necessary to allow for more accurate radiotherapy treatment planning. This will result in higher doses to the tumour whilst sparing healthy tissue. This project involved human trials, where the radiation therapy patient's kidneys were CT scanned under three different conditions; whilst free breathing (FB), breath-hold at normal tidal inspiration (BHIN), and breath-hold at normal tidal expiration (BHEX). The magnitude of motion was measured by recording the outline of the kidney from a Beam's Eye View (BEV). The centre of mass of this 2D shape was calculated for each set using "ImageJ" software and the magnitude of movement determined from the change in the centroid's coordinates between the BHIN and BHEX scans. The movement ranged from, for the left and right kidneys, 4-46mm and 2-44mm in the superior/inferior (axial) plane, 1-21mm and 2- 16mm in the anterior/posterior (coronal) plane, and 0-6mm and 0-8mm in the lateral/medial (sagittal) plane. From exhale to inhale, the kidneys tended to move inferiorly, anteriorly and laterally. A standard radiotherapy plan, designed to treat the para-aortics with opposed lateral fields was performed on the free breathing (planning) CT set. The field size and arrangement was set up using the same parameters for each subject. The prescription was to deliver 45 Gray in 25 fractions. This field arrangement and prescription was then copied over to the breath hold CT sets, and the dosimetric differences were compared using Dose Volume Histograms (DVH). The point of comparison for the three sets was recorded as the percentage volume of kidney receiving less than or equal to 10 Gray. The QUASAR respiratory motion phantom was used with the range of motion determined from the human study. The phantom was imaged, planned and treated with a linear accelerator with dose determined by film. The effect of the motion was measured by the change in the penumbra of the film and compared to the penumbra from the treatment planning system.

Taguchi optimized synthesis of graphene films by copper catalyzed ethanol decomposition

Taguchi method is for the first time applied to optimize the synthesis of graphene films by copper-catalyzed decomposition of ethanol. In order to find the most appropriate experimental conditions for the realization of thin high-grade films, six experiments suitably designed and performed. The influence of temperature (1000–1070 °C) and synthesis duration (1–30 min) and hydrogen flow (0–100 sccm) on the number of graphene layers and defect density in the graphitic lattice was ranked by monitoring the intensity of the 2D- and D-bands relative to the G-band in the Raman spectra. After critical examination and adjusting of the conditions predicted to give optimal results, a continuous film consisting of 2–4 nearly defect-free graphene layers was obtained.

Enabling aircraft emergency landings using active visual site detection

The ability to automate forced landings in an emergency such as engine failure is an essential ability to improve the safety of Unmanned Aerial Vehicles operating in General Aviation airspace. By using active vision to detect safe landing zones below the aircraft, the reliability and safety of such systems is vastly improved by gathering up-to-the-minute information about the ground environment. This paper presents the Site Detection System, a methodology utilising a downward facing camera to analyse the ground environment in both 2D and 3D, detect safe landing sites and characterise them according to size, shape, slope and nearby obstacles. A methodology is presented showing the fusion of landing site detection from 2D imagery with a coarse Digital Elevation Map and dense 3D reconstructions using INS-aided Structure-from-Motion to improve accuracy. Results are presented from an experimental flight showing the precision/recall of landing sites in comparison to a hand-classified ground truth, and improved performance with the integration of 3D analysis from visual Structure-from-Motion.

An evaluation of different features and learning models for anomalous event detection

The huge amount of CCTV footage available makes it very burdensome to process these videos manually through human operators. This has made automated processing of video footage through computer vision technologies necessary. During the past several years, there has been a large effort to detect abnormal activities through computer vision techniques. Typically, the problem is formulated as a novelty detection task where the system is trained on normal data and is required to detect events which do not fit the learned ‘normal’ model. There is no precise and exact definition for an abnormal activity; it is dependent on the context of the scene. Hence there is a requirement for different feature sets to detect different kinds of abnormal activities. In this work we evaluate the performance of different state of the art features to detect the presence of the abnormal objects in the scene. These include optical flow vectors to detect motion related anomalies, textures of optical flow and image textures to detect the presence of abnormal objects. These extracted features in different combinations are modeled using different state of the art models such as Gaussian mixture model(GMM) and Semi- 2D Hidden Markov model(HMM) to analyse the performances. Further we apply perspective normalization to the extracted features to compensate for perspective distortion due to the distance between the camera and objects of consideration. The proposed approach is evaluated using the publicly available UCSD datasets and we demonstrate improved performance compared to other state of the art methods.

Numerical simulations of particle deposition in metal foam heat exchangers

Australia is a high-potential country for geothermal power with reserves currently estimated in the tens of millions of petajoules, enough to power the nation for at least 1000 years at current usage. However, these resources are mainly located in isolated arid regions where water is scarce. Therefore, wet cooling systems for geothermal plants in Australia are the least attractive solution and thus air-cooled heat exchangers are preferred. In order to increase the efficiency of such heat exchangers, metal foams have been used. One issue raised by this solution is the fouling caused by dust deposition. In this case, the heat transfer characteristics of the metal foam heat exchanger can dramatically deteriorate. Exploring the particle deposition property in the metal foam exchanger becomes crucial. This paper is a numerical investigation aimed to address this issue. Two dimensional (2D) numerical simulations of a standard one-row tube bundle wrapped with metal foam in cross-flow are performed and highlight preferential particle deposition areas.

Numerical simulations of particle deposition in metal foam heat exchangers

Australia is a high potential country for geothermal power with reserves currently estimated in the tens of millions of petajoules, enough to power the nation for at least 1000 years at current usage.However, these resources are mainly located in isolated arid regions where water is scarce. Therefore, wet cooling systems for geothermal plants in Australia are the least attractive solution and thus air-cooled heat exchangers are preferred. In order to increase the efficiency of such heat exchangers, metal foams have been used. One issue raised by this solution is the fouling caused by dust deposition. In this case, the heat transfer characteristics of the metal foam heat exchanger can dramatically deteriorate. Exploring the particle deposition property in the metal foam exchanger becomes crucial. This paper is a numerical investigation aimed to address this issue. Two-dimensional(2D numerical simulations of a standard one-row tube bundle wrapped with metal foam in cross-flow are performed and highlight preferential particle deposition areas.