Open Source Field Operation And Manipulation Research Articles

Three-dimensional numerical investigation of vortex-induced vibration of a circular cylinder at low Reynolds number

Three-dimensional numerical investigations of vortex-induced vibration (VIV) of flexibly mounted rigid circular cylinder have been carried out using Open Source Field Operation And Manipulation (OpenFOAM). The cylinder is allowed to oscillate in both streamwise as well as cross-stream directions. The cylinder response and wake transition are studied by varying the mass damping parameter ( m∗ζ , 0.01 and 0.1) and Reynolds number (Re). The computations are conducted for a wide range of reduced velocities ( Ur ), ranging from 3 to 8, at a fixed value of Re equal to 350. Furthermore, the effect of Re, varied from 250 to 400, on VIV is illustrated for a fixed value of Ur = 4.5. The response of the cylinder in terms of streamwise and cross-stream displacements is examined in and outside the lock-in regime. The cylinder trajectories are found to be influenced by the change in Ur , m∗ζ and Re. The dynamic characteristics are also disclosed with the change in Ur , m∗ζ and Re. In order to capture the effect of Re on the wake transition of the circular cylinder undergoing VIV, the value of Re is varied from 250 to 400. Mode A instability is observed to be delayed due to VIV in comparison to that of the stationary cylinder.

Direct numerical simulations on oscillating flow past surface-mounted finite-height circular cylinder

In this work, a surface-mounted circular cylinder with a fixed aspect ratio (ratio of height of the cylinder to its diameter) of 5 is subjected to a non-zero mean oscillating flow with a range of frequencies and amplitudes. Three-dimensional direct numerical simulations are then conducted on this finite-height cylinder. The mass and momentum equations are resolved using the finite volume-based Open Source Field Operation and Manipulation (OpenFOAM). A fixed Reynolds number Re=UoD/ν of 250 is used in this study, which is defined based on mean velocity at the inlet ( Uo ) and cylinder diameter (D). Here ν is the kinematic viscosity of the working fluid. Non-dimensional velocity oscillation amplitude ( A∗=a/Uo ) is varied from 0.1 to 0.3, while the non-dimensional oscillation frequency ( f∗=f/fo ) takes the values of 0.33, 0.5, 1, 2, and 3. Here a and f are the dimensional oscillation amplitude and frequency, respectively and fo is the vortex shedding frequency corresponding to a uniform flow at Re = 250. The three-dimensional vortex structures, presented with the help of iso-Q surfaces, show that the oscillating flow changes the size and shape of the hairpin-shaped vortices. Wake is found to be synchronized with the oscillation frequency at f* = 2 for each value of the A* and results in the maximum lift force on the cylinder. Hilbert Huang transformation analysis of the transverse velocity signals at a specific point in the wake reveals that the wake is more complex and aperiodic in nature for f* values of 0.33, 0.5, and 1, whereas it is periodic for f* = 2 and 3. In order to further disclose the nonlinearity associated with the oscillating flow, the degree of stationarity is discussed corresponding to each value of A* and f*. Dynamic mode decomposition is exploited to obtain information about the coherent vortical structures and their spatial and temporal behavior in the wake with a change in the value of f*. Effects of A* and f* on the dynamic characteristics are also investigated.

Implementation of a level-set-based volume penalization method for solving fluid flows around bluff bodies

A volume penalization-based immersed boundary technique is developed and thoroughly validated for fluid flow problems, specifically flow over bluff bodies. The proposed algorithm has been implemented in an open source field operation and manipulation (OpenFOAM), a computational fluid dynamics solver. The immersed boundary method offers the advantage of inserting a complex solid object inside a Cartesian grid system, and therefore, the governing equations can be applied to such a simpler grid arrangement. For capturing the fluid–solid interface more accurately, the grid is refined near the solid surface using topoSetDict and refineMeshDict utilities in OpenFOAM. In order to avoid any numerical oscillation and to compute the gradients accurately near the interface, the present volume penalization method (VPM) is integrated with a signed distance function, which is also referred to as a level-set function. Benchmark problems, such as flows around a cylinder and a sphere, are considered and thoroughly validated with the results available in the literature. For the flow over a stationary cylinder, the Reynolds number is varied so that it covers from a steady two-dimensional flow to an unsteady three-dimensional flow. The capability of the present solver has been further verified by considering the flow past a vibrating cylinder in the cross-stream direction. In addition, a flow over a sphere, which is inherently three-dimensional due to its geometrical shape, is validated in both steady and unsteady regimes. The results obtained by the present VPM show good agreement with those obtained by a body-fitted grid using the same numerical scheme as that of the VPM, and also with those reported in the literature. The present results indicate that the VPM-based immersed boundary technique can be widely applicable to scientific and engineering problems involving flow past stationary and moving bluff bodies of arbitrary geometry.

