Search results

The reduced problem of the Navier–Stokes and the continuity equations, in two-dimensional Cartesian coordinates with Eulerian description, for incompressible non-Newtonian fluids, is considered. The Ladyzhenskaya model, with a non-linear velocity dependent stress tensor is adopted, and leads to the governing equation of interest. The reduction is based on a self-similar transformation as demonstrated in existing literature, for two spatial variables and one time variable, resulting in an ODE defined on a semi-infinite domain. In our search for classical solutions, existence and uniqueness will be determined depending on the signs of two parameters with physical interpretation in the equation. Illustrations are included to highlight some of the main results.

This study aims to simulate the dendritic growth in Stokes flow by iteratively coupling a domain and boundary type meshless method.

A preconditioned phase-field model for dendritic solidification of a pure supercooled melt is solved by the strong-form space-time adaptive approach based on dynamic quadtree domain decomposition. The domain-type space discretisation relies on monomial augmented polyharmonic splines interpolation. The forward Euler scheme is used for time evolution. The boundary-type meshless method solves the Stokes flow around the dendrite based on the collocation of the moving and fixed flow boundaries with the regularised Stokes flow fundamental solution. Both approaches are iteratively coupled at the moving solid–liquid interface. The solution procedure ensures computationally efficient and accurate calculations. The novel approach is numerically implemented for a 2D case.

The solution procedure reflects the advantages of both meshless methods. Domain one is not sensitive to the dendrite orientation and boundary one reduces the dimensionality of the flow field solution. The procedure results agree well with the reference results obtained by the classical numerical methods. Directions for selecting the appropriate free parameters which yield the highest accuracy and computational efficiency are presented.

A combination of boundary- and domain-type meshless methods is used to simulate dendritic solidification with the influence of fluid flow efficiently.

This paper aims to analyze the pressure distribution of rectangular aerostatic thrust bearing with a single air supply inlet using the complex potential theory and conformal mapping.

The Möbius transform is used to map the interior of a rectangle onto the interior of a unit circle, from which the pressure distribution and load carrying capacity are obtained. The calculation results are verified by finite difference method.

The constructed Möbius formula is very effective for the performance characteristics researches for the rectangular thrust bearing with a single air supply inlet. In addition, it is also noted that to obtain the optimized load carrying capacity, the square thrust bearing can be adopted.

The Möbius transform is found suitable to describe the pressure distribution of the rectangular thrust bearing with a single air supply inlet.

The paper targets on providing new experimental data for validation of the well-established mathematical models within the framework of the lattice Boltzmann method (LBM), which are applied to problems of casting processes in complex mould cavities.

An experimental campaign aiming at the free-surface flow within a system of narrow channels is designed and executed under well-controlled laboratory conditions. An in-house lattice Boltzmann solver is implemented. Its algorithm is described in detail and its performance is tested thoroughly using both the newly recorded experimental data and well-known analytical benchmark tests.

The benchmark tests prove the ability of the implemented algorithm to provide a reliable solution when the surface tension effects become dominant. The convergence of the implemented method is assessed. The two new experimentally studied problems are resolved well by simulations using a coarse computational grid.

A detailed set of original experimental data for validation of computational schemes for simulations of free-surface gravity-driven flow within a system of narrow channels is presented.

The goal of this paper is to give a comprehensive and short review on how to compute the first- and second-order topological derivatives and potentially higher-order topological derivatives for partial differential equation (PDE) constrained shape functionals.

The authors employ the adjoint and averaged adjoint variable within the Lagrangian framework and compare three different adjoint-based methods to compute higher-order topological derivatives. To illustrate the methodology proposed in this paper, the authors then apply the methods to a linear elasticity model.

The authors compute the first- and second-order topological derivatives of the linear elasticity model for various shape functionals in dimension two and three using Amstutz' method, the averaged adjoint method and Delfour's method.

In contrast to other contributions regarding this subject, the authors not only compute the first- and second-order topological derivatives, but additionally give some insight on various methods and compare their applicability and efficiency with respect to the underlying problem formulation.

The purpose of the paper is to predict the aerodynamic performance of a complete scale model H-Darrieus vertical axis wind turbine (VAWT) with end plates at different operating conditions. This paper aims at understanding the flow physics around a model VAWT for three different tip speed ratios corresponding to three different flow regimes.

This study achieves a first three-dimensional hybrid lattice Boltzmann method/very large eddy simulation (LBM-VLES) model for a complete scaled model VAWT with end plates and mast using the solver PowerFLOW. The power curve predicted from the numerical simulations is compared with the experimental data collected at Erlangen University. This study highlights the complexity of the turbulent flow features that are seen at three different operational regimes of the turbine using instantaneous flow structures, mean velocity, pressure iso-contours, blade loading and skin friction plots.

The power curve predicted using the LBM-VLES approach and setup provides a good overall match with the experimental power curve, with the peak and drop after the operational point being captured. Variable turbulent flow structures are seen over the azimuthal revolution that depends on the tip speed ratio (TSR). Significant dynamic stall structures are seen in the upwind phase and at the end of the downwind phase of rotation in the deep stall regime. Strong blade wake interactions and turbulent flow structures are seen inside the rotor at higher TSRs.

