Search

In this paper, a novel surrogate model for shape-parametrized vehicle drag force prediction is proposed. It is assumed that only a limited dataset of high-fidelity CFD results is available, typically less than ten high-fidelity CFD solutions for different shape samples. The idea is to take advantage not only of the drag coefficients but also physical fields such as velocity, pressure, and kinetic energy evaluated on a cutting plane in the wake of the vehicle and perpendicular to the road. This additional “augmented” information provides a more accurate and robust prediction of the drag force compared to a standard surface response methodology. As a first step, an original reparametrization of the shape based on combination coefficients of shape principal components is proposed, leading to a low-dimensional representation of the shape space. The second step consists in determining principal components of the x-direction momentum flux through a cutting plane behind the car. The final step is to find the mapping between the reduced shape description and the momentum flux formula to achieve an accurate drag estimation. The resulting surrogate model is a space-parameter separated representation with shape principal component coefficients and spatial modes dedicated to drag-force evaluation. The algorithm can deal with shapes of variable mesh by using an optimal transport procedure that interpolates the fields on a shared reference mesh. The Machine Learning algorithm is challenged on a car concept with a three-dimensional shape design space. With only two well-chosen samples, the numerical algorithm is able to return a drag surrogate model with reasonable uniform error over the validation dataset. An incremental learning approach involving additional high-fidelity computations is also proposed. The leading algorithm is shown to improve the model accuracy. The study also shows the sensitivity of the results with respect to the initial experimental design. As feedback, we discuss and suggest what appear to be the correct choices of experimental designs for the best results.

Water–air trans-media unmanned vehicle is a kind of aircraft, which can freely fly in the air, sail in the water and pass through free surface. For trans-media aircrafts, the development process from air–surface to air–underwater and from submarine-launched drive to autonomous drive is investigated. By analyzing the characteristic of manta ray, flying fish and existing aircraft, this paper proposes a new water–air trans-media unmanned vehicle with flat dish-airfoil-shaped main body and telescopic NACA-type wing. Then the numerical method to calculate the lift and drag forces is established and validated by the results of classic NACA cases. On this basis, the flow field around the new vehicle is numerically simulated, and its lift–drag performances in different media (air and water) and ground effect are analyzed, comparing it with a model inspired by the Blackwing Unmanned Aerial Vehicle (UAV). The findings illustrate the superior performance of the new vehicle in terms of lift and drag forces, offering an innovative design framework for water–air trans-media UAV applications.

Assessing drag force and Reynolds stresses in turbulent flows is crucial for evaluating the stability and longevity of hydraulic structures. Yet, this task is challenging due to the complex nature of turbulent flows. To address this, physical models are often employed. Nonetheless, this practice is associated with difficulties, especially in the case of high sampling frequency where the inherent randomness of velocity fluctuations becomes mixed with the measurement noise. This study introduces a stochastic approach, which aims to mitigate bias from measurement errors and provide a probabilistic estimate of extreme stress values. To accomplish this, a simple experimental setup with a hydraulic jump was employed to acquire long-duration velocity measurements. Subsequently, a modified first-order autoregressive model was applied through ensemble simulations, demonstrating the benefits of the stochastic approach. The analysis highlights its effectiveness in estimating the uncertainty of extreme events frequency and minimizing the bias induced by the noise in the high-magnitude velocity measurements and by the limited length of observations. These findings contribute to advancing our understanding of turbulent flow analysis and have implications for the design and assessment of hydraulic structures.

This review examines studies aimed at the organization of energy (non-mechanical) control of high-speed flow/flight using spatially multi-component plasma structures and combined energy deposition. The review covers selected works on the experimental acquisition and numerical modeling of multi-component plasma structures and the use of sets of actuators based on plasma of such a spatial type for the purposes of control of shock wave/bow shock wave–energy source interaction, as well as control of shock wave–boundary layer interaction. A series of works on repetitive multiple laser pulse plasma structures is also analyzed from the point of view of examining shock wave/bow shock wave–boundary layer interaction. Self-sustained theoretical models for laser dual-pulse, multi-mode laser pulses, and self-sustained glow discharge are also considered. Separate sections are devoted to high-speed flow control using combined physical phenomena and numerical prediction of flow control possibilities using thermal longitudinally layered plasma structures. The wide possibilities for organization and applying spatially multi-component structured plasma for the purposes of high-speed flow control are demonstrated.