Hydrodynamic characteristics and particle tracking of 90° lateral intakes at an inclined river slope

Lateral intakes are common in rivers. The pump efficiency and sediment deposition are determined by the local hydrodynamic characteristics and mainstream division width. The hydraulic characteristics of lateral withdrawal from inclined river slopes at different intake elevations should be investigated. Meanwhile, the division width exhibits significant vertical non-uniformity at an inclined river slope, which should be clarified. Hence, a three-dimensional (3-D) hydrodynamic and particle-tracking model was developed with the Open Source Field Operation and Manipulation (OpenFOAM), and the model was validated with physical model tests for 90° lateral withdrawal from an inclined side bank. The flow fields, withdrawal sources, and division widths were investigated with different intake bottom elevations, withdrawal discharges, and main channel velocities. This study showed that under inclined side bank conditions, water entered the intake at an oblique angle, causing significant 3-D spiral flows in the intake rather than two-dimensional closed recirculation. A lower withdrawal discharge, a lower bottom elevation of the intake, or a higher main channel velocity could further strengthen this phenomenon. The average division width and turbulent kinetic energy were smaller under inclined side bank conditions than under vertical bank conditions. With a low intake bottom elevation, a low withdrawal discharge, or a high main channel velocity, the sources of lateral withdrawal were in similar ranges near the local inclined bank in the vertical direction. Under inclined slope conditions, sediment deposition near the intake entrance could be reduced, compared to that under vertical slope conditions. The results provide hydrodynamic and sediment references for engineering designs for natural rivers with inclined terrains.

CFD SIMULATION AND VALIDATION FOR MIXING VENTILATION SCALED-DOWN EMPTY AIRCRAFT CABIN USING OPENFOAM

An investigation into the spread of the COVID-19 virus within a confined space including an aircraft cabin is essential in order to find out the mechanism. However, this is time-consuming and limited in scope, so a computational fluid dynamics (CFD) simulation is used instead. Therefore, a prior study and an appropriate choice of turbulence model are required before the simulation. The main objective of this study is to validate and evaluate the results predicted by the Open Source Field Operation and Manipulation (OpenFOAM) software through comparison with the experimental data from the literature which was conducted using particle image velocimetry (PIV) measurement. Three different Reynolds-averaged Navier-Stokes turbulence models were selected; Re-normalisation Group k - ɛ (RNG), Realizable k - ɛ (RLZ) and Low Reynold Number (LRN) to simulate a mixing ventilation system of a scaled-down model of empty aircraft cabin. In the RNG and LRN model cases, a fairly large circulation flows were observed on the right and left sides of the model representing the passenger area. The results were also evaluated quantitatively using the factor of two of observations (FAC2) for the velocity components and turbulent kinetic energy (TKE) with root mean square error (RMSE) for the former and normalised mean square errors (NMSE) for the latter. The simulation results showed that RNG and LRN are capable of predicting the flow field well. However, for TKE prediction LRN performed better than RNG which concluded that LRN is the suitable turbulence model in simulating flow fields in investigated case.

NUMERICAL SIMULATION OF A FLUID FLOW IN THE CYLINDRICAL DUCT CONTAINING TWO DIAPHRAGMS WITH ORIFICES OF DIFFERENT DIAMETERS

