Solving the Navier-Stokes equations using direct numerical simulation (DNS) is computationally impractical, especially at high Reynolds numbers. Recent technological advances in supercomputing have paved the way for Large Eddy Simulations (LES) to circumvent this problem by resolving large scale turbulence motions and modeling only the small (subgrid) scales. However, LES modeling still requires advanced knowledge of the turbulence and LES models are currently very simplistic. Because of this, there has been considerable interest in hybrid turbulence models, which can perform either Reynolds Averaged Navier-Stokes (RANS) modeling or Large Eddy Simulation (LES). The self-adapting model presented is fundamentally different from prior LES models and these current hybrid models in that it achieves a completely natural evolution from RANS to LES to (with enough mesh resolution) DNS. A modified k/e model and a Reynolds stress transport model is implemented in this manner and is compared to DNS data of isotropic decaying turbulence. The results indicate that this modeling approach is practical and efficient. In addition, this approach is extensible and not restricted to a particular (RANS) transport equation.

The present work addresses the fundamental question involving the modeling of subgrid-scale turbulence as a function of resolved field. A new-nonlinear model has been developed from the constitutive equation of subgrid stresses extending the Reynolds stress model proposed by Warsi. The time scale is expressed in terms of subgrid scale kinetic energy as opposed to strain rate tensor. Effort has been made to identify the terms appearing in the modeled subgrid stresses with "Reynolds term", "Leonard's term" and "cross term". The physical nature of these terms can be best understood from the triadic interactions in wave number space. Understanding these three terms leads to decouple the complex nature of the subgrid stresses. Modeling of these terms separately helps to capture the physics of the problem accurately. The turbulent field is assumed to be isotropic and Kolmogrov's hypothesis is used. The model coefficients are expressed as universal constants for Gaussian filter so as to satisfy the dissipation criteria in inertial subrange. Further dissipation term is assumed to be isotropic and equilibrium condition is used. Although the definition of the subgrid stress terms becomes less clear and separate for smooth filter, an attempt has been made to compare the stress terms with the exact definition obtained for sharp cut-off filter. An estimate of the backscatter of energy can be obtained from the Eddy-Damped Quasi Normal Markovian (EDQNM) theory. The model coefficients thus obtained are tested with results of plain homogeneous shear layer. The model results have been compared with the mixed-nonlinear model and Smagorinsky model. A priori test shows that new-nonlinear model has a good correlation with Smagorinsky model, which in turn has good correlation with experimental results, and has the behavior of the mixed-nonlinear model. The above model has been used for solving two-dimensional flow over backward facing step as a test case. The numerical model solves the vertically hydrostatic boundary layer equation. The top boundary is assumed to be a free surface. Terrain following coordinate system has been used. Because of the non-negativity of the subgrid scale dissipation term i.e. backscatter of energy, the nature of the solution is stochastic. The deterministic solution is obtained by clipping the dissipation term. The results are compared with the experimental data of Kim et al. Good agreement with the experimental data is obtained for the velocity profile and SGS kinetic energy. The reattachment point obtained is at 5.2h (h is the step height), which is less compared to 6h as suggested by other authors. This discrepancy may be due to the assumptions involved in the equations, which is being solved. The model is further extended for the diffusion of scalar variables and to include the buoyancy effect. It is implemented to explore the hydrostatic flow over three dimensional elliptical mountain ridges, where Boussinesq approximation is used for variable density. The flow characteristics have been studied for the various aspect ratios of the mountain and Froude?s number (Nh/U) based on Brunt-Vaisala frequency (N). The phenomenon of upstream blocking and Lee-vortices generation has been studied.

Prediction of the near wall region is crucial to the accuracy of turbulent flow computational fluid dynamics (CFD) simulation. However, sufficient near-wall resolution is often prohibitive for high Reynolds number flows with complex geometries, due to high memory and processing requirements. A common approach in these cases is to use wall functions to bridge the region from the first grid node to the wall. This thesis presents an alternative method that relaxes the near wall resolution requirement by solving one dimensional transport equations for velocity and turbulence across a locally defined subgrid contained within wall adjacent grid cells. The addition of the subgrid allows for wall adjacent primary grid sizes to vary arbitrarily from low-Re model sizing (y+ ~ 1) to wall function sizing without significant loss of accuracy or increase in computational cost.

京都大学 / 新制・課程博士 / 博士(工学) / 甲第23852号 / 工博第4939号 / 新制||工||1771(附属図書館) / 京都大学大学院工学研究科社会基盤工学専攻 / (主査)教授 森 信人, 教授 平石 哲也, 准教授 志村 智也 / 学位規則第4条第1項該当 / Doctor of Philosophy (Engineering) / Kyoto University / DFAM