The force exerted on a cylindrical intruder as it penetrates a granular medium was analyzed utilizing both experiments and the discrete element method (DEM). In this work, a series of penetration experiments were performed, considering cylindrical intruders with different nose shapes. We found that the drag force of the intruder with a hemispherical nose is close to that of those with conical noses with apex angles of 53° and 90°. The drag force of the blunt-nosed intruder is bigger; the drag force of the conical-nosed intruder with an apex angle of 37° is the smallest. We studied the interplay between the drag force on an intruder with a hemispherical nose and key variables—the penetration velocity (V), penetrator’s diameter (di), and friction coefficient (μ). From this analysis, two piecewise functions were derived: one for the average drag force versus the penetration velocity, and the other for the scaled drag force versus the friction coefficient. Furthermore, the average drag force per contact point, Fa/P, can be succinctly represented by two linear relationships: Fa/P = 0.232μ + 0.015(N) for μ<0.9, and Fa/P = 0.225(N) for μ≥0.9.

A numerical simulation of the Herschel–Bulkley laminar steady state shear flow around a stationary particle located on a sedimentation layer was carried out. The surface of the sedimentation layer was formed by hemispheres of the same radius as the particle. The drag force, lift force, and torque values were obtained in the following ranges: shear Reynolds numbers for a particle ReSH=2–200, corresponding to laminar flow; power law index n=0.6–1.0; and Bingham number Bn=0–10. A significant difference in the forces and torque acting on a particle in shear flow in comparison to the case of a smooth wall is shown. It is shown that the drag coefficient is on average 6% higher compared to a smooth wall for a Newtonian fluid but decreases with the increase in non-Newtonian properties. At the edge values of n=0.6 and Bn=10, the drag is on average 25% lower compared to the smooth wall. For a Newtonian fluid, the lift coefficient is on average 30% higher compared to a smooth wall. It also decreases with the increase in non-Newtonian properties of the fluid, but at the edge values of n=0.6 and Bn=10, it is on average only 3% lower compared to the smooth wall. Approximation functions for the drag, lift force, and torque coefficient are constructed. A reduction in the drag force and lifting force leads to an increase in critical stresses (Shields number) on the wall on average by 10% for incipient motion (rolling) and by 12% for particle detachment from the sedimentation bed.

A new mathematical model of an aerosol fibrous filter, composed of a variety of nano- and microfibers, is developed. The combination of nano- and microfibers in a mixed-type filter provides a higher overall quality factor compared with filters with monodisperse fibers. In this paper, we propose a mathematical model of the flow of an incompressible fluid in a porous region consisting of a set of cylinders of various diameters in the range of nano- and micrometers to describe a mixed-type aerosol filter. The flow domain is a rectangular periodic cell with one microfiber and many nanofibers. The motion of the carrier medium is described by the boundary value problem in Stokes flow approximation with the no-slip boundary condition for microfibers and the slip condition for nanofibers. The boundary element method taking into account the slip and non-slip conditions is developed. The calculated velocity field, streamlines, vorticity distribution, and drag of separate fibers and the entire periodic cell are presented. Numerical results for the drag force of the porous medium of a mixed-type filter for the various ratios of mass proportion of nano- and microfibers, porosity, and filtration velocity are presented. The obtained results are compared with the analytical formulas based on the approximate theory of filtration of bimodal filters and with known experimental data. It is shown that with an increase in the mass fraction of nanofibers, the total drag force of the cell increases, while the relative contribution of nanofibers to the total drag force tends toward the value that is less than unity. An approximate analytical formula for the drag coefficient of a mixed aerosol filter is derived. The developed flow model and analytical formulas allow for estimating the aerodynamic drag of a mixed filter composed by nano- and microfibers.

Traditional liquid phase methods for growing single crystals are static growth methods, which include seed crystal sedimentation growth and seed crystal clamping growth using seed crystal holders. Single crystals grown via seed crystal sedimentation often have a flat and elongated shape, and the region in contact with the bottom of the container is restricted during growth, resulting in significant defects. Similarly, the seed crystal clamping growth method cannot avoid contact with external objects, leading to abnormal growth at the contact points and along the direction of the seed crystal holder, also resulting in certain defects. Both of these growth methods require processes, such as cutting and grinding, to remove defects, resulting in resource waste. To address the shortcomings of the static liquid phase single-crystal preparation mentioned above, this study successfully designed a dynamic liquid phase method for single crystal growth, which achieved the successful suspension of seed crystals in the mother solution and the growth of high-quality, large-sized single crystals, avoiding contact with the walls and the bottom of the container during the crystal growth process. Based on the dynamic liquid phase single crystal growth apparatus mentioned above, stable and dynamic liquid phase preparation was successfully achieved, ranging from seed crystals with a diameter of approximately 5 mm to single crystals with a diameter of approximately 20 mm, by controlling the cooling rate and adjusting the solution flow rate.