The flow of a viscous incompressible fluid in a cylindrical duct with two serial diaphragms with orifices of different diameters was studied based on the numerical solution of the unsteady Navier-Stokes equations. The solution algorithm was based on the finite volume method using second-order accurate difference schemes in both space and time. The TVD (Total-Variation Diminishing) form of a central-difference scheme with a flux limiter was used for the interpolation of the convective terms. The combined evaluation of the velocity and pressure fields was carried out using the PISO (Pressure Implicit Split Operator) procedure. The problem was solved using OpenFOAM (Open-source Field Operation And Manipulation) open-source toolkit libraries using the computing power of the cluster supercomputer of the V.M. Glushkov Institute of Cybernetics of the National Academy of Sciences of Ukraine. It was shown that within a certain range of the diameter ratio of the orifices of the diaphragms, a circulation movement is established in the cavity between the diaphragms. The boundary layer breaks off from the surface of the first diaphragm and forms an annular shear layer. When approaching the second diaphragm, a sequence of ring vortices is formed in it, which interact with the surface of the diaphragm and lead to the emergence of tonal sound. When the ratio of the orifice diameter of the second diaphragm to the first decreases, the share of the jet’s kinetic energy, which participates in the circulation in the middle of the cavity between the diaphragms, increases. As a result, the amplitude of velocity fluctuations in the orifice of the second diaphragm decreases. When the ratio of the diameters of the orifices increases, the share of energy involved in the circulation decreases, and when a certain value of the ratio is reached, the interaction between the vortices in the shear layer and the surface of the diaphragm does not occur. As a result, the excitation of the tonal sound stops. The Strouhal number of self-oscillations practically does not change when the diameter ratio of the orifices changes. During the calculations, two different modes of self-oscillation were obtained, which is consistent with previous works.

Visualization and Parametric Study on Vortex Shedding Suppression of Cylindrical Structures in Offshore Engineering Using Large Eddy Simulation

Cylindrical structures are widely used in offshore and marine engineering, but they may suffer from vortex-induced vibration under the influence of ocean or wave currents, which can lead to severe fatigue damage. In this study, we applied the open-source software Open-Source Field Operation and Manipulation (OpenFOAM) to investigate the characteristics of fluid flow around offshore cylindrical structures, taking into account the effect of helical strake parameters, such as pitch and strake number. The aim of this study is to explore the possibility of suppressing vortex shedding with different helical strake parameters. Numerical simulation results demonstrated that attaching a helical strake to the bare cylinder destroyed vortex shedding in offshore cylindrical structures. The vortex visualization showed that the helical strake destroyed the three-dimensional vortex structures. Moreover, the lift coefficient data showed that the vibration frequency of the cylinder decreased after attaching the helical strake, indicating that the vortex-induced vibrations on the wake flow tended to fade. The results suggest that the helical strake is a promising option for suppressing the wake vortex shedding of cylindrical structures in offshore engineering.

Effect of size and spacing on the wake characteristics of two spheres placed in tandem

In the present study, the flow past two spheres placed in a tandem arrangement is investigated numerically using open source field operation and manipulation (OpenFOAM) at a fixed Reynolds number (Re) of 300, where the Reynolds number is defined based on the diameter (D) of the downstream sphere (DS) and the freestream velocity at the inlet. Keeping size of the DS constant, the diameter of the upstream sphere (US) has been varied, so that the diameter ratio, DR (ratio of the diameters upstream and downstream spheres), takes the values of 0.4, 0.6, 0.8, 1.0, and 1.5. The spacing between the spheres (S) is also varied from 1D to 5D. Iso-Q surfaces show that both upstream and downstream wakes undergo various transitions with changes in the values of DR and S. For US, transition from steady symmetric to planar symmetric occurs at DR = 0.6 and S = 1, which corresponds to a local Re of 180. For DR = 0.8, steady to unsteady transition occurs at S = 2, whereas for all other values of S, the wake remains steady. For DR > 0.8 and S > 1, both upstream and downstream wakes are found to be unsteady. Hilbert analysis revealed that unsteady wakes are periodic for DR values of 0.4, 0.8, and 1.0 and are aperiodic for DR = 1.5. In addition, the extent of wake nonlinearity has also been quantified in terms of degree of stationarity. The drag force on both the spheres increases with an increase in spacing between the spheres.

Numerical study of rotating convection with bi-directional temperature gradients

