Analyses of Gas Stratification Erosion by a Vertical Jet in Presence of an Obstacle Using the GOTHIC Code

Author(s):  
Michele Andreani ◽  
Ralf Kapulla

The GOTHIC code was validated using three experiments carried out in the PANDA facility in the framework of the OECD/NEA project HYMERES. These tests addressed the mixing of an initially stratified atmosphere by means of a vertical jet in the presence of on obstacle (circular plate). This paper reports on the simulations of three experiments, and complementary, quasi-steady state tests without stratification, where the flow structure above the impingement plate could be observed by means of PIV velocity measurements in a region larger than that considered in the transient experiments. Moreover, simulations of similar tests without obstacle conducted during the OECD/SETH-2 are also discussed. The reference, best-estimate model used for the analyses of the three experiments with different flow rates and initial and pressure boundary conditions was built on the base of a multi-step approach. This was based on mesh and modelling sensitivity studies mostly performed for the complementary tests, to assess the capability to represent the flow structure produced by the jet-plate interaction with different meshes around the plate. Generally, the results show that the use of a coarse mesh and the standard k-ε turbulence model permits a reasonable representation of the erosion process, but with a systematic over prediction of the mixing time. The results with the reference model were more accurate for two experiments with two flow rates and same initial conditions and all complementary tests. For the third test with different initial and boundary conditions, however, poor results were obtained with the reference model, which could only be improved by further refining the mesh. These results indicate that a model “qualified” for certain conditions could be inadequate for other cases, and sensitivity studies are necessary for the specific conditions considered in the analyses.

2020 ◽  
Vol 6 (2) ◽  
Author(s):  
Michele Andreani ◽  
Ralf Kapulla ◽  
Stephan Kelm ◽  
Domenico Paladino ◽  
Sidharth Paranjape

Abstract The GOTHIC code was validated using three experiments carried out in the PANDA facility in the framework of the OECD/NEA HYMERES project. These tests addressed the mixing of an initially stratified atmosphere by means of a vertical jet in the presence of on obstacle (circular plate). This paper reports on the simulations of three experiments, and complementary, quasi-steady-state tests without stratification, where the flow structure above the impingement plate could be observed by means of particle image velocimetry (PIV) velocity measurements in a region larger than that considered in the transient experiments. Moreover, simulations of similar tests without obstacle conducted during the OECD/SETH-2 project are also discussed. The reference, best-estimate model used for the analyses of the three experiments with different flowrates and initial and pressure boundary conditions was built on the base of a multistep approach. This was based on mesh and modeling sensitivity studies mostly performed for the complementary tests, to assess the capability to represent the flow structure produced by the jet–plate interaction with different meshes around the plate. Generally, the results show that the use of a coarse mesh and the standard k–ε turbulence model permits a reasonable representation of the erosion process, but with a systematic over prediction of the mixing time. The results with the reference model were more accurate for two experiments with two flowrates and same initial conditions and all complementary tests. For the third test with different initial and boundary conditions, however, poor results were obtained with the reference model, which could only be improved by further refining the mesh. These results indicate that a model “qualified” for certain conditions could be inadequate for other cases, and sensitivity studies are necessary for the specific conditions considered in the analyses.


2020 ◽  
Vol 4 (41) ◽  
pp. 57-62
Author(s):  
SHAVKAT KLYCHEV ◽  
◽  
BAKHRAMOV SAGDULLA ◽  
VALERIY KHARCHENKO ◽  
VLADIMIR PANCHENKO ◽  
...  

There are needed energy (heat) accumulators to increase the efficiency of solar installations, including solar collectors (water heaters, air heaters, dryers). One of the tasks of designing heat accumulators is to ensure its minimal heat loss. The article considers the problem of determining the distribution of temperatures and heat losses by convection and radiation of the heat insulation-accumulating body (water) system for a ball heat accumulator under symmetric boundary conditions. The problem is solved numerically according to the program developed on the basis of the proposed «gap method». (Research purpose) The research purpose is in determining heat losses by convection and radiation of a two-layer ball heat accumulator with symmetric boundary conditions. (Materials and methods) Authors used the Fourier heat equation for spherical bodies. The article presents the determined boundary and initial conditions for bodies and their surfaces. (Results and discussion) The thickness of the insulation and the volume of the heat accumulator affect the dynamics and values of heat loss. The effect of increasing the thickness of the thermal insulation decreases with increasing its thickness, starting with a certain volume of the heat accumulator or with R > 0.3 meters, the heat losses change almost linearly over time. The dynamics of heat loss decreases with increasing shelf life, but the losses remain large. (Conclusions) Authors have developed a method and program for numerical calculation of heat loss and temperature over time in a spherical two-layer heat accumulator with symmetric boundary conditions, taking into account both falling and intrinsic radiation. The proposed method allows to unify the boundary conditions between contacting bodies.


