International Journal for Numerical Methods in Fluids最新文献

In flapping insect wings, veins support flexible wing membranes such that the wings form feathering and cambering motions passively from large elastic deformations. These motions are essentially important in unsteady aerodynamics of insect flapping flight. Hence, the underlying mechanism of this phenomenon is an important issue in studies on insect flight. Systematic parametric studies on strong coupling between a model wing describing these elastic deformations and the surrounding fluid, which is a direct formulation of this phenomenon, will be effective for solving this issue. The purpose of this study is to develop a robust numerical framework for these systematic parametric studies. The proposed framework consists of two novel numerical methods: (1) A fully parallelized solution method using both algebraic splitting and semi-implicit scheme for monolithic fluid–structure interaction (FSI) equation systems, which is numerically stable for a wide range of properties such as solid-to-fluid mass ratios and large body motions, and large elastic deformations. (2) A structural mechanics model for insect flapping wings using pixel modeling (pixel model wing), which is combined with explicit node-positioning to reduce computational costs significantly in controlling fluid meshes. The validity of the proposed framework is demonstrated for some benchmark problems and a dynamically scaled model incorporating actual insect data. Finally, from a parametric study for the pixel model wing flapped in fluid with a wide range of solid-to-fluid mass ratios, we find a FSI mechanism of feathering and cambering motions in flapping insect wings.

The study of multilayered shallow water equations has developed from a consideration of immiscible layers as a vertical mesh to the layer bounds as imaginary extremes for vertical integration of the flow equations. In the current work, a quasi three-dimensional flow model has been developed with the consideration of the spatiotemporal flexibility/variability of the pervious vertical discretization/layer ratios. Thus, in principle, vertical layering offers a nonuniform grid with a temporal variation. The system of equations thus formulated comprises a conservative part and the appended source/sink terms. These source/sink terms pertain to the inter-layer interactions such as mass/momenta transfer and interfacial stress, which have been treated in a novel implicit form alongwith the subgrid-scale eddy-viscosity for interlayer shear. They are integrated into the system through different physical considerations so as to arrive at a well-balanced numerical scheme in a regular finite volume grid. The model has been validated through the standard test-cases highlighting the conservation properties and the model's adaptability to uniform and nonuniform vertical meshes alongwith the spatiotemporal transitions of layer ratios, with a specific interest in limiting cases of wet/dry fronts. The increase in layer ratios tends to nearly replicate the full-scale model results in experimental scenarios at a lesser computational overhead.

A proper estimation of flow resistance coefficient of river is essential for precise simulations of river hydraulics. In addition to the cross-sectional geometry and hydraulic parameters, the alignment of the channel affects the flow resistance coefficient in case of meandering rivers. In the present study, a rigorous field study of 131 km along the Barak River was conducted to assess the influence of meandering on the flow resistance coefficient. The values of flow resistance co-efficient were calculated using Chezy and Manning's equations with measured field data and the values from both are compared. However, the variation in the flow resistance co-efficient along the channel calculated from Manning's equation is significantly less as it does not consider the undulation and meandering. Using these field data, an artificial neural network (ANN) model has been developed to predict the cross-sectional averaged flow resistance for meandering river. The model considered the influence of relative curvature, depth of flow, bed particle size, Froude number and Reynolds number including water temperature for accurate predictions of flow resistance coefficient. The ANN model was tested and validated using 237 field data sample. The values of the statistical parameters indicate a very good fit to the training dataset with coefficient of determination (R²) = 0.9566 for training and good fit for testing with R² = 0.8131. The developed ANN model has been compared with other model with the same data set to check its applicability.

Particle-resolved simulations of turbulent particle-laden flows provide a powerful research tool to explore detailed flow physics at all scales. However, efficient particle-resolved simulations for wall-bounded turbulent particle-laden flows remain a challenging task. In this article, we develop a simulation approach for a turbulent channel flow laden with finite-size particles on a nonuniform mesh by combining the discrete unified gas kinetic scheme (DUGKS) and the immersed boundary method (IBM). The standard discrete delta function was modified according to reproducible kernel particle method to take into account mesh non-uniformity and correctly conserve force moments. Simulation results based on uniform and nonuniform meshes are compared to validate and examine the accuracy of the nonuniform mesh DUGKS-IBM. Finally, the computational performance of the nonuniform mesh DUGKS-IBM is discussed.

