A direct-forcing fictitious domain (DF/FD) method for the simulation of particulate flows is reported. The new method is a non-Lagrange-multiplier version of our previous DLM/FD code and is obtained ...by employing a discrete
δ-function in the form of bi(tri-) function to transfer explicitly quantities between the Eulerian and Lagrangian nodes, as in the immersed boundary method. Due to the use of the collocation-point approach for the rigidity constraint and the integration over the particle domain, the Lagrangian nodes are retracted a little from the particle boundary. Our method in case of a prescribed velocity on the boundary is verified via the comparison to the benchmark results on the flow over a fixed cylinder in a wide channel and to our spectral-element results for a channel with the width of four cylinder diameters. We then verify our new method for the case of the particulate flows through various typical flow situations, including the sedimentation of a circular particle in a vertical channel, the sedimentation of a sphere in a vertical pipe, the inertial migration of a sphere in a circular Poiseuille flow, the behavior of a neutrally-buoyant sphere in Couette flow, and the rotation of a prolate spheroid in Couette flow. The accuracy and robustness of the new method are fully demonstrated, in particular for the case of relatively low Reynolds numbers and the neutrally-buoyant case.
Full text
Available for:
GEOZS, IJS, IMTLJ, KILJ, KISLJ, NUK, OILJ, PNG, SAZU, SBCE, SBJE, UL, UM, UPCLJ, UPUK
A fictitious domain method is used to perform fully resolved numerical simulations of particle-laden turbulent flow in a horizontal channel. The effects of large particles of diameter 0.05 and 0.1 ...times the channel height on the turbulence statistics and structures are investigated for different settling coefficients and volume fractions (0.79 %–7.08 %) for the channel Reynolds number being 5000. The results indicate the following. (a) When the particle sedimentation effect is negligible (i.e. neutrally buoyant), the presence of particles decreases the maximum r.m.s. of streamwise velocity fluctuation near the wall by weakening the intensity of the large-scale streamwise vortices, while increasing the r.m.s. of the streamwise fluctuating velocity in the region very close to the wall and in the centre region. On the other hand, the particles increase the r.m.s. of transverse and spanwise fluctuating velocities in the near-wall region by inducing the small-scale vortices. (b) When the particle settling effect is so substantial that most particles settle onto the bottom wall and form a particle sediment layer (SL), the SL plays the role of a rough wall and parts of the vortex structures shedding from the SL ascend into the core region and substantially increase the turbulence intensity there. (c) When the particle settling effect is moderate, the effects of particles on the turbulence are a combination of the former two situations, and the Shields number is a good parameter for measuring the particle settling effects (i.e. the particle concentration distribution in the transverse direction). The average velocities of the particle are smaller in the lower half-channel and larger in the upper half-channel compared to the local fluid velocities in the presence of gravity effects. The effects of the smaller particles on the turbulence are found to be stronger at the same particle volume fractions.
Volumetric measurement for the non-cavitating tip vortex in the near field of an elliptical hydrofoil is conducted using tomographic particle image velocimetry, which provides a fully ...three-dimensional diagnose of the vortex formation and development. The wandering motion and flow properties of the near-field tip vortex under different incident angles and Reynolds numbers are investigated in detail. Unlike in the far field, the wandering motion in the near field is mainly subject to the local flow unsteadiness rather than the flow condition. By the “re-centered” post-processing, the deviations introduced by the wandering motion can be technically corrected, and more accurate vortex properties can be thus obtained. In the near field, a turning point of the vortex center trajectory is detected, the position of which is basically independent of the flow condition. By investigating the local flow properties, it is found that this turning point is the position where the tip vortex completely leaves the trailing edge of hydrofoil and enters the wake region. At this turning point, the external supply to the vortex core starts to be restricted, and the vortex circulation reaches a rather constant value. Further according to the local flow properties, the development process of the near-field tip vortex can be divided into three stages: vortex-attached stage, vortex-lifting stage and vortex-detached stage, which are found to be closely relevant to the hydrofoil configuration.
Full text
Available for:
EMUNI, FIS, FZAB, GEOZS, GIS, IJS, IMTLJ, KILJ, KISLJ, MFDPS, NLZOH, NUK, OILJ, PNG, SAZU, SBCE, SBJE, SBMB, SBNM, UKNU, UL, UM, UPUK, VKSCE, ZAGLJ
The minimum-time path for intercepting a moving target with a prescribed impact angle is studied in the paper. The candidate paths from Pontryagin’s maximum principle are synthesized, so that each ...candidate is related to a zero of a real-valued function. It is found that the real-valued functions or their first-order derivatives can be converted to polynomials of at most fourth degree. As a result, each candidate path can be computed within a constant time by embedding a standard polynomial solver into the typical bisection method. The control strategy along the shortest candidate eventually gives rise to the time-optimal guidance law. Finally, the developments of the paper are illustrated and verified by three numerical examples.
Full text
Available for:
GEOZS, IJS, IMTLJ, KILJ, KISLJ, NLZOH, NUK, OILJ, PNG, SAZU, SBCE, SBJE, UILJ, UL, UM, UPCLJ, UPUK, ZAGLJ, ZRSKP
The distributed-Lagrange-multiplier/fictitious-domain (DLM/FD) method of Glowinski et al. R. Glowinski, T.-W. Pan, T.I. Hesla, D.D. Joseph, A distributed Lagrange multiplier/fictitious domain method ...for particulate flows, Int. J. Multiphase Flow 25 (1999) 755–794 is extended to deal with heat transfer in particulate flows in two dimensions. The Boussinesq approximation is employed for the coupling between the flow and temperature fields. The fluid-flow equations are solved with the finite-difference projection method on a half-staggered grid. In our operator splitting scheme, the Lagrange multipliers at the previous time level are kept in the fluid equations, and the new Lagrange multipliers for the rigid-body motion constraint and the Dirichlet temperature boundary condition are determined from the reduced saddle-point problem, whereas a very simple scheme based on the fully explicit computation of the Lagrange multiplier is proposed for the problem in which the solid heat conduction inside the particle boundary is also considered. Our code for the case of fixed temperature on the immersed boundary is verified by comparing favorably our results on the natural convection driven by a hot cylinder eccentrically placed in a square box and on the sedimentation of a cold circular particle in a vertical channel to the data in the literature. The code for the case of freely varying temperature on the boundaries of freely moving particles is applied to analyze the motion of a catalyst particle in a box and in particular the heat conductivities of nanofluids and sheared non-colloidal suspensions, respectively. Our preliminary computational results support the argument that the micro-heat-convection in the fluids is primarily responsible for the unusually high heat conductivity of nanofluids. It is shown that the Peclet number plays a negative role in the diffusion-related heat conductivity of a sheared non-colloidal suspension, whereas the Reynolds number does the opposite.
Full text
Available for:
GEOZS, IJS, IMTLJ, KILJ, KISLJ, NUK, OILJ, PNG, SAZU, SBCE, SBJE, UL, UM, UPCLJ, UPUK
8.
Improvements on the virtual obstacle method Zheng, Yuan; Shao, Xueming; Chen, Zheng ...
International journal of advanced robotic systems,
03/2020, Volume:
17, Issue:
2
Journal Article
Peer reviewed
Open access
While the artificial potential field has been widely employed to design path planning algorithms, it is well-known that artificial potential field-based algorithms suffer a severe problem that a ...robot may sink into a local minimum point. To address such problems, a virtual obstacle method has been developed in the literature. However, a robot may be blocked by virtual obstacles generated during performing the virtual obstacle method if the environments are complex. In this article, an improved virtual obstacle method for local path planning is designed via proposing a new minimum criterion, a new switching condition, and a new exploration force. All the three new contributions allow to overcome the drawbacks of the artificial potential field-based algorithms and the virtual obstacle method. As a consequence, feasible collision-free paths can be found in complex environments, as illustrated by final numerical simulations.
Variable-sweep wings have large shape-changing capabilities and wide flight envelops, which are considered as one of the most promising directions for intelligent morphing UAVs. Aerodynamic ...investigations always focus on several static states in the varying sweep process, which ignore the unsteady aerodynamic characteristics. However, deviations to static aerodynamic forces are inevitably caused by dynamic sweep motion. In this work, first, unsteady aerodynamic characteristics on a typical variable-sweep UAV with large aspect ratio were analyzed. Then, deep mechanism of unsteady aerodynamic characteristics in the varying sweep process was studied. Finally, numerical simulation method integrated with structured moving overset grids was applied to solve the unsteady fluid of varying sweep process. The simulation results of a sweep forward-backward circle show a distinct dynamic hysteresis loop surrounding the static data for the aerodynamic forces. Compared with the static lift coefficients , at the same sweep angles, dynamic lift coefficient in sweep forward process are all smaller, while dynamic sweep backward lift coefficient are all larger. In addition, dynamic deviations to static lift coefficient are positively related with the varying sweep speeds. Mechanism study on the unsteady aerodynamic characteristics indicates that three key factors lead to the dynamic hysteresis loop in varying sweep process. They are the effects of additional velocity caused by varying sweep motion, the effects of flow hysteresis and viscosity. The additional velocity induced by sweep motion affects the transversal flow direction along the wing and the effective angle of attack at the airfoil profile. The physical properties of flow, the hysteresis and viscosity affect the unsteady aerodynamic characteristics by flow separation and induced vortexes.
Full text
Available for:
GEOZS, IJS, IMTLJ, KILJ, KISLJ, NLZOH, NUK, OILJ, PNG, SAZU, SBCE, SBJE, UILJ, UL, UM, UPCLJ, UPUK, ZAGLJ, ZRSKP
Hummingbirds are flapping winged creatures with unique flight mechanisms. Their flight pattern is more similar to insects than other birds. Because their flight pattern provides a large lift force at ...a very small scale, hummingbirds can remain hovering while flapping. This feature is of high research value. In order to understand the high-lift mechanism of hummingbirds' wings, in this study a kinematic model is established based on hummingbirds' hovering and flapping process, and wing models imitating the wing of a hummingbird are designed with different aspect ratios. Therefore, with the help of computational fluid dynamics methods, the effect of aspect ratio changes on the aerodynamic characteristics of hummingbirds' hovering and flapping are explored in this study. Through two different quantitative analysis methods, the results of lift coefficient and drag coefficient show completely opposite trends. Therefore, lift-drag ratio is introduced to better evaluate aerodynamic characteristics under different aspect ratios, and it is found that the lift-drag ratio reaches a higher value when AR = 4. A similar conclusion is also reached following research on the power factor, which shows that the biomimetic hummingbird wing with AR = 4 has better aerodynamic characteristics. Furthermore, the study of the pressure nephogram and vortices diagram in the flapping process are examined, leading to elucidation of the effect of aspect ratio on the flow field around hummingbirds' wings and how these effects ultimately lead to changes in the aerodynamic characteristics of the birds' wings.