Single-point measurements from towers in cities cannot properly quantify the impact of all terms in the turbulent kinetic energy (TKE) budget and are often not representative of horizontally-averaged ...quantities over the entire urban domain. A series of large-eddy simulations (LES) is here performed to quantify the relevance of non-measurable terms, and to explore the spatial variability of the flow field over and within an urban geometry in the city of Basel, Switzerland. The domain has been chosen to be centered around a tower where single-point turbulence measurements at six heights are available. Buildings are represented through a discrete-forcing immersed boundary method and are based on detailed real geometries from a surveying dataset. The local model results at the tower location compare well against measurements under near-neutral stability conditions and for the two prevailing wind directions chosen for the analysis. This confirms that LES in conjunction with the immersed boundary condition is a valuable model to study turbulence and dispersion within a real urban roughness sublayer (RSL). The simulations confirm that mean velocity profiles in the RSL are characterized by an inflection point
z
γ
located above the average building height
z
h
. TKE in the RSL is primarily produced above
z
γ
, and turbulence is transported down into the urban canopy layer. Pressure transport is found to be significant in the very-near-wall regions. Further, spatial variations of time-averaged variables and non-measurable dispersive terms are important in the RSL above a real urban surface and should therefore be considered in future urban canopy parametrization developments.
Two multi-thread based parallel implementations of the lattice Boltzmann method for non-uniform grids on different hardware platforms are compared in this paper: a multi-core CPU implementation and ...an implementation on General Purpose Graphics Processing Units (GPGPU). Both codes employ second order accurate compact interpolation at the interfaces, coupling grids of different resolutions. Since the compact interpolation technique is both simple and accurate, it produces almost no computational overhead as compared to the lattice Boltzmann method for uniform grids in terms of node updates per second. To the best of our knowledge, the current paper presents the first study on multi-core parallelization of the lattice Boltzmann method with inhomogeneous grid spacing and nested time stepping for both CPUs and GPUs.
In this contribution a numerical study of a turbulent jet flow is presented. The simulation results of two different variants of the Lattice Boltzmann method (LBM) are compared. The first is the ...well-established D3Q19 MRT model extended by a Smagorinsky Large Eddy Simulation (LES) model. The second is the D3Q27 Factorized Cascaded Lattice Boltzmann (FCLB) model without any additional explicit turbulence model. For this model no studies of turbulent flow with high resolution on nonuniform grids existed so far. The underlying computational procedure uses a time nested refinement technique and a grid with more than a billion DOF. The simulations were conducted with the parallel multi physics solver VirtualFluids. It is shown that both models are feasible for the present flow case, but the FCLB outperforms the traditional approach in some aspects.
A very efficient implementation of a lattice Boltzmann (LB) kernel in 3D on a graphical processing unit using the compute unified device architecture interface developed by nVIDIA is presented. By ...exploiting the explicit parallelism offered by the graphics hardware, we obtain an efficiency gain of up to two orders of magnitude with respect to the computational performance of a PC. A non-trivial example shows the performance of the LB implementation, which is based on a D3Q13 model that is described in detail.
In this work we use two numerical methods which rely only on the geometry and material parameters to predict capillary hysteresis in a porous material. The first numerical method is a morphological ...pore network (MPN) model, where structural elements are inserted into the imaged pore space to quantify the local capillary forces. Then, based on an invasion–percolation mechanism, the fluid distribution is computed. The second numerical method is a lattice-Boltzmann (LB) approach which solves the coupled Navier–Stokes equations for both fluid phases and describes the dynamics of the fluid/fluid interface. We have developed an optimized version of the model proposed in Tölke J, Freudiger S, Krafczyk M. An adaptive scheme for LBE multiphase flow simulations on hierarchical grids, Comput. Fluids 2006;35:820–30 for the type of flow problems encountered in this work. A detailed description of the model and an extensive validation of different multiphase test cases have been carried out. We investigated pendular rings in a sphere packing, static and dynamic capillary bundle models and the residual saturation in a sphere packing.
A sample of 15
mm in diameter filled with sand particles ranging from 100 to 500
μm was scanned using X-rays from a synchrotron source with a spatial resolution of 11
μm. Based on this geometry we computed the primary drainage, the first imbibition and the secondary drainage branch of the hysteresis loop using both approaches. For the LB approach, we investigated the dependence of the hysteresis loop on the speed of the drainage and the imbibition process. Furthermore we carried out a sensitivity analysis by simulating the hysteretic effect in several subcubes of the whole geometry with extremal characteristic properties. The predicted hysteretic water retention curves were compared to the results of laboratory experiments using inverse modeling based on the Richards equation.
A good agreement for the hysteresis loop between the LB and MPN model has been obtained. The primary and secondary drainage of the hysteresis loop of the LB and MPN model compare very well, and also the experimental results fit well with a slight offset of 10% in the amplitude. Differences for the first imbibition have been observed, but also large differences between two different experimental runs have been observed.
In this work, the suitability of the lattice Boltzmann method is evaluated for the simulation of subcritical turbulent flows around a sphere. Special measures are taken to reduce the computational ...cost without sacrificing the accuracy of the method. A large eddy simulation turbulence model is employed to allow efficient simulation of resolved flow structures on non-uniform computational meshes. In the vicinity of solid walls, where the flow is governed by the presence of a thin boundary layer, local grid-refinement is employed in order to capture the fine structures of the flow. In the test case considered, reference values for the drag force in the Reynolds number range from 2000 to 10
000 and for the surface pressure distribution and the angle of separation at a Reynolds number of 10
000 could be quantitatively reproduced. A parallel efficiency of 80% was obtained on an Opteron cluster.
We propose a scheme to reconstruct the flow field up to second order in a cell composed of four lattices nodes in 2 dimensions. The information contained in the higher order moments of the ...distribution functions is used to construct an interpolation scheme of second order for the velocity field.
Evaporation is an important process in many natural and technical systems, such as the unsaturated zone of the subsurface or microchannel evaporators. For the understanding and prediction of the ...involved processes, numerical simulations of multiphase flow and transport processes are an important tool. In order to achieve an accurate, physically based description of kinetic interphase mass and heat transfer occurring during evaporation, the numerical model has to account for the interfacial areas between phases. A recently developed model for two‐phase flow in porous media is able to account for the involved processes by using interfacial areas explicitly as parameters in the model. The crucial issue, however, is the determination of the relationships between specific interfacial areas, capillary pressure, and saturation in this paper, we present a multiphase lattice Boltzmann model, which allows us to determine these relationships. On the basis of the scanned geometry of a natural porous medium, the relationships between specific interfacial areas, capillary pressure, and saturation are determined. To the best of our knowledge, this is the first time that fluid‐solid specific interfacial area relationships have been obtained from pore‐scale data. Using these functions, we present the results of macroscale simulations of an evaporator device and of drying in a porous medium.
Key Points
Obtaining constitutive relationships from micro scale LB simulations
Thermodynamically consistent method including interfacial areas
Modeling mass and heat tranfer kinetics across interfaces