In this work, we report results from three-dimensional numerical simulations of convection in a rotating cylindrical annulus with mean temperature gradients present along the radial and vertical directions. The simulations have been performed in OpenFOAM (Open-source Field Operation And Manipulation). The bi-directional temperature gradient is achieved by imposing uniform heat flux in a thin circular strip at the outer periphery of the bottom surface, and a uniform temperature at the entire inner cylindrical surface. The study has been carried out for a particular Rayleigh number, Ra = 4.76×108, and a range of Taylor numbers, Ta = 6.5×108–2.7×109. The Ra value is chosen so that some qualitative comparison could be done with available experimental data. A laminar flow solver has been used to perform the computations in a rotating frame of reference. Complex empirical orthogonal function analysis has been carried out for the velocity and temperature fields. The contour maps of the radial velocity field in the r−θ plane clearly show the existence of mode 4 and 4 AV (amplitude vacillation) baroclinic waves at low values of Ta, which break down into eddies as Ta increases. The mean temperature field in the r–z plane supports the existence of Convective Columnar Plumes (CCP) on the outer edge, which interacts with the baroclinic waves, aiding in the transport of heat in the system. The fluctuations in CCP disintegrate into smaller-scale structures at the highest Ta. Quantification of turbulent fluxes shows the effectiveness of the heat transport at Ta∼O(108), as compared to Ta∼O(109). Our results clearly demonstrate the role of waves, plumes, and eddies on heat transport in the presence of bi-directional temperature gradients.

Actuator line model using simplified force calculation methods

Abstract. To simulate transient wind turbine wake interaction problems using limited wind turbine data, two new variants of the actuator line technique are proposed in which the rotor blade forces are computed locally using generic load data. The proposed models, which are extensions of the actuator disk force models proposed by Navarro Diaz et al. (2019a) and Sørensen et al. (2020), only demand thrust and power coefficients and the tip speed ratio as input parameters. In the paper the analogy between the actuator disk model (ADM) and the actuator line model (ALM) is shown, and from this a simple methodology to implement local forces in the ALM without the need for knowledge of blade geometry and local airfoil data is derived. Two simplified variants of ALMs are proposed, an analytical one based on Sørensen et al. (2020) and a numerical one based on Navarro Diaz et al. (2019a). The proposed models are compared to the ADM using analogous data, as well as to the classical ALM based on blade element theory, which provides more detailed force distributions by using airfoil data. To evaluate the local force calculation, the analysis of a partial-wake interaction case between two wind turbines is carried out for a uniform laminar inflow and for a turbulent neutral atmospheric boundary layer inflow. The computations are performed using the large eddy simulation facility in Open Source Field Operation and Manipulation (OpenFOAM), including Simulator for Wind Farm Applications (SOWFA) libraries and the reference National Renewable Energy Laboratory (NREL) 5 MW wind turbine as the test case. In the single-turbine case, computed normal and tangential force distributions along the blade showed a very good agreement between the employed models. The two new ALMs exhibited the same distribution as the ALM based on geometry and airfoil data, with minor differences due to the particular tip correction needed in the ALM. For the challenging partially impacted wake case, both the analytical and the numerical approaches manage to correctly capture the force distribution at the different regions of the rotor area, with, however, a consistent overestimation of the normal force outside the wake and an underestimation inside the wake. The analytical approach shows a slightly better performance in wake impact cases compared to the numerical one. As expected, the ALMs gave a much more detailed prediction of the higher-frequency power output fluctuations than the ADM. These promising findings open the possibility to simulate commercial wind farms in transient inflows using the ALM without having to get access to actual wind turbine and airfoil data, which in most cases are restricted due to confidentiality.

Improving Wind Speed Forecasting for Urban Air Mobility Using Coupled Simulations

