Search results

The main purpose of this paper is to develop two efficient and accurate numerical analytical methods for engineering computation of natural sloshing frequencies and modes i the case of truncated circular conical tanks.

The numerical‐analytical methods are based on a Ritz Treftz variational scheme with two distinct analytical harmonic functional bases.

Comparative numerical analysis detects the limit of applicability of variational methods in terms of the semi‐apex angle and the ratio between radii of the mean free surface and the circular bottom. The limits are caused by different analytical properties of the employed functional bases. However, parallel use of two or more bases makes it possible to give an accurate approximation of the lower natural frequencies for relevant tanks. For V‐shaped tanks, dependencies of the lowest natural frequency versus the semi‐apex angle and the liquid depth are described.

The methods provide the natural sloshing frequencies for V‐shaped tanks that are valuable for designing elevated containers in seismic areas. Approximate natural modes can be used in derivations of nonlinear modal systems, which describe a resonant coupling with structural vibrations.

Although variational methods have been widely used for computing the natural sloshing frequencies, this paper presents their application for truncated conical tanks for the first time. An original point is the use of two distinct functional bases.

The purpose of this paper is to investigate different baffles on mitigating liquid sloshing in a rectangular tank due to a horizontal excitation and to find out the optimal selection of sloshing mitigation for practical applications.

The numerical study is conducted by using a proven improved smoothed particle hydrodynamics (SPH), which is convenient in tracking free surfaces and capable of obtaining smooth and correct pressure field.

Liquid sloshing effects in a rectangular tank with vertical middle baffles, horizontal baffles, T-shape baffles and porous baffles are investigated together with those without any baffles. It is found that the existence of baffles can mitigate sloshing effects and the mitigation performance depends on the shape, structure and location of the baffles. Considering the balance of sloshing mitigation performance and the complexity in structure and design, the I shaped and T shaped baffles can be good choices to mitigate sloshing effects.

The presented methodology and findings can be helpful in practical engineering applications, especially in ocean engineering and problems with large sloshing effects.

The SPH method is a meshfree, Lagrangian particle method, and therefore it is an attractive approach for modeling liquid sloshing with material interfaces, free surfaces and moving boundaries. In most previous literature, only simple baffles are investigated. In this paper, more complicated baffles are investigated, which can be helpful in practical applications and engineering designs.

The average residence time of liquid flowing over the surface of a rotating cone was determined numerically. The development and propagation of the free surface flow was simulated using the volume of fluid (VOF) method. The numerical simulations were validated using laboratory experiments using soy‐oil as a model liquid, and approximate analytical solutions of the simplified governing equations. The numerical simulations revealed the importance of the cone rotation frequencies and the minor influence of the cone angles on the residence times. Higher liquid throughputs produced smaller residence times. As expected, an increasing cone size results in proportionally higher residence times. Furthermore, it was established that even for small cones with a characteristic diameter of, e.g. less than 1m, relatively high (∼1 kg/s) throughputs of liquid are possible. It appears that the combination of the decreasing layer thickness and the increasing size of the numerical grid cells with increasing radial cone coordinate hampers the numerical simulation of this system.

Currently, human hearts destined for transplantation can be used for 4.5 hours which is often insufficient to test the heart, the purpose of this paper is to find a compatible recipient and transport the heart to larger distances. Cooling systems with simultaneous internal and external liquid cooling were numerically simulated as a method to extend the usable life of human hearts.

Coolant was pumped inside major veins and through the cardiac chambers and also between the heart and cooling container walls. In Case 1, two inlets and two outlets on the container walls steadily circulated the coolant. In the Case 2, an additional inlet was specified on the container wall thus creating a steady jet impinging one of the thickest parts of the heart. Laminar internal flow and turbulent external flow were used in both cases. Unsteady periodic inlet velocities at two frequencies were applied in Case 3 and Case 4 that had four inlets and four outlets on walls with turbulent flows used for internal and external circulations.