Evaluating the hydrodynamic force fluctuations acting on each structure in a group of subsea objects of different cross-section shapes, sizes and relative positions represents a challenge due to the sensitivity of the vortex shedding process, especially for a variety of sheared flows. The present study uses the numerical 2D computational fluid dynamics model to estimate the flow-induced forces on a group of small circular and D-shaped cylinders in the linear and parabolic sheared flow, which are placed in proximity to a larger structure of the squared cross-section. This allows us to evaluate loads, which are affected by the presence of subsea equipment located on the seabed. The average Reynolds number of the considered linear flow profile is 3900, while the parabolic flow profile has the maximum Reynolds number of 3900. The k-ω SST turbulence model is used for simulations. The work demonstrates the effect of the cross-sectional shape of smaller cylinders on hydrodynamic coefficients, explores the effect from the spacing in between the structures and highlights differences between loads in the linearly sheared and parabolic flow. The results obtained show that the presence of the squared cylinder notably influences the mean drag coefficient on the first cylinder, for both circular and D-shaped cylinders. The parabolic sheared flow profile in this series leads to the highest mean drag and the highest amplitudes of the fluctuating drag and lift coefficients.

The calculation of the drag force is a fundamental requirement in the design of any submerged system intended for marine exploration. The calculation can be performed by analytic analysis, numerical modeling, or by a direct calculation performed in a designated testing facility. However, for complex structures and especially those with a non-rigid design, the analytic and numerical analyses are not sufficiently accurate, while the direct calculation is a costly operation. In this paper, we propose a simple approach for how to calculate the drag coefficient in-situ. Aimed specifically at the complex case of elastic objects whose modeling via Computer-Aided Design (CAD) is challenging, our approach evaluates the relation between the object’s speed at steady-state and its mass to extract the drag coefficient in any desired direction, the hydro-static force, and, when relevant, also the thruster’s force. We demonstrate our approach for the special case of a highly complex elastic-shaped floater that profiles the water column. The analysis of two such floaters in two different sea environments shows accurate evaluation results and supports our claim for robustness. In particular, the simplicity of the approach makes it appealing for any arbitrary shaped object.

The bullet head plays a principal role in the modern enlargement of an efficient bullet. A bullet’s main design parameters depend upon the lift and drag forces acting on the head. The factors in a bullet’s shape design that affect bullets’ lift and drag forces are essential in aerodynamics, especially in ballistics. Therefore, the effect of wind on the lift and drag forces acting on the bullet, and the role of the bullet head to allow the bullet to travel efficiently through the wind, need to be investigated. This work discusses the parameters that affect the lift and drag force on the bullet. Simulations are performed in Ansys Fluent by varying the key parameters of the bullet head, i.e., the length and angle of attack, while keeping the air velocity at 5.2 m/s. The simulation outcome shows that the size of the bullet and the angle of attack are important factors related to the drag force. Therefore, this work predicts the inspection of a bullet under distinct wind conditions. An evaluation is performed to scrutinize the effect of design factors on the system execution of the bullet and its constructive flight path. It is concluded that when increasing the length of the bullet and its angle of attack (AOA), the drag force and lift forces increase drastically, contributing to the inefficiency of the bullet’s accuracy and penetrating power. A new design is also proposed in which the drag forces are reduced to the minimum.