There are two big challenges which restrict the extensive application of fully coupled geomechanics-reservoir modeling. The first challenge is computational effort. Consider a 3-D simulation combining pressure and heat diffusion, elastoplastic mechanical response, and saturation changes; each node has at least 5 degrees of freedom, each leading to a separate equation. Furthermore, regions of large p, T and σ′ gradients require small-scale discretization for accurate solutions, greatly increasing the number of equations. When the rock mass surrounding the reservoir region is included, it is represented by many elements or nodes. These factors mean that accurate analysis of realistic 3-D problems is challenging, and will so remain as we seek to solve larger and larger coupled problems involving nonlinear responses. To overcome the first challenge, the displacement discontinuity method is introduced wherein a large-scale 3-D case is divided into a reservoir region where Δp, ΔT and non-linear effects are critical and analyzed using FEM, and an outside region in which the reservoir is encased where Δp and ΔT effects are inconsequential and the rock may be treated as elastic, analyzed with a 3D displacement discontinuity formulation. This scheme leads to a tremendous reduction in the degrees of freedom, yet allows for reasonably rigorous incorporation of the reactions of the surrounding rock. The second challenge arises from some forms of numerical instability. There are actually two types of sharp gradients implied in the transient advection-diffusion problem: one is caused by the high Peclet numbers, the other by the sharp gradient which appears during the small time steps due to the transient solution. The way to eliminate the spurious oscillations is different when the sharp gradients are induced by the transient evolution than when they are produced by the advective terms, and existing literature focuses mainly on eliminating the spurious spatial temperature oscillations caused by advection-dominated flow. To overcome the second challenge, numerical instability sources are addressed by introducing a new stabilized finite element method, the subgrid scale/gradient subgrid scale (SGS/GSGS) method.

There are two big challenges which restrict the extensive application of fully coupled geomechanics-reservoir modeling. The first challenge is computational effort. Consider a 3-D simulation combining pressure and heat diffusion, elastoplastic mechanical response, and saturation changes; each node has at least 5 degrees of freedom, each leading to a separate equation. Furthermore, regions of large p, T and σ′ gradients require small-scale discretization for accurate solutions, greatly increasing the number of equations. When the rock mass surrounding the reservoir region is included, it is represented by many elements or nodes. These factors mean that accurate analysis of realistic 3-D problems is challenging, and will so remain as we seek to solve larger and larger coupled problems involving nonlinear responses. To overcome the first challenge, the displacement discontinuity method is introduced wherein a large-scale 3-D case is divided into a reservoir region where Δp, ΔT and non-linear effects are critical and analyzed using FEM, and an outside region in which the reservoir is encased where Δp and ΔT effects are inconsequential and the rock may be treated as elastic, analyzed with a 3D displacement discontinuity formulation. This scheme leads to a tremendous reduction in the degrees of freedom, yet allows for reasonably rigorous incorporation of the reactions of the surrounding rock. The second challenge arises from some forms of numerical instability. There are actually two types of sharp gradients implied in the transient advection-diffusion problem: one is caused by the high Peclet numbers, the other by the sharp gradient which appears during the small time steps due to the transient solution. The way to eliminate the spurious oscillations is different when the sharp gradients are induced by the transient evolution than when they are produced by the advective terms, and existing literature focuses mainly on eliminating the spurious spatial temperature oscillations caused by advection-dominated flow. To overcome the second challenge, numerical instability sources are addressed by introducing a new stabilized finite element method, the subgrid scale/gradient subgrid scale (SGS/GSGS) method.

A low-order discrete dynamical system (DDS) for finite-rate chemistry of H2-air combustion is derived in 3D. Fourier series with a single wavevector are employed to represent dependent variables of subgrid-scale (SGS) behaviors for applications to large-eddy simulation (LES). A Galerkin approximation is applied to the governing equations for comprising the DDS. Regime maps are employed to aid qualitative determination of useful values for bifurcation parameters of the DDS. Both isotropic and anisotropic assumptions are employed when constructing regime maps and studying bifurcation parameters sequences. For H2-air reactions, two reduced chemical mechanisms are studied via the DDS. As input to the DDS, physical quantities from experimental turbulent flow are used. Numerical solutions consisting of time series of velocities, species mass fractions, temperature, and the sum of mass fractions are analyzed. Numerical solutions are compared with experimental data at selected spatial locations within the experimental flame to check whether this model is suitable for an entire flame field. The comparisons show the DDS can mimic turbulent combustion behaviors in a qualitative sense, and the time-averaged computed results of some species are quantitatively close to experimental data.