2021 ◽  
Vol 2021 (1) ◽  
Author(s):  
Ji Lin ◽  
Yuhui Zhang ◽  
Chein-Shan Liu

AbstractFor nonlinear third-order three-point boundary value problems (BVPs), we develop two algorithms to find solutions, which automatically satisfy the specified three-point boundary conditions. We construct a boundary shape function (BSF), which is designed to automatically satisfy the boundary conditions and can be employed to develop new algorithms by assigning two different roles of free function in the BSF. In the first algorithm, we let the free functions be complete functions and the BSFs be the new bases of the solution, which not only satisfy the boundary conditions automatically, but also can be used to find solution by a collocation technique. In the second algorithm, we let the BSF be the solution of the BVP and the free function be another new variable, such that we can transform the BVP to a corresponding initial value problem for the new variable, whose initial conditions are given arbitrarily and terminal values are determined by iterations; hence, we can quickly find very accurate solution of nonlinear third-order three-point BVP through a few iterations. Numerical examples confirm the performance of the new algorithms.


1965 ◽  
Vol 16 (2) ◽  
pp. 122-138 ◽  
Author(s):  
L. F. Henderson

SummaryThe intake is conceived as a device that generates a shock system and a subsonic patch within a supersonic field. The object of the paper is to obtain a qualitative description of this flow structure. The analytical procedure is a simple but fairly powerful mapping technique. Several examples are discussed and these show how the flow on the intake alters under the influence of gradually changing boundary conditions.


2008 ◽  
Vol 2008 ◽  
pp. 1-7 ◽  
Author(s):  
Mantas Povilaitis ◽  
Egidijus Urbonavičius

An issue of the stratified atmospheres in the containments of nuclear power plants is still unresolved; different experiments are performed in the test facilities like TOSQAN and MISTRA. MASPn experiments belong to the spray benchmark, initiated in the containment atmosphere mixing work package of the SARNET network. The benchmark consisted of MASP0, MASP1 and MASP2 experiments. Only the measured depressurisation rates during MASPn were available for the comparison with calculations. When the analysis was performed, the boundary conditions were not clearly defined therefore most of the attention was concentrated on MASP0 simulation in order to develop the nodalisation scheme and define the initial and boundary conditions. After achieving acceptable agreement with measured depressurisation rate, simulations of MASP1 and MASP2 experiments were performed to check the influence of sprays. The paper presents developed nodalisation scheme of MISTRA for the COCOSYS code and the results of analyses. In the performed analyses, several parameters were considered: initial conditions, loss coefficient of the junctions, initial gradients of temperature and steam volume fraction, and characteristic length of structures. Parametric analysis shows that in the simulation the heat losses through the external walls behind the lower condenser installed in the MISTRA facility determine the long-term depressurisation rate.


Author(s):  
Zhihang Song ◽  
Bruce T. Murray ◽  
Bahgat Sammakia

The integration of a simulation-based Artificial Neural Network (ANN) with a Genetic Algorithm (GA) has been explored as a real-time design tool for data center thermal management. The computation time for the ANN-GA approach is significantly smaller compared to a fully CFD-based optimization methodology for predicting data center operating conditions. However, difficulties remain when applying the ANN model for predicting operating conditions for configurations outside of the geometry used for the training set. One potential remedy is to partition the room layout into a finite number of characteristic zones, for which the ANN-GA model readily applies. Here, a multiple hot aisle/cold aisle data center configuration was analyzed using the commercial software FloTHERM. The CFD results are used to characterize the flow rates at the inter-zonal partitions. Based on specific reduced subsets of desired treatment quantities from the CFD results, such as CRAC and server rack air flow rates, the approach was applied for two different CRAC configurations and various levels of CRAC and server rack flow rates. Utilizing the compact inter-zonal boundary conditions, good agreement for the airflow and temperature distributions is achieved between predictions from the CFD computations for the entire room configuration and the reduced order zone-level model for different operating conditions and room layouts.


2018 ◽  
Vol 18 (20) ◽  
pp. 14813-14835 ◽  
Author(s):  
Liza I. Díaz-Isaac ◽  
Thomas Lauvaux ◽  
Kenneth J. Davis