The computational cost and time for such high-fidelity simulations using the LBM-VLES remains expensive. Each simulation requires around a week using supercomputing facilities. Further studies need to be performed to improve analytical VAWT models using inputs/calibration from high fidelity simulation databases. As a future work, the impact of turbulent and nonuniform inflow conditions that are more representative of a typical urban environment also needs to be investigated.

The LBM methodology is shown to be a reliable approach for VAWT power prediction. Dynamic stall and blade wake interactions reduce the aerodynamic performance of a VAWT. An ideal operation close to the peak of the power curve should be favored based on the local wind resource, as this point exhibits a smoother variation of forces improving operational performance. The 3D flow features also exhibit a significant wake asymmetry that could impact the optimal layout of VAWT clusters to increase their power density. The present work also highlights the importance of 3D simulations of the complete model including the support structures such as end plates and mast.

Accurate predictions of power performance for Darrieus VAWTs could help in better siting of wind turbines thus improving return of investment and reducing levelized cost of energy. It could promote the development of onsite electricity generation, especially for industrial sites/urban areas and renew interest for VAWT wind farms.

A first high-fidelity simulation of a complete VAWT with end plates and supporting structures has been performed using the LBM approach and compared with experimental data. The 3D flow physics has been analyzed at different operating regimes of the turbine. These physical insights and prediction capabilities of this approach could be useful for commercial VAWT manufacturers.

Current methods for flow field reconstruction mainly rely on data-driven algorithms which require an immense amount of experimental or field-measured data. Physics-informed neural network (PINN), which was proposed to encode physical laws into neural networks, is a less data-demanding approach for flow field reconstruction. However, when the fluid physics is complex, it is tricky to obtain accurate solutions under the PINN framework. This study aims to propose a physics-based data-driven approach for time-averaged flow field reconstruction which can overcome the hurdles of the above methods.

A multifidelity strategy leveraging PINN and a nonlinear information fusion (NIF) algorithm is proposed. Plentiful low-fidelity data are generated from the predictions of a PINN which is constructed purely using Reynold-averaged Navier–Stokes equations, while sparse high-fidelity data are obtained by field or experimental measurements. The NIF algorithm is performed to elicit a multifidelity model, which blends the nonlinear cross-correlation information between low- and high-fidelity data.

Two experimental cases are used to verify the capability and efficacy of the proposed strategy through comparison with other widely used strategies. It is revealed that the missing flow information within the whole computational domain can be favorably recovered by the proposed multifidelity strategy with use of sparse measurement/experimental data. The elicited multifidelity model inherits the underlying physics inherent in low-fidelity PINN predictions and rectifies the low-fidelity predictions over the whole computational domain. The proposed strategy is much superior to other contrastive strategies in terms of the accuracy of reconstruction.

In this study, a physics-informed data-driven strategy for time-averaged flow field reconstruction is proposed which extends the applicability of the PINN framework. In addition, embedding physical laws when training the multifidelity model leads to less data demand for model development compared to purely data-driven methods for flow field reconstruction.

The purpose of this paper is to present a fast and bare bones implementation of a numerical method for quickly simulating turbulent thermal flows on GPUs. The work also validates earlier research showing that the lattice Boltzmann method (LBM) method is suitable for complex thermal flows.

A dual lattice hydrodynamic (D3Q27) thermal (D3Q7) multiple-relaxation time LBM model capable of thermal DNS calculations is implemented in CUDA.

The model has the same computational performance compared to earlier publications of similar LBM solvers. The solver is validated against three benchmark cases for turbulent thermal flow with available data and is shown to be in excellent agreement.

The combination of a D3Q27 and D3Q7 stencil for a multiple relaxation time -LBM has, to the authors’ knowledge, not been used for simulations of thermal flows. The code is made available in a public repository under a free license.

This study aims to analyse the non-linear losses of a porous media (stack) composed by parallel plates and inserted in a resonator tube in oscillatory flows by proposing numerical correlations between pressure gradient and velocity.

The numerical correlations origin from computational fluid dynamics simulations, conducted at the microscopic scale, in which three fluid channels representing the porous media are taken into account. More specifically, for a specific frequency and stack porosity, the oscillating pressure input is varied, and the velocity and the pressure-drop are post-processed in the frequency domain (Fast Fourier Transform analysis).

It emerges that the viscous component of pressure drop follows a quadratic trend with respect to velocity inside the stack, while the inertial component is linear also at high-velocity regimes. Furthermore, the non-linear coefficient b of the correlation ax + bx² (related to the Forchheimer coefficient) is discovered to be dependent on frequency. The largest value of the b is found at low frequencies as the fluid particle displacement is comparable to the stack length. Furthermore, the lower the porosity the higher the Forchheimer term because the velocity gradients at the stack geometrical discontinuities are more pronounced.

The main novelty of this work is that, for the first time, non-linear losses of a parallel plate stack are investigated from a macroscopic point of view and summarised into a non-linear correlation, similar to the steady-state and well-known Darcy–Forchheimer law. The main difference is that it considers the frequency dependence of both Darcy and Forchheimer terms. The results can be used to enhance the analysis and design of thermoacoustic devices, which use the kind of stacks studied in the present work.