In this study, the aerodynamic characteristics of a vertical-axis wind turbine blade coupled with a high-lift device, such as the Gurney flap at the trailing edge, are investigated. For numerical analysis, the force element theory is used to understand how the Gurney flap influences the force evolution of the lift-type vertical-axis wind turbine. This study shows that the lift and drag can be respectively approximated into four elements, which are induced by volume vorticity, rotational velocity, angular acceleration and surface friction of the flow around the blades. Based on the perspective of the force element theory, the present simulation provides a clear picture of how the Gurney flap influences the formation of the aerodynamic force elements during a rotational cycle for a vertical-axis wind turbine. Simulation results show that the contributions mainly result from the surface vorticities, the rotational acceleration of the airfoil, and the acceleration of the surface.

The objective of this work is to computationally investigate the impact of an incidence-tolerant rotor blade concept on gas turbine engine performance under off-design conditions. When a gas turbine operates at an off-design condition such as hover flight or takeoff, large-scale flow separation can occur around turbine blades, which causes performance degradation, excessive noise, and critical loss of operability. To alleviate this shortcoming, a novel concept which articulates the rotating turbine blades simultaneous with the stator vanes is explored. We use a finite-element-based moving-domain computational fluid dynamics (CFD) framework to model a single high-pressure turbine stage. The rotor speeds investigated range from 100% down to 50% of the designed condition of 44,700 rpm. This study explores the limits of rotor blade articulation angles and determines the maximal performance benefits in terms of turbine output power and adiabatic efficiency. The results show articulating rotor blades can achieve an efficiency gain of 10% at off-design conditions thereby providing critical leap-ahead design capabilities for the U.S. Army Future Vertical Lift (FVL) program.

In glaze ice conditions, beads on the surface usually grow to form roughness elements through coalescence, finally resulting in enhancement of local collection efficiency. However, the effects of roughness elements due to freezing of beads are not reflected on the local collection efficiency in CFD icing simulations. This is problematic for predicting the resultant ice shape, which may lead to inaccurate aerodynamic performance and load distribution. The aim of this study is to propose a macroscopic icing model which can reflect bead microscopic phenomena using the Eulerian approach. To this end, a correction was made for collection efficiency by introducing a novel parameter - the effective impinging angle- which is the angle to calculate the local collection efficiency depending on the physical state of surface. It is assumed that the parameter related to the contact angle represents the state of beads. The computational icing analysis of airfoil was performed using the proposed model both in the rime condition and glaze conditions. The results show that the icing characteristics in the feather region is captured with enhanced accuracy in both conditions.

Residual deformation and failure are two critical issues in powder bed fusion (PBF) additive manufacturing (AM) of metal products. Residual deformation caused by the non-uniform residual stress distribution dramatically affects the quality of AM product and can result in catastrophic failure in operation. Therefore, the development of an effective numerical approach to predict residual deformation and failure characteristics of AM product is always a major concern in industrial applications.

In this paper, a numerical approach in predicting residual distortion, stress and failure in AM products is presented. The conventional inherent strain method used in welding process is modified to consider the specific characteristic of AM process, such as the influences of reheating and scanning pattern. This approach consists of three simulation steps including a detailed process simulation in small-scale, a onetime static mechanical finite element analysis in part-scale, and a material failure analysis. First, the inherent strains are calculated from a thermo-mechanical process simulation in small-scale, which considers AM process parameters, such as laser power, scanning speed and path. The physical state in deposited materials including powder, liquid and solid states are taken into account in the simulation by specifying the solidus and liquidus temperature and corresponding material properties. Then the inherent strains are applied layer by layer to the part-scale simulation, where the residual distortion and stress can be predicted efficiently. Finally, a Lagrange particle method is utilized to study the failure characteristics of AM products. Numerical examples are studied to investigate the effectiveness and applicability of present approach.

To understand the lubrication-dominated permeation through a membrane, numerical simulations of permeation through a moving corrugated permeable membrane is carried out with a fully validated numerical method. Through comparisons between the numerical results and the results of an asymptotic analysis of permeate flux (under an infinitesimal permeability condition) using Reynolds lubrication equation, the effect of permeation on lubrication and its inverse effect (i.e., the dependence of permeation on lubrication) are discussed. The linear and non-linear dependences of the relaxation of the lubrication pressure due to membrane permeation are identified. The effect of the tangential component of the permeate flux is evaluated by a linear analysis, and the limitation of Reynolds-type lubrication is discussed.

The quest for drag reduction is driven by environmental concerns, in general, and the need to reduce fuel consumption in transport applications, in particular. Turbulent friction is especially important in civil aviation, accounting for over 50% of the total drag in cruise. In this context, spatially and/or temporally varying in-plane wall motion, while undoubtedly difficult to implement in practice, has attracted major interest, because of the large drag-reduction margins it yields. It is also a forcing method that is of fundamental interest, as it provokes intriguing interactions between the spanwise Stokes layer induced by the wall motion and the near-wall turbulence-regeneration mechanisms. This article provides a relatively brief, ‘entry-level’, review of research in this area, principally over the past two decades. While far from being exhaustive, the review conveys a reasonably detailed picture of some major physical issues as well as of the outcome of the most important computational and experimental studies. Particular emphasis is placed on the question of how results obtained in idealised laboratory conditions and by simulation at relatively low Reynolds-number values pertain to high values typical of high-speed transport.