Computational results show that the proposed cooling systems are able to reduce the heart temperature from +37°C to almost uniform +5°C within 25 min of cooling, thus reducing its metabolic rate of decay by 95 percent. Calculated combined thermal and hydrodynamic stresses were below the allowable threshold. Unsteady flows did not make any noticeable difference in the speed of cooling and uniformity of temperature field.

This is the pioneering numerical study of conjugate convective cooling schemes capable of cooling organs much faster and more uniformly than currently practiced.

With the valveless stabilised liquid oxygen converter provided for the Tornado (MRCA), Normalair Garrett have produced the most advanced converter to date. Its advantages over the other systems discussed and demonstrated by H. A. S. Hamlin, Project leader—Life support systems.

Hydrothermal waves represent the preferred mode of instability of the so-called Marangoni flow for a wide range of liquids and conditions. The related features in classical rectangular containers have attracted much attention over recent years owing to the relevance of these oscillatory modes to several techniques used for the production of single crystals of semiconductor or oxide materials. Control or a proper knowledge of convective instabilities in these systems is an essential topic from a material/product properties saving standpoint. The purpose of this study is to improve our understanding of these phenomena in less ordinary circumstances.

This short paper reports on a numerical model developed to inquire specifically about the role played by sudden changes in the available cross-section of the shallow cavity hosting the liquid. Although accounting for the spanwise dimension would be necessary to derive quantitative results, the approach is based on the assumption of two-dimensional flow, which, for high-Pr fluids, is believed to retain the essence of the involved physical processes.

Results are presented for the case of a fluid with Pr = 15 filling an open container with a single backward-facing or forward-facing step on the bottom wall or with an obstruction located in the centre. It is shown that the presence of steps in the considered geometry can lead to a variety of situations with significant changes in the local spectral content of the flow and even flow stabilization in certain circumstances. The role of thermal boundary conditions is assessed by considering separately adiabatic and conducting conditions for the bottom wall.

Although a plethora of studies have been appearing over recent years motivated, completely or in part, by a quest to identify new means to mitigate these instabilities and produce accordingly single crystals of higher quality for the industry, unfortunately, most of these research works were focusing on very simple geometries. In the present paper, the causality and interdependence among all the kinematic and thermal effects mentioned above is discussed.

This study aims to focus on understanding how different jet angles and Reynolds numbers influence the phase change materials’ (PCMs) melting process and their capacity to store energy. This approach is intended to offer novel insights into enhancing thermal energy storage systems, particularly for applications where heat transfer efficiency and energy storage are critical.

The research involved an experimental and numerical analysis of PCM with a melting temperature range of 22 °C–26°C under various conditions. Three different jet angles (45°, 90° and 135°) and two container angles (45° and 90°) were tested. Additionally, two different Reynolds numbers (2,235 and 4,470) were used to explore the effects of jet outlet velocities on PCM melting behaviour. The study used a circular container and analysed the melting process using the hot air inclined jet impingement (HAIJI) method.

The obtained results showed that the average temperature for the last time step at Ф = 90° and Re = 4,470 is 6.26% higher for Ф = 135° and 14.23% higher for Ф = 90° compared with the 45° jet angle. It is also observed that the jet angle, especially for Ф = 90°, is a much more important factor in energy storage than the Reynolds number. In other words, the jet angle can be used as a passive control parameter for energy storage.

This study offers a novel perspective on the effective storage of waste heat transferred with air, such as exhaust gases. It provides valuable insights into the role of jet inclination angles and Reynolds numbers in optimizing the melting and energy storage performance of PCMs, which can be crucial for enhancing the efficiency of thermal energy storage systems.

The purpose of this paper is to analyze the liquid sloshing behaviors in two-dimensional tanks with various porous baffles under the external excitation.

Adopting the finite element method (FEM) and control variable method to study the impacts of the height, length, number, location, shape, porous-effect parameter of the porous baffle, the external load frequency and the shape of the tank on the liquid sloshing response.

