-
Assessment of flamelet manifolds for turbulent flame-wall interactions in Large-Eddy Simulations
Authors:
Yujuan Luo,
Matthias Steinhausen,
Driss Kaddar,
Christian Hasse,
Federica Ferraro
Abstract:
A turbulent side-wall quenching (SWQ) flame in a fully developed channel flow is studied using Large-Eddy Simulation (LES) with a tabulated chemistry approach. Three different flamelet manifolds with increasing levels of complexity are applied: the Flamelet-Generated Manifold (FGM) considering varying enthalpy levels, the Quenching Flamelet-Generated Manifold (QFM), and the recently proposed Quenc…
▽ More
A turbulent side-wall quenching (SWQ) flame in a fully developed channel flow is studied using Large-Eddy Simulation (LES) with a tabulated chemistry approach. Three different flamelet manifolds with increasing levels of complexity are applied: the Flamelet-Generated Manifold (FGM) considering varying enthalpy levels, the Quenching Flamelet-Generated Manifold (QFM), and the recently proposed Quenching Flamelet-Generated Manifold with Exhaust Gas Recirculation (QFM-EGR), with the purpose being to assess their capability to predict turbulent flame-wall interactions (FWIs), which are highly relevant to numerical simulations of real devices such as gas turbines and internal combustion engines.
The accuracy of the three manifolds is evaluated and compared a posteriori, using the data from a previously published flame-resolved simulation with detailed chemistry for reference. For LES with the FGM, the main characteristics such as the mean flow field, temperature, and major species can be captured well, while notable deviations from the reference results are observed for the near-wall region, especially for pollutant species such as \ce{CO}. In accordance with the findings from laminar FWI, improvement is also observed in the simulation with QFM under turbulent flow conditions. Although LES with the QFM-EGR shows a similar performance in the prediction of mean quantities as LES with QFM, it presents significantly better agreement with the reference data regarding instantaneous thermo-chemical states near the quenching point. This indicates the necessity to take into account the mixing effects in the flamelet manifold to correctly capture the flame-vortex interaction near the flame tip in turbulent configurations. Based on the findings from this study, suitable flamelet manifolds can be chosen depending on the aspects of interest in practical applications.
△ Less
Submitted 15 August, 2023;
originally announced August 2023.
-
Application of dense neural networks for manifold-based modeling of flame-wall interactions
Authors:
Julian Bissantz,
Jeremy Karpowski,
Matthias Steinhausen,
Yujuan Luo,
Federica Ferraro,
Arne Scholtissek,
Christian Hasse,
Luc Vervisch
Abstract:
Artifical neural networks (ANNs) are universal approximators capable of learning any correlation between arbitrary input data with corresponding outputs, which can also be exploited to represent a low-dimensional chemistry manifold in the field of combustion. In this work, a procedure is developed to simulate a premixed methane-air flame undergoing side-wall quenching utilizing an ANN chemistry ma…
▽ More
Artifical neural networks (ANNs) are universal approximators capable of learning any correlation between arbitrary input data with corresponding outputs, which can also be exploited to represent a low-dimensional chemistry manifold in the field of combustion. In this work, a procedure is developed to simulate a premixed methane-air flame undergoing side-wall quenching utilizing an ANN chemistry manifold. In the investigated case, the flame characteristics are governed by two canonical problems: the adiabatic flame propagation in the core flow and the non-adiabatic flame-wall interaction governed by enthalpy losses to the wall. Similar to the tabulation of a Quenching Flamelet-Generated Manifold (QFM), the neural network is trained on a 1D head-on quenching flame database to learn the intrinsic chemistry manifold. The control parameters (i.e. the inputs) of the ANN are identified from thermo-chemical state variables by a sparse principal component analysis (PCA) without using prior knowledge about the flame physics. These input quantities are then transported in the coupled CFD solver and used for manifold access during simulation runtime. The chemical source terms are corrected at the manifold boundaries to ensure boundedness of the thermo-chemical state at all times. Finally, the ANN model is assessed by comparison to simulation results of the 2D side-wall quenching (SWQ) configuration with detailed chemistry and with a flamelet-based manifold (QFM).
△ Less
Submitted 10 August, 2023;
originally announced August 2023.
-
Numerical Investigation of the Local Thermo-Chemical State in a Thermo-Acoustically Unstable Dual Swirl Gas Turbine Model Combustor
Authors:
T. Jeremy P. Karpowski,
Federica Ferraro,
Matthias Steinhausen,
Sebastian Popp,
Christoph M. Arndt,
Christian Kraus,
Henning Bockhorn,
Wolfgang Meier,
Christian Hasse
Abstract:
In this work, the thermo-acoustic instabilities of a gas turbine model combustor, the so-called SFB606 combustor, are numerically investigated using Large Eddy Simulation (LES) combined with tabulated chemistry and Artificial Thickened Flame (ATF) approach. The main focus is a detailed analysis of the thermo-acoustic cycle and the accompanied equivalence ratio oscillations and their associated con…
▽ More
In this work, the thermo-acoustic instabilities of a gas turbine model combustor, the so-called SFB606 combustor, are numerically investigated using Large Eddy Simulation (LES) combined with tabulated chemistry and Artificial Thickened Flame (ATF) approach. The main focus is a detailed analysis of the thermo-acoustic cycle and the accompanied equivalence ratio oscillations and their associated convective time delay. In particular, the variations of the thermo-chemical state and flame characteristics over the thermo-acoustic cycle are investigated. For the operating point flame B ($P_{th}=25\,$kW), the burner exhibits thermo-acoustic instabilities with a dominant frequency of 392Hz, the acoustic eigenmode of the inner air inlet duct. These oscillations are accompanied by an equivalence ratio oscillation, which exhibits a convective time delay between the injection in the inner swirler and the flame zone. Two LES, one adiabatic and one accounting for heat losses at the walls by prescribing the wall temperatures from experimental data and Conjugated Heat Transfer (CHT) simulations, are conducted. Results with the enthalpy-dependent table are found to predict the time-averaged flow field in terms of velocity, major species, and temperature with higher accuracy than in the adiabatic case. Further, they indicate, that heat losses should be accounted for to correctly predict the flame position. Subsequently, the thermo-chemical state variations over the thermo-acoustic cycle for the enthalpy-dependant case are analyzed in detail and compared with experimental data in terms of phase-conditioned averaged profiles and conditional averages. An overall good prediction is observed. The results provide a detailed quantitative analysis of the thermo-acoustic feedback mechanism of this burner.
△ Less
Submitted 10 August, 2023;
originally announced August 2023.
-
Combined effects of heat loss and curvature on turbulent flame-wall interaction in a premixed dimethyl ether/air flame
Authors:
Driss Kaddar,
Matthias Steinhausen,
Thorsten Zirwes,
Henning Bockhorn,
Christian Hasse,
Federica Ferraro
Abstract:
This study investigates the effects of curvature on the local heat release rate and mixture fraction during turbulent flame-wall interaction of a lean dimethyl ether/air flame using a fully resolved simulation with a reduced skeletal chemical reaction mechanism and mixture-averaged transport. The region in which turbulent flame-wall interaction affects the flame is found to be restricted to a wall…
▽ More
This study investigates the effects of curvature on the local heat release rate and mixture fraction during turbulent flame-wall interaction of a lean dimethyl ether/air flame using a fully resolved simulation with a reduced skeletal chemical reaction mechanism and mixture-averaged transport. The region in which turbulent flame-wall interaction affects the flame is found to be restricted to a wall distance less than twice the laminar flame thickness. In regions without heat losses, heat release rate and curvature, as well as mixture fraction and curvature, are negatively correlated, which is in accordance with experimental findings. Flame-wall interaction alters the correlation between heat release rate and curvature. An inversion in the sign of the correlation from negative to positive is observed as the flame starts to experience heat losses to the wall. The correlation between mixture fraction and curvature, however, is unaffected by flame-wall interactions and remains negative. Similarly to experimental findings, the investigated turbulent side-wall quenching flame shows both head-on quenching and side-wall quenching-like behavior. The different quenching events are associated with different curvature values in the near-wall region. Furthermore, for medium heat loss, the correlations between heat release rate and curvature are sensitive to the quenching scenario.
△ Less
Submitted 24 February, 2023;
originally announced February 2023.
-
Flamelet modeling of thermo-diffusively unstable hydrogen-air flames
Authors:
Hannes Böttler,
Haris Lulic,
Matthias Steinhausen,
Xu Wen,
Christian Hasse,
Arne Scholtissek
Abstract:
In order to reduce CO2 emissions, hydrogen combustion has become increasingly relevant for technical applications. In this context, lean H2-air flames show promising features but, among other characteristics, they tend to exhibit thermo-diffusive instabilities. The formation of cellular structures associated with these instabilities leads to an increased flame surface area which further promotes t…
▽ More
In order to reduce CO2 emissions, hydrogen combustion has become increasingly relevant for technical applications. In this context, lean H2-air flames show promising features but, among other characteristics, they tend to exhibit thermo-diffusive instabilities. The formation of cellular structures associated with these instabilities leads to an increased flame surface area which further promotes the flame propagation speed, an important reference quantity for design, control, and safe operation of technical combustors. While many studies have addressed the physical phenomena of intrinsic flame instabilities in the past, there is also a demand to predict such flame characteristics with reduced-order models to allow computationally efficient simulations. In this work, a H2-air spherical expanding flame, which exhibits thermo-diffusive instabilities, is studied with flamelet-based modeling approaches both in a-priori and a-posteriori manner. A recently proposed Flamelet/Progress Variable (FPV) model, with a manifold based on unstretched planar flames, and a novel FPV approach, which takes into account a large curvature variation in the tabulated manifold, are compared to detailed chemistry (DC) calculations. First, both FPV approaches are assessed in terms of an a-priori test with the DC reference dataset. Thereafter, the a-posteriori assessment contains two parts: a linear stability analysis of perturbed planar flames and the simulation of the spherical expanding flame. Both FPV models are systematically analyzed considering global and local flame properties in comparison to the DC reference data. It is shown that the new FPV model, incorporating large curvature variations in the manifold, leads to improved predictions for the microstructure of the corrugated flame front and the formation of cellular structures, while global flame properties are reasonably well reproduced by both models.
△ Less
Submitted 17 January, 2023;
originally announced January 2023.
-
Turbulent flame-wall interaction of premixed flames using Quadrature-based Moment Methods (QbMM) and tabulated chemistry: An a priori analysis
Authors:
M. Steinhausen,
T. Zirwes,
F. Ferraro,
S. Popp,
F. Zhang,
H. Bockhorn,
C. Hasse
Abstract:
Presumed probability density function (PDF) and transported PDF methods are commonly applied to model the turbulence chemistry interaction in turbulent reacting flows. However, little focus has been given to the turbulence chemistry interaction PDF closure for flame-wall interaction. In this study, a quasi-DNS of a turbulent premixed, stoichiometric methane-air flame ignited in a fully developed t…
▽ More
Presumed probability density function (PDF) and transported PDF methods are commonly applied to model the turbulence chemistry interaction in turbulent reacting flows. However, little focus has been given to the turbulence chemistry interaction PDF closure for flame-wall interaction. In this study, a quasi-DNS of a turbulent premixed, stoichiometric methane-air flame ignited in a fully developed turbulent channel flow undergoing side-wall quenching is investigated. The objective of this study is twofold. First, the joint PDF of the progress variable and enthalpy that needs to be accounted for in turbulence chemistry interaction closure models is analyzed in the quasi-DNS configuration, both in the core flow and the near-wall region. Secondly, a transported PDF closure model, based on a Conditional Quadrature Method of Moments approach, and a presumed PDF approach are examined in an a priori analysis using the quasi-DNS as a reference both in the context of Reynolds-Averaged Navier Stokes (RANS) and Large-Eddy Simulations (LESs). The analysis of the joint PDF demonstrates the high complexity of the reactive scalar distribution in the near-wall region. Here, a high correlation between the progress variable and enthalpy is found, where the flame propagation and quenching are present simultaneously. The transported PDF approach presented in this work, based on the Conditional Quadrature Method of Moments, accounts for the moments of the joint PDF of progress variable and enthalpy coupled to a Quenching-Flamelet Generated Manifold. In the a priori analysis both turbulence chemistry interaction PDF closure models show a high accuracy in the core flow. In the near-wall region, however, only the Conditional Quadrature Method of Moments approach is suitable to predict the flame structures.
△ Less
Submitted 7 December, 2022;
originally announced December 2022.
-
Effect of flame retardants on side-wall quenching of partially premixed laminar flames
Authors:
Matthias Steinhausen,
Federica Ferraro,
Max Schneider,
Florian Zentgraf,
Max Greifenstein,
Andreas Dreizler,
Christian Hasse,
Arne Scholtissek
Abstract:
A combined experimental and numerical investigation of partially premixed laminar methane-air flames undergoing side-wall quenching (SWQ) is performed. A well-established SWQ burner is adapted to allow the seeding of the main flow with additional gaseous products issued from a (secondary) wall inlet close to the flame's quenching point. First, the characteristics of the partially premixed flame th…
▽ More
A combined experimental and numerical investigation of partially premixed laminar methane-air flames undergoing side-wall quenching (SWQ) is performed. A well-established SWQ burner is adapted to allow the seeding of the main flow with additional gaseous products issued from a (secondary) wall inlet close to the flame's quenching point. First, the characteristics of the partially premixed flame that quenches at the wall are assessed using planar laser-induced fluorescence measurements of the OH radical, and a corresponding numerical simulation with fully-resolved transport and chemistry is conducted. A boundary layer of enriched mixture is formed at the wall, leading to a reaction zone parallel to the wall for high injection rates from the wall inlet. Subsequently, in a numerical study, the wall inflow is mixed with dimethylmethylphosphonat (DMMP), a phosphor-based flame retardant. The DMMP addition allows the assessment of the combined effects of heat loss and flame retardants on the flame structure during flame-wall interaction. With an increasing amount of DMMP in the injected mixture, the flame stabilizes further away from the wall and shows a decrease in the local heat-release rate. Thereby, the maximum wall heat flux is significantly reduced. That results in a lower thermal load on the quenching wall. The flame structure analysis shows an accumulation of the intermediate species HOPO at the wall similar to the CO accumulation during the quenching of premixed flames without flame retardant addition. The study shows how the structure of a partially premixed flame is influenced by a wall that releases either additional fuel or a mixture of fuel and flame retardant. The insights gained from the canonical configuration can lead to a better understanding of the combined effects of flame retardants and heat losses in near-wall flames.
△ Less
Submitted 7 December, 2022;
originally announced December 2022.
-
Flame-vortex interaction during turbulent side-wall quenching and its implications for flamelet manifolds
Authors:
Matthias Steinhausen,
Thorsten Zirwes,
Federica Ferraro,
Arne Scholtissek,
Henning Bockhorn,
Christian Hasse
Abstract:
In this study, the thermochemical state during turbulent flame-wall interaction of a stoichiometric methane-air flame is investigated using a fully resolved simulation with detailed chemistry. The turbulent side-wall quenching flame shows both head-on quenching and side-wall quenching-like behavior that significantly affects the CO formation in the near-wall region. The detailed insights from the…
▽ More
In this study, the thermochemical state during turbulent flame-wall interaction of a stoichiometric methane-air flame is investigated using a fully resolved simulation with detailed chemistry. The turbulent side-wall quenching flame shows both head-on quenching and side-wall quenching-like behavior that significantly affects the CO formation in the near-wall region. The detailed insights from the simulation are used to evaluate a recently proposed flame (tip) vortex interaction mechanism identified from experiments on turbulent side-wall quenching. It describes the entrainment of burnt gases into the fresh gas mixture near the flame's quenching point. The flame behavior and thermochemical states observed in the simulation are similar to the phenomena observed in the experiments. A novel chemistry manifold is presented that accounts for both the effects of flame dilution due to exhaust gas recirculation in the flame vortex interaction area and enthalpy losses to the wall. The manifold is validated in an a-priori analysis using the simulation results as a reference. The incorporation of exhaust gas recirculation effects in the manifold leads to a significantly increased prediction accuracy in the near-wall regions of flame-vortex interactions.
△ Less
Submitted 7 December, 2022;
originally announced December 2022.
-
Short dissipation times of proto-planetary discs - an artifact of selection effects?
Authors:
Susanne Pfalzner,
Manuel Steinhausen,
Karl Menten
Abstract:
The frequency of discs around young stars, a key parameter for understanding planet formation, is most readily determined in young stellar clusters where many relatively coeval stars are located in close proximity. Observational studies seem to show that the disc frequency decreases rapidly with cluster age with <10% of cluster stars retaining their discs for longer than 2-6 Myr. Given that at lea…
▽ More
The frequency of discs around young stars, a key parameter for understanding planet formation, is most readily determined in young stellar clusters where many relatively coeval stars are located in close proximity. Observational studies seem to show that the disc frequency decreases rapidly with cluster age with <10% of cluster stars retaining their discs for longer than 2-6 Myr. Given that at least half of all stars in the field seem to harbor one or more planets, this would imply extremely fast disc dispersal and rapid planet growth. Here we question the validity of this constraint by demonstrating that the short disc dissipation times inferred to date might have been heavily underestimated by selection effects. Critically, for ages >3Myr only stars that originally populated the densest areas of very populous clusters, which are prone to disc erosion, are actually considered. This tiny sample may not be representative of the majority of stars. In fact, the higher disc fractions in co-moving groups indicate that it is likely that over 30% of all field stars retain their discs well beyond 10 Myr, leaving ample time for planet growth. Equally our solar system, with a likely formation time > 10 Myr, need no longer be an exception but in fact typical of planetary systems.
△ Less
Submitted 3 September, 2014;
originally announced September 2014.
-
The evolutionary tracks of young massive star clusters
Authors:
S. Pfalzner,
G. Parmentier,
M. Steinhausen,
K. Vincke,
K. Menten
Abstract:
Stars mostly form in groups consisting of a few dozen to several ten thousand members. For 30 years, theoretical models provide a basic concept of how such star clusters form and develop: they originate from the gas and dust of collapsing molecular clouds. The conversion from gas to stars being incomplete, the left over gas is expelled, leading to cluster expansion and stars becoming unbound. Obse…
▽ More
Stars mostly form in groups consisting of a few dozen to several ten thousand members. For 30 years, theoretical models provide a basic concept of how such star clusters form and develop: they originate from the gas and dust of collapsing molecular clouds. The conversion from gas to stars being incomplete, the left over gas is expelled, leading to cluster expansion and stars becoming unbound. Observationally, a direct confirmation of this process has proved elusive, which is attributed to the diversity of the properties of forming clusters. Here we take into account that the true cluster masses and sizes are masked, initially by the surface density of the background and later by the still present unbound stars. Based on the recent observational finding that in a given star-forming region the star formation efficiency depends on the local density of the gas, we use an analytical approach combined with \mbox{N-body simulations, to reveal} evolutionary tracks for young massive clusters covering the first 10 Myr. Just like the Hertzsprung-Russell diagram is a measure for the evolution of stars, these tracks provide equivalent information for clusters. Like stars, massive clusters form and develop faster than their lower-mass counterparts, explaining why so few massive cluster progenitors are found.
△ Less
Submitted 26 August, 2014;
originally announced August 2014.
-
Sizes of protoplanetary discs after star-disc encounters
Authors:
Andreas Breslau,
Manuel Steinhausen,
Kirsten Vincke,
Susanne Pfalzner
Abstract:
Most stars do not form in isolation, but as part of a star cluster or association. These young stars are initially surrounded by protoplanetary discs. In these cluster environments tidal interactions with other cluster members can alter the disc properties. Besides the disc frequency, its mass, angular momentum, and energy, in particular the disc's size is prone to being changed by a passing star.…
▽ More
Most stars do not form in isolation, but as part of a star cluster or association. These young stars are initially surrounded by protoplanetary discs. In these cluster environments tidal interactions with other cluster members can alter the disc properties. Besides the disc frequency, its mass, angular momentum, and energy, in particular the disc's size is prone to being changed by a passing star. So far the change in disc size was only investigated for a small number of very specific encounters. Several studies investigated the effect of the cluster environment on the sizes of planetary systems, like our own solar system, based on a generalisation of information from this limited sample. We performed numerical simulations covering the wide parameter space typical for young star clusters, to test the validity of this approach. Here the sizes of discs after encounters are presented, based on a size definition which is comparable to that one used in observational studies. We find that, except for encounters between equal-mass stars, the usually applied estimates are insufficient. They tend to severely overestimate the remaining disc size. We show that the disc size after an encounter can be described by a relatively simple dependence on the periastron distance and the mass ratio of the encounter partners. This knowledge allows, for example, to pin down the types of encounter possibly responsible for the structure of today's solar system.
△ Less
Submitted 31 March, 2014;
originally announced March 2014.
-
Does the mass distribution in discs influence encounter-induced losses in young star clusters?
Authors:
Manuel Steinhausen,
Susanne Pfalzner
Abstract:
One mechanism for the external destruction of protoplanetary discs in young dense clusters is tidal disruption during the flyby of another cluster member. The degree of mass loss in such an encounter depends, among other parameters, on the distribution of the material within the disc. Previous work showed that this is especially so in encounters that truncate large parts of the outer disc. The exp…
▽ More
One mechanism for the external destruction of protoplanetary discs in young dense clusters is tidal disruption during the flyby of another cluster member. The degree of mass loss in such an encounter depends, among other parameters, on the distribution of the material within the disc. Previous work showed that this is especially so in encounters that truncate large parts of the outer disc. The expectation is that the number of completely destroyed discs in a cluster depends also on the mass distribution within the discs. Here we test this hypothesis by determining the influence of encounters on the disc fraction and average disc mass in clusters of various stellar densities for different mass distributions in the discs. This is done by performing Nbody6 simulation of a variety of cluster environments, where we track the encounter dynamics and determine the mass loss due to these encounters for different disc-mass distributions. We find that although the disc mass distribution has a significant impact on the disc losses for specific star-disc encounters, the overall disc frequency generally remains rather unaffected. The reason is that in single encounters the dependence on the mass distribution is strongest if both stars have very different masses. Such encounters are rather infrequent in sparse clusters. In dense clusters such encounters are more common, however, here the disc frequency is largely determined by encounters between low-mass stars such that the overall disc frequency does not change significantly. For tidal disruption the disc destruction in clusters is fairly independent of the actual distribution of the material in the disc. The all determining factor remains the cluster density.
△ Less
Submitted 7 March, 2014;
originally announced March 2014.
-
Disc-mass distribution in star-disc encounters
Authors:
M. Steinhausen,
C. Olczak,
S. Pfalzner
Abstract:
Investigations of stellar encounters in cluster environments have demonstrated their potential influence on the mass and angular momentum of protoplanetary discs around young stars. In this study it is investigated in how far the initial surface density in the disc surrounding a young star influences the outcome of an encounter. Based on a power-law ansatz for the surface density,…
▽ More
Investigations of stellar encounters in cluster environments have demonstrated their potential influence on the mass and angular momentum of protoplanetary discs around young stars. In this study it is investigated in how far the initial surface density in the disc surrounding a young star influences the outcome of an encounter. Based on a power-law ansatz for the surface density, $Σ(r) \propto r^{-p}$, a parameter study of star-disc encounters with different initial disc-mass distributions has been performed using N-body simulations. It is demonstrated that the shape of the disc-mass distribution has a significant impact on the quantity of the disc-mass and angular momentum losses in star-disc encounters. Most sensitive are the results where the outer parts of the disc are perturbed by high-mass stars. By contrast, disc-penetrating encounters lead more or less independently of the disc-mass distribution always to large losses. However, maximum losses are generally obtained for initially flat distributed disc material. Based on the parameter study a fit formula is derived, describing the relative mass and angular momentum loss dependent on the initial disc-mass distribution index p. Generally encounters lead to a steepening of the density profile of the disc. The resulting profiles can have a r^{-2}-dependence or even steeper independent of the initial distribution of the disc material. From observations the initial density distribution in discs remains unconstrained, so the here demonstrated strong dependence on the initial density distribution might require a revision of the effect of encounters in young stellar clusters. The steep surface density distributions induced by some encounters might be the prerequisite to form planetary systems similar to our own solar system.
△ Less
Submitted 10 November, 2011;
originally announced November 2011.
-
Accretion bursts in young stars driven by cluster environment
Authors:
S. Pfalzner,
J. Tackenberg,
M. Steinhausen
Abstract:
The standard picture of accretion is a steady flow of matter from the disc onto the young star - a concept which assumes the star-disc system to be completely isolated. However, in a dense cluster environment star-disc systems do interact gravitationally. The aim here is to estimate the encounter-induced accretion rate in an ONC-like environment. Combining simulations of the cluster dynamics wit…
▽ More
The standard picture of accretion is a steady flow of matter from the disc onto the young star - a concept which assumes the star-disc system to be completely isolated. However, in a dense cluster environment star-disc systems do interact gravitationally. The aim here is to estimate the encounter-induced accretion rate in an ONC-like environment. Combining simulations of the cluster dynamics with simulations of the effect of encounters on star-disc systems we determine the likelihood and degree of encounter-triggered accretion processes. We show that accretion bursts triggered by encounters of star-disc systems are common in young dense clusters like the ONC leading in the outburst phase to typical accretion rates of 10^(-7)-10^(-4) M(sun)/yr. Up to a third of stars presently in the Trapezium region accreted at least 1% of their disc mass via this mechanism in the last 1Myr. Accretion of over 6-7% of the disc material can occur in a single encounter. Despite losing their discs quickly, the total percentage of disc matter accreted per star is largest for the massive stars. Supplementing the steady accretion flow there exist episodic periods of high accretion in dense cluster environments. Due to their high accretion rate these processes should be observable even now in some of the low-mass stars in the ONC.
△ Less
Submitted 26 June, 2008;
originally announced June 2008.