In the present study, we use direct numerical simulation to investigate the density-driven convection in a two-dimensional anisotropic heterogeneous porous media associated with significant laminated formation. At first, the heterogeneous porous media are randomly generated to represent laminated structure, in which the horizontal correlation length of permeability field is much longer than the vertical counterpart. Then, a highly accurate pseudo-spectral method and compact finite difference scheme with higher order of accuracy are employed to numerically reproduce the convection flow in the laminated porous media. The results show that the laminated structures restrict interactions among the downward plumes of heavier fluid. The plumes tend to descend more straightly in a laminated porous medium associated with a slower growth rate. As a result, the laminated distribution of permeability is considered having an inhibiting effect on the convection flow.

Existing numerical schemes used to solve the governing equations for compressible flow suffer from dissipation errors which tend to smear out sharp discontinuities. Hybrid schemes show potential improvements in this challenging problem; however, the solution quality of a hybrid scheme heavily depends on the criterion to switch between the different candidate reconstruction functions. This work presents a new type of switching criterion (or selector) using machine learning techniques. The selector is trained with randomly generated samples of continuous and discontinuous data profiles, using the exact solution of the governing equation as a reference. Neural networks and random forests were used as the machine learning frameworks to train the selector, and it was later implemented as the indicator function in a hybrid scheme which includes THINC and WENO-Z as the candidate reconstruction functions. The trained selector has been verified to be effective as a reliable switching criterion in the hybrid scheme, which significantly improves the solution quality for both advection and Euler equations.

Direct numerical simulations have been applied to simulate flows with polymer additives. FENE-P (finite-extensible-nonlinear-elastic-Peterlin) dumbbell model solving for the conformation tensor is adopted to investigate the influence of the polymer on the flowfield. Boundary treatments of the conformation tensor on the flowfield are examined first, where boundary condition based on the linear extrapolation scheme provides more accurate results with second-order accurate error norms. Further simulations of the turbulent channel flow at different Weissenberg numbers are also conducted to investigate the influence on drag reduction. Drag reduction increases in tandem with the increase of Weissenberg number and the increase saturates at Weτ~200, where the drag reduction is close to the maximum drag reduction (MDR) limit. At the regime of y+ > 5, the viscous layer thickens with the increase of the Weissenberg number showing a departure from the traditional log-law profile, and the velocity profiles approach the MDR line at high Weissenberg number. The Reynolds stress decreases in tandem with the increase of Weτ, whereas the levels of laminar stress and polymer stress act adversely. However, as the Weissenberg number increases, the proportion of the laminar stress in the total stress increases, and this contributes to the drag reduction of the polymer flow.

Icing on aircraft can drastically reduce aerodynamic performance and lead to serious accidents. Therefore, prediction of the accreted ice shape and area and its effects on aerodynamic performance is crucial during the design phase of an aircraft. However, numerical simulations based on conventional grid-based methods such as the finite volume method cannot accurately reproduce the complex ice shapes, which involve horn growth, feather growth, air voids, and severe surface roughness. In the present study, instead of the grid-based method, a hybrid grid- and particle-based method was newly proposed and applied to the icing problem on a NACA0012 airfoil. The explicit moving particle semi-implicit method was employed as the particle-based method due to its short computing time. The numerical simulations effectively reproduced feather-shaped ice, air voids, and surface roughness. Finally, by computing the flow around the iced airfoil, it was confirmed that flow separation around the leading edge occurred due to the ice layer, which resulted in a thicker boundary layer and wake and an increase in the drag coefficient of approximately 70% after a residence time of only 60 seconds.

This work examines the μ(I) relation that describes the effective friction coefficient μ of a dense granular flow as a function of flow inertial number I, at the center of a rotating drum from its flow onset to steady state using DEM. We want to see how the internal friction coefficient of an accelerating flow may be predicted so that the associated tangential stress can be estimated with the proper knowledge of the normal stress. Under the three investigated drum speeds (3, 4.5 and 6 rpm), the bulk normal stress, σn(y), is found to be a consistent linear depth profile throughout the flow development with a slope degraded from the hydrostatic value, Ph(y), due to lateral wall friction. With the discovery of a non-constant depth-decaying effective wall friction coefficient, we derive analytically a wall-degradation factor K(h) to give σn(y)= K(h)Ph(y). The depth profile of tangential stress, however, varies in time from a concave shape upon acceleration, τa(y), to a more linear trend at the steady state, τss(y). Hence, the μa-Ia profile (with μa=τ/σn) upon flow acceleration offsets from the steady μss(Iss) relation. A pseudo-steady acceleration modification number, ΔI, is proposed to shift the inertial number in the acceleration phase to I* = Ia+ΔI so that the μa-I* data converge to μss(Iss). This finding shall allow us to predict a transient tangential stress by τa(y) = μss(I*)K(y)Ph(y) using the well-accepted knowledge of steady flow rheology, hydrostatic pressure, and the currently developed wall-degradation factor.

In this research, the thermal transport behavior of the branched carbon nanotube (CNT) with T-junction was investigated using non-equilibrium molecular dynamics simulation. Both symmetric and asymmetric temperature-controlled simulations were imposed to evaluate how the heat flowed inside the branched CNT with three branches of equal length and same chirality. The branch length and strain effects on the heat flow were examined. In addition, the simulated heat flow was compared with the prediction made by conventional thermal circuit calculation based on diffusive phonon transport. The heat was observed to flow straight rather than sideway inside the branched CNT with T-junction under the asymmetric temperature setup; this finding contradicts the conventional thermal circuit calculation. There are two possible explanations for this phenomenon. One is ballistic phonon transport and the other is phonons have different interactions or scattering with the defective atomic configurations at the T-junction. Moreover, the tensile strain could tune the heat flow, a finding that might be useful in thermal management applications.