The rapid electrokinetic patterning (REP) technique has been demonstrated to enable dynamic particle manipulation in biomedical applications. Previous studies on REP have generally considered particles with a size less than 5 μm. In this study, a REP platform was used to manipulate polystyrene particles with a size of 3~11 μm in a microfluidic channel sandwiched between two ITO conductive glass plates. The effects of the synergy force produced by the REP electrothermal vortex on the particle motion were investigated numerically for fixed values of the laser power, AC driving voltage, and AC driving frequency, respectively. The simulation results showed that the particles were subject to a competition effect between the drag force produced by the toroidal vortex, which prompted the particles to recirculate in the bulk flow adjacent to the laser illumination spot on the lower electrode, and the trapping force produced by the particle and electrode interactions, which prompted the particles to aggregate in clusters on the surface of the illuminated spot. The experimental results showed that as the laser power increased, the toroidal flow range over which the particles circulated in the bulk flow increased, while the cluster range over which the particles were trapped on the electrode surface reduced. The results additionally showed that the particle velocity increased with an increasing laser power, particularly for particles with a smaller size. The excitation frequency at which the particles were trapped on the illuminated hot-spot reduced as the particle size increased. The force and velocity of polystyrene particles by the REP toroidal vortex has implications for further investigating the motion behavior at the biological cell level.

In the simulation of dilute gas-solid flows such as those seen in many industrial applications, the Lagrangian Particle Tracking method is used to track packets of individual particles through a converged fluid field. In the tracking of these particles, the most dominant forces acting upon the particles are those of gravity and drag. In order to accurately predict particle motion, the determination of the aforementioned forces become of the upmost importance, and hence an improved drag force formula was developed to incorporate the effects of particle concentration and particle Reynolds number. The present CFD study examines the individual effects of particles located both perpendicular and parallel to the flow direction, as well as the effect of a particle entrain within an infinite matrix of evenly distributed particles. Results show that neighbouring particles perpendicular to the flow (Model 2) have an effect of increasing the drag force at close separation distances, but this becomes negligible between 5–10 particle diameters depending on particle Reynolds number (Rep). When entrained in an infinite line of particles co-aligned with the flow (Model 1), the drag force is remarkably reduced at close separation distances and increases as the distance increases. The results of the infinite matrix of particles (Model 3) show that, although not apparent in the individual model, the effect of side particles is experienced many particle diameters downstream.

Entrapped air in the chambers of a T-girder bridge generates a considerably large vertical force which is harmful to the safety of the bridge superstructure. The contribution of the entrapped air is not considered when calculating flood force in the related specifications, and also it has not been taken into account by researchers. Two-dimensional, scaled-down models at a 1:40 scale are selected as the research object. Analytical and numerical methods are employed to investigate the contribution of entrapped air in T-girder chambers to flood forces. Results show that the entrapped air in the rear chambers could escape easier than that in the front chambers and the compressibility of entrapped air has a small influence on drag force and vertical force (generated by dynamic pressure) coefficients, but it reduces the buoyancy of the T-girder. Considering the entrapped air compressibility and entrapped air escape, the calculation method of the buoyancy of entrapped air is proposed to improve the precision of the existing flood force calculation method.

Masonry arch bridges form an essential part of existing transport infrastructure around the world, including mainland Europe and the northeastern US. Recent extreme flood events highlight that masonry arch bridges spanning watercourses are vulnerable to flood-induced hydrodynamic and debris impact loads. When the flow interacts with the bridge superstructure, with or without discrete floating debris, a complex interaction is observed. This paper presents both experimental and numerical studies to investigate this complex phenomenon, including fluid–structure and structure–structure interactions. A typical single-span masonry arch bridge and large woody debris representing a tree log are investigated. Experimental observations from a scaled hydraulic model, with and without debris in the flow, are first presented for the case where the abutment of the bridge is fully submerged. Next, the capability of the numerical method smoothed particle hydrodynamics (SPH) in simulating the hydrodynamic behaviour and debris impact observed in the experiment is discussed. Following this, both hydrodynamic and debris-induced pressure–time histories on the bridge are obtained using the SPH model. Results reveal that the debris impact leads to a significantly more localised load on the bridge compared to the situation with hydrodynamic load only.

Micromotors have spread widely in order to meet the needs of new applications, including cell operation, drug delivery, biosensing, precise surgery and environmental decontamination, due to their small size, low energy consumption and large propelling power, especially the newly designed multifunctional micromotors that combine many extra shape features in one device. Features such as rod-like receptors, dendritic biosensors and ball-like catalyzing enzymes are added to the outer surface of the tubular micromotor during fabrication to perform their special mission. However, the structural optimization of motion performance is still unclear. The main factor restricting the motion performance of the micromotors is the drag forces. The complex geometry of a micromotor makes its dynamic behavior more complicated in a fluid environment. This study aimed to design the optimum structure of tubular micromotors with minimum drag forces and obtain the magnitude of drag forces considering both the internal and external fluids of the micromotors. By using the computational fluid dynamics software Fluent 18.0 (ANSYS), the drag force and the drag coefficient of different conical micromotors were calculated. Moreover, the influence of the Reynolds numbers Re, the semi-cone angle δ and the ratios ξ and η on the drag coefficient was analyzed. The results show the drag force monotonically increased with Reynolds numbers Re and the ratio η. The extreme point of the drag curve is reached when the semi-cone angle δ is 8° and the ratio ξ is 3.846. This work provides theoretical support and guidance for optimizing the design and development of conical micromotors.

