Abstract: When a partially or completely immersed solid moves in a liquid such as water, it undergoes a force called hydrodynamic drag. Reducing this force has always been the objective of hydrodynamic engineers to make water slide better on submerged bodies. This paper deals with the examination of the different terms composing the analytical solution of the flow over an obstacle embedded at the bottom of a hydraulic channel. We have chosen to use a linear method to study a two-dimensional flow over an obstacle, in order to understand the evolution of the drag. We set the following assumptions: incompressible inviscid fluid, irrotational flow, low obstacle height compared to the water height. Those assumptions allow overcoming the difficulties associated with modelling these waves. We will mathematically formulate the equations that allow the determination of the stream function, and then the free surface equation. A similar method is used to determine the exact analytical solution for an obstacle in the shape of a sinusoidal arch.
Abstract: The earliest theories of sloshing waves and solitary waves based on potential theory idealisations and irrotational flow have been extended to be applicable to more realistic domains. To this end, the computational fluid dynamics (CFD) methods are widely used. Three-dimensional CFD methods such as Navier-Stokes solvers with volume of fluid treatment of the free surface and Navier-Stokes solvers with mappings of the free surface inherently impose high computational expense; therefore, considerable effort has gone into developing depth-averaged approaches. Examples of such approaches include Green–Naghdi (GN) equations. In Cartesian system, GN velocity profile depends on horizontal directions, x-direction and y-direction. The effect of vertical direction (z-direction) is also taken into consideration by applying weighting function in approximation. GN theory considers the effect of vertical acceleration and the consequent non-hydrostatic pressure. Moreover, in GN theory, the flow is rotational. The present study illustrates the application of GN equations to propagation of sloshing waves and solitary waves. For this purpose, GN equations solver is verified for the benchmark tests of Gaussian hump sloshing and solitary wave propagation in shallow basins. Analysis of the free surface sloshing of even harmonic components of an initial Gaussian hump demonstrates that the GN model gives predictions in satisfactory agreement with the linear analytical solutions. Discrepancies between the GN predictions and the linear analytical solutions arise from the effect of wave nonlinearities arising from the wave amplitude itself and wave-wave interactions. Numerically predicted solitary wave propagation indicates that the GN model produces simulations in good agreement with the analytical solution of the linearised wave theory. Comparison between the GN model numerical prediction and the result from perturbation analysis confirms that nonlinear interaction between solitary wave and a solid wall is satisfactorilly modelled. Moreover, solitary wave propagation at an angle to the x-axis and the interaction of solitary waves with each other are conducted to validate the developed model.
Abstract: One of the most perspective methods to produce SoG-Si is refinement via metallurgical route. The most critical part of this route is refinement from boron and phosphorus. Therefore, a new approach could address this problem. We propose an approach of creating surface waves on silicon melt’s surface in order to enlarge its area and accelerate removal of boron via chemical reactions and evaporation of phosphorus. A two dimensional numerical model is created which includes coupling of electromagnetic and fluid dynamic simulations with free surface dynamics. First results show behaviour similar to experimental results from literature.
Abstract: Local irregular topography has a great impact on earthquake ground motion. For scarp topography, using numerical simulation method, the influence extent and scope of the scarp terrain on scarp's upside and downside ground motion are discussed in case of different vertical incident SV waves. The results show that: (1) The amplification factor of scarp's upside region is greater than that of the free surface, while the amplification factor of scarp's downside part is less than that of the free surface; (2) When the slope angle increases, for x component, amplification factors of the scarp upside also increase, while the downside part decrease with it. For z component, both of the upside and downside amplification factors will increase; (3) When the slope angle changes, the influence scope of scarp's downside part is almost unchanged, but for the upside part, it slightly becomes greater with the increase of slope angle; (4) Due to the existence of the scarp, the z component ground motion appears at the surface. Its amplification factor increases for larger slope angle, and the peaks of the surface responses are related with incident waves. However, the input wave has little effects on the x component amplification factors.
Abstract: As a well-known conventional finishing process, the grinding is commonly used to manufacture seal mating surfaces and bearing surfaces, but is also creates twisted surfaces. The machined surfaces by turning or grinding usually have twist structure on the surfaces, which can convey lubricants such as conveyor screw. To avoid this phenomenon, have to use special techniques or machines, for example start-stop turning, tangential turning, ultrasonic protection or special toll geometries. All of these solutions have high cost and difficult usability. In this paper, we describe a system and summarize the results of the experimental research carried out mainly in the field of Magnetic Abrasive Polishing (MAP) and Magnetic Roller Burnishing (MRB). These technologies are simple and also green while able to produce twist-free surfaces. During the tests, C45 normalized steel was used as workpiece material which was machined by simple and Wiper geometrical turning inserts in a CNC turning lathe. After the turning, the MAP and MRB technologies can be used directly to reduce the twist of surfaces. The evaluation was completed by advanced measuring and IT equipment.
Abstract: Reactive rotational molding (RRM) is a process to manufacture hollow plastic parts with reactive material has several advantages compared to conventional roto molding of thermoplastic powders: process cycle time is shorter; raw material is less expensive because polymerization occurs during processing and high-performance polymers may be used such as thermosets, thermoplastics or blends. However, several phenomena occur during this process which makes the optimization of the process quite complex. In this study, we have used a mixture of isocyanate and polyol as a reactive system. The chemical transformation of this system to polyurethane has been studied by thermal analysis and rheology tests. Thanks to these results of the curing process and rheological measurements, the kinetic and rheokinetik of polyurethane was identified. Smoothed Particle Hydrodynamics, a Lagrangian meshless method, was chosen to simulate reactive fluid flow in 2 and 3D configurations of the polyurethane during the process taking into account the chemical, and chemiorehological results obtained experimentally in this study.
Abstract: The present study addresses the problem of ammonia evaporation during filling of a vertical cylindrical tank and the influence of various external factors on the stability of storage by determining the conditions for minimum evaporation. Numerical simulation is carried out by solving the governing equations namely, continuity, momentum, energy, and diffusion of species. The effect of temperature of surrounding air, the filling speed of the reservoir and the temperature of the filling liquid ammonia on the evaporation rate is investigated. Results show that the temperature of the filling liquid has little effect on the liquid ammonia for a short period, which, in fact, is function of the filling speed. The evaporation rate along the free surface of the liquid is non-uniform. The inlet temperature affects the vapor ammonia temperature because of pressure increase. The temperature of the surrounding air affects the temperature of the vapor phase rather than the liquid phase. The maximum of evaporation is reached at the final step of filling. In order to minimize loss of ammonia vapors automatically causing losses in quantity of the liquid stored, it is suggested to ensure the proper insulation for the walls and roof of the reservoir and to increase the filling speed.
Abstract: The damage caused by surge waves generated in water
bodies by flow-like landslides can be very high in terms of human
lives and economic losses. The complicated phenomena occurred in
this highly unsteady process are difficult to model because three
interacting phases: air, water and sediment are involved. The problem
therefore is challenging since the effects of non-Newtonian fluid
describing the rheology of the flow-like landslides, multi-phase flow
and free surface have to be included in the simulation. In this work, the
commercial computational fluid dynamics (CFD) package FLUENT is
used to model the surge waves due to flow-like landslides. The
comparison between the numerical results and experimental data
reported in the literature confirms the accuracy of the method.
Abstract: The aim of this study is to examine, through
experimentation in the laboratory, the supercritical flow in the
presence of an obstacle in a rectangular channel. The supercritical
regime in the whole hydraulic channel is achieved by adding a
convergent. We will observe the influence of the obstacle shape and
dimension on the characteristics of the supercritical flow, mainly the
free-surface elevation and the velocity profile. The velocity
measurements have been conducted with the one dimension laser
anemometry technique.
Abstract: The present study is concerned with the problem of determining the shape of the free surface flow in a hydraulic channel which has an uneven bottom. For the mathematical formulation of the problem, the fluid of the two-dimensional irrotational steady flow in water is assumed inviscid and incompressible. The solutions of the nonlinear problem are obtained by using the usual conformal mapping theory and Hilbert’s technique. An experimental study, for comparing the obtained results, has been conducted in a hydraulic channel (subcritical regime and supercritical regime).
Abstract: Steady three-dimensional and two free surface waves
generated by moving bodies are presented, the flow problem to be
simulated is rich in complexity and poses many modeling challenges
because of the existence of breaking waves around the ship hull, and
because of the interaction of the two-phase flow with the turbulent
boundary layer. The results of several simulations are reported. The
first study was performed for NACA0012 of hydrofoil with different
meshes, this section is analyzed at h/c= 1, 0345 for 2D. In the second
simulation a mathematically defined Wigley hull form is used to
investigate the application of a commercial CFD code in prediction of
the total resistance and its components from tangential and normal
forces on the hull wetted surface. The computed resistance and wave
profiles are used to estimate the coefficient of the total resistance for
Wigley hull advancing in calm water under steady conditions. The
commercial CFD software FLUENT version 12 is used for the
computations in the present study. The calculated grid is established
using the code computer GAMBIT 2.3.26. The shear stress k-ωSST
model is used for turbulence modeling and the volume of fluid
technique is employed to simulate the free-surface motion. The
second order upwind scheme is used for discretizing the convection
terms in the momentum transport equations, the Modified HRIC
scheme for VOF discretization. The results obtained compare well
with the experimental data.
Abstract: Steady three-dimensional and two free surface waves
generated by moving bodies are presented, the flow problem to be
simulated is rich in complexity and poses many modeling challenges
because of the existence of breaking waves around the ship hull, and
because of the interaction of the two-phase flow with the turbulent
boundary layer. The results of several simulations are reported. The
first study was performed for NACA0012 of hydrofoil with different
meshes, this section is analyzed at h/c= 1, 0345 for 2D. In the second
simulation a mathematically defined Wigley hull form is used to
investigate the application of a commercial CFD code in prediction of
the total resistance and its components from tangential and normal
forces on the hull wetted surface. The computed resistance and wave
profiles are used to estimate the coefficient of the total resistance for
Wigley hull advancing in calm water under steady conditions. The
commercial CFD software FLUENT version 12 is used for the
computations in the present study. The calculated grid is established
using the code computer GAMBIT 2.3.26. The shear stress k-ωSST
model is used for turbulence modeling and the volume of fluid
technique is employed to simulate the free-surface motion. The
second order upwind scheme is used for discretizing the convection
terms in the momentum transport equations, the Modified HRIC
scheme for VOF discretization. The results obtained compare well
with the experimental data.
Abstract: Water spray cooling is a technique typically used in
heat treatment and other metallurgical processes where controlled
temperature regimes are required. Water spray cooling is used in
static (without movement) or dynamic (with movement of the steel
plate) regimes. The static regime is notable for the fixed position of
the hot steel plate and fixed spray nozzle. This regime is typical for
quenching systems focused on heat treatment of the steel plate. The
second application of spray cooling is the dynamic regime. The
dynamic regime is notable for its static section cooling system and
moving steel plate. This regime is used in rolling and finishing mills.
The fixed position of cooling sections with nozzles and the
movement of the steel plate produce nonhomogeneous water
distribution on the steel plate. The length of cooling sections and
placement of water nozzles in combination with the nonhomogeneity
of water distribution lead to discontinued or interrupted cooling
conditions. The impact of static and dynamic regimes on cooling
intensity and the heat transfer coefficient during the cooling process
of steel plates is an important issue.
Heat treatment of steel is accompanied by oxide scale growth. The
oxide scale layers can significantly modify the cooling properties and
intensity during the cooling. The combination of static and dynamic
(section) regimes with the variable thickness of the oxide scale layer
on the steel surface impact the final cooling intensity. The study of
the influence of the oxide scale layers with different cooling regimes
was carried out using experimental measurements and numerical
analysis. The experimental measurements compared both types of
cooling regimes and the cooling of scale-free surfaces and oxidized
surfaces. A numerical analysis was prepared to simulate the cooling
process with different conditions of the section and samples with
different oxide scale layers.
Abstract: Cesiumiodide with Na doping (CsI(Na)) solution or melt is easily forming three- dimension dendrites on the free surface. The defects or bobbles form inside the CsI(Na) during the solution or melt solidification. The defects or bobbles can further effect the x-ray path in the CsI(Na) crystal and decrease the scintillation characteristics of CsI(Na). In order to enhance the CsI(Na) scintillated property we made single crystal of CsI(Na) column in the anodic aluminum oxide (AAO) template by hydraulic pressure injection method. It is interesting that when CsI(Na) melt is confined in the small AAO channels, the column grow as stable single column without any dendrites. The high aspect ratio (100~10000) of AAO and nano to sub-micron channel structure which is a suitable template for single of crystal CsI(Na) formation. In this work, a new low-cost approach to fabricate scintillator crystals using anodic aluminum oxide (AAO) rather than Si is reported, which can produce scintillator crystals with a wide range of controllable size to optimize their performance in X-ray detection.
Abstract: A two-dimensional linear wave-body interaction problem can be solved using a desingularized integral method by placing free surface Rankine sources over calm water surface and satisfying boundary conditions at prescribed collocation points on the
calm water surface. A new free-surface Rankine source distribution scheme, determined by the intersection points of free surface and body surface, is developed to reduce numerical computation cost. Associated with this, a new treatment is given to the intersection point. The present scheme results are in good agreement with traditional numerical results and measurements.
Abstract: The simulation of extrusion process is studied widely
in order to both increase products and improve quality, with broad
application in wire coating. The annular tube-tooling extrusion was
set up by a model that is termed as Navier-Stokes equation in
addition to a rheological model of differential form based on singlemode
exponential Phan-Thien/Tanner constitutive equation in a twodimensional
cylindrical coordinate system for predicting the
contraction point of the polymer melt beyond the die. Numerical
solutions are sought through semi-implicit Taylor-Galerkin pressurecorrection
finite element scheme. The investigation was focused on
incompressible creeping flow with long relaxation time in terms of
Weissenberg numbers up to 200. The isothermal case was considered
with surface tension effect on free surface in extrudate flow and no
slip at die wall. The Stream Line Upwind Petrov-Galerkin has been
proposed to stabilize solution. The structure of mesh after die exit
was adjusted following prediction of both top and bottom free
surfaces so as to keep the location of contraction point around one
unit length which is close to experimental results. The simulation of
extrusion process is studied widely in order to both increase products
and improve quality, with broad application in wire coating. The
annular tube-tooling extrusion was set up by a model that is termed
as Navier-Stokes equation in addition to a rheological model of
differential form based on single-mode exponential Phan-
Thien/Tanner constitutive equation in a two-dimensional cylindrical
coordinate system for predicting the contraction point of the polymer
melt beyond the die. Numerical solutions are sought through semiimplicit
Taylor-Galerkin pressure-correction finite element scheme.
The investigation was focused on incompressible creeping flow with
long relaxation time in terms of Weissenberg numbers up to 200. The
isothermal case was considered with surface tension effect on free
surface in extrudate flow and no slip at die wall. The Stream Line
Upwind Petrov-Galerkin has been proposed to stabilize solution. The
structure of mesh after die exit was adjusted following prediction of
both top and bottom free surfaces so as to keep the location of
contraction point around one unit length which is close to
experimental results.
Abstract: The paper discuses the effect of initial stresses on the reflection coefficients of plane waves in a dissipative medium. Basic governing equations are formulated in context of Biot's incremental deformation theory. These governing equations are solved analytically to obtain the dimensional phase velocities of plane waves propagating in plane of symmetry. Closed-form expressions for the reflection coefficients of P and SV waves- incident at the free surface of an initially stressed dissipative medium are obtained. Numerical computations, using these expressions, are carried out for a particular model. Computations made with the results predicted in presence and absence of the initial stresses and the results have been shown graphically. The study shows that the presence of compressive initial stresses increases the velocity of longitudinal wave (P-wave) but diminishes that of transverse wave (SV-wave). Also the numerical results presented indicate that initial stresses and dissipation might affect the reflection coefficients significantly.
Abstract: Displacement measurement was conducted on compact normal and shear specimens made of acrylic homogeneous material subjected to mixed-mode loading by digital image correlation. The intelligent hybrid method proposed by Nishioka et al. was applied to the stress-strain analysis near the crack tip. The accuracy of stress-intensity factor at the free surface was discussed from the viewpoint of both the experiment and 3-D finite element analysis. The surface images before and after deformation were taken by a CMOS camera, and we developed the system which enabled the real time stress analysis based on digital image correlation and inverse problem analysis. The great portion of processing time of this system was spent on displacement analysis. Then, we tried improvement in speed of this portion. In the case of cracked body, it is also possible to evaluate fracture mechanics parameters such as the J integral, the strain energy release rate, and the stress-intensity factor of mixed-mode. The 9-points elliptic paraboloid approximation could not analyze the displacement of submicron order with high accuracy. The analysis accuracy of displacement was improved considerably by introducing the Newton-Raphson method in consideration of deformation of a subset. The stress-intensity factor was evaluated with high accuracy of less than 1% of the error.
Abstract: Axisymmetric vibration of an infinite Pyrocomposite
circular hollow cylinder made of inner and outer pyroelectric layer of
6mm-class bonded together by a Linear Elastic Material with Voids
(LEMV) layer is studied. The exact frequency equation is obtained
for the traction free surfaces with continuity condition at the
interfaces. Numerical results in the form of data and dispersion
curves for the first and second mode of the axisymmetric vibration of
the cylinder BaTio3 / Adhesive / BaTio3 by taking the Adhesive layer
as an existing Carbon Fibre Reinforced Polymer (CFRP) are
compared with a hypothetical LEMV layer with and without voids
and as well with a pyroelectric hollow cylinder. The damping is
analyzed through the imaginary parts of the complex frequencies.
Abstract: The Marangoni convective instability in a horizontal
fluid layer with the insoluble surfactant and nondeformable free
surface is investigated. The surface tension at the free surface is
linearly dependent on the temperature and concentration gradients.
At the bottom surface, the temperature conditions of uniform
temperature and uniform heat flux are considered. By linear stability
theory, the exact analytical solutions for the steady Marangoni
convection are derived and the marginal curves are plotted. The
effects of surfactant or elasticity number, Lewis number and Biot
number on the marginal Marangoni instability are assessed. The
surfactant concentration gradients and the heat transfer mechanism at
the free surface have stabilizing effects while the Lewis number
destabilizes fluid system. The fluid system with uniform temperature
condition at the bottom boundary is more stable than the fluid layer
that is subjected to uniform heat flux at the bottom boundary.