Sensitivity of convection with an endothermic phase change to the form of governing equations, initial conditions, boundary conditions, and equation of state

1994 ◽  
Vol 99 (B8) ◽  
pp. 15919 ◽  
Author(s):  
Joel Ita ◽  
Scott D. King
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 13 (5) ◽  
pp. 2590
Author(s):  
S. A. M. Mehryan ◽  
Kaamran Raahemifar ◽  
Leila Sasani Gargari ◽  
Ahmad Hajjar ◽  
Mohamad El Kadri ◽  
...  

A Nano-Encapsulated Phase-Change Material (NEPCM) suspension is made of nanoparticles containing a Phase Change Material in their core and dispersed in a fluid. These particles can contribute to thermal energy storage and heat transfer by their latent heat of phase change as moving with the host fluid. Thus, such novel nanoliquids are promising for applications in waste heat recovery and thermal energy storage systems. In the present research, the mixed convection of NEPCM suspensions was addressed in a wavy wall cavity containing a rotating solid cylinder. As the nanoparticles move with the liquid, they undergo a phase change and transfer the latent heat. The phase change of nanoparticles was considered as temperature-dependent heat capacity. The governing equations of mass, momentum, and energy conservation were presented as partial differential equations. Then, the governing equations were converted to a non-dimensional form to generalize the solution, and solved by the finite element method. The influence of control parameters such as volume concentration of nanoparticles, fusion temperature of nanoparticles, Stefan number, wall undulations number, and as well as the cylinder size, angular rotation, and thermal conductivities was addressed on the heat transfer in the enclosure. The wall undulation number induces a remarkable change in the Nusselt number. There are optimum fusion temperatures for nanoparticles, which could maximize the heat transfer rate. The increase of the latent heat of nanoparticles (a decline of Stefan number) boosts the heat transfer advantage of employing the phase change particles.


2021 ◽  
Vol 37 ◽  
pp. 346-358
Author(s):  
Fuchun Yang ◽  
Xiaofeng Jiang ◽  
Fuxin Du

Abstract Free vibrations of rotating cylindrical shells with distributed springs were studied. Based on the Flügge shell theory, the governing equations of rotating cylindrical shells with distributed springs were derived under typical boundary conditions. Multicomponent modal functions were used to satisfy the distributed springs around the circumference. The natural responses were analyzed using the Galerkin method. The effects of parameters, rotation speed, stiffness, and ratios of thickness/radius and length/radius, on natural response were also examined.


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.


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.


2006 ◽  
Vol 5-6 ◽  
pp. 407-414 ◽  
Author(s):  
Mohammad Mohammadi Aghdam ◽  
M.R.N. Farahani ◽  
M. Dashty ◽  
S.M. Rezaei Niya

Bending analysis of thick laminated rectangular plates with various boundary conditions is presented using Generalized Differential Quadrature (GDQ) method. Based on the Reissner first order shear deformation theory, the governing equations include a system of eight first order partial differential equations in terms of unknown displacements, forces and moments. Presence of all plate variables in the governing equations provide a simple procedure to satisfy different boundary condition during application of GDQ method to obtain accurate results with relatively small number of grid points even for plates with free edges .Illustrative examples including various combinations of clamped, simply supported and free boundary condition are given to demonstrate the accuracy and convergence of the presented GDQ technique. Results are compared with other analytical and finite element predictions and show reasonably good agreement.


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.


Sign in / Sign up

Export Citation Format

Share Document