The amplitude of the free surface can be reduced effectively when the baffle opening is appropriate. The anti-sway ability of the system increases in pace with the baffle’s height growing. Under the same conditions, the shapes of the baffles have an important effect on improving the anti-sway ability of the system.

As there exist the differences of the velocity potential between each side of the porous baffle, which means that there are two different velocity potentials at a point on the porous baffle, the conventional finite element modeling technologies are not suitable to be applied here. To deal with this problem, the points on the porous baffle are regarded as two nodes with the same coordinate to model and calculate.

The purpose of the present investigation is to compute the circulation flow of a liquid in a closed chamber when the liquid is fired by a gas jet through number of nozzles.

The conservation equations for mass and momentum have been solved in a closed container along with the conservation of volume fraction of the secondary phase in order to take into account the gas phase present in the liquid. The drag force created by the gas on the liquid has been incorporated in the momentum equation as a source term and the resulting equations have been solved numerically using a finite volume technique in an unstructured grid employing a phase coupled pressure linked velocity solver for the pressure correction equation, which is usually known as the Eulerian Scheme for two phase flow solution. An eddy viscosity based k‐ε turbulence model for the mixture was considered to update the fluid viscosity with iterations and capture the turbulence in the overall mixture rather than computing the individual turbulence in both the phases, which was found to be extremely time‐consuming and computationally unstable to some extent.

The model thus developed was tried to predict the circulation flow rate in an experimental setup where air was injected to drive the water in a long U tube setup. The computed circulation flow rate was found to be within 15 percent deviation from the experimentally observed values. The circulation flow rate of water was found to be increasing with the injected airflow rate. After this model validation, circulation flow rate of steel in an industrial size Ruhrstal‐Haraeus (RH)‐degasser was computed by injecting argon into the liquid steel through the up‐leg of the RH vessel. It was found that the circulation flow rate of steel in the RH degasser was increasing when the argon flow was being varied from 800 to 1,600 NL/min, which confirms the industrial findings.

The present computation could not use the energy equation to compute the swelling of the gas bubbles inside the chamber due to huge computing time requirement.

The present computation could compute realistically the circulation flow rate of water in a U tube when fired by a gas jet by using a two‐phase Eulerian model and hence this model can be effectively used for industrial applications where two‐phase flow comes into picture.

The original contribution of the paper is in the use of the state‐of the‐art Eulerian two‐phase flow model to predict circulation flow in an industrial size RH degasser.

The paper aims to study the effect of liquid flow velocity on corrosion behavior of 20# steel at initial stage under (CO2/aqueous solution) gas–liquid two-phase plug flow conditions.

Weight loss, scanning electron microscopy, energy-dispersive X-ray spectroscopy and XPS methods were used in this study.

The corrosion rate increased with the increasing liquid flow velocity at any different corrosion time. The corrosion rate decreased with the extension of corrosion time at the same liquid flow velocity. There was no continuous corrosion products film on the whole pipe wall at any different corrosion time. The macroscopic brown-yellow corrosion products on the pipe wall surface decreased with the increasing liquid flow velocity and the loose floccus corrosion products decreased gradually until these products were transformed into un-continuous needle-like dense products with the increasing liquid velocity. The main elements among the products film were Fe, C and O, and the main phases of products film on the pipe wall were Fe3C, FeCO3, FeOOH and Fe3O4. When the corrosion time was 1 h under different liquid–velocity condition, the thickness of local corrosion products film was from 3.5 to 3.8 µm.

The ion mass transfer model of corrosion process in pipe was put forward under gas–liquid two-phase plug flow condition. The total thickness of diffusion sublayer and turbulence sublayer decreased as well as the turbulence propagation coefficient increased with the increasing liquid velocity, which led to the increasing velocity of ion transfer during corrosion process. This was the fundamental reason for the increase of corrosion rate with the increasing liquid velocity.