An understanding of fish movement behavior in response to flow field variables is important for exploring the hydrodynamic strategies of fish in fish passages. In this paper, bighead carps were taken as an example. The fish movement behavior response to water flow field information by means of estimating the energetic expenditure using an IBM approach in an experimental fishway was investigated. Fish swimming velocity, drag force, and energy expenditure were analyzed in varied flow conditions related to hydraulic variables, including velocity (V), turbulent kinetic energy (TKE), and strain rate (SR). The result indicated that the fish will require more energy in high TKE zones. This study provides a reference for optimizing the design of fish passages and fisheries management. This method can be applied to assess the efficiency of fish bypass structures and conduct fish survival studies.

Numerical simulations of contaminated spherical drops falling through a stagnant liquid at low Reynolds numbers are carried out using the finite difference method. The numerical results are used to describe the behavior of the surfactant concentrations and to understand the surfactant effects on the fluid motions in detail. The predicted interfacial surfactant concentration, Γ, is almost zero for angles, θ, below a certain value (the stagnant-cap angle, θcap), whereas it steeply increases and reaches a large value for θ>θcap (the stagnant-cap region). The increase in the initial surfactant concentration, C0, in the drop enhances the adsorption from the drop to the interface, which results in the increase in Γ and the decrease in θcap. Peaks appear in the predicted Marangoni stresses around θcap, which causes similar peaks in the pressure distribution. The high-pressure spots prevent the fluid motion along the interface, which results in the formation of the stagnant-cap region and the attenuation of the tangential velocity in the continuous phase. The surfactant flux from the bulk to the interface decreases C in the vicinity of the interface for θ<θcap and the weak diffusion cannot compensate for the reduction in C by adsorption, which results in C at the interface smaller than C0. The pattern of the low C region is determined by the advection and does not smear out because of a small diffusive flux.

Experimental tests regarding the M-346 aircraft model made via 3D printing were carried out in order to obtain numerical data and characteristics in the form of graphs of basic aerodynamic forces and coefficients. The tests were carried out for the left side of the airframe model in a clean configuration, without additional suspension equipment; the flight control surfaces and the aerodynamic brake were in neutral positions. Based on the scan of the base model in 1:48 scale using a Nikon Model Maker MMDx laser scanning head, followed by the generation and optimization of some of the airframe elements in SolidWorks software, a test model ready for printing was prepared. Using the MakerBot Print program, the printing parameters were set, and the process itself was completed using a MakerBot Replicator Z18 3D printer. The next step was manual treatment in order to remove the material excess from the melted thermoplastic material, join the elements and appropriately polish the surface of the tested model in order to obtain the desired quality. The test was carried out using a Gunt HM 170 wind tunnel for fixed airflow velocities at variable angles of attack. On this basis, the numerical values of lift force, Pz, and drag force, Px, were obtained; then, the lift force indices, Cz, and drag force indices, Cx, were computed for the steady states, which were for angle α from −12 to 16°. The use of 3D printing contributed to the generation of geometry, which, for research purposes, was scaled down in order to fully use the available measurement space of the wind tunnel. The final stage of the work was to compare the obtained curves of particular characteristics with the literature data.

Offshore installations with pile foundations in shallow water are vulnerable to submarine landslides, which cause serious damage to engineering facilities, loss of life, and loss of money. Due to a shortage of real observation data and the difficulty of reproduction, we lack insight into the interaction behavior between submarine landslides and monopiles. This study capitalized on ANSYS Fluent 20.0 to develop a three-dimensional biphasic (water and slurry) numerical model. This CFD model was used to analyze the interaction between a monopile and submarine landslides at different flow heights. The velocities of submarine landslides were from low to high values. Two modes of interactional forces acting on the monopile are proposed, which are (i) interaction force with peak value and (ii) interaction force without peak value. The influence of flow height and velocity on interaction forces was investigated. Results show that the effect of the flow heights on the interaction force is significant at low velocity stage, while the peak force representing a hazard level of the pile was non-negligible under high flow velocity and low flow height conditions, which should be considered in a future study. The related mechanisms are revealed with a hybrid model considering different components of the force.