Hazardous weather, turbulence, wind, and thermals pose a ubiquitous challenge to Unmanned Aircraft Systems (UAS) and electric-Vertical Take-Off and Landing (e-VTOL) aircrafts, and the safe integration of UAS into urban area requires accurate high-granularity wind data especially during landing and takeoff phases. Two models, namely, Open-Source Field Operation and Manipulation (OpenFOAM) software package and Weather Research and Forecasting (WRF) model, are used in the present study to simulate airflow over Downtown Oklahoma City, Oklahoma, United States. Results show that computational fluid dynamics wind simulation driven by the atmospheric simulation significantly improves the simulated wind speed because the accurate modeling of the buildings affects wind patterns. The evaluation of different simulations against six Micronet stations shows that WRF-CFD numerical evaluation is a reliable method to understand the complicated wind flow within built-up areas. The comparison of wind distributions of simulations at different resolutions shows better description of wind variability and gusts generated by the urban flows. Simulations assuming anisotropy and isotropy of turbulence show small differences in the predicted wind speeds over Downtown Oklahoma City given the stable atmospheric stratification showing that turbulent eddy scales at the evaluation locations are within the inertial subrange and confirming that turbulence is locally isotropic.

OpenFOAM for computational hydrodynamics using finite volume method

Partial differential equations may explain anything from planetary movement to tectonic plate, yet it is notoriously difficult to resolve them. Turbulence is present in nearly all fluid flows, and pure laminar flow is extremely unusual in practice. The Large Eddy Simulation (LES) computational model is employed for the simulation of turbulence flow on a spillway having four inlets with a single outlet. Such flows are observed at hydroelectric power dams all over the world. The fluctuated flows produced a large amount of energy in terms of electricity that costs a very low amount compared to the energy obtained in tidal power sectors. In the production of hydropower energy, the flow simulation is of great interest. This paper focuses on the study of turbulence kinetic energy with the help of a LES model. The spillway considered in this paper contains four inlets and a single outlet. The four inlets will allow more flow which will insert more pressure nearer the outlet. The kinetic energy is computed at the inlets and outlet in the turbulent flow. The fluctuated velocity along with the mean velocity at the inlets and outlet is also computed along with the pressure. The C++-based programming is made, which is simulated in Open-source Field Operation and Manipulation (OpenFOAM). The graphs are presented for a better understanding of readers.

Selection of Characteristic Particle Size of Drilling Cuttings Based on Adsorption-Desorption Properties: Experiment and Simulation

The direct method for measuring gas content based on image processing, that is, comparing the desorption curves of drilling cuttings with the characteristic particle size with the desorption image in the database to determine the gas content of coal seams, can solve the problems of long measurement time and larger measurement errors of the traditional direct method. However, the precondition to realize this method is to define the characteristic particle size of drilling cuttings. Therefore, this paper performed adsorption and desorption experiments under different particle size ranges and dissimilar adsorption equilibrium pressures, solved the gas migration control equation in coal particles, and determined the characteristic particle size from experimental results and simulation effect. The results showed that gas desorption quantity was positively correlated with adsorption equilibrium pressure, while negatively correlated with the size of a particle size range. The greater the adsorption equilibrium pressure, the better the linear fitting relationship with the maximum desorption quantity, and the maximum correlation coefficient could reach 0.9658. Based on Darcy’s law, Langmuir adsorption equation, and mass conservation law, the theoretical control equation of gas migration in coal particles was derived and solved by Open source Field Operation and Manipulation (OpenFOAM). The accuracy of the calculation results of the solver was verified by the derived analytical solution of the gas pressure change of coal particles without considering the adsorption. From the experimental results and simulation results, the characteristic particle size of drilling cuttings was determined to be 0.5-1 mm. The research can provide help for the subsequent development of direct gas content measurement technology based on image processing.

A multi-scale model for day-ahead wind speed forecasting: A case study of the Houhoku wind farm, Japan

Accurate day-ahead wind prediction from numerical weather prediction (NWP) models is crucial for integrating wind power into existing grids and wind farm arrangements. Most wind farms in Japan are scattered in mountainous areas susceptible to strong winds and wind patterns that are highly fluctuating and cannot be captured well with mesoscale NWP models. Aiming to provide reliable day-ahead forecasts in rugged zones with complex topography, we introduce a multi-scale model based on coupling the Weather Research Forecasting (WRF) model and the Open Source Field Operation and Manipulation (OpenFOAM) Computational Fluid Dynamics model. The WRF’s predictions are used as initial/boundary conditions for a microscale simulation using OpenFOAM. This coupled model can overcome the shortcomings when employing each component independently. Numerical simulations and validations are carried out at a wind farm in Houhoku, Japan, during summer and winter. The multi-scale model results show a considerable reduction in the bias, mean absolute error, and root mean square error compared to the WRF-alone forecasts. It also has a superior skill against most post-processing algorithms or artificial intelligence techniques to reduce random forecasting errors. In addition, analysis of computational cost demonstrates the promising capability of this multi-scale model being a routine tool.