Les modèles de prévision opérationnelle du temps résolvent numériquement les équations de la mécanique des fluides en calculant l'évolution de champs (pression, température, humidité, vitesses) définis comme moyennes horizontales à l'échelle des mailles d'une grille (et à différents niveaux verticaux). Les processus d'échelle inférieure à la maille jouent néanmoins un rôle essentiel dans les transferts et les bilans de chaleur, humidité et quantité de mouvement. Les paramétrisations physiques visent à évaluer les termes de source correspondant à ces phénomènes, et apparaissant dans les équations des champs moyens aux points de grille. Lorsque l'on diminue la taille des mailles afin de représenter plus finement l'évolution des phénomènes atmosphériques, certaines hypothèses utilisées dans ces paramétrisations perdent leur validité. Le problème se pose surtout quand la taille des mailles passe en dessous d'une dizaine de kilomètres, se rapprochant de la taille des grands systèmes de nuages convectifs (systèmes orageux, lignes de grain). Ce travail s'inscrit dans le cadre des développements du modèle à mailles fines ARPÈGE ALADIN, utilisé par une douzaine de pays pour l'élaboration de prévisions à courte échéance (jusque 48 heures). Nous décrivons d'abord l'ensemble des paramétrisations physiques du modèle. Suit une analyse détaillée de la paramétrisation actuelle de la convection profonde. Nous présentons également notre contribution personnelle à celle ci, concernant l'entraînement de la quantité de mouvement horizontale dans le nuage convectif. Nous faisons ressortir les principaux points faibles ou hypothèses nécessitant des mailles de grandes dimensions, et dégageons les voies pour de nouveaux développements. Nous approfondissons ensuite deux des aspects sortis de cette discussion: l'usage de variables pronostiques de l'activité convective, et la prise en compte de différences entre l'environnement immédiat du nuage et les valeurs des champs à grande échelle. Ceci nous conduit à la réalisation et la mise en œuvre d'un schéma pronostique de la convection profonde. A ce schéma devraient encore s'ajouter une paramétrisation pronostique des phases condensées suspendues (actuellement en cours de développement par d'autres personnes) et quelques autres améliorations que nous proposons. Des tests de validation et de comportement du schéma pronostique ont été effectués en modèle à aire limitée à différentes résolutions et en modèle global. Dans ce dernier cas l'effet du nouveau schéma sur les bilans globaux est également examiné. Ces expériences apportent un éclairage supplémentaire sur le comportement du schéma convectif et les problèmes de partage entre la schéma de convection profonde et le schéma de précipitation de grande échelle. La présente étude fait donc le point sur le statut actuel des différentes paramétrisations du modèle, et propose des solutions pratiques pour améliorer la qualité de la représentation des phénomènes convectifs. L'utilisation de mailles plus petites que 5 km nécessite enfin de lever l'hypothèse hydrostatique dans les équations de grande échelle, et nous esquissons les raffinements supplémentaires de la paramétrisation possibles dans ce cas.

No description available.

This dissertation consists of two parts, i.e. dynamic approaches for subgrid-scale (SGS) stress modelling for large eddy simulation and advanced assessment of the resolved scale motions related to turbulence geometrical statistics and topologies. The numerical simulations are based on turbulent Couette flow. The first part of the dissertation presents four contributions to the development of dynamic SGS models. The conventional integral type dynamic localization SGS model is in the form of a Fredholm integral equation of the second kind. This model is mathematically consistent, but demanding in computational cost. An efficient solution scheme has been developed to solve the integral system for turbulence with homogeneous dimensions. Current approaches to the dynamic two-parameter mixed model (DMM2) are mathematically inconsistent. As a second contribution, the DMM2 has been optimized and a modelling system of two integral equations has been rigorously obtained. The third contribution relates to the development of a novel dynamic localization procedure for the Smagorinsky model using the functional variational method. A sufficient and necessary condition for localization is obtained and a Picard's integral equation for the model coefficient is deduced. Finally, a new dynamic nonlinear SGS stress model (DNM) based on Speziale's quadratic constitutive relation [J. Fluid Mech., 178, p.459, 1987] is proposed. The DNM allows for a nonlinear anisotropic representation of the SGS stress, and exhibits a significant local stability and flexibility in self-calibration. In the second part, the invariant properties of the resolved velocity gradient tensor are studied using recently developed methodologies, i.e. turbulence geometrical statistics and topology. The study is a posteriori based on the proposed DNM, which is different than most of the current a priori approaches based on experimental or DNS databases. The performance of the DNM is further validated in terms of its capability of simulating advanced geometrical and topological features of resolved scale motions. Phenomenological results include, e.g. the positively skewed resolved enstrophy generation, the alignment between the vorticity and vortex stretching vectors, and the pear-shape joint probability function contour in the tensorial invariant phase plane. The wall anisotropic effect on these results is also examined.