Abstract. Atmospheric transport model errors are one of the main contributors to the uncertainty affecting CO2 inverse flux estimates. In this study, we determine the leading causes of transport errors over the US upper Midwest with a large set of simulations generated with the Weather Research and Forecasting (WRF) mesoscale model. The various WRF simulations are performed using different meteorological driver datasets and physical parameterizations including planetary boundary layer (PBL) schemes, land surface models (LSMs), cumulus parameterizations and microphysics parameterizations. All the different model configurations were coupled to CO2 fluxes and lateral boundary conditions from the CarbonTracker inversion system to simulate atmospheric CO2 mole fractions. PBL height, wind speed, wind direction, and atmospheric CO2 mole fractions are compared to observations during a month in the summer of 2008, and statistical analyses were performed to evaluate the impact of both physics parameterizations and meteorological datasets on these variables. All of the physical parameterizations and the meteorological initial and boundary conditions contribute 3 to 4 ppm to the model-to-model variability in daytime PBL CO2 except for the microphysics parameterization which has a smaller contribution. PBL height varies across ensemble members by 300 to 400 m, and this variability is controlled by the same physics parameterizations. Daily PBL CO2 mole fraction errors are correlated with errors in the PBL height. We show that specific model configurations systematically overestimate or underestimate the PBL height averaged across the region with biases closely correlated with the choice of LSM, PBL scheme, and cumulus parameterization (CP). Domain average PBL wind speed is overestimated in nearly every model configuration. Both planetary boundary layer height (PBLH) and PBL wind speed biases show coherent spatial variations across the Midwest, with PBLH overestimated averaged across configurations by 300–400 m in the west, and PBL winds overestimated by about 1 m s−1 on average in the east. We find model configurations with lower biases averaged across the domain, but no single configuration is optimal across the entire region and for all meteorological variables. We conclude that model ensembles that include multiple physics parameterizations and meteorological initial conditions are likely to be necessary to encompass the atmospheric conditions most important to the transport of CO2 in the PBL, but that construction of such an ensemble will be challenging due to ensemble biases that vary across the region.


Ocean Science ◽  
2017 ◽  
Vol 13 (2) ◽  
pp. 235-257 ◽  
Author(s):  
Reiner Onken

Abstract. The Regional Ocean Modeling System (ROMS) has been employed to explore the sensitivity of the forecast skill of mixed-layer properties to initial conditions, boundary conditions, and vertical mixing parameterisations. The initial and lateral boundary conditions were provided by the Mediterranean Forecasting System (MFS) or by the MERCATOR global ocean circulation model via one-way nesting; the initial conditions were additionally updated through the assimilation of observations. Nowcasts and forecasts from the weather forecast models COSMO-ME and COSMO-IT, partly melded with observations, served as surface boundary conditions. The vertical mixing was parameterised by the GLS (generic length scale) scheme Umlauf and Burchard (2003) in four different set-ups. All ROMS forecasts were validated against the observations which were taken during the REP14-MED survey to the west of Sardinia. Nesting ROMS in MERCATOR and updating the initial conditions through data assimilation provided the best agreement of the predicted mixed-layer properties with the time series from a moored thermistor chain. Further improvement was obtained by the usage of COSMO-ME atmospheric forcing, which was melded with real observations, and by the application of the k-ω vertical mixing scheme with increased vertical eddy diffusivity. The predicted temporal variability of the mixed-layer temperature was reasonably well correlated with the observed variability, while the modelled variability of the mixed-layer depth exhibited only agreement with the observations near the diurnal frequency peak. For the forecasted horizontal variability, reasonable agreement was found with observations from a ScanFish section, but only for the mesoscale wave number band; the observed sub-mesoscale variability was not reproduced by ROMS.


Author(s):  
Ali Pourfard ◽  
Esmaeel Khanmirza ◽  
Reza Madoliat

Simulation of a natural gas network operation is a prerequisite for optimization and control tasks. Treating gas in a transient manner is necessary for accurate simulation of gas networks. However, solving the governing nonlinear partial differential equations of pipe flows is a challenging task. In this paper, a novel approach is proposed based on using an intelligent algorithm called teaching–learning-based optimization. This approach simplifies transient simulation of gas networks with a specified type of boundary conditions. Teaching–learning-based optimization estimates different values for network inlet flow rates. Then by knowing the inlet boundary conditions of the network, the discretized flow equations become linear and the flow equations of each pipe can be solved independently. Thus, the network outlet flow variables can be easily obtained. The differences of obtained and actual network outlet flow rates are considered as a cost function or error. Finally, this intelligent algorithm determines the optimum inlet flow rates at each time level, which minimize the error. The proposed approach is implemented on the in-service gas network. To validate the simulation results, a conventional gradient-based method called trust region dogleg is also used for simulation of the gas network. The comparison of numerical results confirms the accuracy and efficiency of this approach, while it is more computationally efficient. Moreover, the substitution of teaching–learning-based optimization with another powerful intelligent optimization algorithm would not improve the performance of the proposed approach.


Sign in / Sign up

Export Citation Format

Share Document