An optimized Eulerian–Lagrangian method for two-phase flow with coarse particles: Implementation in open-source field operation and manipulation, verification, and validation

The solid–liquid two-phase flow with coarse particles is ubiquitous in natural phenomena and engineering practice, which is characterized by coarse particles, high particle concentration, and large particle size distribution. In this work, the numerical models describing two-phase flows are reviewed, which given that the Eulerian–Lagrangian method is applicable in this work. Then, some modified models are proposed for the situation where the conventional Eulerian–Lagrangian method is not applicable to deal with coarse particles. The continuous phase equations of liquid are solved based on the finite volume method. The pressure implicit with splitting of operators algorithm for solving the Navier–Stokes (N–S) equations of the pseudo-single-phase flow, considering phase fraction and momentum exchange source term, is proposed. The discrete coarse particle is tracked in the Lagrangian method. A virtual mass distribution function is proposed for calculating coarse particle volume fraction. A weighted function method relating to the particle size is given for the interpolation between the Eulerian and Lagrangian fields. The barycentric coordinates are introduced into the particle localization. All the modified models are algorithmically implanted in the open-source field operation and manipulation (OpenFOAM) as a new solver named coarse discrete particle method FOAM (CoarseDPMFoam). Subsequently, the applicability of the numerical simulation method is verified by some typical test cases. The proposed numerical simulation method provides new ideas and methods for the mechanism investigation and engineering application of the two-phase flow with coarse particles.

The research into the propagation law of the shock wave of a gas explosion inside a building

Based on the dissipation rate conservation equations of turbulent kinetic energy in the k-ε turbulence model, a complicated three-dimensional finite element model of a kitchen filled with gas mixture is developed by using the open source field operation and manipulation (OpenFOAM). Two representative kitchens were used to investigate the propagation law of the shock wave of a gas explosion inside a building by considering the key characteristics of the blast shock wave. The influence of some crucial parameters, such as initial conditions and kitchen parameters, on the properties of the blast shock wave is investigated. The basic steps to predict the peak pressure of the blast shock wave are given in consideration of the initial condition and the kitchen whilst the injury effect of the blast shock wave on the humans and animals is evaluated. The research results indicate that the pressure time history and the peak pressure space distribution are greatly influenced by the kitchen design layout. The coupled interaction between the initial temperature and gas volume concentration, especially at the upper and lower explosion limits of the gas, significantly affects the peak pressure. The peak pressure varies significantly with the opening and the buffer; however, it has little relation with the width, length, and height of the kitchen. The proposed method can accurately and effectively predict the peak pressure of the blast shock wave inside buildings. In terms of the peak pressure space distribution of the explosion shock wave, the peak pressure is much higher than the threshold of the killing pressure, which is unsafe for the humans and animals in the building.

A Two-Day Case Study: Comparison of Turbulence Data from an Unmanned Aircraft System with a Model Chain for Complex Terrain

The airborne measurement platform MASC-3 (Multi-Purpose Airborne Sensor Carrier) is used for measurements over a forested escarpment in the Swabian Alps to evaluate the wind field. Data from flight legs between 20 and 200 m above the ground on two consecutive days with uphill (westerly) flow in September 2018 are analyzed. In the lowest 140 m above the ground a speed-up is found with increased turbulence and changes in wind direction directly over the escarpment, whereas in the lowest 20 to 50 m above the ground a deceleration of the flow is measured. Additionally, simulation results from a numerical model chain based on the Weather Research and Forecasting (WRF) model and an OpenFOAM (Open Source Field Operation and Manipulation) model, developed for complex terrain, are compared to the data captured by MASC-3. The models and measurements compare well for the mean wind speed and inclination angle.

Investigation of the pressure vessel lower head potential failure under IVR-ERVC condition during a severe accident scenario in APR1400 reactors

