Weak first- and second-order numerical schemes for stochastic differential equations appearing in Lagrangian two-phase flow modeling

Verfasser / Beitragende:
[Jean-Pierre Minier, Eric Peirano, Sergio Chibbaro]
Ort, Verlag, Jahr:
2003
Enthalten in:
Monte Carlo Methods and Applications, 9/2(2003-04-01), 93-133
Format:
Artikel (online)
ID: 378876503
LEADER caa a22 4500
001 378876503
003 CHVBK
005 20180305123421.0
007 cr unu---uuuuu
008 161128e20030401xx s 000 0 eng
024 7 0 |a 10.1515/156939603322663312  |2 doi 
035 |a (NATIONALLICENCE)gruyter-10.1515/156939603322663312 
245 0 0 |a Weak first- and second-order numerical schemes for stochastic differential equations appearing in Lagrangian two-phase flow modeling  |h [Elektronische Daten]  |c [Jean-Pierre Minier, Eric Peirano, Sergio Chibbaro] 
520 3 |a Weak first- and second-order numerical schemes are developed to integrate the stochastic differential equations that arise in mean-field - pdf methods (Lagrangian stochastic approach) for modeling polydispersed turbulent two-phase flows. These equations present several challenges, the foremost being that the problem is characterized by the presence of different time scales that can lead to stiff equations, when the smallest time-scale is significantly less than the time-step of the simulation. The numerical issues have been detailed by Minier [Monte Carlo Meth. and Appl. 7 295-310, (2000)] and the present paper proposes numerical schemes that satisfy these constraints. This point is really crucial for physical and engineering applications, where various limit cases can be present at the same time in different parts of the domain or at different times. In order to build up the algorithm, the analytical solutions to the equations are first carried out when the coefficients are constant. By freezing the coefficients in the analytical solutions, first and second order unconditionally stable weak schemes are developed. A prediction/ correction method, which is shown to be consistent for the present stochastic model, is used to devise the second-order scheme. A complete numerical investigation is carried out to validate the schemes, having included also a comprehensive study of the different error sources. The final method is demonstrated to have the required stability, accuracy and efficiency. 
540 |a Copyright 2003, Walter de Gruyter 
700 1 |a Minier  |D Jean-Pierre  |u Electricité de France, MFTT, 6 Quati Watier 78400 France.  |4 aut 
700 1 |a Peirano  |D Eric  |u ADEME-Der, 500 Route des Lucioles 06560 Valbonne France.  |4 aut 
700 1 |a Chibbaro  |D Sergio  |u Dipartimento di fisica and INFN, Università di Cagliari, Cittadella Universitaria Monserrato CA, Italy.  |4 aut 
773 0 |t Monte Carlo Methods and Applications  |d Walter de Gruyter  |g 9/2(2003-04-01), 93-133  |x 0929-9629  |q 9:2<93  |1 2003  |2 9  |o mcma 
856 4 0 |u https://doi.org/10.1515/156939603322663312  |q text/html  |z Onlinezugriff via DOI 
908 |D 1  |a research article  |2 jats 
950 |B NATIONALLICENCE  |P 856  |E 40  |u https://doi.org/10.1515/156939603322663312  |q text/html  |z Onlinezugriff via DOI 
950 |B NATIONALLICENCE  |P 700  |E 1-  |a Minier  |D Jean-Pierre  |u Electricité de France, MFTT, 6 Quati Watier 78400 France  |4 aut 
950 |B NATIONALLICENCE  |P 700  |E 1-  |a Peirano  |D Eric  |u ADEME-Der, 500 Route des Lucioles 06560 Valbonne France  |4 aut 
950 |B NATIONALLICENCE  |P 700  |E 1-  |a Chibbaro  |D Sergio  |u Dipartimento di fisica and INFN, Università di Cagliari, Cittadella Universitaria Monserrato CA, Italy  |4 aut 
950 |B NATIONALLICENCE  |P 773  |E 0-  |t Monte Carlo Methods and Applications  |d Walter de Gruyter  |g 9/2(2003-04-01), 93-133  |x 0929-9629  |q 9:2<93  |1 2003  |2 9  |o mcma 
900 7 |b CC0  |u http://creativecommons.org/publicdomain/zero/1.0  |2 nationallicence 
898 |a BK010053  |b XK010053  |c XK010000 
949 |B NATIONALLICENCE  |F NATIONALLICENCE  |b NL-gruyter