Lift force is an important parameter for the performance evaluation of an H-Darrieus wind turbine. The rotational direction of the streamlined force is effective on the performance of the wind turbine. In order to analyze the flow characteristics around the turbine blades in real-time, a numerical analysis using three-dimensional unsteady Reynold-averaged Navier–Stokes equations has been introduced. Experimental data were obtained from a field test facility constructed on an island in South Korea and was introduced to compare the numerical simulation results with measured data. The optimum tip speed ratio (TSR) was investigated via a multi-variable optimization approach and was determined to be 3.5 for the NACA 0015 blade profile. The turbine displays better performance with the maximum power coefficient at the optimum TSR. It is due to the delay in the flow separation from the blade surface and the relatively lower strength of the tip vortices. Furthermore, the ratio between lift and drag forces is also the highest at the optimum TSR, as most of the aerodynamic force is directly converted into lift force. For one rotation of the turbine blade at the optimum TSR, the first quarter of motion produces the highest lift as the static pressure difference is maximum at the leading edge, which helps to generate maximum lift. At a TSR less than the optimum TSR, small-lift generation is dominant, whereas at a higher TSR, large drag production is observed. Both of these lead to lower performance of the turbine. Apart from the TSR, the optimum wind angle of attack is also investigated, and the results are prepared against each TSR.

Gas bubbles are of major importance in most metallurgical processes. They promote chemical reactions, homogenize the melt, or float inclusions. Thus, their dynamics are of crucial interest for the optimization of metallurgical processes. In this work, the state of knowledge of bubble dynamics at the bubble scale in liquid metals is reviewed. Measurement methods, with emphasis on liquid metals, are presented, and difficulties and shortcomings are analyzed. The bubble formation mechanism at nozzles and purging plugs is discussed. The uncertainty regarding the prediction of the bubble size distribution in real processes is demonstrated using the example of the steel casting ladle. Finally, the state of knowledge on bubble deformation and interfacial forces is summarized and the scalability of existing correlations to liquid metals is critically discussed. It is shown that the dynamics of bubbles, especially in liquid metals, are far from understood. While the drag force can be predicted reasonably well, there are large uncertainties regarding the bubble size distribution, deformation, and lift force. In particular, the influence of contaminants, which cannot yet be quantified in real processes, complicates the discussion and the comparability of experimental measurements. Further open questions are discussed and possible solutions are proposed.

This article presents the results of numerical investigation of the influence of the streamlined body temperature on drag force and on the falling velocity in a water channel. The experimental data reflecting the cooling dynamics and body temperature influence on the falling velocity are presented as well. k − ε turbulence model and homogenous heat transfer model were chosen for the numerical 3D simulation. Drag force changes induced by the alteration of the body temperature were investigated. Velocity of the streamlined body under different temperatures of water was investigated experimentally, and the results were compared to the data obtained during the numerical simulation. The increase of the falling velocity and decrease of drag force were found to have been affected by the increase of the body temperature, which had influence on the change of the water parameters (density, phase, etc.) near the surface of the body. Simulation showed that the drag force and a velocity also depended on the water temperature. The drag force of the streamlined body decreased by 32% in comparison to the cold body for the body temperature equal to 150 °C and water temperature close to the saturation temperature (98 °C). Experimentally, it was determined that the velocity of the streamlined body covered by vapor film depended on the falling time and increased by 10–30%. Velocity difference was very small for the cold and hot bodies at the initial moment of the drop; however, it reached 20% and more after 0.3 s of the falling process.

The interaction between a shock wave and solid particles involves complex gas–solid two-phase flow, which is widely used in industrial processes. Theoretical analysis, an experimental test, and simulation were combined to investigate the interaction process between a shock wave and quartz sand particles. The variation of physical parameters of the two phases during the interaction process was considered theoretically. Then, a novel vertical shock tube generator was employed to record the pressure attenuation and dispersion process of solid particles. Finally, the complex gas–solid two-phase flow was simulated based on the computational fluid dynamics method. The results showed that a nonequilibrium state was formed during the interaction process and momentum exchange generated, resulting in a drag force of the shock wave on the particles. The shock intensity obviously attenuated after the shock wave passed through the solid particles, and this part of the energy was work on the solid particles to drive their dispersion. A three-dimensional annular vortex was generated around the solid particles due to the entrainment effect of airflow. Under the shock wave action of 1.47 Ma, the three types of solid particles with average diameters of 2.5, 0.95, and 0.42 mm presented different motion laws. The particles with smaller size were easier to disperse, and the cloud that formed was larger and more uniform.