In the event of a core meltdown in a high-power reactor, the integrity of the reactor pressure vessel is presumably protected by severe accident mitigation systems such as in-vessel retention external reactor vessel cooling (IVR-ERVC). However, in the late phase of the accident, two possible locations on the RPV are prone to failure: the location of the focusing effect and location of in-core instrument penetration. These two potential points of damage in the RPV are investigated in this study. A numerical model for the prediction of the natural convection, melting, and solidification processes for IVR-ERVC is presented. The model is based on the enthalpy-porosity approach with an extension for continuous liquid fraction function. The model is implemented in open-source field operation and manipulation (OpenFOAM) computational fluid dynamic code to produce a new solver which is based on the combination of conjugate heat transfer solver and buoyant-driven natural convection solver and the new solver is validated against the melting Gallium experimental test, in-core instrumentation failure experimental test, and BALI experimental test. This numerical model is applied for the investigation of the RPV rupture at the location of the focusing effect and in-core instrumentation penetrations. Severe ablations of the cladding and the weld materials are observed at a heat load of about ~1800 K which is expected to lead to the ejection of the penetration tubes if the force holding the penetration tube in place is lower than the force exerted by the system pressure. Subsequently, a two-layer IVR configuration is assessed and the integrity of the RPV is found not to be compromised under external reactor vessel cooling. However, in the case of a boiling crisis, the temperature of the ex-vessel wall is expected to rise quickly and this is simulated by increasing the ex-vessel wall temperature. The RPV is found to fail near the beltline due to a phenomenon known as focusing effect when the ex-vessel wall temperature rises above 1200 K.

Effects of High-speed Wind, Humidity, and Temperature on the Generation of a SARS-CoV-2 Aerosol; a Novel Point of View

The new coronavirus (SARS-CoV-2) is rapidly spreading across communities around the world. Respiratory droplet transmission is a common transmission route for many airborne diseases, including novel coronavirus disease (COVID-19). Wearing the face mask prevents respiratory droplet transmission. Both face mask leakage and non-use of the face mask under high-speed wind conditions can increase the risk of SARS-CoV-2 transmission. The respiratory droplets’ behavior during sneezing or coughing (i.e., the size and the distance between droplets) depends on face mask wearing. The respiratory droplets during coughing and sneezing break apart into extremely small respiratory droplets (i.e., cloud of aerosol) upon interaction with the high-speed wind condition. The volume-of-fluid (VOF) method has been used to study the deformation and breakup of a single respiratory droplet placed in high-speed wind flow in the presence of smaller neighboring respiratory droplets. The effect of the diameter of respiratory droplets and the distance between them on the breakup mechanism has been investigated using open-source field operation and manipulation (OpenFOAM) libraries. The results are presented for the dimensionless distance and the normalized diameter of the secondary respiratory droplet in the ranges of 2 to 6 and 0.4 to 0.8, respectively. According to the results of simulations, the spreading of SARS-CoV-2 increases with a decrease in respiratory droplet breakup time. In addition, the effects of ambient relative humidity and temperature on the cloud of respiratory particles were analyzed.

Numerical Simulation of Wave Interaction with Payloads of Different Postures Using OpenFOAM

A three-dimensional numerical wave tank (NWT) is established with Open Source Field Operation and Manipulation (OpenFOAM) software and waves2foam to investigate wave interaction with payloads with different postures in the process of offshore lifting or lowering. Numerical results of regular wave interaction with a vertically suspending cylinder are presented first for validation by comparison with the published data. A series of simulation experiments are carried out, and the forces and the moments exerted by the regular waves on a fixed suspending cylinder payload and a fixed suspending cuboid payload with different postures are presented. It can be concluded from the results that the rotating rectangular payload (cuboid and cylinder) suffers a drastically changed moment when it is initially vertically placed, and the projection area of payload vertical to the force affects the corresponding force. The simulation results also show how the forces and the moments change with different posture angles. With some certain posture, the suspending payload suffers minimum forces and moments. Parametric study for the cuboid payload is done in the case of normal incidence. The influence of the payload’s size and wave parameters on forces and moments are analyzed. All of the numerical simulation results and conclusions provide the fundamentals for further research and safe control of offshore lifting or lowering.