In this paper, a slip-wall model for large eddy simulation (LES) is improved and implemented in an immersed boundary method named the local domain-free discretization (DFD) method. Considering that the matching location may be in the viscous sublayer, the physics-based interpretation of a Robin-type wall closure is complemented. Then, the slip-wall wall model is improved, in which the slip length is redefined and the Robin boundary condition is imposed at the solid wall. The improved slip-wall model is implemented in the local DFD method to evaluated the tangential velocity at an exterior dependent node, and then the requirement on high resolution of boundary layers can be alleviated. The non-equilibrium effects are accounted for by adding an explicit correction to the wall shear stress. In order to validate the present wall-modeled LES/DFD method, a series of turbulent channel flows at various Reynolds numbers, the flow over periodic hills and the flows over a NACA 4412 airfoil at a high Reynolds number are simulated. The predicted results agree well with the referenced experimental data and numerical results. Especially, the results of the separated flow over the airfoil at a near-stall condition demonstrate the performance of the present wall-modeled LES/DFD method for complex flows.

The Reynolds-Averaged Navier–Stokes (RANS) model is the main model in engineering applications today. However, the normal value of the closure coefficient of the RANS turbulence model is determined based on some simple basic flows and may no longer be applicable for complex flows. In this paper, the closure coefficient of shear stress transport (SST) turbulence model is recalibrated by combining Bayesian method and particle swarm optimization algorithm, so as to improve the numerical simulation accuracy of wall pressure in supersonic flow. First, the obtained prior samples were numerically calculated, and the Sobol index of the closure coefficient was calculated by sensitivity analysis method to characterize the sensitivity of the wall pressure to the model parameters. Second, combined with the uncertainty of propagation parameters by non-intrusive polynomial chaos (NIPC). Finally, Bayesian optimization is used to quantify the uncertainty and obtain the maximum likelihood function estimation and optimal parameters. The results show that the maximum relative error of wall pressure predicted by the SST turbulence model decreases from 29.71% to 9.00%, and the average relative error decreases from 9.86% to 3.67% through the parameter calibration of Bayesian optimization method. In addition, the system evaluated the calibration effect of three criteria, and the calibration results parameters under the three criteria were all better than the calculated results of the nominal values. Meanwhile, the velocity profile and density profile of the flow field were also analyzed. Finally, the same calibration method was applied to the supersonic hollow cylinder and BSL (Baseline) turbulence model, and the same calibration results were obtained, which verified the universality of the calibration method.

We present a novel and simple yet intuitive approach to the stabilization problem for the numerically solved Euler equations with gravity source term relying on a low-order nodal Discontinuous Galerkin Method (DGM). Instead of assuming isothermal or polytropic solutions, we only take a hydrostatic balance as a given property of the flow and use the hydrostatic equation to calculate a hydrostatic pressure reconstruction that replaces the gravity source term. We compare two environments that both solve the Euler equations using the DGM: deal.II and StormFlash. We utilize StormFlash as it allows for the use of the novel stabilization method. Without stabilization, StormFlash does not yield results that resemble correct physical behavior while the results with stabilization for StormFlash, as well as deal.II model the fluid flow more accurately. Convergence rates for deal.II do not match the expected order while the convergence rates for StormFlash with the stabilization scheme (with the exceptions for the L$��{}_{�2�}�$ errors for momentum) meet the expectation. The results from StormFlash with stabilization also fit reference solutions from the literature much better than those from deal.II. We conclude that this novel scheme is a low cost approach to stabilize the Euler equations while not limiting the flow in any way other than it being in hydrostatic balance.

This paper addresses the challenges in studying the interaction between high-intensity sound waves and large-deformable hyperelastic solids, which are characterized by nonlinearities of the hyperelastic material, the finite-amplitude acoustic wave, and the large-deformable fluid–solid interface. An implicit coupling method is proposed for predicting nonlinear structural-acoustic responses of the large-deformable hyperelastic solid submerged in a compressible viscous fluid of infinite extent. An arbitrary Lagrangian–Eulerian (ALE) formulation based on an unsplit complex-frequency-shifted perfectly matched layer method is developed for long-time simulation of the nonlinear acoustic wave propagation without exhibiting long-time instabilities. The solid and acoustic fluid domains are discretized using the finite element method, and two different options of staggered implicit coupling procedures for nonlinear structural-acoustic interactions are developed. Theoretical formulations for stability analysis of the implicit methods are provided. The accuracy, robustness, and convergence properties of the proposed methods are evaluated by a benchmark problem, that is, a hyperelastic rod interacting with finite-amplitude acoustic waves. The numerical results substantiate that the present methods are able to provide long-time steady-state solutions for a nonlinear coupled hyperelastic solid and viscous acoustic fluid system without numerical constraints of small time step sizes and long-time instabilities. The methods are applied to investigate nonlinear dynamic behaviors of coupled hyperelastic elliptical ring and acoustic fluid systems. Physical insights into 2:1 and 4:2:1 internal resonances of the hyperelastic elliptical ring and period-doubling bifurcations of the structural and acoustic responses of the system are provided.