CFD-DEM (computational fluid dynamic-discrete element method) is a promising approach for simulating fluid–solid flows in fluidized beds. This approach generally under-predicts the granular temperature due to the use of drag models for the average drag force. This work develops a simple model to improve the granular temperature through increasing the drag force fluctuations on the particles. The increased drag force fluctuations are designed to match those obtained from PR-DNSs (particle-resolved direct numerical simulations). The impacts of the present model on the granular temperatures are demonstrated by posteriori tests. The posteriori tests of tri-periodic gas–solid flows show that simulations with the present model can obtain transient as well as steady-state granular temperature correctly. Moreover, the posteriori tests of fluidized beds indicated that the present model could significantly improve the granular temperature for the homogenous or slightly inhomogeneous systems, while it showed negligible improvement on the granular temperature for the significantly inhomogeneous systems.

Herein, we compare the drag area estimated using unsteady Reynolds-averaged Navier-Stokes (URANS), using the γ−ReΘ transitional shear stress transport (SST) k−ω (SSTLM) turbulence model with published experimental measurements of a static full-scale cyclist mannequin in an open test section wind tunnel, with the left leg fully extended. The turbulence model employs a local empirical correlation based upon a classical Blasius boundary layer behavior to predict flow transition. For a given mesh density, we aim to improve drag area estimation by modifying the empirical correlation coefficient to better capture actual boundary layer transition location around the arms and legs, to facilitate computationally economical cyclist simulations. Large Eddy Simulation (LES), in conjunction with experimental wake data in the vicinity of the arms and legs, is used to assess boundary layer shape factors, which are related to the empirical coefficient. Overall, the drag area predicted by LES is within 3.7% of the measured results, while the original SSTLM is within 7.8%. By tuning the correlation coefficient, the drag area error is improved to 6.0% at no additional computational cost. The tuning was relatively coarse, and was only considered for the appendages. In other regions, the original SSTLM coefficient seems to perform better, suggesting that local coefficient selection may lead to further improvements in results over the currently employed global value.

In Japan, the marine aquaculture net cage has an important role in farming pacific bluefin tuna farming in oceans, and the design of the net cage needs to ensure robustness against hostile oceanic conditions. Accordingly, this study focuses on the drag forces and the cage volume of the net cage, and on their variations induced by different design parameters (netting solidity ratio, netting height, and bottom weight). A series of parametric studies on drag force and deformation of the net cage was conducted using a numerical simulation model. Accordingly, the contribution of each parameter to the drag and volume was analyzed using a generalized additive model. The results indicate that the bottom weight had the highest contribution to the holding ratio of the cage volume, whereas the netting height had the highest contribution to the drag coefficient of the net cage. Finally, a fast prediction model was created by a backpropagation (BP) neural network model and was examined for the accurate prediction of the objective variables.

In the context of dam breaks, tsunami, and flash floods, it is paramount to quantify the time-history of forces by the rapidly transient flow to vertical structures and the characteristics of the induced flow patterns. To resemble on-land tsunami-induced flow, a free-surface-piercing structure is exposed to long leading depression waves in a tsunami flume where long waves run up and down a 1:40 smooth and impermeable sloping beach after its generation by a volume-driven wave maker. The structure and its surrounding were monitored with force transducers, pressure gauges and cameras. Preparatory steady-state experiments were accomplished to determine the drag force coefficient of the square cylinder at various water depths. The flow during wave run-up and draw-down acting on the structure resulted in distinct flow pattern which were characteristic for the type of flow-structure interaction. Besides bow wave propagating upstream, a standing or partially-standing wave was observed in front of the structure together with a wake formation downstream, while a von Kármán vortex street developed during the deceleration phase of the flow motion and during draw-down. Force measurements indicated a sudden increase in the stream-wise total force starting with the arrival of the flow front during initial run-up. Lateral velocities showed significant oscillations in correlation with the von Kármán vortex street development. A comparison of the total measured base force with the analytically-calculated share of the drag force revealed that forces were prevailingly drag-dominated.

In underwater towing operations, the drag force and vertical offset angle of towropes are important considerations when choosing and setting up towing equipment. The aim of this paper is to study the variation in drag force, vertical offset angle, resistance, and attitude for towing operations with a view to optimizing these operations. An underwater experiment was conducted using a 1:8 scale physical model of a subsea module. A comprehensive series of viscous Computational Fluid Dynamics (CFD) simulations were carried out based on Reynolds-averaged Navier−Stokes equations for uniform velocity towing. The results of the simulation were compared with experimental data and showed good agreement. Numerical results of the vorticity field and streamlines at the towing speeds were presented to analyze the distribution of vortexes and flow patterns. The resistance components were analyzed based on the numerical result. It was found that the lateral direction was a better direction for towing operations because of the smaller drag force, resistance, and offset angle. Similar patterns and locations of streamlines and vortexes were present in both the longitudinal and lateral directions, the total resistance coefficient decreases at a Reynolds number greater than that of a cylinder.

Adequate design of pipelines used for oil, gas, water, and wastewater transmission is essential not only for their proper operation but particularly to avoid failure and the possible extreme consequences. This is even more drastic in nearshore environments, where pipelines are potentially exposed to extreme hydrodynamic events, such as tsunami- or storm-surge-induced inundation. The American Society of Civil Engineers (ASCE), in its ASCE7 Chapter 6 on Tsunami Loads and Effects which is the new standard for tsunami impacts and loading, specifically stresses the need to study loads on pipelines located in tsunami-prone areas. To address this issue, this study is the first of its kind to investigate loading on pipelines due to tsunami-like bores. A comprehensive program of physical model experiments was conducted in the Dam-Break Hydraulic Flume at the University of Ottawa, Canada. The tests simulated on-land tsunami flow inundation propagating over a coastal plain. This allowed to record and investigate the hydrodynamic forces exerted on the pipe due to the tsunami-like, dam-break waves. Different pipe configurations, as well as various flow conditions, were tested to investigate their influence on exerted forces and moments. The goal of this study was to propose, based on the results of this study, resistance and lift coefficients which could be used for the design of pipelines located in tsunami-prone areas. The values of the resistance and lift coefficients investigated were found to be in the range of 1 < C R < 3.5 and 0.5 ≤ C L < 3 , respectively. To that end, the study provides an upper envelope of resistance and lift coefficients over a wide range of Froude numbers for design purposes.

A numerical simulation of the laminar flow field and convection⁻diffusion mass transfer in a regular system of parallel fully absorbing fibers for the range of Reynolds numbers up to Re = 300 is performed. An isolated row of equidistant circular fibers arranged normally to the external flow is considered as the simplest model for a hollow-fiber membrane contactor. The drag forces acting on the fibers with dependence on Re and on the ratio of the fiber diameter to the distance between the fiber axes, as well as the fiber Sherwood number versus Re and the Schmidt number, Sc, are calculated. A nonlinear regression formula is proposed for calculating the fiber drag force versus Re in a wide range of the interfiber distances. It is shown that the Natanson formula for the fiber Sherwood number as a function of the fiber drag force, Re, and Sc, which was originally derived in the limit of high Peclet numbers, is applicable for small and intermediate Reynolds numbers; intermediate and large Peclet numbers, where Pe = Re × Sc; and for sparse and moderately dense rows of fibers.

In past decades, considerable advances have been achieved in micro and nanomotors. Particular attention has been given to self-propelled catalytic micromotors, which have been widely used in cell separation, drug delivery, microsurgery, lithography and environmental remediation. Fast moving, long life micromotors appear regularly, however it seems there are no solutions yet that thoroughly clarify the hydrodynamic behavior of catalytic micromotors moving in fluid. Dynamic behavior of this kind of micromotors is mainly determined by the driving force and drag force acting on the micromotors. Based on the hydromechanics theory, a hydrodynamic model is established to predict the drag force for a conical micromotor immersed in the flow field. By using the computational fluid dynamics software Fluent 18.0 (ANSYS), the drag force and the drag coefficient of different conical micromotors are calculated. A mathematical model was proposed to describe the relationship among Reynolds numbers Re, the ratio λ, the semi-cone angle δ and the drag coefficient Cd of the micromotors. This work provides theoretical support and reference for optimizing the design and development of conical micromotors.