Venusian Applications of 3D Convection Modeling
NASA Technical Reports Server (NTRS)
Bonaccorso, Timary Annie
2011-01-01
This study models mantle convection on Venus using the 'cubed sphere' code OEDIPUS, which models one-sixth of the planet in spherical geometry. We are attempting to balance internal heating, bottom mantle viscosity, and temperature difference across Venus' mantle, in order to create a realistic model that matches with current planetary observations. We also have begun to run both lower and upper mantle simulations to determine whether layered (as opposed to whole-mantle) convection might produce more efficient heat transfer, as well as to model coronae formation in the upper mantle. Upper mantle simulations are completed using OEDIPUS' Cartesian counterpart, JOCASTA. This summer's central question has been how to define a mantle plume. Traditionally, we have defined a hot plume the region with temperature at or above 40% of the difference between the maximum and horizontally averaged temperature, and a cold plume as the region with 40% of the difference between the minimum and average temperature. For less viscous cases (1020 Pa?s), the plumes generated by that definition lacked vigor, displaying buoyancies 1/100th of those found in previous, higher viscosity simulations (1021 Pa?s). As the mantle plumes with large buoyancy flux are most likely to produce topographic uplift and volcanism, the low viscosity cases' plumes may not produce observable deformation. In an effort to eliminate the smallest plumes, we experimented with different lower bound parameters and temperature percentages.
Towards implementing plate tectonics in 3D mantle convection simulations
NASA Astrophysics Data System (ADS)
Bollada, Peter; Davies, Huw
2010-05-01
One of the great challenges in numerical mantle convection simulations is to achieve models that naturally develop plate tectonic like behaviour at the surface. In this work we are looking to achieve such models by investigating the set of models where a single consistent rheology is used for the whole model. We have started by investigating a viscoelastic rheology, related to the Oldroyd-B model from the field of polymers. The goal will be to have the parameter that controls the relaxation between elastic and viscous behaviour to depend upon temperature, pressure and strain-rate. With an appropriate choice of this dependence we have, on the near surface, high viscous/elastic regions interfaced with lower, pure viscous, regions of high strain-rate; while it also becomes more viscous at depth in the interior. In this way we hope to obtain plate like behaviour at the surface which naturally progresses to viscous convective behaviour in the interior. We have started to implement this model in the established mantle 3D finite element spherical mantle convection code TERRA (Baumgardner, 1984). Some parts of the model have been implemented as a force (to be combined with the gravitational body force) on the right hand side. The work has required us to develop and code in TERRA: (i) methods to overcome the continuity problem of the stress field stemming from the fact that the velocity field is represented by linear finite elements; (ii) new operators to handle stress and its gradients; (iii) methods to analyse plate-like behaviour at the surface (iv) the necessary functional dependence of viscosity and elastic relaxation time on temperature, strain-rate and pressure We will present the background to the work, its implementation and results.
3-D Spherical Mantle Convection with Radial Basis Functions
NASA Astrophysics Data System (ADS)
Flyer, N.; Wright, G. B.; Yuen, D.
2008-12-01
In the past 25 years a wide variety of numerical methods, such as finite-difference, finite-volume , finite- elements, and pseudospectral methods have been employed to study the problem of 3-D mantle convection. All have specialized strengths but also serious weaknesses. The first three methods are generally considered low-order and can involve high algorithmic complexity (as in triangular elements). Spectrally accurate methods do not practically allow for local mesh refinement and often involve cumbersome algebra. Here, we introduce a new grid/mesh-free approach using radial basis functions (RBFs). It has the advantage of being spectrally accurate for arbitrary node layouts in multi-dimensions with extreme algorithmic simplicity, and naturally permits local node refinement. It has been shown for shallow-water equations and vortex flows that RBFs outperform other numerical methods in the sense that they obtain a much higher accuracy for the same spatial resolution while being able to take unusually large time steps. One virtue of the RBF scheme is the ability to use a simple Cartesian geometry while implementing the required boundary conditions for the temperature, velocity and stresses on a spherical surface of both the outer( planetary surface ) and inner shell ( core-mantle boundary ). The velocity and stress components are expressed in terms of the scalar potential approach (Zebib and Schubert, 1982) and the other remaining variable is the perturbed temperature field. We have studied the problem from the onset of convection to a modest nonlinear regime.
Convective instability in sedimentation: 3-D numerical study
NASA Astrophysics Data System (ADS)
Yu, Xiao; Hsu, Tian-Jian; Balachandar, S.
2014-11-01
To provide a probable explanation on the field observed rapid sedimentation process near river mouths, we investigate the convective sedimentation in stably stratified saltwater using 3-D numerical simulations. Guided by the linear stability analysis, this study focuses on the nonlinear interactions of several mechanisms, which lead to various sediment finger patterns, and the effective settling velocity for sediment ranging from clay (single-particle settling velocity V0 = 0.0036 and 0.0144 mm/s, or particle diameter d = 2 and 4 μm) to silt (V0 = 0.36 mm/s, or d = 20 μm). For very fine sediment with V0 = 0.0036 mm/s, the convective instability is dominated by double diffusion, characterized by millimeter-scale fingers. Gravitational settling slightly increases the growth rate; however, it has notable effect on the downward development of vertical mixing shortly after the sediment interface migrates below the salt interface. For sediment with V0 = 0.0144 mm/s, Rayleigh-Taylor instabilities become dominant before double-diffusive modes grow sufficiently large. Centimeter-scale and highly asymmetric sediment fingers are obtained due to nonlinear interactions between different modes. For sediment with V0 = 0.36 mm/s, Rayleigh-Taylor mechanism dominates and the resulting centimeter-scale sediment fingers show a plume-like structure. The flow pattern is similar to that without ambient salt stratification. Rapid sedimentation with effective settling velocity on the order of 1 cm/s is likely driven by convective sedimentation for sediment with V0 greater than 0.1 mm/s at concentration greater than 10-20 g/L.
3D Simulations of methane convective storms on Titan's atmosphere
NASA Astrophysics Data System (ADS)
Hueso, R.; Sánchez-Lavega, A.
2005-08-01
The arrival of the Cassini/Huygens mission to Titan has opened an unprecedented opportunity to study the atmosphere of this satellite. Under the pressure-temperature conditions on Titan, methane, a large atmospheric component amounting perhaps to a 3-5% of the atmosphere, is close to its triple point, potentially playing a similar role as water on Earth. The Huygens probe has shown a terrain shaped by erosion of probably liquid origin, suggestive of past rain. On the other hand, Voyager IRIS spectroscopic observations of Titan imply a saturated atmosphere of methane (amounting perhaps to 150 covered by methane clouds, if we think on Earth meteorology. However, observations from Earth and Cassini have shown that clouds are localized, transient and fast evolving, in particular in the South Pole (currently in its summer season). This might imply a lack of widespread presence on Titan of nuclei where methane could initiate condensation and particle growth with subsequent precipitation. We investigate different scenarios of moist convective storms on Titan using a complete 3D atmospheric model that incorporates a full microphysics treatment required to study cloud formation processes under a saturated atmosphere with low concentration of condensation nuclei. We study local convective development under a variety of atmospheric conditions: sub-saturation, super-saturation, abundances of condensation nuclei fall, condensation nuclei lifted from the ground or gently falling from the stratosphere. We show that under the appropriate circumstances, precipitation rates comparable to typical tropical storms on Earth can be found. Acknowledgements: This work has been funded by Spanish MCYT PNAYA2003-03216, fondos FEDER and Grupos UPV 15946/2004. R. Hueso acknowledges a post-doc fellowship from Gobierno Vasco.
3d Modelling of Convective Flow In The Rhine Graben
NASA Astrophysics Data System (ADS)
Bächler, D.; Kohl, T.; Rybach, L.
Detailed investigations of the temperature distribution in the Rhine Graben indi- cate regular pattern of thermal anomalies following major north-south striking faults. These anomalies remain unexplained by conventional Rhine Graben studies based on 2D east-west striking sections. First analytical solutions for convective flow in vertical faults are applied for a clearly observable anomalous temperature pattern along ma- jor Rhine Graben faults. By these calculations the fault height, fault aperture, minimal fault permeability and time to convective onset is derived from the observed distances. Since analytical solutions are limited to simple model geometries further improvement was achieved by numerical model simulations, which allow to assume more com- plex initial and boundary conditions. Using the finite volume code TOUGH2 series of anomalies following the same fault were simulated by a 3D numerical model. Fo- cussing on the predominant north-south permeability structure the model consists of a vertical north-south striking fault and surrounding matrix. The fault geometries are based on the analytically predicted fault geometries (aperture=200m, height=3500m) and on the observed temperatures. Comparison of simulation results with observed temperatures shows that the fault is situated between 500 to 600m and 4200m. The fault permeability is taken as 5*10-13m2 and the fluid velocity in the fault is calcu- lated as 10-9 to 10-10 m/s. These results indicate the importance of our considerations since mass flux is much higher in the faults than across them. The minimal age of the anomaly is considered to be 77'000 years, since steady state is reached after this time span. The study proves that the observed temperature anomaly pattern along the gamma fault at Landau can be explained by north-south striking convection systems within fault zones. Similar situations have been found at Soultz. This may be a hint on a general feature of the major north-south striking
Computational study of 3-D Benard convection with gravitational modulation
NASA Technical Reports Server (NTRS)
Biringen, S.; Peltier, L. J.
1989-01-01
In this numerical study the effects of a modulated gravitational field on three-dimensional Rayleigh-Benard convection with heating from above or from below is investigated. The full, nonlinear, time-dependent, Boussinesq Navier-Stokes equations and the energy equation are solved by a semiimplicit, pseudo-spectral procedure. This study has been motivated by the need to better understand the effects of vibration (G-Jitter) on fluids systems especially in the low gravity environment.
Numerical simulation of 3-D Benard convection with gravitational modulation
NASA Technical Reports Server (NTRS)
Biringen, S.; Peltier, L. J.
1990-01-01
In this numerical study, randomly and sinusoidally modulated gravitational fields imposed on three-dimensional Rayleigh-Benard convection are investigated in an effort to understand the effects of vibration (G-Jitter) on fluid systems. The time-dependent, Navier-Stokes equations and the energy equation with Boussinesq approximations are solved by a semi-implicit, pseudospectral procedure. An analysis of energy balances indicates that with increasing modulation amplitude, transition from synchronous to relaxation oscillation goes through the subharmonic response. Random modulations are found to be less stabilizing than sinusoidal and are shown to impose three-dimensionality on the flow for some parameter ranges both at terrestrial and zero base gravity conditions.
3-D simulations of solar and stellar convection and magnetoconvection
NASA Astrophysics Data System (ADS)
Nordlund, Å.; Stein, R. F.
1990-05-01
We present the key components of a 3-D code designed for simulating the hydrodynamics and magnetohydrodynamics of stellar atmospheres and envelopes. Some particular properties of the code are: (1) the ability to handle strong stratification (extensive simulations with bottom/top pressure ratios of 3×104 have been performed, and simulations with pressure ratios of 5×106 are being initiated); (2) a detailed treatment of the radiating surface; (3) a functional form of the subgrid-scale diffusion designed to minimize the influence on resolved motions; (4) boundary conditions open to flows. The top boundary allows the transmission of short period waves, while the bottom boundary condition was designed to enforce a displacement node for radial pressure modes.
Heat distribution by natural convection
Balcomb, J.D.
1985-01-01
Natural convection can provide adequate heat distribution in many situtations that arise in buildings. This is appropriate, for example, in passive solar buildings where some rooms tend to be more strongly solar heated than others or to reduce the number of heating units required in a building. Natural airflow and heat transport through doorways and other internal building apertures is predictable and can be accounted for in the design. The nature of natural convection is described, and a design chart is presented appropriate to a simple, single-doorway situation. Natural convective loops that can occur in buildings are described and a few design guidelines are presented.
NASA Astrophysics Data System (ADS)
Liu, X.; Zhong, S.
2013-12-01
Classic mantle dynamic models for the Earth's geoid are mostly based on whole mantle convection and constrain that the upper mantle is significantly weaker than the lower mantle. Whole mantle convection models with such mantle viscosity structure have successfully explained the long-wavelength structure in the mantle. However, with increasing consensus on the existence of chemically distinct piles above the core mantle boundary (CMB) (also known as large low shear velocity provinces or LLSVPs), questions arise as to what extent the chemical piles influence the Earth's geoid and long-wavelength mantle convection. Some recent studies suggested that the chemical piles have a controlling effect on the Earth's degree two mantle structure, geoid, and true polar wander, although the chemical piles are estimated to be of small volume (~2% of the whole mantle) by seismic studies. We have formulated dynamically consistent 3D mantle convection models using CitcomS and studied how the chemical piles above CMB influence the long-wavelength convective structure and geoid. The models have free slip boundary conditions and temperature dependent viscosity. By comparing with purely thermal convection models, we found that the long wavelength convective structure is not sensitive to the presence of the chemical piles. By determining the geoid from the buoyance of a certain layer of the mantle, we found that for both purely thermal and thermochemical convection, the geoid is mostly contributed by the upper part of the mantle, with ~80% geoid explained by the buoyancy in the upper half of the mantle. In purely thermal convection, the contribution to the geoid from the bottom layer of the mantle always has the same sign with the total geoid (a bottom ~ 600 km thick layer gives ~3.5% of the total geoid). However, in the thermochemical convection, the bottom layer with overall negatively buoyant chemical piles gives rise to the geoid that has opposite sign with the total geoid and has a
Similarities between 2D and 3D convection for large Prandtl number
NASA Astrophysics Data System (ADS)
Pandey, Ambrish; Verma, Mahendra K.; Chatterjee, Anando G.; Dutta, Biplab
2016-06-01
Using direct numerical simulations of Rayleigh-B\\'{e}nard convection (RBC), we perform a comparative study of the spectra and fluxes of energy and entropy, and the scaling of large-scale quantities for large and infinite Prandtl numbers in two (2D) and three (3D) dimensions. We observe close similarities between the 2D and 3D RBC, in particular the kinetic energy spectrum $E_u(k) \\sim k^{-13/3}$, and the entropy spectrum exhibits a dual branch with a dominant $k^{-2}$ spectrum. We showed that the dominant Fourier modes in the 2D and 3D flows are very close. Consequently, the 3D RBC is quasi two-dimensional, which is the reason for the similarities between the 2D and 3D RBC for large- and infinite Prandtl numbers.
High Rayleigh Number 3-D Spherical Mantle Convection with Radial Basis Functions
NASA Astrophysics Data System (ADS)
Flyer, N.; Yuen (3), G. Wright, D.
2009-04-01
In the last quarter of a century many numerical methods, such as finite-differences, finite-volume, their yin-yang variants, finite-elements and pseudo-spectral methods have been used to study the problem of 3-D spherical convection. All have their respective strengths but also serious weaknesses, such as low-order and can involve high algorithmic complexity, as in triangular elements. Spectrally accurate methods do not practically allow for local mesh refinement and often involve cumbersome algebra. We have recently introduced a new grid/mesh-free approach, using radial basis functions ( RBFs) . It has the advantage of being spectrally accurate for arbitrary node layouts in multi-dimensions with extreme algorithmic simplicity, and allows naturally node-refinement. One virtue of the RBF scheme is the ability to use a simple Cartesian geometry while implementing the required boundary conditions for the temperature, velocity and stresses on a spherical surface of both the outer( planetary surface ) and inner shell ( core-mantle boundary ). The velocity and stress components are expressed in terms of the scalar potential approach and the other remaining variable is the perturbed temperature field. We have studied the problem from the weakly nonlinear to a moderately nonlinear regime involving a Rayleigh number, about 1000 times super-critical. Both purely basal and partially internal -heating cases have been considered
High Rayleigh Number 3-D Spherical Mantle Convection with Radial Basis Functions
NASA Astrophysics Data System (ADS)
Flyer, N.; Wright, G.; Yuen, D. A.
2009-04-01
In the last quarter of a century many numerical methods, such as finite-differences, finite-volume, their yin-yang variants, finite-elements and pseudo-spectral methods have been used to study the problem of 3-D spherical convection. All have their respective strengths but also serious weaknesses, such as low-order and can involve high algorithmic complexity, as in triangular elements. Spectrally accurate methods do not practically allow for local mesh refinement and often involve cumbersome algebra. We have recently introduced a new grid/mesh-free approach, using radial basis functions (RBFs). It has the advantage of being spectrally accurate for arbitrary node layouts in multi-dimensions with extreme algorithmic simplicity, and allows naturally node-refinement. One virtue of the RBF scheme is the ability to use a simple Cartesian geometry while implementing the required boundary conditions for the temperature, velocity and stresses on a spherical surface of both the outer(planetary surface) and inner shell (core-mantle boundary). The velocity and stress components are expressed in terms of the scalar potential approach and the other remaining variable is the perturbed temperature field. We have studied the problem from the weakly onlinear to a moderately nonlinear regime involving a Rayleigh number, about 1000 times super-critical. Both purely basal and partially internal-heating cases have been considered.
New Era in 3-D Modeling of Convection and Magnetic Dynamos in Stellar Envelopes and Cores
NASA Astrophysics Data System (ADS)
Toomre, J.; Augustson, K. C.; Brown, B. P.; Browning, M. K.; Brun, A. S.; Featherstone, N. A.; Miesch, M. S.
2012-09-01
The recent advances in asteroseismology and spectropolarimetry are beginning to provide estimates of differential rotation and magnetic structures for a range of F and G-type stars possessing convective envelopes, and in A-type stars with convective cores. It is essential to complement such observational work with theoretical studies based on 3-D simulations of highly turbulent convection coupled to rotation, shear and magnetic fields in full spherical geometries. We have so employed the anelastic spherical harmonic (ASH) code, which deals with compressible magnetohydrodynamics (MHD) in spherical shells, to examine the manner in which the global-scale convection can establish differential rotation and meridional circulations under current solar rotation rates, and these make good contact with helioseismic findings. For younger G stars rotating 3 to 5 times faster than the current Sun, the convection establishes ever stronger angular velocity contrasts between their fast equators and slow poles, and these are accompanied by prominent latitudinal temperature contrasts as well. Turning to MHD simulation of magnetic dynamo action within these younger G stars, the resulting magnetism involves wreaths of strong toroidal magnetic fields (up to 50 to 100 kG strengths) in the bulk of the convection zone, typically of opposite polarity in the northern and southern hemispheres. These fields can persist for long intervals despite being pummeled by the fast convective downflows, but they can also exhibit field reversals and cycles. Turning to shallower convective envelopes in the more luminous F-type stars that range in mass from 1.2 to 1.4 solar masses and for various rotation rates, we find that the convection can again establish solar-like differential rotation profiles with a fast equator and slow poles, but the opposite is achieved at the slower rotation rates. The F stars are also capable of building strong magnetic fields, often as wreaths, through dynamo action. We also
NASA Astrophysics Data System (ADS)
Wright, G. B.; Barnett, G. A.; Yuen, D. A.
2009-12-01
, 533, 1984. Isosurfaces of the temperature field from a 3-D mantle convection simulation at Rayleigh number 10**7 during the transition from a purely conductive state to a double-layer convection state. Simulation was performed using the compact fourth order finite difference scheme at a resolution of 200-by-200-by-100 (length-by-width-by-height).
Fast and Robust Sixth Order Multigrid Computation for 3D Convection Diffusion Equation.
Wang, Yin; Zhang, Jun
2010-10-15
We present a sixth order explicit compact finite difference scheme to solve the three dimensional (3D) convection diffusion equation. We first use multiscale multigrid method to solve the linear systems arising from a 19-point fourth order discretization scheme to compute the fourth order solutions on both the coarse grid and the fine grid. Then an operator based interpolation scheme combined with an extrapolation technique is used to approximate the sixth order accurate solution on the fine grid. Since the multigrid method using a standard point relaxation smoother may fail to achieve the optimal grid independent convergence rate for solving convection diffusion equation with a high Reynolds number, we implement the plane relaxation smoother in the multigrid solver to achieve better grid independency. Supporting numerical results are presented to demonstrate the efficiency and accuracy of the sixth order compact scheme (SOC), compared with the previously published fourth order compact scheme (FOC). PMID:21151737
3D coexisting modes of thermal convection in the faulted Lower Yarmouk Gorge
NASA Astrophysics Data System (ADS)
Magri, Fabien; Inbar, Nimrod; Möller, Peter; Raggad, Marwan; Rödiger, Tino; Rosenthal, Eliyahu; Siebert, Christian
2016-04-01
Numerical investigations of 3D modes of large-scale convection in faulted aquifers are presented with the aim to infer possible transport mechanisms supporting the formation of thermal springs in the Lower Yarmouk Gorge (LYG), at the border between Israel and Jordan. The transient finite elements models are based on a geological model of the LYG that introduces more realistic structural features of the basin, compared to previous existing models of the area (Magri et al., submitted). The sensitivity analysis of the fault permeability showed that faults cross-cutting the main regional flow direction allow groundwater to be driven laterally by convective forces within the fault planes. Therein thermal waters can either discharge along the fault traces or exit the fault through adjacent permeable aquifers. The location of springs can migrate with time, is not strictly constrained to the damage zones and reflects the interplay between the wavelength of the multicellular regime in the fault zone and the regional flow toward discharge areas in the lowlands. The results presented here suggest that in the LYG case, crossing flow paths result from the coexistence of fault convection, that can develop for example along NE-SW oriented faults within the Gorge, and additional flow fields that can be induced either by topography N-S gradients, e.g. perpendicular to the major axe of the Gorge, or by local thermal convection in permeable aquifers below Eocene aquiclude. The sensitivity analysis is consistent with the analytical solutions based on viscous-dependent Rayleigh theory. It indicates that in the LYG coexisting transport processes likely occur at fault hydraulic conductivity ranging between 2.3e-7 m/s and 9.3e- 7 m/s (i.e. 7 m/yr and 30 m/yr). The LYG numerical example and the associated Rayleigh analysis can be applied to study the onset of thermal convection and resulting flow patterns of any fractured hydrothermal basin. References Magri F, Möller S, Inbar N, M
NASA Astrophysics Data System (ADS)
Yanagisawa, Takatoshi; Yamagishi, Yasuko; Hamano, Yozo; Stegman, Dave R.; Suetsugu, Daisuke; Bina, Craig; Inoue, Toru; Wiens, Douglas; Jellinek, Mark
2010-11-01
Seismic tomography reveals the natural mode of convection in the Earth is whole mantle with subducted slabs clearly seen as continuous features into the lower mantle. However, simultaneously existing alongside these deep slabs are stagnant slabs which are, if only temporarily, trapped in the upper mantle. Previous numerical models of mantle convection have observed a range of behavior for slabs in the transition zone depending on viscosity stratification and mineral phase transitions, but typically only exhibit flat-lying slabs when mantle convection is layered or trench migration is imposed. We use 3-D spherical models of mantle convection which range up to Earth-like conditions in Rayleigh number to systematically investigate three effects on mantle dynamics: (1) the mineral phase transitions, (2) a strongly temperature-dependent viscosity with plastic yielding at shallow depth, and (3) a viscosity increase in the lower mantle. First a regime diagram is constructed for isoviscous models over a wide range of Rayleigh number and Clapeyron slope for which the convective mode is determined. It agrees very well with previous results from 2-D simulations by Christensen and Yuen (1985), suggesting present-day Earth is in the intermittent convection mode rather than layered or strictly whole mantle. Two calculations at Earth-like conditions (Ra and RaH = 2 í 107 and 5 í 108, respectively) which include effects (2) and (3) are produced with and without the effect of the mineral phase transitions. The first calculation (without the phase transition) successfully produces plate-like behavior with a long wavelength structure and surface heat flow similar to Earth's value. While the observed convective flow pattern in the lower mantle is broader compared to isoviscous models, it basically shows the behavior of whole mantle convection, and does not exhibit any slab flattening at the viscosity increase at 660 km depth. The second calculation which includes the phase
NASA Astrophysics Data System (ADS)
Oh, J.; Jiang, X.; Waliser, D. E.; Moncrieff, M. W.; Johnson, R. H.
2013-12-01
As one of the most prominent tropical atmospheric variability modes, the Madden-Julian Oscillation (MJO) exerts profound influences on global weather and climate, and serves as a critical predictability source for extend-range forecast. While credible representation of the MJO still represents a great challenge for current general circulation models (GCMs), previous studies on the vertical structure of the MJO have largely focused on collective impacts from multi-scale convective systems on thermodynamic properties of the MJO. Most recently, limited observational studies and idealized modeling work suggested that convective momentum transport (CMT) could also play an important role in interpreting the observed MJO features. In this study, the 3D CMT structure associated with the MJO is examined by analyzing model output from three recent high-quality reanalysis systems, including NOAA's Climate Forecast System Reanalysis (CFSR), NASA's Modern Era Retrospective-analysis for Research and Applications (MERRA), and ECMWF-the Year of Tropical Convection (YOTC) reanalysis. Consistent with previous cloud-resolving model study, a well-organized three-layer vertical structure in the CMT associated with the MJO is also discerned based on reanalyses. The result suggests that CMT tends to intensify the MJO circulation, particularly in the lower troposphere. Relative roles of meso-scale systems (MCS) and synoptic waves in contributing the total CMT profiles of the MJO will also be explored. Differences in CMT profiles in these several reanalysis models will be discussed.
3D modelling of coupled mass and heat transfer of a convection-oven roasting process.
Feyissa, Aberham Hailu; Gernaey, Krist V; Adler-Nissen, Jens
2013-04-01
A 3D mathematical model of coupled heat and mass transfer describing oven roasting of meat has been developed from first principles. The proposed mechanism for the mass transfer of water is modified and based on a critical literature review of the effect of heat on meat. The model equations are based on a conservation of mass and energy, coupled through Darcy's equations of porous media - the water flow is mainly pressure-driven. The developed model together with theoretical and experimental assessments were used to explain the heat and water transport and the effect of the change in microstructure (permeability, water binding capacity and elastic modulus) that occur during the meat roasting process. The developed coupled partial differential equations were solved by using COMSOL Multiphysics®3.5 and state variables are predicted as functions of both position and time. The proposed mechanism was partially validated by experiments in a convection oven where temperatures were measured online. PMID:23305831
Natural convection in porous media
Prasad, V.; Hussain, N.A.
1986-01-01
This book presents the papers given at a conference on free convection in porous materials. Topics considered at the conference included heat transfer, nonlinear temperature profiles and magnetic fields, boundary conditions, concentrated heat sources in stratified porous media, free convective flow in a cavity, heat flux, laminar mixed convection flow, and the onset of convection in a porous medium with internal heat generation and downward flow.
Interactive Visualization of 3-D Mantle Convection Extended Through AJAX Applications
NASA Astrophysics Data System (ADS)
McLane, J. C.; Czech, W.; Yuen, D.; Greensky, J.; Knox, M. R.
2008-12-01
We have designed a new software system for real-time interactive visualization of results taken directly from large-scale simulations of 3-D mantle convection and other large-scale simulations. This approach allows for intense visualization sessions for a couple of hours as opposed to storing massive amounts of data in a storage system. Our data sets consist of 3-D data for volume rendering with over 10 million unknowns at each timestep. Large scale visualization on a display wall holding around 13 million pixels has already been accomplished with extension to hand-held devices, such as the OQO and Nokia N800 and recently the iPHONE. We are developing web-based software in Java to extend the use of this system across long distances. The software is aimed at creating an interactive and functional application capable of running on multiple browsers by taking advantage of two AJAX-enabled web frameworks: Echo2 and Google Web Toolkit. The software runs in two modes allowing for a user to control an interactive session or observe a session controlled by another user. Modular build of the system allows for components to be swapped out for new components so that other forms of visualization could be accommodated such as Molecular Dynamics in mineral physics or 2-D data sets from lithospheric regional models.
Natural convection between concentric spheres
NASA Technical Reports Server (NTRS)
Garg, Vijay K.
1992-01-01
A finite-difference solution for steady natural convective flow in a concentric spherical annulus with isothermal walls has been obtained. The stream function-vorticity formulation of the equations of motion for the unsteady axisymmetric flow is used; interest lying in the final steady solution. Forward differences are used for the time derivatives and second-order central differences for the space derivatives. The alternating direction implicit method is used for solution of the discretization equations. Local one-dimensional grid adaptation is used to resolve the steep gradients in some regions of the flow at large Rayleigh numbers. The break-up into multi-cellular flow is found at high Rayleigh numbers for air and water, and at significantly low Rayleigh numbers for liquid metals. Excellent agreement with previous experimental and numerical data is obtained.
Modelling natural convection of fluid in cuvette
NASA Astrophysics Data System (ADS)
Kucher, D.; Manukhin, B.; Andreeva, O.; Chivilikhin, S.
2014-09-01
Convection is a process of transfer liquid from a hot region to a cool region. This phenomenon is involved in many physical processes. The main characteristic of convection is a temperature field. Modelling of convection allows to get the information about temperature field at any time of process. In this paper the results of modelling natural convection of fluid in cuvette are presented. All results are approved by experimental data. For modelling the process of natural convection Navier-Stokes equations under Boussinesq approximation were used. An experimental setup based on digital holographic interferometry was developedin order to make an experiment. The results for three stadiums of convection, such as: jet initiation, initial jet formation, jet development with formation of mushroom-shaped convective stream, are presented.
NASA Astrophysics Data System (ADS)
Antonini Alves, Thiago; Santos, Paulo H. D.; Barbur, Murilo A.
2015-09-01
In this research, the temperatures of threedimensional (3D) protruding heaters mounted on a conductive substrate in a horizontal rectangular channel with laminar airflow are related to the independent power dissipation in each heater by using a matrix G + with invariant coefficients, which are dimensionless. These coefficients are defined in this study as the conjugate influence coefficients ( g +) caused by the forced convection- conduction nature of the heaters' cooling process. The temperature increase of each heater in the channel is quantified to clearly identify the contributions attributed to the self-heating and power dissipation in the other heaters (both upstream and downstream). The conjugate coefficients are invariant with the heat generation rate in the array of heaters when assuming a defined geometry, invariable fluid and flow rate, and constant substrate and heater conductivities. The results are numerically obtained by considering three 3D protruding heaters on a twodimensional (2D) array by ANSYS/Fluent™ 15.0 software. The conservation equations are solved by a coupled procedure within a single calculation domain comprising of solid and fluid regions and by considering a steady state laminar airflow with constant properties. Some examples are shown, indicating the effects of substrate thermal conductivity and Reynolds number on conjugate influence coefficients.
NASA Astrophysics Data System (ADS)
Crameri, Fabio; Tackley, Paul J.
2014-07-01
We present temporally evolving 3-D global mantle convection models with single-sided subduction and a free surface in both 3-D Cartesian and fully spherical geometry. Special focus is given to the spontaneous development of three-dimensional structures at the surface and in the upper mantle. We find that an arcuate shape is the natural form for trenches and slabs. Cartesian models are used first to study the dynamic evolution of subduction zones, spreading ridges, and interconnected transform features. These experiments highlight the strong variation of spontaneously developing, arcuate slab curvature and subduction polarity along the trench strike. The spontaneous development of spreading ridges leads to lateral offsets between separated segments that are characterized by normal transform motion. Spherical models then allow insights into the evolution of plate tectonics on a sphere. Investigated are the spontaneous evolution of slab geometry, trench motion, and subduction-induced mantle flow. Two new dynamical features are discovered: "back-slab spiral flow" and "slab tunneling." 2014. American Geophysical Union. All Rights Reserved.
Natural convection in low-g environments
NASA Technical Reports Server (NTRS)
Grodzka, P. G.; Bannister, T. C.
1974-01-01
The present state of knowledge in the area of low-g natural convection is reviewed, taking into account a number of experiments conducted during the Apollo 14, 16, and 17 space flights. Convections due to steady low-g accelerations are considered. Steady g-levels result from spacecraft rotation, gravity gradients, solar wind, and solar pressure. Varying g-levels are produced by engine burns, attitude control maneuvers, and onboard vibrations from machinery or astronaut movement. Thermoacoustic convection in a low-g environment is discussed together with g-jitter convection, surface tension-driven convection, electrohydrodynamics under low-g conditions, phase change convection, and approaches for the control and the utilization of convection in space.
NASA Astrophysics Data System (ADS)
Magri, Fabien; Möller, Sebastian; Inbar, Nimrod; Siebert, Christian; Möller, Peter; Rosenthal, Eliyahu; Kühn, Michael
2015-04-01
It has been shown that thermal convection in faults can also occur for subcritical Rayleigh conditions. This type of convection develops after a certain period and is referred to as "delayed convection" (Murphy, 1979). The delay in the onset is due to the heat exchange between the damage zone and the surrounding units that adds a thermal buffer along the fault walls. Few numerical studies investigated delayed thermal convection in fractured zones, despite it has the potential to transport energy and minerals over large spatial scales (Tournier, 2000). Here 3D numerical simulations of thermally driven flow in faults are presented in order to investigate the impact of delayed convection on deep fluid processes at basin-scale. The Tiberias Basin (TB), in the Jordan Rift Valley, serves as study area. The TB is characterized by upsurge of deep-seated hot waters along the faulted shores of Lake Tiberias and high temperature gradient that can locally reach 46 °C/km, as in the Lower Yarmouk Gorge (LYG). 3D simulations show that buoyant flow ascend in permeable faults which hydraulic conductivity is estimated to vary between 30 m/yr and 140 m/yr. Delayed convection starts respectively at 46 and 200 kyrs and generate temperature anomalies in agreement with observations. It turned out that delayed convective cells are transient. Cellular patterns that initially develop in permeable units surrounding the faults can trigger convection also within the fault plane. The combination of these two convective modes lead to helicoidal-like flow patterns. This complex flow can explain the location of springs along different fault traces of the TB. Besides being of importance for understanding the hydrogeological processes of the TB (Magri et al., 2015), the presented simulations provide a scenario illustrating fault-induced 3D cells that could develop in any geothermal system. References Magri, F., Inbar, N., Siebert, C., Rosenthal, E., Guttman, J., Möller, P., 2015. Transient
Convection and chemistry effects in CVD: A 3-D analysis for silicon deposition
NASA Technical Reports Server (NTRS)
Gokoglu, S. A.; Kuczmarski, M. A.; Tsui, P.; Chait, A.
1989-01-01
The computational fluid dynamics code FLUENT has been adopted to simulate the entire rectangular-channel-like (3-D) geometry of an experimental CVD reactor designed for Si deposition. The code incorporated the effects of both homogeneous (gas phase) and heterogeneous (surface) chemistry with finite reaction rates of important species existing in silane dissociation. The experiments were designed to elucidate the effects of gravitationally-induced buoyancy-driven convection flows on the quality of the grown Si films. This goal is accomplished by contrasting the results obtained from a carrier gas mixture of H2/Ar with the ones obtained from the same molar mixture ratio of H2/He, without any accompanying change in the chemistry. Computationally, these cases are simulated in the terrestrial gravitational field and in the absence of gravity. The numerical results compare favorably with experiments. Powerful computational tools provide invaluable insights into the complex physicochemical phenomena taking place in CVD reactors. Such information is essential for the improved design and optimization of future CVD reactors.
Coupled three-dimensional conduction and natural convection heat transfer
NASA Astrophysics Data System (ADS)
Tolpadi, Anil Kumar
1987-09-01
A numerical and experimental investigation of three-dimensional natural convection heat transfer coupled with conduction was performed. This general problem is of great importance because of its widespread applicability in areas such as compact natural convection heat exchangers, cooling of electronic equipment, and porous media flows. The determination of flow patterns and heat transfer coefficients in such situations is necessary because of its practical use in various industries. A vectorized finite difference code was developed for the Cray-2 supercomputer which has the capability of simulating a wide class of three-dimensional coupled conduction-convection problems. This program numerically solves the transient form of the complete laminar Navier-Stokes equations of motion using the vorticity-vector potential methods. Using this program, numerical solutions were obtained for 3-D natural convection from a horizontal isothermal heat exchanger tube with an attached circular cooling fin array. Experiments were performed to measure three-dimensional temperature fields using Mach-Zehnder interferometry. Software was developed to digitize and process fringe patterns and inversion algorithms used to compute the 3-D temperature field.
Natural convection in nonvertical wells
Denbow, D.A.; Murphy, H.D.; McEligot, D.M.
1985-01-01
Convective instabilities and the shapes of the ensuing convection cells were experimentally studied for nonvertical wellbores. Steady-state temperature distributions were measured for three inclination angles over a wide range of heating rates to demonstrate the effects of drilling angle and Rayleigh number. In addition, velocities were estimated by measuring the time-of-flight of tracers formed by the Thymol blue technique. 8 refs., 6 figs.
Interactive Visualization and Monitoring of Large-Scale 3-D Mantle Convection Runs
NASA Astrophysics Data System (ADS)
Damon, M.; Yuen, D.; Kameyama, M.; Knox, M.; Porter, D.; Sevre, E. O.; Woodward, P.
2007-12-01
With the imminent arrival of petascale computing in the United States by 2011, new strategies for visualizing and monitoring high-resolution numerical simulations on massively parallel computers are needed to overcome the extreme data and resource requirements. We have employed a visualization system consisting of 14 powerful Dell workstations, each with a multi-terabyte disk, connected via a high-speed network with a bandwidth on the order of a few gigabits per second to a locally situated massively parallel system with approximately 2,000 processing elements. This system has been constructed at the Laboratory of Computational Sciences and Engineering at the University of Minnesota. Near real-time interactive analysis of 3-D mantle convection using around 10 million grid points has been carried out using a client-server application capable of streaming gigabytes of simulated data to a remote Powerwall with 13 million pixels. Concurrently, we have constructed a web-portal that allows a user to monitor the same run at home or in a hotel room, using a laptop. In our case, interactive computing takes on the meaning of performing such runs for a limited duration of time, say 1 to 2 hours. This calls for a balance between grid resolution and the number of processing elements required to provide the level of interactivity needed to achieve one to a few frames per second. Our mode of operation represents a new paradigm in numerical modeling that supports a trend toward both real-time visualization and monitoring of high-resolution models and a consequent reduction in storage of raw output data, since the interactive periods are by definition short. Using this interactive strategy periodically we can facilitate long heroic runs extending over a few days.
3D recovery of human gaze in natural environments
NASA Astrophysics Data System (ADS)
Paletta, Lucas; Santner, Katrin; Fritz, Gerald; Mayer, Heinz
2013-01-01
The estimation of human attention has recently been addressed in the context of human robot interaction. Today, joint work spaces already exist and challenge cooperating systems to jointly focus on common objects, scenes and work niches. With the advent of Google glasses and increasingly affordable wearable eye-tracking, monitoring of human attention will soon become ubiquitous. The presented work describes for the first time a method for the estimation of human fixations in 3D environments that does not require any artificial landmarks in the field of view and enables attention mapping in 3D models. It enables full 3D recovery of the human view frustum and the gaze pointer in a previously acquired 3D model of the environment in real time. The study on the precision of this method reports a mean projection error ≈1.1 cm and a mean angle error ≈0.6° within the chosen 3D model - the precision does not go below the one of the technical instrument (≈1°). This innovative methodology will open new opportunities for joint attention studies as well as for bringing new potential into automated processing for human factors technologies.
3D printing of natural organic materials by photochemistry
NASA Astrophysics Data System (ADS)
Da Silva Gonçalves, Joyce Laura; Valandro, Silvano Rodrigo; Wu, Hsiu-Fen; Lee, Yi-Hsiung; Mettra, Bastien; Monnereau, Cyrille; Schmitt Cavalheiro, Carla Cristina; Pawlicka, Agnieszka; Focsan, Monica; Lin, Chih-Lang; Baldeck, Patrice L.
2016-03-01
In previous works, we have used two-photon induced photochemistry to fabricate 3D microstructures based on proteins, anti-bodies, and enzymes for different types of bio-applications. Among them, we can cite collagen lines to guide the movement of living cells, peptide modified GFP biosensing pads to detect Gram positive bacteria, anti-body pads to determine the type of red blood cells, and trypsin columns in a microfluidic channel to obtain a real time biochemical micro-reactor. In this paper, we report for the first time on two-photon 3D microfabrication of DNA material. We also present our preliminary results on using a commercial 3D printer based on a video projector to polymerize slicing layers of gelatine-objects.
NASA Technical Reports Server (NTRS)
Shie, Chung-Lin; Tao, Wei-Kuo; Simpson, Joanne
2003-01-01
The 1999 Kwajalein Atoll field experiment (KWAJEX), one of several major TRMM (Tropical Rainfall Measuring Mission) field experiments, has successfully obtained a wealth of information and observation data on tropical convective systems over the western Central Pacific region. In this paper, clouds and convective systems that developed during three active periods (Aug 7-12, Aug 17-21, and Aug 29-Sep 13) around Kwajalein Atoll site are simulated using both 2D and 3D Goddard Cumulus Ensemble (GCE) models. Based on numerical results, the clouds and cloud systems are generally unorganized and short lived. These features are validated by radar observations that support the model results. Both the 2D and 3D simulated rainfall amounts and their stratiform contribution as well as the heat, water vapor, and moist static energy budgets are examined for the three convective episodes. Rainfall amounts are quantitatively similar between the two simulations, but the stratiform contribution is considerably larger in the 2D simulation. Regardless of dimension, fo all three cases, the large-scale forcing and net condensation are the two major physical processes that account for the evolution of the budgets with surface latent heat flux and net radiation solar and long-wave radiation)being secondary processes. Quantitative budget differences between 2D and 3D as well as between various episodes will be detailed.Morover, simulated radar signatures and Q1/Q2 fields from the three simulations are compared to each other and with radar and sounding observations.
NASA Astrophysics Data System (ADS)
Passos, D.; Charbonneau, P.; Miesch, M. S.
2016-04-01
The solar meridional circulation is a "slow", large scale flow that transports magnetic field and plasma throughout the convection zone in the (r,θ) plane and plays a crucial role in controlling the magnetic cycle solutions presented by flux transport dynamo models. Observations indicate that this flow speed varies in anti-phase with the solar cycle at the solar surface. A possible explanation for the source of this variation is based on the fact that inflows into active regions alter the global surface pattern of the meridional circulation. In this work we examine the meridional circulation profile that emerges from a 3D global simulation of the solar convection zone, and its associated dynamics. We find that at the bottom of the convection zone, in the region where the toroidal magnetic field accumulates, the meridional circulation is highly modulated through the action of a magnetic torques and thus provides evidence for a new mechanism to explain the observed cyclic variations.
Onset of Natural Convection in Saline Aquifers
NASA Astrophysics Data System (ADS)
Riaz, A.
2013-05-01
Sequestration of carbon dioxide in saline aquifers has emerged as the preferred method of permanently storing CO2 in the subsurface. In order to be successful over geologic time scales, sequestration in saline aquifers relies upon enhanced dissolution of CO2 in brine by natural convection. In this talk we review the progress made thus far towards the modeling and prediction of the onset time for natural convection that occurs due to an unstable stratification of aqueous CO2. We show how the onset of natural convection is connected to a preceding event of the onset of instability with respect to small amplitude perturbations that originate within the aqueous boundary layer. Our analysis indicates that the onset time for instability is uncertain within an initial transient period where perturbation growth depends on the specific form of the initial condition. A constrained adjoint based optimization is employed to determine the upper bound and the mean of perturbation growth. With the help of a weakly nonlinear analysis, we show that the time at which convection initiates is associated with fixed perturbation amplitude. The influence of permeability heterogeneity is studied with this approach. For certain permeability structures, the marginal stability curve bifurcates to form multiple stable and unstable zones in the space of the perturbation wavenumber and time. The transition toward bifurcation governs the behavior of the most dangerous mode in the linear regime and determines the route to the onset of natural convection.
Rockfall assessment of natural hazards by 3D-simulation
NASA Astrophysics Data System (ADS)
Tobler, D.; Graf, K.; Krummenacher, B.
2009-04-01
Keywords: 3D rockfall modelling, protection measures, case study, ROFMOD 4.1 model In mountain areas many residential areas as well as important lifelines are exposed to potential rockfall events. The assessment of the hazard potential in a complex morphology with changing surface parameters and inhomogeneous forestation, call for an in-depth analysis of the physical processes involved. Therefore the different problems concerning the assessment of rockfall hazards have to be solved adequately. Moreover, due to the actual discussion about cost efficiency the planning of protective measures calls for (more) detailed information about the intensity and probability of expected rockfall incidents along a typical trajectory or for a given area. Today simulation models like GEOTEST's - ROFMOD 4.1 (2D and 3D) are important tools for hazard assessment. This simulation model produces profiles and maps with intensity and range and allows to calculate decisive parameters for the dimensioning optimized positioning of protection measures. Within the impact analysis the relevant scenarios according to their intensities are determined using a 3D rockfall model. It calculates the trajectories of rock-bodies with realistic forms based on physical laws from defined starting points within the mapped rock area as well as steep slopes with loose blocks. The movement is of any chosen particle with defined axis length ratio, rounding and specific weight is strongly influenced by the spatial resolution of the digital terrain model (DTM as well as the interaction with the underground and possible tree-impacts). Therefore, the following parameters are implemented in the rockfall model: geometry of blocks, underground plasticity (value of attenuation), surface roughness as well as forest quality (amount of trees per area (ha) and trunk diameter). A very important aspect is the geometry of the blocks. For realistic results only a simulation with real block shapes will be sucessful. For these
A hybrid radial basis function-pseudospectral method for thermal convection in a 3-D spherical shell
NASA Astrophysics Data System (ADS)
Wright, G. B.; Flyer, N.; Yuen, D. A.
2010-07-01
A novel hybrid spectral method that combines radial basis function (RBF) and Chebyshev pseudospectral methods in a "2 + 1" approach is presented for numerically simulating thermal convection in a 3-D spherical shell. This is the first study to apply RBFs to a full 3-D physical model in spherical geometry. In addition to being spectrally accurate, RBFs are not defined in terms of any surface-based coordinate system such as spherical coordinates. As a result, when used in the lateral directions, as in this study, they completely circumvent the pole issue with the further advantage that nodes can be "scattered" over the surface of a sphere. In the radial direction, Chebyshev polynomials are used, which are also spectrally accurate and provide the necessary clustering near the boundaries to resolve boundary layers. Applications of this new hybrid methodology are given to the problem of convection in the Earth's mantle, which is modeled by a Boussinesq fluid at infinite Prandtl number. To see whether this numerical technique warrants further investigation, the study limits itself to an isoviscous mantle. Benchmark comparisons are presented with other currently used mantle convection codes for Rayleigh number (Ra) 7 × 103 and 105. Results from a Ra = 106 simulation are also given. The algorithmic simplicity of the code (mostly due to RBFs) allows it to be written in less than 400 lines of MATLAB and run on a single workstation. We find that our method is very competitive with those currently used in the literature.
Solar Hot Water Heating by Natural Convection.
ERIC Educational Resources Information Center
Noble, Richard D.
1983-01-01
Presents an undergraduate laboratory experiment in which a solar collector is used to heat water for domestic use. The working fluid is moved by natural convection so no pumps are required. Experimental apparatus is simple in design and operation so that data can be collected quickly and easily. (Author/JN)
A Simple Classroom Demonstration of Natural Convection
ERIC Educational Resources Information Center
Wheeler, Dean R.
2005-01-01
This article explains a simple way to demonstrate natural convection, such as from a lit candle, in the classroom using an overhead projector. The demonstration is based on the principle of schlieren imaging, commonly used to visualize variations in density for gas flows.
NASA Astrophysics Data System (ADS)
Stemmer, K.; Harder, H.; Hansen, U.
2004-12-01
The style of convection in planetary mantles is presumably dominated by the strong dependence of the viscosity of the mantle material on temperature and pressure. While several efforts have been undertaken in cartesian geometry to investigate convection in media with strong temperature dependent viscosity, spherical models are still in their infancy and still limited to modest parameters. Spectral approaches are usually employed for spherical convection models which do not allow to take into account lateral variations, like temperature dependent viscosity. We have developed a scheme, based on a finite volume discretization, to treat convection in a spherical shell with strong temperature dependent viscosity. Our approach has been particularly tailored to run efficiently on parallel computers. The spherical shell is topologically divided into six cubes. The equations are formulated in primitive variables, and are treated in the cartesian cubes. In order to ensure mass conservation a SIMPLER pressure correction procedure is applied and to handle strong viscosity variations up to Δ η =106 and high Rayleigh-numbers up to Ra=108 the pressure correction algorithm is combined with a pressure weighted interpolation method to satisfy the incompressibility condition and to avoid oscillations. We study thermal convection in a basal and mixed-mode heated shell with stress free and isothermal boundary conditions, as a function of the Rayleigh-number and viscosity contrast. Besides the temperature dependence we have further explored the effects of pressure on the viscosity. As a general result we observe the existence of three regimes (mobile, sluggish and stagnant lid), characterized by the type of surface motion. Laterally averaged depth-profiles of velocity, temperature and viscosity exhibit significant deviations from the isoviscous case. As compared to cartesian geometries, convection in a spherical shell possesses strong memory for the initial state. At strong
Stellar models with mixing length and T(τ) relations calibrated on 3D convection simulations
NASA Astrophysics Data System (ADS)
Salaris, Maurizio; Cassisi, Santi
2015-05-01
The calculation of the thermal stratification in the superadiabatic layers of stellar models with convective envelopes is a long-standing problem of stellar astrophysics, and has a major impact on predicted observational properties such as radius and effective temperature. The mixing length theory, almost universally used to model the superadiabatic convective layers, contains one free parameter to be calibrated (αml) whose value controls the resulting effective temperature. Here we present the first self-consistent stellar evolution models calculated by employing the atmospheric temperature stratification, Rosseland opacities, and calibrated variable αml (dependent on effective temperature and surface gravity) from a recently published large suite of three-dimensional radiation hydrodynamics simulations of stellar convective envelopes and atmospheres for solar stellar composition. From our calculations (with the same composition of the radiation hydrodynamics simulations), we find that the effective temperatures of models with the hydro-calibrated variable αml (that ranges between ~1.6 and ~2.0 in the parameter space covered by the simulations) present only minor differences, by at most ~30-50 K, compared to models calculated at constant solar αml (equal to 1.76, as obtained from the same simulations). The depth of the convective regions is essentially the same in both cases. We also analyzed the role played by the hydro-calibrated T(τ) relationships in determining the evolution of the model effective temperatures, when compared to alternative T(τ) relationships often used in stellar model computations. The choice of the T(τ) can have a larger impact than the use of a variable αml compared to a constant solar value. We found that the solar semi-empirical T(τ) by Vernazza et al. (1981, ApJS, 45, 635) provides stellar model effective temperatures that agree quite well with the results with the hydro-calibrated relationships.
Heterogeneous nanofluids: natural convection heat transfer enhancement
2011-01-01
Convective heat transfer using different nanofluid types is investigated. The domain is differentially heated and nanofluids are treated as heterogeneous mixtures with weak solutal diffusivity and possible Soret separation. Owing to the pronounced Soret effect of these materials in combination with a considerable solutal expansion, the resulting solutal buoyancy forces could be significant and interact with the initial thermal convection. A modified formulation taking into account the thermal conductivity, viscosity versus nanofluids type and concentration and the spatial heterogeneous concentration induced by the Soret effect is presented. The obtained results, by solving numerically the full governing equations, are found to be in good agreement with the developed solution based on the scale analysis approach. The resulting convective flows are found to be dependent on the local particle concentration φ and the corresponding solutal to thermal buoyancy ratio N. The induced nanofluid heterogeneity showed a significant heat transfer modification. The heat transfer in natural convection increases with nanoparticle concentration but remains less than the enhancement previously underlined in forced convection case. PMID:21711755
Heterogeneous nanofluids: natural convection heat transfer enhancement.
Oueslati, Fakhreddine Segni; Bennacer, Rachid
2011-01-01
Convective heat transfer using different nanofluid types is investigated. The domain is differentially heated and nanofluids are treated as heterogeneous mixtures with weak solutal diffusivity and possible Soret separation. Owing to the pronounced Soret effect of these materials in combination with a considerable solutal expansion, the resulting solutal buoyancy forces could be significant and interact with the initial thermal convection. A modified formulation taking into account the thermal conductivity, viscosity versus nanofluids type and concentration and the spatial heterogeneous concentration induced by the Soret effect is presented. The obtained results, by solving numerically the full governing equations, are found to be in good agreement with the developed solution based on the scale analysis approach. The resulting convective flows are found to be dependent on the local particle concentration φ and the corresponding solutal to thermal buoyancy ratio N. The induced nanofluid heterogeneity showed a significant heat transfer modification. The heat transfer in natural convection increases with nanoparticle concentration but remains less than the enhancement previously underlined in forced convection case. PMID:21711755
Heterogeneous nanofluids: natural convection heat transfer enhancement
NASA Astrophysics Data System (ADS)
Oueslati, Fakhreddine Segni; Bennacer, Rachid
2011-12-01
Convective heat transfer using different nanofluid types is investigated. The domain is differentially heated and nanofluids are treated as heterogeneous mixtures with weak solutal diffusivity and possible Soret separation. Owing to the pronounced Soret effect of these materials in combination with a considerable solutal expansion, the resulting solutal buoyancy forces could be significant and interact with the initial thermal convection. A modified formulation taking into account the thermal conductivity, viscosity versus nanofluids type and concentration and the spatial heterogeneous concentration induced by the Soret effect is presented. The obtained results, by solving numerically the full governing equations, are found to be in good agreement with the developed solution based on the scale analysis approach. The resulting convective flows are found to be dependent on the local particle concentration φ and the corresponding solutal to thermal buoyancy ratio N. The induced nanofluid heterogeneity showed a significant heat transfer modification. The heat transfer in natural convection increases with nanoparticle concentration but remains less than the enhancement previously underlined in forced convection case.
The Phenix ultimate natural convection test
Gauthe, P.; Pialla, D.; Tenchine, D.; Vasile, A.; Rochwerger, D.
2012-07-01
The French sodium cooled fast reactor Phenix was shut down in 2009 after 35 years of operation. Before decommissioning, a final set of tests were performed by the CEA during 9 months. Several topics were involved such as thermal hydraulics, core physics and fuel behaviour. Among these ultimate experiments, two thermal hydraulic tests were performed: an asymmetrical test consisting in a trip of one secondary pump and a natural convection test in the primary circuit. Recognizing the unique opportunity offered by these Phenix ultimate tests, IAEA decided in 2007 to launch a Coordinated Research Project (CRP) devoted to benchmarking analyses with system codes on the Phenix natural convection test. One objective of the natural convection test in Phenix reactor is the assessment of the CATHARE system code for safety studies on future and advanced sodium cooled fast reactors. The aim of this paper is to describe this test, which was performed on June 22-23, 2009, and the associated benchmark specifications for the CRP work. The paper reminds briefly the Phenix reactor with the main physical parameters and the instrumentation used during the natural convection test. After that, the test scenario is described: - initial state at a power of 120 MWth, - test beginning resulting from a manual dry out of the two steam generators, - manual scram, - manual trip on the three primary pumps without back-up by pony motors, - setting and development of natural convection in the primary circuit, in a first phase without significant heat sink in the secondary circuits and in a second phase with significant heat sink in the secondary circuits, by opening the casing of steam generators to create an efficient heat sink, by air natural circulation in the steam generators casing. The benchmark case ends after this second phase, which corresponds to the experimental test duration of nearly 7 hours. The paper presents also the benchmark specifications data supplied by the CEA to all
Bifurcations and unfoldings in natural convection
Decker, W.J.; Dorning, J.
1996-12-31
Extensive numerical studies of bifurcations and unfoldings have been carried out for two important problems in natural convection. These are (a) the Rayleigh-Benard convection (RBC) problem-a rectangular cavity, with insulated sidewalls, heated at constant uniform temperature along the bottom and cooled at constant uniform temperature along the top; and (b) the volumetric heating convection (VHC) problem - a rectangular cavity, with insulated sidewalls and bottom, heated by a constant uniform volumetric heat source and cooled at constant uniform temperature along the top. The information available in the literature on RBC was used to evaluate and justify the approximations made in the current research, which has shed additional light on nonlinear phenomena in RBC and led to new basic information on the bifurcations and unfoldings that occur in VHC for which there were essentially no previous results available. Both problems arise in many important technological and scientific contexts, including reactor safety analysis and meteorological phenomena. In particular, VHC is relevant to the development of an understanding of the natural convective motion driven by the radioactive decay heat in the molten core mixture (corium) in the core catcher following a hypothetical reactor core meltdown accident and of that which occurs in the atmosphere due to the deposition of radiant solar energy. The calculations were done using newly developed versions of the nodal integral method (NIM) for steady-state flows in conjunction with extended system methods for numerical bifurcation analysis for the saddle-node and pitchfork bifurcation computations.
NASA Astrophysics Data System (ADS)
Magic, Z.; Weiss, A.; Asplund, M.
2015-01-01
Aims: We investigate the relation between 1D atmosphere models that rely on the mixing length theory and models based on full 3D radiative hydrodynamic (RHD) calculations to describe convection in the envelopes of late-type stars. Methods: The adiabatic entropy value of the deep convection zone, sbot, and the entropy jump, Δs, determined from the 3D RHD models, were matched with the mixing length parameter, αMLT, from 1D hydrostatic atmosphere models with identical microphysics (opacities and equation-of-state). We also derived the mass mixing length parameter, αm, and the vertical correlation length of the vertical velocity, C[vz,vz], directly from the 3D hydrodynamical simulations of stellar subsurface convection. Results: The calibrated mixing length parameter for the Sun is α๏MLT (Sbot) = 1.98. . For different stellar parameters, αMLT varies systematically in the range of 1.7 - 2.4. In particular, αMLT decreases towards higher effective temperature, lower surface gravity and higher metallicity. We find equivalent results for α๏MLT (ΔS). In addition, we find a tight correlation between the mixing length parameter and the inverse entropy jump. We derive an analytical expression from the hydrodynamic mean-field equations that motivates the relation to the mass mixing length parameter, αm, and find that it qualitatively shows a similar variation with stellar parameter (between 1.6 and 2.4) with the solar value of α๏m = 1.83.. The vertical correlation length scaled with the pressure scale height yields 1.71 for the Sun, but only displays a small systematic variation with stellar parameters, the correlation length slightly increases with Teff. Conclusions: We derive mixing length parameters for various stellar parameters that can be used to replace a constant value. Within any convective envelope, αm and related quantities vary strongly. Our results will help to replace a constant αMLT. Appendices are available in electronic form at http
Self-propulsion via natural convection
NASA Astrophysics Data System (ADS)
Ardekani, Arezoo; Mercier, Matthieu; Allshouse, Michael; Peacock, Thomas
2014-11-01
Natural convection of a fluid due to a heated or cooled boundary has been studied within a myriad of different contexts due to the prevalence of the phenomenon in environmental systems such as glaciers, katabatic winds, or magmatic chambers; and in engineered problems like natural ventilation of buildings, or cooling of electronic components. It has, however, hitherto gone unrecognized that boundary-induced natural convection can propel immersed objects. We experimentally investigate the motion of a wedge-shaped object, immersed within a two-layer fluid system, due to a heated surface. The wedge resides at the interface between the two fluid layers of different density, and its concomitant motion provides the first demonstration of the phenomenon of propulsion via boundary-induced natural convection. Established theoretical and numerical models are used to rationalize the propulsion speed by virtue of balancing the propulsion force against the appropriate drag force. We successfully verified the influence of various fluid and heat parameters on the predicted speed. now at IMFT (Institut de Mécanique des Fluides de Toulouse).
MERIDIONAL CIRCULATION DYNAMICS FROM 3D MAGNETOHYDRODYNAMIC GLOBAL SIMULATIONS OF SOLAR CONVECTION
Passos, Dário; Charbonneau, Paul; Miesch, Mark
2015-02-10
The form of solar meridional circulation is a very important ingredient for mean field flux transport dynamo models. However, a shroud of mystery still surrounds this large-scale flow, given that its measurement using current helioseismic techniques is challenging. In this work, we use results from three-dimensional global simulations of solar convection to infer the dynamical behavior of the established meridional circulation. We make a direct comparison between the meridional circulation that arises in these simulations and the latest observations. Based on our results, we argue that there should be an equatorward flow at the base of the convection zone at mid-latitudes, below the current maximum depth helioseismic measures can probe (0.75 R{sub ⊙}). We also provide physical arguments to justify this behavior. The simulations indicate that the meridional circulation undergoes substantial changes in morphology as the magnetic cycle unfolds. We close by discussing the importance of these dynamical changes for current methods of observation which involve long averaging periods of helioseismic data. Also noteworthy is the fact that these topological changes indicate a rich interaction between magnetic fields and plasma flows, which challenges the ubiquitous kinematic approach used in the vast majority of mean field dynamo simulations.
Zemskova, Varvara; Garaud, Pascale; Deal, Morgan; Vauclair, Sylvie
2014-11-10
Iron-rich layers are known to form in the stellar subsurface through a combination of gravitational settling and radiative levitation. Their presence, nature, and detailed structure can affect the excitation process of various stellar pulsation modes and must therefore be modeled carefully in order to better interpret Kepler asteroseismic data. In this paper, we study the interplay between atomic diffusion and fingering convection in A-type stars, as well as its role in the establishment and evolution of iron accumulation layers. To do so, we use a combination of three-dimensional idealized numerical simulations of fingering convection (which neglect radiative transfer and complex opacity effects) and one-dimensional realistic stellar models. Using the three-dimensional simulations, we first validate the mixing prescription for fingering convection recently proposed by Brown et al. (within the scope of the aforementioned approximation) and identify what system parameters (total mass of iron, iron diffusivity, thermal diffusivity, etc.) play a role in the overall evolution of the layer. We then implement the Brown et al. prescription in the Toulouse-Geneva Evolution Code to study the evolution of the iron abundance profile beneath the stellar surface. We find, as first discussed by Théado et al., that when the concurrent settling of helium is ignored, this accumulation rapidly causes an inversion in the mean molecular weight profile, which then drives fingering convection. The latter mixes iron with the surrounding material very efficiently, and the resulting iron layer is very weak. However, taking helium settling into account partially stabilizes the iron profile against fingering convection, and a large iron overabundance can accumulate. The opacity also increases significantly as a result, and in some cases it ultimately triggers dynamical convection. The direct effects of radiative acceleration on the dynamics of fingering convection (especially in the
Natural convective heat transfer from square cylinder
NASA Astrophysics Data System (ADS)
Novomestský, Marcel; Smatanová, Helena; Kapjor, Andrej
2016-06-01
This article is concerned with natural convective heat transfer from square cylinder mounted on a plane adiabatic base, the cylinders having an exposed cylinder surface according to different horizontal angle. The cylinder receives heat from a radiating heater which results in a buoyant flow. There are many industrial applications, including refrigeration, ventilation and the cooling of electrical components, for which the present study may be applicable
Turbulent Convection: Is 2D a good proxy of 3D?
NASA Technical Reports Server (NTRS)
Canuto, V. M.
2000-01-01
Several authors have recently carried out 2D simulations of turbulent convection for both solar and massive stars. Fitting the 2D results with the MLT, they obtain that alpha(sub MLT) greater than 1 specifically, 1.4 less than alpha(sub MLT) less than 1.8. The authors further suggest that this methodology could be used to calibrate the MLT used in stellar evolutionary codes. We suggest the opposite viewpoint: the 2D results show that MLT is internally inconsistent because the resulting alpha(sub MLT) greater than 1 violates the MLT basic assumption that alpha(sub MLT) less than 1. When the 2D results are fitted with the CM model, alpha(sub CMT) less than 1, in accord with the basic tenet of the model. On the other hand, since both MLT and CM are local models, they should be replaced by the next generation of non-local, time dependent turbulence models which we discuss in some detail.
Natural-convection promoter for geothermal wells
Allis, R.G.; James, R.
1980-09-01
Many geothermal wells stand with relatively cold water overlying hot water. If a pipe is inserted into such a well, natural convection will occur and hot water will flow to the top of the well. The convection-promoting pipe enables domestic wells which would normally require the use of a downhole pump or airlift (with attendant environmental problems of fluid disposal) to be satisfactorily operated with a downhole heat exchanger. In potentially powerful steam-water wells which are difficult to discharge, a pipe positioned beneath the water level should raise wellhead pressure to the point where spontaneous discharge is possible. In both cases, the permeability and temperature of the feed zones are the limiting factors for the heat output of the well.
Studies of heat source driven natural convection
NASA Technical Reports Server (NTRS)
Kulacki, F. A.; Nagle, M. E.; Cassen, P.
1974-01-01
Natural convection energy transport in a horizontal layer of internally heated fluid with a zero heat flux lower boundary, and an isothermal upper boundary, has been studied. Quantitative information on the time-mean temperature distribution and the fluctuating component of temperature about the mean temperature in steady turbulent convection are obtained from a small thermocouple inserted into the layer through the upper bounding plate. Data are also presented on the development of temperature at several vertical positions when the layer is subject to both a sudden increase and to a sudden decrease in power input. For changes of power input from zero to a value corresponding to a Rayleigh number much greater than the critical linear stability theory value, a slight hysteresis in temperature profiles near the upper boundary is observed between the heat-up and cool-down modes.
NASA Astrophysics Data System (ADS)
Al-Momar, S. A.; Deierling, W.; Williams, J. K.; Hoffman, E. G.
2014-12-01
Convectively induced turbulence (CIT) is commonly listed as a cause or factor in weather-related commercial aviation accidents. In-cloud CIT is generated in part by shears between convective updrafts and downdrafts. Total lightning is also dependent on a robust updraft and the resulting storm electrification. The relationship between total lightning and turbulence could prove useful in operational aviation settings with the use of future measurements from the geostationary lightning mapper (GLM) onboard the GOES-R satellite. Providing nearly hemispheric coverage of total lightning, the GLM could help identify CIT in otherwise data-sparse locations. For a severe thunderstorm case on 7 June 2012 in northeast Colorado, in-cloud eddy dissipation rate estimates from the NCAR/NEXRAD Turbulence Detection Algorithm were compared with cloud electrification data from the Colorado Lightning Mapping Array and radar products from the Denver, Colorado WSR-88D. These comparisons showed that high concentrations of very high frequency (VHF) source densities emitted by lightning occurred near and downstream of the storm's convective core. Severe turbulence was also shown to occur near this area, extending near the melting level of the storm and spreading upward and outward. Additionally, increases/decreases in VHF sources and turbulence volumes occurred within a few minutes of each other; although, light turbulence was shown to increase near one storm's dissipation. This may be due to increased shear from the now downdraft dominate storm. The 3D wind field from this case, obtained by either a dual-Doppler or a Variational Doppler Radar Assimilation System (VDRAS) analysis, will also be examined to further study the relationships between total lightning and thunderstorm kinematics. If these results prove to be robust, lightning may serve as a strong indicator of the location of moderate or greater turbulence.
On natural solutal convection in magnetic fluids
NASA Astrophysics Data System (ADS)
Ivanov, A. S.; Pshenichnikov, A. F.
2015-09-01
An experiment was carried out to investigate natural solutal convection in a magnetic fluid caused by non-homogeneous initial distribution of colloidal particles in a vertical Hele-Shaw cell. For experiment, we used a dilute magnetic fluid of the "magnetite-kerosene-oleic acid" type. The initial distribution of particles was formed by magnetophoresis of the drop-like aggregates and their sedimentation on the surface of the diamagnetic disk located in the center of the cell. Application of the magnetic field on the system led to the onset of the Rayleigh-Taylor instability and formation of descending convective jets. The velocity of the flow at the front of descending jets was measured for different values of cell thickness (up to 0.18 mm) and strength of the magnetic field generating the drop-like aggregates (up to 21 kA/m). The solutal Rayleigh numbers varied in the range Ra = 50-105. It was shown that the intensity of the convective flow characterized by the Reynolds number Re, increases with the Rayleigh number according to the power law: Re = 1.16 × 10-5Ra0.86.
Laminar natural convection under nonuniform gravity.
NASA Technical Reports Server (NTRS)
Lienhard, J.; Eichhorn, R.; Dhir, V.
1972-01-01
Laminar natural convection is analyzed for cases in which gravity varies with the distance from the leading edge of an isothermal plate. The study includes situations in which gravity varies by virtue of the varying slope of a surface. A general integral solution method which includes certain known integral solutions as special cases is developed to account for arbitrary position-dependence of gravity. A series method of solution is also developed for the full equations. Although it is more cumbersome it provides verification of the integral method.
The power spectrum of solar convection flows from high-resolution observations and 3D simulations
NASA Astrophysics Data System (ADS)
Yelles Chaouche, L.; Moreno-Insertis, F.; Bonet, J. A.
2014-03-01
Context. Understanding solar surface magnetoconvection requires the study of the Fourier spectra of the velocity fields. Nowadays, observations are available that resolve very small spatial scales, well into the subgranular range, almost reaching the scales routinely resolved in numerical magnetoconvection simulations. Comparison of numerical and observational data at present can provide an assessment of the validity of the observational proxies. Aims: Our aims are: (1) to obtain Fourier spectra for the photospheric velocity fields using the spectropolarimetric observations with the highest spatial resolution so far (~120 km), thus reaching for the first time spatial scales well into the subgranular range; (2) to calculate corresponding Fourier spectra from realistic 3D numerical simulations of magnetoconvection and carry out a proper comparison with their observational counterparts considering the residual instrumental degradation in the observational data; and (3) to test the observational proxies on the basis of the numerical data alone, by comparing the actual velocity field in the simulations with synthetic observations obtained from the numerical boxes. Methods: (a) For the observations, data from the SUNRISE/IMaX spectropolarimeter are used. (b) For the simulations, we use four series of runs obtained with the STAGGER code for different average signed vertical magnetic field values (0, 50, 100, and 200 G). Spectral line profiles are synthesized from the numerical boxes for the same line observed by IMaX (Fe I 5250.2 Å) and degraded to match the performance of the IMaX instrument. Proxies for the velocity field are obtained via Dopplergrams (vertical component) and local correlation tracking (LCT, for the horizontal component). Fourier power spectra are calculated and a comparison between the synthetic and observational data sets carried out. (c) For the internal comparison of the numerical data, velocity values on constant optical depth surfaces are used
NASA Astrophysics Data System (ADS)
Wittkowski, M.; Chiavassa, A.; Freytag, B.; Scholz, M.; Höfner, S.; Karovicova, I.; Whitelock, P. A.
2016-03-01
Aims: We aim at comparing spectro-interferometric observations of Mira variable asymptotic giant branch (AGB) stars with the latest 1D dynamic model atmospheres based on self-excited pulsation models (CODEX models) and with 3D dynamic model atmospheres including pulsation and convection (CO5BOLD models) to better understand the processes that extend the molecular atmosphere to radii where dust can form. Methods: We obtained a total of 20 near-infrared K-band spectro-interferometric snapshot observations of the Mira variables o Cet, R Leo, R Aqr, X Hya, W Vel, and R Cnc with a spectral resolution of about 1500. We compared observed flux and visibility spectra with predictions by CODEX 1D dynamic model atmospheres and with azimuthally averaged intensities based on CO5BOLD 3D dynamic model atmospheres. Results: Our visibility data confirm the presence of spatially extended molecular atmospheres located above the continuum radii with large-scale inhomogeneities or clumps that contribute a few percent of the total flux. The detailed structure of the inhomogeneities or clumps show a variability on time scales of 3 months and above. Both modeling attempts provided satisfactory fits to our data. In particular, they are both consistent with the observed decrease in the visibility function at molecular bands of water vapor and CO, indicating a spatially extended molecular atmosphere. Observational variability phases are mostly consistent with those of the best-fit CODEX models, except for near-maximum phases, where data are better described by near-minimum models. Rosseland angular diameters derived from the model fits are broadly consistent between those based on the 1D and the 3D models and with earlier observations. We derived fundamental parameters including absolute radii, effective temperatures, and luminosities for our sources. Conclusions: Our results provide a first observational support for theoretical results that shocks induced by convection and pulsation in the
Conjugate natural convection between horizontal eccentric cylinders
NASA Astrophysics Data System (ADS)
Nasiri, Davood; Dehghan, Ali Akbar; Hadian, Mohammad Reza
2016-06-01
This study involved the numerical investigation of conjugate natural convection between two horizontal eccentric cylinders. Both cylinders were considered to be isothermal with only the inner cylinder having a finite wall thickness. The momentum and energy equations were discretized using finite volume method and solved by employing SIMPLER algorithm. Numerical results were presented for various solid-fluid conductivity ratios (KR) and various values of eccentricities in different thickness of inner cylinder wall and also for different angular positions of inner cylinder. From the results, it was observed that in an eccentric case, and for KR < 10, an increase in thickness of inner cylinder wall resulted in a decrease in the average equivalent conductivity coefficient (overline{{K_{eq} }} ); however, a KR > 10 value caused an increase in overline{{K_{eq} }} . It was also concluded that in any angular position of inner cylinder, the value of overline{{K_{eq} }} increased with increase in the eccentricity.
NASA Astrophysics Data System (ADS)
Magri, F.; Inbar, N.; Raggad, M.; Möller, S.; Siebert, C.; Möller, P.; Kuehn, M.
2014-12-01
Lake Kinneret (Lake Tiberias or Sea of Galilee) is the most important freshwater reservoir in the Northern Jordan Valley. Simulations that couple fluid flow, heat and mass transport are built to understand the mechanisms responsible for the salinization of this important resource. Here the effects of permeability distribution on 2D and 3D convective patterns are compared. 2D simulations indicate that thermal brine in Haon and some springs in the Yamourk Gorge (YG) are the result of mixed convection, i.e. the interaction between the regional flow from the bordering heights and thermally-driven flow (Magri et al., 2014). Calibration of the calculated temperature profiles suggests that the faults in Haon and the YG provides paths for ascending hot waters, whereas the fault in the Golan recirculates water between 1 and 2 km depths. At higher depths, faults induce 2D layered convection in the surrounding units. The 2D assumption for a faulted basin can oversimplify the system, and the conclusions might not be fully correct. The 3D results also point to mixed convection as the main mechanism for the thermal anomalies. However, in 3D the convective structures are more complex allowing for longer flow paths and residence times. In the fault planes, hydrothermal convection develops in a finger regime enhancing inflow and outflow of heat in the system. Hot springs can form locally at the surface along the fault trace. By contrast, the layered cells extending from the faults into the surrounding sediments are preserved and are similar to those simulated in 2D. The results are consistent with the theory from Zhao et al. (2003), which predicts that 2D and 3D patterns have the same probability to develop given the permeability and temperature ranges encountered in geothermal fields. The 3D approach has to be preferred to the 2D in order to capture all patterns of convective flow, particularly in the case of planar high permeability regions such as faults. Magri, F., et al., 2014
Natural and Enhanced Attenuation of Chlorinated Solvents Using RT3D
Johnson, Christian D.; Truex, Michael J.; Clement, T P.
2006-07-25
RT3D (Reactive Transport in 3-Dimensions) is a reactive transport code that can be applied to model solute fate and transport for many different purposes. This document specifically addresses application of RT3D for modeling related to evaluation and implementation of Monitored Natural Attenuation (MNA). Selection of MNA as a remedy requires an evaluation process to demonstrate that MNA will meet the remediation goals. The U.S. EPA, through the Office of Solid Waste and Emergency Response (OSWER) Directive 9200.4?17P, provides the regulatory context for the evaluation and implementation of MNA. In a complementary fashion, the context for using fate and transport modeling as part of MNA evaluation is described in the EPA?s technical protocol for chlorinated solvent MNA, the Scenarios Evaluation Tool for Chlorinated Solvent MNA, and in this document. The intent of this document is to describe (1) the context for applying RT3D for chlorinated solvent MNA and (2) the attenuation processes represented in RT3D, (3) dechlorination reactions that may occur, and (4) the general approach for using RT3D reaction modules (including a summary of the RT3D reaction modules that are available) to model fate and transport of chlorinated solvents as part of MNA or for combinations of MNA and selected types of active remediation.
Natural and forced convection during solidification
NASA Astrophysics Data System (ADS)
Neufeld, Jerome A.
The following work marries theoretical and experimental approaches to study the interaction of an external shear flow with a solidifying porous medium. The porous medium, a dendritic 'mushy layer', is created when a super-eutectic binary alloy is cooled leading to solid crystals bathed in an interstitial fluid which is compositionally enriched. This compositional enrichment leads to natural buoyant instabilities in the solidifying porous medium coupled with instabilities in the adjoining liquid layer. Theoretically, the effect of an external shear flow on the convective instabilities inherent to this mushy layer is investigated using a linear stability analysis. The external flow is coupled to advective perturbations in the liquid and to flow in the mush through a perturbed mush-liquid interface. A complete numerical solution of the stability of the system is performed and a critical porous medium Rayleigh number is found which is a function of both the external flow speed and the wavenumber of the interfacial perturbations. By neglecting the effects of buoyancy in the liquid and solving only for the pressure perturbations on the corrugated mush-liquid interface induced by the external flow, a reduced model is constructed and solved analytically. These theoretical results are compared with experimental observations obtained in a laboratory flume in which an ammonium-chloride solution is solidified from below at a constant rate. The experimental results reveal that at flow speeds above critical, convection is forced within the mush leading to a series of zero solid fraction tesselations aligned perpendicular to the applied shear flow. The results of the experiments compare favorably to the linear stability analysis.
Natural Convection in Enclosed Porous or Fluid Media
ERIC Educational Resources Information Center
Saatdjian, Esteban; Lesage, François; Mota, José Paulo B.
2014-01-01
In Saatdjian, E., Lesage, F., and Mota, J.P.B, "Transport Phenomena Projects: A Method to Learn and to Innovate, Natural Convection Between Porous, Horizontal Cylinders," "Chemical Engineering Education," 47(1), 59-64, (2013), the numerical solution of natural convection between two porous, concentric, impermeable cylinders was…
fVisiOn: glasses-free tabletop 3D display to provide virtual 3D media naturally alongside real media
NASA Astrophysics Data System (ADS)
Yoshida, Shunsuke
2012-06-01
A novel glasses-free tabletop 3D display, named fVisiOn, floats virtual 3D objects on an empty, flat, tabletop surface and enables multiple viewers to observe raised 3D images from any angle at 360° Our glasses-free 3D image reproduction method employs a combination of an optical device and an array of projectors and produces continuous horizontal parallax in the direction of a circular path located above the table. The optical device shapes a hollow cone and works as an anisotropic diffuser. The circularly arranged projectors cast numerous rays into the optical device. Each ray represents a particular ray that passes a corresponding point on a virtual object's surface and orients toward a viewing area around the table. At any viewpoint on the ring-shaped viewing area, both eyes collect fractional images from different projectors, and all the viewers around the table can perceive the scene as 3D from their perspectives because the images include binocular disparity. The entire principle is installed beneath the table, so the tabletop area remains clear. No ordinary tabletop activities are disturbed. Many people can naturally share the 3D images displayed together with real objects on the table. In our latest prototype, we employed a handmade optical device and an array of over 100 tiny projectors. This configuration reproduces static and animated 3D scenes for a 130° viewing area and allows 5-cm-tall virtual characters to play soccer and dance on the table.
Natural convection around the human head.
Clark, R P; Toy, N
1975-01-01
1. Factors determining the convective flow patterns around the human head in 'still' conditions are discussed in relation to body posture. 2. The flow patterns have been visualized using a schlieren optical system which reveals that the head has a thicker 'insulating' layer of convecting air in the erect posture than in the supine position. 3. Local convective and radiative heat transfer measurements from the head have been using surface calorimeters. These results are seen to be closely related to the thickness of the convective boundary layer flows. 4. The total convective and radiative heat loss from the head of a subject in the erect and supine position has been evaluated from the local measurements. For the head of the supine subject the heat loss was found to be 30% more than when the subject was standing. PMID:1142118
Combining volumetric edge display and multiview display for expression of natural 3D images
NASA Astrophysics Data System (ADS)
Yasui, Ryota; Matsuda, Isamu; Kakeya, Hideki
2006-02-01
In the present paper the authors present a novel stereoscopic display method combining volumetric edge display technology and multiview display technology to realize presentation of natural 3D images where the viewers do not suffer from contradiction between binocular convergence and focal accommodation of the eyes, which causes eyestrain and sickness. We adopt volumetric display method only for edge drawing, while we adopt stereoscopic approach for flat areas of the image. Since focal accommodation of our eyes is affected only by the edge part of the image, natural focal accommodation can be induced if the edges of the 3D image are drawn on the proper depth. The conventional stereo-matching technique can give us robust depth values of the pixels which constitute noticeable edges. Also occlusion and gloss of the objects can be roughly expressed with the proposed method since we use stereoscopic approach for the flat area. We can attain a system where many users can view natural 3D objects at the consistent position and posture at the same time in this system. A simple optometric experiment using a refractometer suggests that the proposed method can give us 3-D images without contradiction between binocular convergence and focal accommodation.
Modeling Images of Natural 3D Surfaces: Overview and Potential Applications
NASA Technical Reports Server (NTRS)
Jalobeanu, Andre; Kuehnel, Frank; Stutz, John
2004-01-01
Generative models of natural images have long been used in computer vision. However, since they only describe the of 2D scenes, they fail to capture all the properties of the underlying 3D world. Even though such models are sufficient for many vision tasks a 3D scene model is when it comes to inferring a 3D object or its characteristics. In this paper, we present such a generative model, incorporating both a multiscale surface prior model for surface geometry and reflectance, and an image formation process model based on realistic rendering, the computation of the posterior model parameter densities, and on the critical aspects of the rendering. We also how to efficiently invert the model within a Bayesian framework. We present a few potential applications, such as asteroid modeling and Planetary topography recovery, illustrated by promising results on real images.
Suppression of Natural Convection in a Thermoacoustic Pulse Tube Refrigerator
NASA Astrophysics Data System (ADS)
Han, Jun-Qing; Liu, Qiu-Sheng
2013-05-01
The effects of gravity on the efficiency of thermoacoustic engines are investigated theoretically and experimentally, especially for thermoacoustic pulse tube refrigerators. The significant effects of gravity are found to be due to the presence of natural convection in the thermoacoustic pulse tube when the hot side of the tube is lower than the cold side. This kind of natural convection influences and reduces the efficiency of the thermoacoustic working system. Thus, how to suppress this natural convection becomes important for increasing the efficiency of thermoacoustic engines. Unlike the method of inserting a silk screen in a pulse tube, the present study uses a numerical simulation method to research the natural convection in pulse tubes, and we try to change the shape of the pulse tube to suppress this convection.
On Unsteady Natural Convection Between Spherical Shells
NASA Astrophysics Data System (ADS)
Feldman, Yuri; Colonius, Tim
2011-11-01
Natural convection between two concentric spheres is investigated with three-dimensional numerical simulations. Buoyancy is achieved by preserving a temperature difference between the internal hotter and the external colder boundaries of the spherical shell. The numerical simulations were performed for the two basic configurations characterized by external to internal radius ratios of 1.2 and 1.5. Slightly supercritical laminar regimes characterized by the Rayleigh numbers of order Ra ~ O(104-105) were simulated by utilizing a Direct Numerical Simulation (DNS) approach while a Large Eddy Simulation (LES) was used for investigation of turbulent regimes for Ra ~ O (108-109) . We discuss the topological characteristics of the both laminar and turbulent flows. One of the possible scenarios of steady-unsteady transition is proposed as well. Implications of the results for the design of a double-walled Montgolfiere aerobot for the exploration of Titan's atmosphere are discussed. Research supported by Jet Propulsion Laboratory with Dr. Jeffrey Hall as monitor.
Using natural versus artificial stimuli to perform calibration for 3D gaze tracking
NASA Astrophysics Data System (ADS)
Maggia, Christophe; Guyader, Nathalie; Guérin-Dugué, Anne
2013-03-01
The presented study tests which type of stereoscopic image, natural or artificial, is more adapted to perform efficient and reliable calibration in order to track the gaze of observers in 3D space using classical 2D eye tracker. We measured the horizontal disparities, i.e. the difference between the x coordinates of the two eyes obtained using a 2D eye tracker. This disparity was recorded for each observer and for several target positions he had to fixate. Target positions were equally distributed in the 3D space, some on the screen (with a null disparity), some behind the screen (uncrossed disparity) and others in front of the screen (crossed disparity). We tested different regression models (linear and non linear) to explain either the true disparity or the depth with the measured disparity. Models were tested and compared on their prediction error for new targets at new positions. First of all, we found that we obtained more reliable disparities measures when using natural stereoscopic images rather than artificial. Second, we found that overall a non-linear model was more efficient. Finally, we discuss the fact that our results were observer dependent, with variability's between the observer's behavior when looking at 3D stimuli. Because of this variability, we proposed to compute observer specific model to accurately predict their gaze position when exploring 3D stimuli.
NASA Astrophysics Data System (ADS)
Crameri, Fabio; Tackley, Paul
2014-05-01
The work presented aims at a better understanding of plate tectonics, a crucial dynamical feature within the global framework of mantle convection. Special focus is given to the interaction of subduction-related mantle flow and surface topography. Thereby, the application of a numerical model with two key functional requirements is essential: an evolution over a long time period to naturally model mantle flow and a physically correct topography calculation. The global mantle convection model presented in Crameri et al. (2012a) satisfies both of these requirements. First, it is efficiently calculated by the finite-volume code Stag-YY (e.g., Tackley 2008) using a multi-grid method on a fully staggered grid. Second, it applies the sticky-air method (Matsumoto and Tomoda 1983; Schmeling et al, 2008) and thus approximates a free surface when the sticky-air parameters are chosen carefully (Crameri et al., 2012b). This leads to dynamically self-consistent mantle convection with realistic, single-sided subduction. New insights are thus gained into the interplay of obliquely sinking plates, toroidal mantle flow and the arcuate shape of slabs and trenches. Numerous two-dimensional experiments provide optimal parameter setups that are applied to three-dimensional models in Cartesian and fully spherical geometries. Features observed and characterised in the latter experiments give important insight into the strongly variable behaviour of subduction zones along their strike. This includes (i) the spontaneous development of arcuate trench geometry, (ii) regional subduction polarity reversals and slab tearing, and the newly discovered features (iii) 'slab tunnelling' and (iv) 'back-slab spiral flow'. Overall, this study demonstrates the strong interaction between surface topography and mantle currents and highlights the variability of subduction zones and their individual segments. REFERENCES Crameri, F., P. J. Tackley, I. Meilick, T. V. Gerya, and B. J. P. Kaus (2012a), A free
Transient natural convection in heated inclined tubes
NASA Astrophysics Data System (ADS)
McEligot, Donald M.; Denbow, David A.; Murphy, Hugh D.
1990-05-01
To simulate natural convection flow patterns in directionally drilled wellbores, experiments and analyses were conducted for a circular tube with length-to-diameter (L/D) ratio of 36 at angles of 0, 20, and 35 degrees from the vertical. The tube was heated at the bottom and cooled at the top, and the insulation was adjusted so that approximately one- to two-thirds of the power dissipated was transferred through the tube wall to the surroundings. An aqueous solution of polyvinyl alcohol was employed as the working fluid in order to obtain low Rayleigh numbers corresponding to conditions in geothermal wellbores. Results were primarily qualitative but were useful in providing insight into the phenomena occurring. Steady-state temperature distributions were measured for the three orientations and for several heating rates to demonstrate the effects of tube angle and Rayleigh number. Transient measurements of the temperature distribution were obtained during cooling from a higher temperature without a heat source to calibrate the heat losses. With the electrical heat source, temporal data were taken during heating to examine the approach to steady state. Quasi-steady flow conditions were approached rapidly, but the overall time constant of the apparatus was of the order of one-third of a day. Predictions with the three-dimensional TEMPEST code were first tested by comparison with simple conduction analyses. Comparison with actual data showed good agreement of the predicted temperature levels for the maximum inclination, 35 degrees, and slightly poorer agreement for the other limit, a vertical tube. Trends of temperature level and Nusselt number with heating rate or Rayleigh number were reasonable, but the predicted variation of the end Nusselt number versus inclination was in the opposite direction from the experiment.
Transient natural convection in heated inclined tubes
McEligot, D.M. . Oceanic Div.); Denbow, D.A. ); Murphy, H.D. )
1990-05-01
To simulate natural convection flow patterns in directionally drilled wellbores, experiments and analyses were conducted for a circular tube with length-to-diameter (L/D) ratio of 36 at angles of 0{degree}, 20{degree}, and 35{degree} from the vertical. The tube was heated at the bottom and cooled at the top, and the insulation was adjusted so that approximately one- to two-thirds of the power dissipated was transferred through the tube wall to the surroundings. An aqueous solution of polyvinyl alcohol was employed as the working fluid in order to obtain low Rayleigh numbers corresponding to conditions in geothermal wellbores. Results were primarily qualitative but were useful in providing insight into the phenomena occurring. Steady-state temperature distributions were measured for the three orientations and for several heating rates to demonstrate the effects of tube angle and Rayleigh number. transient measurements of the temperature distribution were obtained during cooling from a higher temperature without a heat source to calibrate the heat losses. With the electrical heat source, temporal data were taken during heating to examine the approach to steady state. Quasi-steady flow conditions were approached rapidly, but the overall time constant of the apparatus was of the order of one-third of a day. Predictions with the three-dimensional TEMPEST code were first tested by comparison with simple conduction analyses. Comparison with actual data showed good agreement of the predicted temperature levels for the maximum inclination, 35{degree}, and slightly poorer agreement for the other limit, a vertical tube. Trends of temperature level and Nusselt number with heating rate or Rayleigh number were reasonable, but the predicted variation of the end Nusselt number versus inclination was in the opposite direction from the experiment. 75 refs., 20 figs., 8 tabs.
NASA Technical Reports Server (NTRS)
Bune, Andris V.; Gillies, Donald C.; Lehoczky, Sandor L.
1998-01-01
Numerical simulation of the HgCdTe growth by the vertical Bridgman method was performed using FIDAP finite element code. Double-diffusive melt convection is analyzed, as the primary factor at controls inhomogeneity of the solidified material. Temperature and concentration fields in the model are also coupled via material properties, such as thermal and solutal expansion coefficients with the dependence on both temperature and concentration, and melting temperature evaluation from pseudobinary CdTe-HgTe phase diagram. Experimental measurements were used to obtain temperature boundary conditions. Parametric study of the melt convection dependence on the gravity conditions was undertaken. It was found, that the maximum convection velocity in the melt can be reduced under certain conditions. Optimal conditions to obtain a near flat solidified interface are discussed. The predicted interface shape is in agreement with one obtained experimentally by quenching. The results of 3-D calculations are compared with previous 2- D findings. A video film featuring 3-D melt convection will be presented.
Zhang, Yimeng; Li, Xiong; Samonds, Jason M; Lee, Tai Sing
2016-03-01
Bayesian theory has provided a compelling conceptualization for perceptual inference in the brain. Central to Bayesian inference is the notion of statistical priors. To understand the neural mechanisms of Bayesian inference, we need to understand the neural representation of statistical regularities in the natural environment. In this paper, we investigated empirically how statistical regularities in natural 3D scenes are represented in the functional connectivity of disparity-tuned neurons in the primary visual cortex of primates. We applied a Boltzmann machine model to learn from 3D natural scenes, and found that the units in the model exhibited cooperative and competitive interactions, forming a "disparity association field", analogous to the contour association field. The cooperative and competitive interactions in the disparity association field are consistent with constraints of computational models for stereo matching. In addition, we simulated neurophysiological experiments on the model, and found the results to be consistent with neurophysiological data in terms of the functional connectivity measurements between disparity-tuned neurons in the macaque primary visual cortex. These findings demonstrate that there is a relationship between the functional connectivity observed in the visual cortex and the statistics of natural scenes. They also suggest that the Boltzmann machine can be a viable model for conceptualizing computations in the visual cortex and, as such, can be used to predict neural circuits in the visual cortex from natural scene statistics. PMID:26712581
Numerical prediction of natural convection in square partitioned enclosures
Kelkar, K.M. ); Patankar, S.V. . Dept. of Mechanical Engineering)
1990-01-01
This paper provides a detailed study of flow and heat transfer phenomena in partitioned enclosures that is useful in understanding the more complex processes that occur in natural convection flows in passive solar heated buildings, solar collectors, and other applications. Two-dimensional natural convection flows in square enclosures with partitions are analyzed for laminar flow. Side walls are assumed to be isothermal, while the top and bottom walls are adiabatic.
Analysis of natural convection in a low gravity environment
NASA Technical Reports Server (NTRS)
Mattor, Ethan E.; Durgin, William W.; Bloznalis, Peter; Schoenberg, Richard
1992-01-01
Natural convection inside a spherical container was studied experimentally with two apparatuses at low buoyancy levels. The data generated by these experiments, plotted nondimensionally as the Nusselt versus Rayleigh numbers, give correlations for Rayleigh numbers between 1000 and 10 exp 8, a range previously untested. These results show that natural convection has significant effects at a Rayleigh number of 1000 and higher, although the behavior of the Nusselt number as the conduction limit is approached is still unknown for a spherical geometry.
NASA Astrophysics Data System (ADS)
Rozel, Antoine; Golabek, Gregor; Tackley, Paul
2014-05-01
Thermodynamically consistent models of single phase grain size evolution have been proposed in the past years [Austin and Evans (2007), Ricard and Bercovici (2009), Rozel et al. (2011), Rozel (2012)]. In a recently updated version [Bercovici and Ricard (2012), PEPI], the mechanics of two-phase grain aggregates has been formulated following the same physical approach. Several non-linear mechanisms such as dynamic recrystallization or Zener pinning are now available in a single non-equilibrium formulation of grain size distributions evolution. The self-consistent generation of localized plate boundaries is predicted in [Bercovici and Ricard (2012), EPSL] using this model, but it has not been tested in a dynamically consistent way. We propose the first set of three-dimensional numerical simulations of mantle convection incorporating this formalism using the finite volume code StagYY [Tackley (2008)]. First, we detail how the model is numerically implemented. Pressure and velocity fields are solved on a staggered grid using a SIMPLER-like method. Multigrid W-cycles and extra coarse-grid relaxations are employed to enhance the convergence of Stokes and continuity equations. The grain size is stored on a large number of tracers advected through the computational domain, which prevent numerical diffusion and allows a high resolution in the shear zones developing in the lithosphere. We also describe the physical formalism itself and propose the set of free parameters of the model. Normal growth, dynamic recrystallization and phase transitions all have a strong effect on the average grain size. We use a visco-plastic rheology in which the viscous strain rate is obtained by summation of dislocation, diffusion and grain boundary sliding creep. Second, we describe the 3D grain size distribution in the mantle and in the lithosphere. We characterize in which conditions plate margins can form, mainly investigating grain growth, recrystallization and rheology related parameters
NASA Astrophysics Data System (ADS)
Fargier, Yannick; Dore, Ludovic; Antoine, Raphael; Palma Lopes, Sérgio; Fauchard, Cyrille
2016-04-01
The extraction of subsurface materials is a key element for the economy of a nation. However, natural degradation of underground quarries is a major issue from an economic and public safety point of view. Consequently, the quarries stakeholders require relevant tools to define hazards associated to these structures. Safety assessment methods of underground quarries are recent and mainly based on rock physical properties. This kind of method leads to a certain homogeneity assumption of pillar internal properties that can cause an underestimation of the risk. Electrical Resistivity Imaging (ERI) is a widely used method that possesses two advantages to overcome this limitation. The first is to provide a qualitative understanding for the detection and monitoring of anomalies in the pillar body (e.g. faults). The second is to provide a quantitative description of the electrical resistivity distribution inside the pillar. This quantitative description can be interpreted with constitutive laws to help decision support (water content decreases the mechanical resistance of a chalk). However, conventional 2D and 3D Imaging techniques are usually applied to flat surface surveys or to surfaces with moderate topography. A 3D inversion of more complex media (case of the pillar) requires a full consideration of the geometry that was never taken into account before. The Photogrammetric technique presents a cost effective solution to obtain an accurate description of the external geometry of a complex media. However, this method has never been fully coupled with a geophysical method to enhance/improve the inversion process. Consequently we developed a complete procedure showing that photogrammetric and ERI tools can be efficiently combined to assess a complex 3D structure. This procedure includes in a first part a photogrammetric survey, a processing stage with an open source software and a post-processing stage finalizing a 3D surface model. The second part necessitates the
Nature versus nurture in shallow convection
NASA Astrophysics Data System (ADS)
Romps, D. M.; Kuang, Z.
2009-12-01
We use tracers in a large-eddy simulation of shallow convection to show that stochastic entrainment, not cloud-base properties, determine the fate of convecting parcels. The tracers are used to diagnose the correlations between a parcel's state above the cloud base and both the parcel's state at the cloud base and its entrainment history. We find that the correlation with the cloud-base state goes to zero a few hundred meters above the cloud base. On the other hand, correlations between a parcel's state and its net entrainment are large. Evidence is found that the entrainment events may be described as a stochastic Poisson process. We construct a parcel model with stochastic entrainment that is able to replicate flux profiles and, more importantly, the observed variability. Turning off cloud-base variability has little effect on the results, which suggests that stochastic mass-flux models may be initialized with a single set of properties. The success of the stochastic parcel model suggests that it holds promise as the framework for a convective parameterization.
A decoupled monolithic projection method for natural convection problems
NASA Astrophysics Data System (ADS)
Pan, Xiaomin; Kim, Kyoungyoun; Lee, Changhoon; Choi, Jung-Il
2016-06-01
We propose an efficient monolithic numerical procedure based on a projection method for solving natural convection problems. In the present monolithic method, the buoyancy, linear diffusion, and nonlinear convection terms are implicitly advanced by applying the Crank-Nicolson scheme in time. To avoid an otherwise inevitable iterative procedure in solving the monolithic discretized system, we use a linearization of the nonlinear convection terms and approximate block lower-upper (LU) decompositions along with approximate factorization. Numerical simulations demonstrate that the proposed method is more stable and computationally efficient than other semi-implicit methods, preserving temporal second-order accuracy.
NATURAL CONVECTION IN PASSIVE SOLAR BUILDINGS: EXPERIMENTS, ANALYSIS AND RESULTS
Gadgil, A.; Bauman, F.; Kammerud, R.
1981-04-01
Computer programs have been developed to numerically simulate natural convection in two- and three-dimensional room geometries. The programs have been validated using published data from the literature, results from a full-scale experiment performed at the Massachusetts Institute of Technology, and results from a small-scale experiment performed at LBL. One of the computer programs has been used to study the influence of natural convection on the thermal performance of a single zone in a direct-gain passive solar building. It is found that the convective heat transfer coefficients between the air and the enclosure surfaces can be substantially different from the values assumed in the standard building energy analysis methods, and can exhibit significant variations across a given surface. This study implies that the building heating loads calculated by standard building energy analysis methods may have substantial errors as a result of their use of common assumptions regarding the convection processes which occur in an enclosure.
Natural convection heat transfer within horizontal spent nuclear fuel assemblies
Canaan, R.E.
1995-12-01
Natural convection heat transfer is experimentally investigated in an enclosed horizontal rod bundle, which characterizes a spent nuclear fuel assembly during dry storage and/or transport conditions. The basic test section consists of a square array of sixty-four stainless steel tubular heaters enclosed within a water-cooled rectangular copper heat exchanger. The heaters are supplied with a uniform power generation per unit length while the surrounding enclosure is maintained at a uniform temperature. The test section resides within a vacuum/pressure chamber in order to subject the assembly to a range of pressure statepoints and various backfill gases. The objective of this experimental study is to obtain convection correlations which can be used in order to easily incorporate convective effects into analytical models of horizontal spent fuel systems, and also to investigate the physical nature of natural convection in enclosed horizontal rod bundles in general. The resulting data consist of: (1) measured temperatures within the assembly as a function of power, pressure, and backfill gas; (2) the relative radiative contribution for the range of observed temperatures; (3) correlations of convective Nusselt number and Rayleigh number for the rod bundle as a whole; and (4) correlations of convective Nusselt number as a function of Rayleigh number for individual rods within the array.
NASA Astrophysics Data System (ADS)
Glišović, Petar; Forte, Alessandro M.
2016-06-01
The 3-D distribution of buoyancy in the convecting mantle drives a suite of convection-related manifestations. Although seismic tomography is providing increasingly resolved images of the present-day mantle heterogeneity, the distribution of mantle density variations in the geological past is unknown, and, by implication, this is true for the convection-related observables. The one major exception is tectonic plate motions, since geologic data are available to estimate their history and they currently provide the only available constraints on the evolution of 3-D mantle buoyancy in the past. We developed a new back-and-forth iterative method for time-reversed convection modeling with a procedure for matching plate velocity data at different instants in the past. The crucial aspect of this reconstruction methodology is to ensure that at all times plates are driven by buoyancy forces in the mantle and not vice versa. Employing tomography-based retrodictions over the Cenozoic, we estimate the global amplitude of the following observables: dynamic surface topography, the core-mantle boundary ellipticity, the free-air gravity anomalies, and the global divergence rates of tectonic plates. One of the major benefits of the new data assimilation method is the stable recovery of much shorter wavelength changes in heterogeneity than was possible in our previous work. We now resolve what appears to be two-stage subduction of the Farallon plate under the western U.S. and a deeply rooted East African Plume that is active under the Ethiopian volcanic fields during the Early Eocene.
NASA Astrophysics Data System (ADS)
Geenen, T.; Heister, T.; Van Den Berg, A. P.; Jacobs, M.; Bangerth, W.
2011-12-01
We present high resolution 3D results of the complex mineral phase distribution in the transition zone obtained by numerical modelling of mantle convection. We extend the work by [Jacobs and van den Berg, 2011] to 3D and illustrate the efficiency of adaptive mesh refinement for capturing the complex spatial distribution and sharp phase transitions as predicted by their model. The underlying thermodynamical model is based on lattice dynamics which allows to predict thermophysical properties and seismic wave speeds for the applied magnesium-endmember olivine-pyroxene mineralogical model. The use of 3D geometry allows more realistic prediction of phase distribution and seismic wave speeds resulting from 3D flow processes involving the Earth's transition zone and more significant comparisons with interpretations from seismic tomography and seismic reflectivity studies aimed at the transition zone. Model results are generated with a recently developed geodynamics modeling application based on dealII (www.dealii.org). We extended this model to incorporate both a general thermodynamic model, represented by P,T space tabulated thermophysical properties, and a solution strategy that allows for compressible flow. When modeling compressible flow in the so called truncated anelastic approximation framework we have to adapt the solver strategy that has been proven by several authors to be highly efficient for incompressible flow to incorporate an extra term in the continuity equation. We present several possible solution strategies and discuss their implication in terms of robustness and computational efficiency.
Capturing natural-colour 3D models of insects for species discovery and diagnostics.
Nguyen, Chuong V; Lovell, David R; Adcock, Matt; La Salle, John
2014-01-01
Collections of biological specimens are fundamental to scientific understanding and characterization of natural diversity-past, present and future. This paper presents a system for liberating useful information from physical collections by bringing specimens into the digital domain so they can be more readily shared, analyzed, annotated and compared. It focuses on insects and is strongly motivated by the desire to accelerate and augment current practices in insect taxonomy which predominantly use text, 2D diagrams and images to describe and characterize species. While these traditional kinds of descriptions are informative and useful, they cannot cover insect specimens "from all angles" and precious specimens are still exchanged between researchers and collections for this reason. Furthermore, insects can be complex in structure and pose many challenges to computer vision systems. We present a new prototype for a practical, cost-effective system of off-the-shelf components to acquire natural-colour 3D models of insects from around 3 mm to 30 mm in length. ("Natural-colour" is used to contrast with "false-colour", i.e., colour generated from, or applied to, gray-scale data post-acquisition.) Colour images are captured from different angles and focal depths using a digital single lens reflex (DSLR) camera rig and two-axis turntable. These 2D images are processed into 3D reconstructions using software based on a visual hull algorithm. The resulting models are compact (around 10 megabytes), afford excellent optical resolution, and can be readily embedded into documents and web pages, as well as viewed on mobile devices. The system is portable, safe, relatively affordable, and complements the sort of volumetric data that can be acquired by computed tomography. This system provides a new way to augment the description and documentation of insect species holotypes, reducing the need to handle or ship specimens. It opens up new opportunities to collect data for research
Capturing Natural-Colour 3D Models of Insects for Species Discovery and Diagnostics
Nguyen, Chuong V.; Lovell, David R.; Adcock, Matt; La Salle, John
2014-01-01
Collections of biological specimens are fundamental to scientific understanding and characterization of natural diversity—past, present and future. This paper presents a system for liberating useful information from physical collections by bringing specimens into the digital domain so they can be more readily shared, analyzed, annotated and compared. It focuses on insects and is strongly motivated by the desire to accelerate and augment current practices in insect taxonomy which predominantly use text, 2D diagrams and images to describe and characterize species. While these traditional kinds of descriptions are informative and useful, they cannot cover insect specimens “from all angles” and precious specimens are still exchanged between researchers and collections for this reason. Furthermore, insects can be complex in structure and pose many challenges to computer vision systems. We present a new prototype for a practical, cost-effective system of off-the-shelf components to acquire natural-colour 3D models of insects from around 3 mm to 30 mm in length. (“Natural-colour” is used to contrast with “false-colour”, i.e., colour generated from, or applied to, gray-scale data post-acquisition.) Colour images are captured from different angles and focal depths using a digital single lens reflex (DSLR) camera rig and two-axis turntable. These 2D images are processed into 3D reconstructions using software based on a visual hull algorithm. The resulting models are compact (around 10 megabytes), afford excellent optical resolution, and can be readily embedded into documents and web pages, as well as viewed on mobile devices. The system is portable, safe, relatively affordable, and complements the sort of volumetric data that can be acquired by computed tomography. This system provides a new way to augment the description and documentation of insect species holotypes, reducing the need to handle or ship specimens. It opens up new opportunities to collect data
Thermally induced natural convection effects in Yucca Mountain drifts.
Webb, Stephen W; Francis, Nicholas D; Dunn, Sandra Dalvit; Itamura, Michael T; James, Darryl L
2003-01-01
Thermally induced natural convection from the heat produced by emplaced waste packages is an important heat and mass transfer mechanism within the Yucca Mountain Project (YMP) drifts. Various models for analyzing natural convection have been employed. The equivalent porous medium approach using Darcy's law has been used in many YMP applications. However, this approach has questionable fidelity, especially for turbulent flow conditions. Computational fluid dynamics (CFD), which is based on the fundamental Navier-Stokes equations, is currently being evaluated as a technique to calculate thermally induced natural convection in YMP. Data-model comparisons for turbulent flow conditions show good agreement of CFD predictions with existing experiments including YMP-specific data. PMID:12714318
Nature of stress accommodation in sheared granular material: Insights from 3D numerical modeling
NASA Astrophysics Data System (ADS)
Mair, Karen; Hazzard, James F.
2007-07-01
Active faults often contain distinct accumulations of granular wear material. During shear, this granular material accommodates stress and strain in a heterogeneous manner that may influence fault stability. We present new work to visualize the nature of contact force distributions during 3D granular shear. Our 3D discrete numerical models consist of granular layers subjected to normal loading and direct shear, where gouge particles are simulated by individual spheres interacting at points of contact according to simple laws. During shear, we observe the transient microscopic processes and resulting macroscopic mechanical behavior that emerge from interactions of thousands of particles. We track particle translations and contact forces to determine the nature of internal stress accommodation with accumulated slip for different initial configurations. We view model outputs using novel 3D visualization techniques. Our results highlight the prevalence of transient directed contact force networks that preferentially transmit enhanced stresses across our granular layers. We demonstrate that particle size distribution (psd) controls the nature of the force networks. Models having a narrow (i.e. relatively uniform) psd exhibit discrete pipe-like force clusters with a dominant and focussed orientation oblique to but in the plane of shear. Wider psd models (e.g. power law size distributions D = 2.6) also show a directed contact force network oblique to shear but enjoy a wider range of orientations and show more out-of-plane linkages perpendicular to shear. Macroscopic friction level, is insensitive to these distinct force network morphologies, however, force network evolution appears to be linked to fluctuations in macroscopic friction. Our results are consistent with predictions, based on recent laboratory observations, that force network morphologies are sensitive to grain characteristics such as particle size distribution of a sheared granular layer. Our numerical
An optimal sensing strategy for recognition and localization of 3-D natural quadric objects
NASA Technical Reports Server (NTRS)
Lee, Sukhan; Hahn, Hernsoo
1991-01-01
An optimal sensing strategy for an optical proximity sensor system engaged in the recognition and localization of 3-D natural quadric objects is presented. The optimal sensing strategy consists of the selection of an optimal beam orientation and the determination of an optimal probing plane that compose an optimal data collection operation known as an optimal probing. The decision of an optimal probing is based on the measure of discrimination power of a cluster of surfaces on a multiple interpretation image (MII), where the measure of discrimination power is defined in terms of a utility function computing the expected number of interpretations that can be pruned out by a probing. An object representation suitable for active sensing based on a surface description vector (SDV) distribution graph and hierarchical tables is presented. Experimental results are shown.
NASA Technical Reports Server (NTRS)
Zhang, Jun; Ge, Lixin; Kouatchou, Jules
2000-01-01
A new fourth order compact difference scheme for the three dimensional convection diffusion equation with variable coefficients is presented. The novelty of this new difference scheme is that it Only requires 15 grid points and that it can be decoupled with two colors. The entire computational grid can be updated in two parallel subsweeps with the Gauss-Seidel type iterative method. This is compared with the known 19 point fourth order compact differenCe scheme which requires four colors to decouple the computational grid. Numerical results, with multigrid methods implemented on a shared memory parallel computer, are presented to compare the 15 point and the 19 point fourth order compact schemes.
Numerical modeling of crystal growth on a centrifuge for unstable natural convection configurations
NASA Technical Reports Server (NTRS)
Ramachandran, N.; Downey, J. P.; Curreri, P. A.; Jones, J. C.
1993-01-01
The fluid mechanics associated with crystal growth processes on centrifuges is modeled using 2D and 3D models. Two-dimensional calculations show that flow bifurcations exist in such crystal growth configurations where the ampoule is oriented in the same direction as the resultant gravity vector and a temperature gradient is imposed on the melt. A scaling analysis is formulated to predict the flow transition point from the natural convection dominated regime to the Coriolis force dominated regime. Results of 3D calculations are presented for two thermal configurations of the crystal growth cell: top heated and bottom heated with respect to the centrifugal acceleration. In the top heated configuration, a substantial reduction in the convection intensity within the melt can be attained by centrifuge operations, and close to steady diffusion-limited thermal conditions can be achieved over a narrow range of the imposed microgravity level. In the bottom heated configuration the Coriolis force has a stabilizing effect on fluid motion by delaying the onset of unsteady convection.
Numerical Solution of Natural Convection in Eccentric Annuli
Pepper, D.W.
2001-09-18
The governing equations for transient natural convection in eccentric annular space are solved with two high-order accurate numerical algorithms. The equation set is transformed into bipolar coordinates and split into two one-dimensional equations: finite elements are used in the direction normal to the cylinder surfaces; the pseudospectral technique is used in the azimuthal direction. This report discusses those equations.
Natural convection of a magnetic fluid in a cubic enclosure
NASA Astrophysics Data System (ADS)
Kikura, H.; Sawada, T.; Tanahashi, T.
1993-04-01
Laminar natural convection heat transfer of a magnetic fluid in a cubic enclosure is examined experimentally. Wall-temperature distributions are visualized by thermosensitive liquid crystal sheets. The effect of the magnetic field on the transient temperature distributions, and the local and averaged Nusselt numbers are discussed.
Visual navigation of the UAVs on the basis of 3D natural landmarks
NASA Astrophysics Data System (ADS)
Karpenko, Simon; Konovalenko, Ivan; Miller, Alexander; Miller, Boris; Nikolaev, Dmitry
2015-12-01
This work considers the tracking of the UAV (unmanned aviation vehicle) on the basis of onboard observations of natural landmarks including azimuth and elevation angles. It is assumed that UAV's cameras are able to capture the angular position of reference points and to measure the angles of the sight line. Such measurements involve the real position of UAV in implicit form, and therefore some of nonlinear filters such as Extended Kalman filter (EKF) or others must be used in order to implement these measurements for UAV control. Recently it was shown that modified pseudomeasurement method may be used to control UAV on the basis of the observation of reference points assigned along the UAV path in advance. However, the use of such set of points needs the cumbersome recognition procedure with the huge volume of on-board memory. The natural landmarks serving as such reference points which may be determined on-line can significantly reduce the on-board memory and the computational difficulties. The principal difference of this work is the usage of the 3D reference points coordinates which permits to determine the position of the UAV more precisely and thereby to guide along the path with higher accuracy which is extremely important for successful performance of the autonomous missions. The article suggests the new RANSAC for ISOMETRY algorithm and the use of recently developed estimation and control algorithms for tracking of given reference path under external perturbation and noised angular measurements.
CamMedNP: Building the Cameroonian 3D structural natural products database for virtual screening
2013-01-01
Background Computer-aided drug design (CADD) often involves virtual screening (VS) of large compound datasets and the availability of such is vital for drug discovery protocols. We present CamMedNP - a new database beginning with more than 2,500 compounds of natural origin, along with some of their derivatives which were obtained through hemisynthesis. These are pure compounds which have been previously isolated and characterized using modern spectroscopic methods and published by several research teams spread across Cameroon. Description In the present study, 224 distinct medicinal plant species belonging to 55 plant families from the Cameroonian flora have been considered. About 80 % of these have been previously published and/or referenced in internationally recognized journals. For each compound, the optimized 3D structure, drug-like properties, plant source, collection site and currently known biological activities are given, as well as literature references. We have evaluated the “drug-likeness” of this database using Lipinski’s “Rule of Five”. A diversity analysis has been carried out in comparison with the ChemBridge diverse database. Conclusion CamMedNP could be highly useful for database screening and natural product lead generation programs. PMID:23590173
Two- and three-dimensional natural and mixed convection simulation using modular zonal models
Wurtz, E.; Nataf, J.M.; Winkelmann, F.
1996-07-01
We demonstrate the use of the zonal model approach, which is a simplified method for calculating natural and mixed convection in rooms. Zonal models use a coarse grid and use balance equations, state equations, hydrostatic pressure drop equations and power law equations of the form {ital m} = {ital C}{Delta}{sup {ital n}}. The advantage of the zonal approach and its modular implementation are discussed. The zonal model resolution of nonlinear equation systems is demonstrated for three cases: a 2-D room, a 3-D room and a pair of 3-D rooms separated by a partition with an opening. A sensitivity analysis with respect to physical parameters and grid coarseness is presented. Results are compared to computational fluid dynamics (CFD) calculations and experimental data.
On the convective-absolute nature of river bedform instabilities
NASA Astrophysics Data System (ADS)
Vesipa, Riccardo; Camporeale, Carlo; Ridolfi, Luca; Chomaz, Jean Marc
2014-12-01
River dunes and antidunes are induced by the morphological instability of stream-sediment boundary. Such bedforms raise a number of subtle theoretical questions and are crucial for many engineering and environmental problems. Despite their importance, the absolute/convective nature of the instability has never been addressed. The present work fills this gap as we demonstrate, by the cusp map method, that dune instability is convective for all values of the physical control parameters, while the antidune instability exhibits both behaviors. These theoretical predictions explain some previous experimental and numerical observations and are important to correctly plan flume experiments, numerical simulations, paleo-hydraulic reconstructions, and river works.
NASA Astrophysics Data System (ADS)
Wright, G.; Flyer, N.; Yuen, D. A.; Monnereau, M.; Zhang, S.; Wang, S. M.
2009-05-01
Many numerical methods, such as finite-differences, finite-volume, their yin-yang variants, finite-elements and spectral methods have been employed to study 3-D mantle convection. All have their own strengths, but also serious weaknesses. Spectrally accurate methods do not practically allow for node refinement and often involve cumbersome algebra while finite difference, volume, or element methods are generally low-order, adding excessive numerical diffusion to the model. For the 3-D mantle convection problem, we have introduced a new mesh-free approach, using radial basis functions (RBF). This method has the advantage of being algorithmic simple, spectrally accurate for arbitrary node layouts in multi-dimensions and naturally allows for node-refinement. One virtue of the RBF scheme allows the user to use a simple Cartesian geometry, while implementing the required boundary conditions for the temperature, velocities and stress components on a spherical surface at both the planetary surface and the core-mantle boundary. We have studied time- dependent mantle convection, using both a RBF-pseudospectral code and a code which uses spherical- harmonics in the angular direction and second-order finite volume in the radial direction. We have employed a third code , which uses spherical harmonics and higher-order finite-difference method a la Fornberg in the radial coordinate.We first focus on the onset of time-dependence at Rayleigh number Ra of 70,000. We follow the development of stronger time-dependence to a Ra of one million, using high enough resolution with 120 to 200 points in the radial direction and 128 to 256 spherical harmonics.
NASA Astrophysics Data System (ADS)
Ketkova, L. A.
2015-09-01
3D laser ultramicroscopy (3D LUM) is intended specially for determining the concentration and size distribution of submicron inclusions in the bulk samples of high-purity materials for visible and IR fiber optics. In this work the 3D LUM technique is shown to be able to identify the nature of individual inclusions detected. The measurement of the light scattered by an inclusion at a varied probe beam wavelength and polarization and at a varied scattered light collection angle makes it possible to determine the inclusion refractive index. The 3D LUM possibilities are illustrated by the example of studying the inclusion nature in the As2S3 glass samples prepared by the direct synthesis from elements in a quartz container at elevated temperatures.
NASA Astrophysics Data System (ADS)
Tackley, P. J.
2013-12-01
StagYY is a well-established code for modelling mantle convection in 3D spherical geometry (Tackley, PEPI 2008), incorporating several physical complexities such as compressibility, phase transitions, compositional variations, strongly temperature-dependent, non-linear rheology, tracers to track composition, continents, partial melting and melt migration. It uses a finite volume discretization (primitive variables on a staggered grid) on the yin-yang spherical grid (minimum overlap version). Geometric multigrid is used for simultaneous solution of the Stokes and mass conservation equations. Here, parallelization using MPI is discussed, and performance and scaling of the current StagYY version on up to 4096 cores on grids of up to 768x2304x512x2 cells (1.8 billion, corresponding to 7.2 billion unknowns) is demonstrated. Complexities related to scaling further to 100,000s to millions of cores are discussed together with possible solutions and performance projections.
Mimicking Natural Laminar to Turbulent Flow Transition: A Systematic CFD Study Using PAB3D
NASA Technical Reports Server (NTRS)
Pao, S. Paul; Abdol-Hamid, Khaled S.
2005-01-01
For applied aerodynamic computations using a general purpose Navier-Stokes code, the common practice of treating laminar to turbulent flow transition over a non-slip surface is somewhat arbitrary by either treating the entire flow as turbulent or forcing the flow to undergo transition at given trip locations in the computational domain. In this study, the possibility of using the PAB3D code, standard k-epsilon turbulence model, and the Girimaji explicit algebraic stresses model to mimic natural laminar to turbulent flow transition was explored. The sensitivity of flow transition with respect to two limiters in the standard k-epsilon turbulence model was examined using a flat plate and a 6:1 aspect ratio prolate spheroid for our computations. For the flat plate, a systematic dependence of transition Reynolds number on background turbulence intensity was found. For the prolate spheroid, the transition patterns in the three-dimensional boundary layer at different flow conditions were sensitive to the free stream turbulence viscosity limit, the reference Reynolds number and the angle of attack, but not to background turbulence intensity below a certain threshold value. The computed results showed encouraging agreements with the experimental measurements at the corresponding geometry and flow conditions.
Microscopic magnetic nature of K2NiF4-type 3d transition metal oxides
NASA Astrophysics Data System (ADS)
Sugiyama, J.; Nozaki, H.; Umegaki, I.; Higemoto, W.; Ansaldo, E. J.; Brewer, J. H.; Sakurai, H.; Kao, T.-H.; Yang, H.-D.; Månsson, M.
2014-12-01
In order to elucidate the magnetic nature of K2NiF4-type 3d transition metal oxides, we have measured μ+SR spectra for Sr2VO4, LaSrVO4, and Sr2CrO4 using powder samples. ZF- and wTF-μ+SR measurements propose that Sr2VO4 enters into the static antiferromagnetic (AF) order phase below 8 K. In addition, TF-μ+SR measurements evidence that the transition at 105 K is not magnetic but structural and/or electronic in origin. For LaSrVO4, static long-range order has not been observed down to 20 K, while, as T decreases from 145 K, wTF asymmetry starts to decrease below 60 K, suggesting the appearance and evolution of localized magnetic moments below 60 K. For Sr2CrO4, by contrast, both ZF- and wTF-μ+SR have confirmed the presence of antiferromagnetic order below 117 K, as predicted in the χ(T) curve.
Drift natural convection and seepage at the Yucca Mountain repository
NASA Astrophysics Data System (ADS)
Halecky, Nicholaus Eugene
The decay heat from radioactive waste that is to be disposed in the once proposed geologic repository at Yucca Mountain (YM) will significantly influence the moisture conditions in the fractured rock near emplacement tunnels (drifts). Additionally, large-scale convective cells will form in the open-air drifts and will serve as an important mechanism for the transport of vaporized pore water from the fractured rock, from the hot drift center to the cool drift end. Such convective processes would also impact drift seepage, as evaporation could reduce the build up of liquid water at the tunnel wall. Characterizing and understanding these liquid water and vapor transport processes is critical for evaluating the performance of the repository, in terms of water- induced canister corrosion and subsequent radionuclide containment. To study such processes, we previously developed and applied an enhanced version of TOUGH2 that solves for natural convection in the drift. We then used the results from this previous study as a time-dependent boundary condition in a high-resolution seepage model, allowing for a computationally efficient means for simulating these processes. The results from the seepage model show that cases with strong natural convection effects are expected to improve the performance of the repository, since smaller relative humidity values, with reduced local seepage, form a more desirable waste package environment.
NASA Astrophysics Data System (ADS)
Ma, Qisheng; Wu, Sheng; Tang, Yongchun
2008-11-01
Formation of the Carbon-13 ( 13C) and deuterium (D) doubly-substituted methane isotopologues ( 13CH 3D) in natural gases is studied utilizing both first-principle quantum mechanism molecular calculation and direct FTIR laboratorial measurements of specifically synthesized high isotope concentration methane gas. For 13CH 3D, the symmetrically breathing mode A0 emerges as IR-detectable attributed to the molecular symmetry lowering to C3v from Td of the non-isotopic methane (CH 4), along with a large vibrational frequency shift from ˜3000 to ˜2250 cm -1. Our studies also indicate that the concentration of 13CH 3D is dependent on the environmental temperature through isotope exchanges among methane isotopologues; and the Gibbs' Free Energy difference due to Quantum Mechanics Zero-Point vibrational motions has the major contribution to this temperature dependency. Potential geologic applications of the 13CH 3D measurement to natural gas exploration and assessments are also discussed. In order to detect the 13CH 3D concentration change of each 50 °C in the natural gas system, a 10 -9 resolution is desirable. Such a measurement could provide important add-on information to distinguish natural gas origin and distribution.
Temporal response of laser power standards with natural convective cooling.
Xu, Tao; Gan, Haiyong; Yu, Jing; Zang, Erjun
2016-01-25
Laser power detectors with natural convective cooling are convenient to use and hence widely applicable in a power range below 150 W. However, the temporal response characteristics of the laser power detectors need to be studied in detail for accurate measurement. The temporal response based on the absolute laser power standards with natural convective cooling is studied through theoretical analysis, numerical simulations, and experimental verifications. Our results show that the response deviates from a single exponential function and that an ultimate response balance is difficult to achieve because the temperature rise of the heat sink leads to continuous increase of the response. To determine the measurement values, an equal time reading method is proposed and validated by the laser power calibrations. PMID:26832477
Influence of geometry on natural convection in buildings
White, M.D.; Winn, C.B.; Jones, G.F.; Balcomb, J.D.
1985-01-01
Strong free convection airflows occur within passive solar buildings resulting from elevated temperatures of surfaces irradiated by solar energy compared with the cooler surfaces not receiving radiation. The geometry of a building has a large influence on the directions and magnitudes of natural airflows, and thus heat transfer between zones. This investigation has utilized a variety of reduced-scale building configurations to study the effects of geometry on natural convection heat transfer. Similarity between the reduced-scale model and a full-scale passive solar building is achieved by having similar geometries and by replacing air with Freon-12 gas as the model's working fluid. Filling the model with Freon-12 gas results in similarity in Prandtl numbers and Rayleigh numbers based on temperature differences in the range from 10/sup 9/ to 10/sup 11/. Results from four geometries are described with an emphasis placed on the effects of heat loss on zone temperature stratification shifts.
Verification of a numerical simulation technique for natural convection
Gadgil, A.; Bauman, F.; Altmayer, E.; Kammerud, R.C.
1983-03-01
The present paper describes a verification of CONVEC2 for single-zone geometries by comparison with the results of two natural convection experiments performed in small-scale rectangular enclosures. These experiments were selected because of the high Rayleigh numbers obtained and the small heat loss through the insulated surfaces. Comparisons are presented for (1) heat transfer rates, (2) fluid temperature profiles, and (3) surface heat flux distributions.
Prandtl Number Dependent Natural Convection with Internal Heat Sources
Kang Hee Lee; Seung Dong Lee; Kune Y. Suh; Joy L. Rempe; Fan-Bill Cheung; Sang B. Kim
2004-06-01
Natural convection plays an important role in determining the thermal load from debris accumulated in the reactor vessel lower head during a severe accident. Recently, attention is being paid to the feasibility of external vessel flooding as a severe accident management strategy and to the phenomena affecting the success path for retaining the molten core material inside the vessel. The heat transfer inside the molten core material can be characterized by the strong buoyancy-induced flows resulting from internal heating due to decay of fission products. The thermo-fluid dynamic characteristics of such flow depend strongly on the thermal boundary conditions. The spatial and temporal variation of heat flux on the pool wall boundaries and the pool superheat are mainly characterized by the natural convection flow inside the molten pool. In general, the natural convection heat transfer phenomena involving the internal heat generation are represented by the modified Rayleigh number (Ra’), which quantifies the internal heat source and hence the strength of the buoyancy force. In this study, tests were conducted in a rectangular section 250 mm high, 500 mm long and 160 mm wide. Twenty-four T-type thermocouples were installed in the test section to measure temperatures. Four T-type thermocouples were used to measure the boundary temperatures. The thermocouples were placed in designated locations after calibration. A direct heating method was adopted in this test to simulate the uniform heat generation. The experiments covered a range of Ra' between 1.5x106 and 7.42x1015 and the Prandtl number (Pr) between 0.7 and 6.5. Tests were conducted with water and air as simulant. The upper and lower boundary conditions were maintained uniform. The results demonstrated feasibility of the direct heating method to simulate uniform volumetric heat generation. Particular attentions were paid to the effect of Pr on natural convection heat transfer within the rectangular pool.
Uncovering the true nature of deformation microstructures using 3D analysis methods
NASA Astrophysics Data System (ADS)
Ferry, M.; Quadir, M. Z.; Afrin, N.; Xu, W.; Loeb, A.; Soe, B.; McMahon, C.; George, C.; Bassman, L.
2015-08-01
Three-dimensional electron backscatter diffraction (3D EBSD) has emerged as a powerful technique for generating 3D crystallographic information in reasonably large volumes of a microstructure. The technique uses a focused ion beam (FIB) as a high precision serial sectioning device for generating consecutive ion milled surfaces of a material, with each milled surface subsequently mapped by EBSD. The successive EBSD maps are combined using a suitable post-processing method to generate a crystallographic volume of the microstructure. The first part of this paper shows the usefulness of 3D EBSD for understanding the origin of various structural features associated with the plastic deformation of metals. The second part describes a new method for automatically identifying the various types of low and high angle boundaries found in deformed and annealed metals, particularly those associated with grains exhibiting subtle and gradual variations in orientation. We have adapted a 2D image segmentation technique, fast multiscale clustering, to 3D EBSD data using a novel variance function to accommodate quaternion data. This adaptation is capable of segmenting based on subtle and gradual variation as well as on sharp boundaries within the data. We demonstrate the excellent capabilities of this technique with application to 3D EBSD data sets generated from a range of cold rolled and annealed metals described in the paper.
Topological analysis of a mixing flow generated by natural convection
NASA Astrophysics Data System (ADS)
Contreras, Pablo Sebastián; de la Cruz, Luis Miguel; Ramos, Eduardo
2016-01-01
We use topological tools to describe the natural convective motion and the Lagrangian trajectories of a flow generated by stepwise, alternating heating and cooling protocol of opposite vertical walls of a cubic container. The working fluid considered is Newtonian and the system is in presence of the acceleration of gravity but the nonlinear terms are neglected, i.e., we study the piece-wise steady and linear problem. For this convective mixing flow, we identify invariant surfaces formed by the Lagrangian orbits of massless tracers that are topologically equivalent to spherical shells and period-1 lines with elliptic and hyperbolic segments that are located on symmetry planes. We describe the previous features as functions of the Rayleigh number in the range 3 × 104 ≤ Ra ≤ 5 × 105. We show that this system shares properties with other systems with non-toroidal invariant surfaces.
NASA Astrophysics Data System (ADS)
Tabary, P.; Bousquet, O.; Sénési, S.; Josse, P.
2009-09-01
Airports are recognized to become critical areas in the future given the expected doubling in air traffic by 2020. The increased density of aircrafts in the airport airspaces calls for improved systems and products to monitor in real-time potential hazards and thus meet the airport objectives in terms of safety and throughput. Among all meteorological hazards, convection is certainly the most impacting one. We describe here some innovative radar products that have recently been developed and tested at Météo France around the Paris airports. Those products rely on the French Doppler radar network consisting today of 24 elements with some of them being polarimetric. Reflectivity and Doppler volumetric data are concentrated from all 24 radar sites in real-time at the central level (Toulouse) where 3D Cartesian mosaics covering the entire French territory (i.e. a typical 1,000 by 1,000 km² area) are elaborated. The innovation with respect to what has been done previously is that the three components of the wind are retrieved by operational combination of the radial velocities. The final product, available in real-time every 15 minutes with a spatial resolution of 2.5 km horizontally and 500 m vertically, is a 3D grid giving the interpolated reflectivity and wind field (u, v and w) values. The 2.5 km resolution, arising from the fact that the retrieval is carried out every 15 minutes from radars typically spaced apart by 150 km, is not sufficient for airport airspace monitoring but is valuable for en-route monitoring. Its extension to the entire European space is foreseen. To address the specific needs in the airport areas, a downscaling technique has been proposed to merge the above-mentioned low-resolution 3D wind and reflectivity fields with the high resolution (5 minutes and 1 km²) 2D imagery of the Trappes radar that is the one that covers the Paris airports. The merging approach is based on the assumption that the Vertical Profile of Reflectivity (i.e. the
An experimental investigation of a natural convection solar air loop
Mastrullo, R.; Mazzei, P.; Vanoli, R.
1983-12-01
The interest that has been shown in the use of solar energy to heat dwellings following the ''passive'' design criteria does not correspond to the development of accurate theoretical and experimental analysis. This is particularly true for natural circulation solar air heaters. A significant application of these components is wall panel to complement south-facing windows in supplying solar heat directly to buildings. This idea, formerly suggested by Trombe et al., leads to various realizations, one of which was theoretically investigated by present authors. A convective loop panel consists of a glass layer and a black absorber that is backed by insulation. In the configuration shown the air flows in the channel in front of the absorber and the deflecting panel allows cool air to settle to the bottom of the U channel, preventing reverse thermocirculation during night or very low insolation periods. Since thermocirculation is the primary mode of heat transfer for the solar air heaters, the definition of an accurate convection model for the channel is essential for performance predictions. Studies on this subject - free convection between asymmetrically heated vertical planes - deal mainly with theoretical solutions for laminar flow, with the two usual boundary conditions. As the heat transfer process in the solar air loop cannot be expected to follow this model, there is the need of extensive experimental investigation.
Instabilities of Natural Convection in a Periodically Heated Layer
NASA Astrophysics Data System (ADS)
Hossain, M. Z.; Floryan, Jerzy M.
2013-11-01
Natural convection in a horizontal layer subject to a spatially periodic heating along the lower wall has been investigated. The heating produces sinusoidal temperature variations characterized by the wave number α and the Rayleigh number Rap. The primary response has the form of stationary rolls with axis orthogonal to the heating wave vector. For large α convection is limited to a thin layer adjacent to the lower wall with a uniform conduction above it. Linear stability was used to determine conditions leading to a secondary convection. Two mechanisms of instability have been identified. For α = 0(1), the parametric resonance dominates and leads to the pattern of instability that is locked-in with the pattern of the heating according to the relation δcr = α /2, where δcr denotes the component of the critical disturbance wave vector parallel to the heating wave vector. The second mechanism, Rayleigh-Bénard (RB) mechanism, dominates for large α. Competition between these mechanisms gives rise to non-commensurable states and appearance of soliton lattices, to the formation of distorted transverse rolls, and to the appearance of the wave vector component in the direction perpendicular to the forcing direction.
Santo, Vítor E; Popa, Elena G; Mano, João F; Gomes, Manuela E; Reis, Rui L
2015-06-01
The role of Platelet Lysates (PLs) as a source of growth factors (GFs) and as main element of three-dimensional (3D) hydrogels has been previously described. However, the resulting hydrogels usually suffer from high degree of contraction, limiting their usefulness. This work describes the development of a stable biomimetic 3D hydrogel structure based on PLs, through the spontaneous assembling of a high concentration of chitosan-chondroitin sulfate nanoparticles (CH/CS NPs) with PLs loaded by adsorption. The interactions between the NPs and the lysates resemble the ones observed in the extracellular matrix (ECM) native environment between glycosaminoglycans and ECM proteins. In vitro release studies were carried out focusing on the quantification of PDGF-BB and TGF-β1 GFs. Human adipose derived stem cells (hASCs) were entrapped in these 3D hydrogels and cultured in vitro under chondrogenic stimulus, in order to assess their potential use for cartilage regeneration. Histological, immunohistological and gene expression analysis demonstrated that the PL-assembled constructs entrapping hASCs exhibited results similar to the positive control (hASCS cultured in pellets), concerning the levels of collagen II expression and immunolocalization of collagen type I and II and aggrecan. Moreover, the deposition of new cartilage ECM was detected by alcian blue and safranin-O positive stainings. This work demonstrates the potential of PLs to act simultaneously as a source/carrier of GFs and as a 3D structure of support, through the application of a "bottom-up" approach involving the assembly of NPs, resulting in an enriched construct for cartilage regeneration applications. PMID:25795623
Thermally optimum spacing of vertical, natural convection cooled, parallel plates
NASA Astrophysics Data System (ADS)
Bar-Cohen, A.; Rohsenow, W. M.
Vertical two-dimensional channels formed by parallel plates or fins are a frequently encountered configuration in natural convection cooling in air of electronic equipment. In connection with the complexity of heat dissipation in vertical parallel plate arrays, little theoretical effort is devoted to thermal optimization of the relevant packaging configurations. The present investigation is concerned with the establishment of an analytical structure for analyses of such arrays, giving attention to useful relations for heat distribution patterns. The limiting relations for fully-developed laminar flow, in a symmetric isothermal or isoflux channel as well as in a channel with an insulated wall, are derived by use of a straightforward integral formulation.
NASA Astrophysics Data System (ADS)
Britz, Dieter
Convection has long been coupled with electrochemistry, and the name hydrodynamic voltammetry has become standard. In electroanalytical chemistry we mainly seek reproducible conditions. These are almost always attained by systems in which a steady convective state is achieved, although not always. Thus, the once popular dropping mercury electrode (see texts such as [74, 257]) has convection around it, but is never in steady state; it might be called a reproducible periodic dynamic state.
Particle filter based on thermophoretic deposition from natural convection flow
Sasse, A.G.B.M.; Nazaroff, W.W. ); Gadgil, A.J. )
1994-04-01
We present an analysis of particle migration in a natural convection flow between parallel plates and within the annulus of concentric tubes. The flow channel is vertically oriented with one surface maintained at a higher temperature than the other. Particle migration is dominated by advection in the vertical direction and thermophoresis in the horizontal direction. From scale analysis it is demonstrated that particles are completely removed from air flowing through the channel if its length exceeds L[sub c] = (b[sup 4]g/24K[nu][sup 2]), where b is the width of the channel, g is the acceleration of gravity, K is a thermophoretic coefficient of order 0.5, and [nu] is the kinematic viscosity of air. Precise predictions of particle removal efficiency as a function of system parameters are obtained by numerical solution of the governing equations. Based on the model results, it appears feasible to develop a practical filter for removing smoke particles from a smoldering cigarette in an ashtray by using natural convection in combination with thermophoresis. 22 refs., 8 figs., 1 tab.
Effect of enclosure shape on natural convection velocities
NASA Technical Reports Server (NTRS)
Robertson, S. J.; Nicholson, L. A.
1982-01-01
A numerical analysis was performed to compare natural convection velocities in two dimensional enclosures of various shape. The following shapes were investigated: circle, square, horizontal and upright 2 x 1 aspect ratio rectangles, horizontal and upright half circles, diamond. In all cases, the length scale in the various dimensionless parameters, such as Rayleigh number, is defined as the diameter of the equal area circle. Natural convection velocities were calculated for Rayleigh numbers of 1000 and 5000 with the temperature difference taken to be across (1) the maximum horizontal dimension, (2) the median horizontal line (line through centroid) and (3) the horizontal distance such that the temperature gradient is the same for shapes of equal area. For the class of shapes including the square, upright half circle and upright rectangle, the computed velocities were found to agree very closely with that of the equal area circle when the temperature difference is taken to be across the maximum horizontal dimension (condition (a)). The velocities for the horizontal rectangle and half circle were found to be approximately one half that of the equal area circle for the same condition. Better overall agreement among all shapes was obtained by setting the temperature difference across a distance such that the temperature gradients were equal for shapes of equal area.
NASA Astrophysics Data System (ADS)
Brodu, N.; Lague, D.
2012-04-01
3D point clouds derived from Terrestrial laser scanner (TLS) and photogrammetry are now frequently used in geomorphology to achieve greater precision and completeness in surveying natural environments than what was feasible a few years ago. Yet, scientific exploitation of these large and complex 3D data sets remains difficult and would benefit from automated classification procedures that could pre-process the raw point cloud data. Typical examples of applications are the separation of vegetation from ground or cliff outcrops, the distinction between fresh rock surfaces and rockfall, the classification of flat or rippled bed, and more generally the classification of 3D surfaces according to their morphology directly in the native point cloud data organization rather than after a sometime cumbersome meshing or gridding phase. Yet developing such classification procedures remains difficult because of the 3D nature of the data generated from ground based systems (as opposed to the 2.5D nature of aerial lidar data) and the heterogeneity and complexity of natural surfaces. We present a new software suite (CANUPO) that can classify raw point clouds in 3D based on a new geometrical measure: the multi-scale dimensionality. This method exploits the multi-resolution characteristics high-resolution datasets covering scales ranging from a few centimeters to hundred of meters. The dimensionality characterizes the local 3D organization of the point cloud within spheres centered on the measured points and varies from being 1D (points set along a line), 2D (points forming a plane) to the full 3D volume. By varying the diameter of the sphere, we track how the local cloud geometry behaves across scales (typically ranging from 5 cm to 1 m). We present the technique and illustrate its efficiency on two examples : separating riparian vegetation from ground, and classifying a steep mountain stream as vegetation, rock, gravel or water surface. In these two cases, separating the
Porous media flow problems: Natural convection and non-Newtonian
NASA Astrophysics Data System (ADS)
Walker, K. L.
1980-03-01
Natural convection of a Newtonian fluid and one dimensional flow of a nonNewtonian fluid are studied. Convection in a rectangular porous cavity driven by heating in the horizontal is analyzed by a number of different techniques which yield a fairly complete description of the two dimensional solutions. The solutions are governed by two dimensionless parameters: the Darcy-Rayleigh number R and cavity aspect ratio A. The flow behavior of a dilute solution of polyacrylamide in corn syrup flowing through porous media is also studied. Measurement of the pressure drop and flow rate are made for the solution flowing through a packed bed of glass beads. At low velocities the pressure drop as a function of velocity is the same as that for a Newtonian fluid of equal viscosity. At high flow rates the nonNewtonian fluid exhibited significantly higher pressure drops than a Newtonian fluid. Careful rheological measurements of the fluid are made using a Weissenberg rheogoniometer. From measurements of the dynamic viscosity shear it is determined that elastic effects are negligible. It is believed that the increased pressure gradients are caused by nonlinear viscous effects resulting from the extensional components of the flow.
Special session: computational predictability of natural convection flows in enclosures
Christon, M A; Gresho, P M; Sutton, S B
2000-08-14
Modern thermal design practices often rely on a ''predictive'' simulation capability--although predictability is rarely quantified and often difficult to confidently achieve in practice. The computational predictability of natural convection in enclosures is a significant issue for many industrial thermal design problems. One example of this is the design for mitigation of optical distortion due to buoyancy-driven flow in large-scale laser systems. In many instances the sensitivity of buoyancy-driven enclosure flows can be linked to the presence of multiple bifurcation points that yield laminar thermal convective processes that transition from steady to various modes of unsteady flow. This behavior is brought to light by a problem as ''simple'' as a differentially-heated tall rectangular cavity (8:1 height/width aspect ratio) filled with a Boussinesq fluid with Pr = 0.71--which defines, at least partially, the focus of this special session. For our purposes, the differentially-heated cavity provides a virtual fluid dynamics laboratory.
Morphogenesis and mechanostabilization of complex natural and 3D printed shapes
Tiwary, Chandra Sekhar; Kishore, Sharan; Sarkar, Suman; Mahapatra, Debiprosad Roy; Ajayan, Pulickel M.; Chattopadhyay, Kamanio
2015-01-01
The natural selection and the evolutionary optimization of complex shapes in nature are closely related to their functions. Mechanostabilization of shape of biological structure via morphogenesis has several beautiful examples. With the help of simple mechanics-based modeling and experiments, we show an important causality between natural shape selection as evolutionary outcome and the mechanostabilization of seashells. The effect of biological growth on the mechanostabilization process is identified with examples of two natural shapes of seashells, one having a diametrically converging localization of stresses and the other having a helicoidally concentric localization of stresses. We demonstrate how the evolved shape enables predictable protection of soft body parts of the species. The effect of bioavailability of natural material is found to be a secondary factor compared to shape selectivity, where material microstructure only acts as a constraint to evolutionary optimization. This is confirmed by comparing the mechanostabilization behavior of three-dimensionally printed synthetic polymer structural shapes with that of natural seashells consisting of ceramic and protein. This study also highlights interesting possibilities in achieving a new design of structures made of ordinary materials which have bio-inspired optimization objectives. PMID:26601170
Adjoint optimization of natural convection problems: differentially heated cavity
NASA Astrophysics Data System (ADS)
Saglietti, Clio; Schlatter, Philipp; Monokrousos, Antonios; Henningson, Dan S.
2016-06-01
Optimization of natural convection-driven flows may provide significant improvements to the performance of cooling devices, but a theoretical investigation of such flows has been rarely done. The present paper illustrates an efficient gradient-based optimization method for analyzing such systems. We consider numerically the natural convection-driven flow in a differentially heated cavity with three Prandtl numbers (Pr=0.15{-}7 ) at super-critical conditions. All results and implementations were done with the spectral element code Nek5000. The flow is analyzed using linear direct and adjoint computations about a nonlinear base flow, extracting in particular optimal initial conditions using power iteration and the solution of the full adjoint direct eigenproblem. The cost function for both temperature and velocity is based on the kinetic energy and the concept of entransy, which yields a quadratic functional. Results are presented as a function of Prandtl number, time horizons and weights between kinetic energy and entransy. In particular, it is shown that the maximum transient growth is achieved at time horizons on the order of 5 time units for all cases, whereas for larger time horizons the adjoint mode is recovered as optimal initial condition. For smaller time horizons, the influence of the weights leads either to a concentric temperature distribution or to an initial condition pattern that opposes the mean shear and grows according to the Orr mechanism. For specific cases, it could also been shown that the computation of optimal initial conditions leads to a degenerate problem, with a potential loss of symmetry. In these situations, it turns out that any initial condition lying in a specific span of the eigenfunctions will yield exactly the same transient amplification. As a consequence, the power iteration converges very slowly and fails to extract all possible optimal initial conditions. According to the authors' knowledge, this behavior is illustrated here
From 2D to 3D--a New Dimension for Modelling the Effect of Natural Products on Human Tissue.
Wrzesinski, Krzysztof; Fey, Stephen J
2015-01-01
Natural products, or their synthetic derivatives are a treasure trove to find potential candidates for novel drugs for human treatment. The selection of diamonds from the huge pile of worthless stone is a critical--and difficult--stage in the discovery pipeline. Of all the factors to be considered, perhaps the most important, is that the compound should have the desired effect on the tissue in vivo. Since it is not possible (or ethical) to test all compounds in vivo one must preselect using a surrogate assay system. While animal models have the advantage of being holistic and current 3D culture systems are reductionistic, they at least can be constructed from human cell types. In this review we will consider some of the evidence demonstrating that cells grown in 3D cultures have physiological performances that mimic functions seen in human tissues significantly better than cells grown using classical 2D culture systems. We will discuss advantages and disadvantages of these new culture technologies and highlight theoretical reasons for the differences. 3D cell culture technologies are more labour intensive than 2D culture systems and therefore their introduction is a trade-off between the value of obtaining data that is more relevant to the human condition against their through-put. It is already clear that future in vitro 3D systems will become more complex, using multiple cell types to more faithfully represent a particular tissue or even organ system. And one thing is sure - the diamonds are not easy to find! PMID:26429710
Natural Convection and Boiling for Cooling SRP Reactors During Loss of Circulation Conditions
Buckner, M.R.
2001-06-26
This study investigated natural convection and boiling as a means of cooling SRP reactors in the event of a loss of circulation accident. These studies show that single phase natural convection cooling of SRP reactors in shutdown conditions with the present piping geometry is probably not feasible.
Three-dimensional, transient natural convection in inclined wellbores
McEligot, D.M. . Oceanic Div.); Denbow, D.A. ); Murphy, H.D. )
1990-01-01
The occurrence of natural conduction in a wellbore can affect geothermal gradient measurements and heat flow estimates. In the Hot Dry Rock geothermal concept, the wellbores are purposely inclined in the deep regions to enhance heat production. To simulate natural convection flow patterns in directionally drilled wellbores, experiments and analyses were conducted for a circular tube with length to diameter (L/D) ratio of 36 at angles of 0{degrees}, 20{degrees}, and 35{degrees} from the vertical. The tube was heated at the bottom and cooled at the top, and the insulation was adjusted so that approximately one- to two-thirds of the power dissipated was transferred through the tube wall to the surroundings. An aqueous solution of polyvinyl alcohol was employed as the working fluid in order to obtain low Rayleigh numbers corresponding to conditions in geothermal wellbores. Temperature distributions were measured for the three orientations and for several heating rates to demonstrate the effects of tube angle and Rayleigh number. Comparison with measurements showed good agreement of the predicted temperature levels for the maximum inclination and slightly poorer agreement for the other limit, a vertical tube. 50 refs., 9 figs.
Torrance, K.E.; Catton, I.
1980-01-01
Natural convection in low aspect ratio rectangular enclosures is considered along with three-dimensional convection within rectangular boxes, natural convection flow visualization in irradiated water cooled by air flow over the surface, free convection in vertical slots, the stratification in natural convection in vertical enclosures, the flow structure with natural convection in inclined air-filled enclosures, and natural convection across tilted, rectangular enclosures of small aspect ratio. Attention is given to the effect of wall conduction and radiation on natural convection in a vertical slot with uniform heat generation of the heated wall, a numerical study of thermal insulation enclosure, free convection in a piston-cylinder enclosure with sinusoidal piston motion, natural convection heat transfer between bodies and their spherical enclosure, an experimental study of the steady natural convection in a horizontal annulus with irregular boundaries, three-dimensional natural convection in a porous medium between concentric inclined cylinders, a numerical solution for natural convection in concentric spherical annuli, and heat transfer by natural convection in porous media between two concentric spheres.
Zhang, S.Q.; Jioa, D.; Meng, Y.F.; Fan, Y.
1997-08-01
Three kinds of reservoir cores (limestone, sandstone, and shale with natural fractures) were used to study the effect of morphology of fracture surfaces on stress sensitivity. The cores, obtained from the reservoirs with depths of 2170 to 2300 m, have fractures which are mated on a large scale, but unmated on a fine scale. A specially designed photoelectric scanner with a computer was used to describe the topography of the fracture surfaces. Then, theoretical analysis of the fracture closure was carried out based on the fracture topography generated. The scanning results show that the asperity has almost normal distributions for all three types of samples. For the tested samples, the fracture closure predicted by the elastic-contact theory is different from the laboratory measurements because plastic deformation of the aspirates plays an important role under the testing range of normal stresses. In this work, the traditionally used elastic-contact theory has been modified to better predict the stress sensitivity of reservoir fractures. Analysis shows that the standard deviation of the probability density function of asperity distribution has a great effect on the fracture closure rate.
Evolving Views on the Scale and Nature of Mantle Convection
NASA Astrophysics Data System (ADS)
van der Hilst, R. D.
2014-12-01
Since seminal studies of transition zone discontinuities in the 1960ies and the advent of seismic tomography a decade later much progress has been made with the understanding of the scale and nature of mantle convection. First order questions remain, however, about the fluxes between the canonical upper and lower parts of Earth's mantle and the origin and nature of deep mantle heterogeneity. The first generation of tomographic models depicted fast shear wave propagation in the lowermost mantle beneath the circum-Pacific subduction zones and large low shear wspeed anomalies beneath Africa and the central Pacific. In P-wave models these structures are less apparent, and the anomalous Vp/Vs ratios and related variables are suggestive of chemical heterogeneity. Later tomographic studies revealed the pattern of subducted oceanic lithosphere in more detail and discovered that some slabs sink deep into the lower mantle whereas others remain, at least temporarily, in the transition zone. The complex flow trajectories and the evidence for compositional heterogeneity render simple end-member models of strict layering or unobstructed mantle flow untenable. Various seismic imaging methods have been used to map with increasing precision the variations in depth to the major mantle discontinuities, and also these results are not fully consistent with expectations for simple convection models. In addition, renewed scrutiny with more data and better methods suggest that the models of phase transitions around 410 and 660 km depth in the olivine component of a pyrolitic mantle composition are oversimplifications. Indeed, interfaces are also found at other depths, and many exceptions to the expected anti correlation of the interface topographies have been reported. Some of these observations can be explained with experimental and computational studies of the mineralogy and phase chemistry of deep mantle assemblages, but with such studies still restricted to fairly simple bulk
Clarke, S A; Choi, S Y; McKechnie, Melanie; Burke, G; Dunne, N; Walker, G; Cunningham, E; Buchanan, F
2016-02-01
Bone tissue engineering may provide an alternative to autograft, however scaffold optimisation is required to maximize bone ingrowth. In designing scaffolds, pore architecture is important and there is evidence that cells prefer a degree of non-uniformity. The aim of this study was to compare scaffolds derived from a natural porous marine sponge (Spongia agaricina) with unique architecture to those derived from a synthetic polyurethane foam. Hydroxyapatite scaffolds of 1 cm(3) were prepared via ceramic infiltration of a marine sponge and a polyurethane (PU) foam. Human foetal osteoblasts (hFOB) were seeded at 1 × 10(5) cells/scaffold for up to 14 days. Cytotoxicity, cell number, morphology and differentiation were investigated. PU-derived scaffolds had 84-91% porosity and 99.99% pore interconnectivity. In comparison marine sponge-derived scaffolds had 56-61% porosity and 99.9% pore interconnectivity. hFOB studies showed that a greater number of cells were found on marine sponge-derived scaffolds at than on the PU scaffold but there was no significant difference in cell differentiation. X-ray diffraction and inductively coupled plasma mass spectrometry showed that Si ions were released from the marine-derived scaffold. In summary, three dimensional porous constructs have been manufactured that support cell attachment, proliferation and differentiation but significantly more cells were seen on marine-derived scaffolds. This could be due both to the chemistry and pore architecture of the scaffolds with an additional biological stimulus from presence of Si ions. Further in vivo tests in orthotopic models are required but this marine-derived scaffold shows promise for applications in bone tissue engineering. PMID:26704539
Kim, Kwangtaek; Kim, Joongrock; Choi, Jaesung; Kim, Junghyun; Lee, Sangyoun
2015-01-01
Vision-based hand gesture interactions are natural and intuitive when interacting with computers, since we naturally exploit gestures to communicate with other people. However, it is agreed that users suffer from discomfort and fatigue when using gesture-controlled interfaces, due to the lack of physical feedback. To solve the problem, we propose a novel complete solution of a hand gesture control system employing immersive tactile feedback to the user's hand. For this goal, we first developed a fast and accurate hand-tracking algorithm with a Kinect sensor using the proposed MLBP (modified local binary pattern) that can efficiently analyze 3D shapes in depth images. The superiority of our tracking method was verified in terms of tracking accuracy and speed by comparing with existing methods, Natural Interaction Technology for End-user (NITE), 3D Hand Tracker and CamShift. As the second step, a new tactile feedback technology with a piezoelectric actuator has been developed and integrated into the developed hand tracking algorithm, including the DTW (dynamic time warping) gesture recognition algorithm for a complete solution of an immersive gesture control system. The quantitative and qualitative evaluations of the integrated system were conducted with human subjects, and the results demonstrate that our gesture control with tactile feedback is a promising technology compared to a vision-based gesture control system that has typically no feedback for the user's gesture inputs. Our study provides researchers and designers with informative guidelines to develop more natural gesture control systems or immersive user interfaces with haptic feedback. PMID:25580901
Kim, Kwangtaek; Kim, Joongrock; Choi, Jaesung; Kim, Junghyun; Lee, Sangyoun
2015-01-01
Vision-based hand gesture interactions are natural and intuitive when interacting with computers, since we naturally exploit gestures to communicate with other people. However, it is agreed that users suffer from discomfort and fatigue when using gesture-controlled interfaces, due to the lack of physical feedback. To solve the problem, we propose a novel complete solution of a hand gesture control system employing immersive tactile feedback to the user's hand. For this goal, we first developed a fast and accurate hand-tracking algorithm with a Kinect sensor using the proposed MLBP (modified local binary pattern) that can efficiently analyze 3D shapes in depth images. The superiority of our tracking method was verified in terms of tracking accuracy and speed by comparing with existing methods, Natural Interaction Technology for End-user (NITE), 3D Hand Tracker and CamShift. As the second step, a new tactile feedback technology with a piezoelectric actuator has been developed and integrated into the developed hand tracking algorithm, including the DTW (dynamic time warping) gesture recognition algorithm for a complete solution of an immersive gesture control system. The quantitative and qualitative evaluations of the integrated system were conducted with human subjects, and the results demonstrate that our gesture control with tactile feedback is a promising technology compared to a vision-based gesture control system that has typically no feedback for the user's gesture inputs. Our study provides researchers and designers with informative guidelines to develop more natural gesture control systems or immersive user interfaces with haptic feedback. PMID:25580901
Experimental analysis of natural convection within a thermosyphon
Clarksean, R.
1993-09-01
The heat transfer characteristics of a thermosyphon designed to passively cool cylindrical heat sources are experimentally studied. The analysis is based on recognizing the physics of the flow within different regions of the thermosyphon to develop empirical heat transfer correlations. The basic system consists of three concentric cylinders, with an outer channel between the outer two cylinders, and an inner channel between the inner two cylinders. Tests were conducted. with two different process material container diameters, representing the inner cylinder, and several different power levels. The experimentally determined local and average Nu numbers for the inner channel are in good agreement with previous work for natural convection between vertical parallel plates, one uniformly heated and the other thermally insulated. The implication is that the heat transfer off of each surface is independent of the adjacent surface for sufficiently high Ra numbers. The heat transfer is independent because of limited interaction between the boundary layers at sufficiently high Ra numbers. As a result of the limited interaction, the maximum temperature within the system remained constant, or decreased slightly when the radii of the inner cylinders increased for the same amount of heat removal.
Conjugate natural convection flow over a vertical surface with radiation
NASA Astrophysics Data System (ADS)
Siddiqa, Sadia; Hossain, Md. Anwar; Gorla, Rama Subba Reddy
2016-06-01
Numerical study of conjugate natural convection flow over a finite vertical surface with radiation is reported in this article. Rosseland diffusion approximation is used to express the radiative heat flux term. The governing boundary-layer equations are made dimensionless by means of a suitable form of non-similarity transformation. These equations are obtained in three regimes: (1) upstream (when ξ → 0), (2) downstream (when ξ → ∞ ) and (3) entire regime and are solved numerically. The solutions in the upstream and downstream regimes are obtained via shooting method whereas two-point implicit finite difference method is used to get the solutions for the entire regime. It is seen that asymptotic solutions give accurate results when compared with the numerical solution of the entire regime. The results indicate that the flow field and the temperature distributions are greatly influenced by thermal radiation parameter , R_d, surface temperature parameter, θ _w and Prandtl number Pr. It is established from the analysis that recirculation occurs in the flow specifically for R_d=1.5.
Analysis of Phenix natural convection test with the TRACE code
Chenu, A.; Mikityuk, K.; Chawla, R.
2012-07-01
Experimental data from the Natural Convection (NC) test performed in the Phenix reactor prior to its final shutdown have been used to further validate the single-phase sodium flow modeling in TRACE. The experimental data for the benchmark have been shared by the CEA in the frame of a Coordinated Research Project (CRP), initiated by the IAEA Technical Working Group on Fast Reactors (TWG-FR). This paper presents a complete TRACE model of the Phenix primary circuit developed for the analysis. Steady-state calculations at nominal (350 MWth) and reduced (120 MWth) power are compared to the experimental data for the validation of the model. We presents results from the 'blind' comparison, i.e. the comparison of the test results with those computed prior to the communication of the experimental data, so-called 'pre-test' results. 'Post-test' results, calculated from a model improved on the basis of the discrepancies identified from the blind comparison, are also presented. The analysis highlights the need to accurately simulate the reactor structures, since these define the thermal inertia of the system during the first phase of the transient. Furthermore, it shows the limitations of computed 1D-results when applied to the simulation of highly-stratified temperature fields. Nevertheless, the simulated reactor behavior and temperatures are found to match very well with the experimental data after the first two hours and, in general, the TRACE blind predictions may be considered as having been quite satisfactory. (authors)
Natural convection of ferrofluids in partially heated square enclosures
NASA Astrophysics Data System (ADS)
Selimefendigil, Fatih; Öztop, Hakan F.; Al-Salem, Khaled
2014-12-01
In this study, natural convection of ferrofluid in a partially heated square cavity is numerically investigated. The heater is located to the left vertical wall and the right vertical wall is kept at constant temperature lower than that of the heater. Other walls of the square enclosure are assumed to be adiabatic. Finite element method is utilized to solve the governing equations. The influence of the Rayleigh number (104≤Ra≤5×105), heater location (0.25H≤yh≤0.75H), strength of the magnetic dipole (0≤γ≤2), horizontal and vertical location of the magnetic dipole (-2H≤a≤-0.5H, 0.2H≤b≤0.8H) on the fluid flow and heat transfer characteristics are investigated. It is observed that different velocity components within the square cavity are sensitive to the magnetic dipole source strength and its position. The length and size of the recirculation zones adjacent to the heater can be controlled with magnetic dipole strength. Averaged heat transfer increases with decreasing values of horizontal position of the magnetic dipole source. Averaged heat transfer value increases from middle towards both ends of the vertical wall when the vertical location of the dipole source is varied. When the heater location is changed, a symmetrical behavior in the averaged heat transfer plot is observed and the minimum value of the averaged heat transfer is attained when the heater is located at the mid of vertical wall.
Transient natural convection of cold water in a vertical channel
NASA Astrophysics Data System (ADS)
Chiba, Ryoichi
2016-05-01
The two-dimensional differential transform method (DTM) is applied to analyse the transient natural convection of cold water in a vertical channel. The cold water gives rise to a density variation with temperature that may not be linearized. The vertical channel is composed of doubly infinite parallel plates, one of which has a constant prescribed temperature and the other of which is insulated. Considering the temperature-dependent viscosity and thermal conductivity of the water, approximate analytical (series) solutions for the temperature and flow velocity are derived. The transformed functions included in the solutions are obtained through a simple recursive procedure. Numerical computation is performed for the entire range of water temperature conditions around the temperature at the density extremum point, i.e. 4°C. Numerical results illustrate the effects of the temperature-dependent properties on the transient temperature and flow velocity profiles, volumetric flow rate, and skin friction. The DTM is a powerful tool for solving nonlinear transient problems as well as steady problems.
NASA Technical Reports Server (NTRS)
Bune, Andris V.; Gillies, Donald C.; Lehoczky, Sandor L.
1998-01-01
A numerical calculation for a non-dilute alloy solidification was performed using the FIDAP finite element code. For low growth velocities plane front solidification occurs. The location and the shape of the interface was determined using melting temperatures from the HgCdTe liquidus curve. The low thermal conductivity of the solid HgCdTe causes thermal short circuit through the ampoule walls, resulting in curved isotherms in the vicinity of the interface. Double-diffusive convection in the melt is caused by radial temperature gradients and by material density inversion with temperature. Cooling from below and the rejection at the solid-melt interface of the heavier HgTe-rich solute each tend to reduce convection. Because of these complicating factors dimensional rather then non-dimensional modeling was performed. Estimates of convection contributions for various gravity conditions was performed parametrically. For gravity levels higher then 1 0 -7 of earth's gravity it was found that the maximum convection velocity is extremely sensitive to gravity vector orientation and can be reduced at least by factor of 50% for precise orientation of the ampoule in the microgravity environment. The predicted interface shape is in agreement with one obtained experimentally by quenching. The results of 3-D modeling are compared with previous 2-D finding. A video film featuring melt convection will be presented.
Effects of finiteness on the thermo-fluid-dynamics of natural convection above horizontal plates
NASA Astrophysics Data System (ADS)
Guha, Abhijit; Sengupta, Sayantan
2016-06-01
A rigorous and systematic computational and theoretical study, the first of its kind, for the laminar natural convective flow above rectangular horizontal surfaces of various aspect ratios ϕ (from 1 to ∞) is presented. Two-dimensional computational fluid dynamic (CFD) simulations (for ϕ → ∞) and three-dimensional CFD simulations (for 1 ≤ ϕ < ∞) are performed to establish and elucidate the role of finiteness of the horizontal planform on the thermo-fluid-dynamics of natural convection. Great care is taken here to ensure grid independence and domain independence of the presented solutions. The results of the CFD simulations are compared with experimental data and similarity theory to understand how the existing simplified results fit, in the appropriate limiting cases, with the complex three-dimensional solutions revealed here. The present computational study establishes the region of a high-aspect-ratio planform over which the results of the similarity theory are approximately valid, the extent of this region depending on the Grashof number. There is, however, a region near the edge of the plate and another region near the centre of the plate (where a plume forms) in which the similarity theory results do not apply. The sizes of these non-compliance zones decrease as the Grashof number is increased. The present study also shows that the similarity velocity profile is not strictly obtained at any location over the plate because of the entrainment effect of the central plume. The 3-D CFD simulations of the present paper are coordinated to clearly reveal the separate and combined effects of three important aspects of finiteness: the presence of leading edges, the presence of planform centre, and the presence of physical corners in the planform. It is realised that the finiteness due to the presence of physical corners in the planform arises only for a finite value of ϕ in the case of 3-D CFD simulations (and not in 2-D CFD simulations or similarity theory
Role of natural convection in the dissolution of sessile droplets
NASA Astrophysics Data System (ADS)
Dietrich, Erik; Wildeman, Sander; Visser, Claas Willem; Hofhuis, Kevin; Kooij, E. Stefan; Zandvliet, Harold J. W.; Lohse, Detlef
2016-05-01
The dissolution process of small (initial (equivalent) radius $R_0 < 1$ mm) long-chain alcohol (of various types) sessile droplets in water is studied, disentangling diffusive and convective contributions. The latter can arise for high solubilities of the alcohol, as the density of the alcohol-water mixture is then considerably less as that of pure water, giving rise to buoyancy driven convection. The convective flow around the droplets is measured, using micro-particle image velocimetry ($\\mu$PIV) and the schlieren technique. When nondimensionalizing the system, we fnd a universal $Sh\\sim Ra^1/4$ scaling relation for all alcohols (of different solubilities) and all droplets in the convective regime. Here Sh is the Sherwood number (dimensionless mass flux) and Ra the Rayleigh number (dimensionless density difference between clean and alcohol-saturated water). This scaling implies the scaling relation $\\tau_c \\sim R^5/4$ of the convective dissolution time $\\tau_c$, which is found to agree with experimental data. We show that in the convective regime the plume Reynolds number (the dimensionless velocity) of the detaching alcohol-saturated plume follows $Re_p \\sim Sc^-1 Ra^5/8$, which is confirmed by the $\\mu$PIV data. Here, Sc is the Schmidt number. The convective regime exists when $Ra > Ra_t$, where $Ra_t = 12$ is the transition Ra-number as extracted from the data. For $Ra < Ra_t$ and smaller, convective transport is progressively overtaken by diffusion and the above scaling relations break down.
Joosik Yoo; Jun Young Choi; Moonuhn Kim . Dept. of Mechanical Engineering)
1994-01-01
Two-dimensional natural convection of a fluid of low Prandtl number (Pr = 0.02) in an annulus between two concentric horizontal cylinders is numerically investigated in a wide range of gap widths. For low Grashof numbers, a steady unicellular convection is obtained. Above a transition Grashof number that depends on the gap width, a steady bicellular flow occurs. With further increase of the Grashof number, steady or time-periodic multicellular convection occurs, and finally, complex unsteady convective flow appears. A plot is presented that predicts the type of flow patterns for various combination of gap widths and Grashof numbers.
Experimental and numerical study of turbulent natural convection in an open cubic cavity
NASA Astrophysics Data System (ADS)
Maytorena, V. M.; Piña-Ortiz, A.; Hinojosa, J. F.
2015-09-01
Study of natural convection in an open cubic cavity with side length of 1 m is presented. The experimental setup was built with the air as the heat transfer fluid. The vertical wall opposite to the aperture is subjected to uniform heat flux condition with the four different heat flux values in the range 55-333 W/m2, whereas the remaining walls were kept thermally insulated. The temperature at discrete locations inside the cavity was obtained which followed evaluation of heat transfer coefficient and Rayleigh number. The thermal and flow analysis in 3-D was based on the standard k-ɛ turbulence model and implemented using CFD software Fluent 6.3. The spatial distribution for temperature, velocity and turbulent viscosity are determined and analyzed in the perspective of experimental observations. The experimentally determined ranges of Rayleigh number, Nusselt number and heat transfer coefficient are 1.66 × 1011-7.1 × 1011, 185.94-243.31 and 4.88-6.83 W/m2 K, respectively. The observed maximum difference between the experimental and numerical values for heat transfer coefficient and Nusselt number are 10.8 and 14 % respectively.
Theoretical analysis of solar-driven natural convection energy conversion systems
Jacobs, E.W.; Lasier, D.D.
1984-01-01
This report presents a theoretical study of solar-powered natural convection tower (chimney) performance. Both heated and cooled towers are analyzed; the latter uses evaporating water as the cooling mechanism. The results, which are applicable to any open-cycle configuration, show that the ideal conversion efficiencies of both heated and cooled natural convection towers are linear functions of height. The performance of a heated tower in an adiabatic atmosphere ideally approaches the Carnot efficiency limit of approx. = 3.4%/km (1.0%/1000 ft). Including water pumping requirements, the ideal limit to cooled tower performance is approx. = 2.75%/km (0.85%/1000 ft). Ambient atmospheric conditions such as vertical temperature gradient (lapse rate) and relative humidity can have significantly adverse effects on natural convection tower performance. The combined effects of lapse rate and ambient relative humidity are especially important to cooled natural convection towers.
NASA Technical Reports Server (NTRS)
Ukanwa, A. O.; Stermole, F. J.; Golden, J. O.
1972-01-01
Natural convection effects in phase change thermal control devices were studied. A mathematical model was developed to evaluate natural convection effects in a phase change test cell undergoing solidification. Although natural convection effects are minimized in flight spacecraft, all phase change devices are ground tested. The mathematical approach to the problem was to first develop a transient two-dimensional conduction heat transfer model for the solidification of a normal paraffin of finite geometry. Next, a transient two-dimensional model was developed for the solidification of the same paraffin by a combined conduction-natural-convection heat transfer model. Throughout the study, n-hexadecane (n-C16H34) was used as the phase-change material in both the theoretical and the experimental work. The models were based on the transient two-dimensional finite difference solutions of the energy, continuity, and momentum equations.
Barthold, W.P.
1984-08-01
The scope of work is to summarize inherent safety advantages that are unique to the use of a carbide based fuel system and to summarize the technical issues regarding natural convection flow in LMFBR cores. As discussed in this report, carbide fuel provides the designer with far greater flexibility than oxide fuel. Carbide fuel systems can be designed to eliminate major accident initiators. They turn quantitative advantages into a qualitative advantage. The author proposed to LANL a series of core design and component concepts that would greatly enhance the safety of carbide over oxide systems. This report cites a series of safety advantages which potentially exist for a carbide fuel system. Natural convection issues have not been given much attention in the past. Only during the last few years has this issue been addressed in some detail. Despite claims to the contrary by some of the LMR contractors, the author does not think that the natural convection phenomena is fully understood. Some of the approximations made in natural convection transient analyses have probably a greater impact on calculated transient temperatures than the effects under investigation. Only integral in-pile experimental data and single assembly out-of-pile detailed data are available for comparisons with analytical models and correlations. Especially for derated cores, the natural convection capability of a LMR should be far superior to that of a LWR. The author ranks the natural convection capability of the LMR as the most important inherent safety feature.
Study of plasma natural convection induced by electron beam in atmosphere [
Deng, Yongfeng Han, Xianwei; Tan, Yonghua
2014-06-15
Using high-energy electron beams to ionize air is an effective way to produce a large-size plasma in the atmosphere. In particular, with a steady-state high power generator, some unique phenomena can be achieved, including natural convection of the plasma. The characteristics of this convection are studied both experimentally and numerically. The results show that an asymmetrical temperature field develops with magnitudes that vary from 295 K to 389 K at a pressure of 100 Torr. Natural convection is greatly enhanced under 760 Torr. Nevertheless, plasma transport is negligible in this convection flow field and only the plasma core tends to move upward. Parameter analysis is performed to discern influencing factors on this phenomenon. The beam current, reflecting the Rayleigh number Ra effect, correlates with convection intensity, which indicates that energy deposition is the underlying key factor in determining such convections. Finally, natural convection is concluded to be an intrinsic property of the electron beam when focused into dense air, and can be achieved by carefully adjusting equipment operations parameters.
Nature, theory and modelling of geophysical convective planetary boundary layers
NASA Astrophysics Data System (ADS)
Zilitinkevich, Sergej
2015-04-01
Geophysical convective planetary boundary layers (CPBLs) are still poorly reproduced in oceanographic, hydrological and meteorological models. Besides the mean flow and usual shear-generated turbulence, CPBLs involve two types of motion disregarded in conventional theories: 'anarchy turbulence' comprised of the buoyancy-driven plumes, merging to form larger plumes instead of breaking down, as postulated in conventional theory (Zilitinkevich, 1973), large-scale organised structures fed by the potential energy of unstable stratification through inverse energy transfer in convective turbulence (and performing non-local transports irrespective of mean gradients of transporting properties). C-PBLs are strongly mixed and go on growing as long as the boundary layer remains unstable. Penetration of the mixed layer into the weakly turbulent, stably stratified free flow causes turbulent transports through the CPBL outer boundary. The proposed theory, taking into account the above listed features of CPBL, is based on the following recent developments: prognostic CPBL-depth equation in combination with diagnostic algorithm for turbulence fluxes at the CPBL inner and outer boundaries (Zilitinkevich, 1991, 2012, 2013; Zilitinkevich et al., 2006, 2012), deterministic model of self-organised convective structures combined with statistical turbulence-closure model of turbulence in the CPBL core (Zilitinkevich, 2013). It is demonstrated that the overall vertical transports are performed mostly by turbulence in the surface layer and entrainment layer (at the CPBL inner and outer boundaries) and mostly by organised structures in the CPBL core (Hellsten and Zilitinkevich, 2013). Principal difference between structural and turbulent mixing plays an important role in a number of practical problems: transport and dispersion of admixtures, microphysics of fogs and clouds, etc. The surface-layer turbulence in atmospheric and marine CPBLs is strongly enhanced by the velocity shears in
Numerical simulation of natural convection in a sessile liquid droplet
NASA Astrophysics Data System (ADS)
Bartashevich, M. V.; Marchuk, I. V.; Kabov, O. A.
2012-06-01
Heat transfer in a sessile liquid droplet was studied with numerical methods. A computer code was developed for solving the problem of convection in an axisymmetric hemispherical droplet and in a spherical layer as well. The problem of establishing an equilibrium state in a droplet was solved using several variables: temperature, stream function, and vorticity. Simulation was performed for droplets of water, ethyl alcohol, and model liquids. Variable parameters: intensity of heat transfer from droplet surface, Rayleigh and Marangoni dimensionless criteria, and the characteristic temperature difference. It was revealed that the curve of convective flow intensity versus heat transfer intensity at droplet surface has a maximum. A dual-vortex structure was obtained in a stationary hemispherical profile of liquid droplet for the case of close values for thermocapillary and thermogravitational forces. Either thermocapillary or thermogravitational vortex might be dominating phenomena in the flow structure.
NASA Astrophysics Data System (ADS)
Zhang, Xu; Bao, Jian-Wen; Chen, Baode
2016-04-01
This presentation highlights a study in which a series of dry convective boundary layer (CBL) simulations are carried out using a generalized 3-dimensional (3-D) TKE-based parameterization scheme of sub-grid turbulent mixing in the Weather Research and Forecasting (WRF) model. The simulated characteristics of dry CBL are analyzed for the purpose of evaluating this scheme in comparison with a commonly-used scheme for sub-grid turbulent mixing in NWP models (i.e., the Mellor-Yamada 1.5-order TKE scheme). The same surface layer scheme is used in all the simulations so that only the sensitivity of the WRF model to different parameterizations of the sub-grid turbulent mixing above the surface layer is examined. The effect of horizontal grid resolution on the simulated CBL is also examined by running the model with grid sizes of 200, 400 m, 600 m, 1 km and 3 km. We will first compare the characteristics of the simulated CBL using the two schemes with the WRF LES dataset. We will then illustrate the importance of including the non-local component in the vertical buoyancy specification in the 3-D TKE-based scheme. Finally, comparing the results from the simulations against coarse-grained WRF LES dataset, we will show the feasibility and advantage of replacing conventional planetary boundary layer parameterization schemes with a scale-aware 3-D TKE-based scheme in the WRF model.
Decellularized kidney in the presence of chondroitin sulfate as a natural 3D scaffold for stem cells
Rafighdoust, Alireza; Shahri, Nasser Mahdavi; Baharara, Javad
2015-01-01
Objective(s): Use of biological scaffolds and automating the cells directing process with materials such as growth factors and glycosaminoglycans (GAGs) in a certain path may have beneficial effects in tissue engineering and regenerative medicine in future. In this research, chondroitin sulfate sodium was used for impregnation of the scaffolds. It is a critical component in extracellular matrix and plays an important role in signaling pathway; however, little is known about its role within mammalian development and cell linage specification. Materials and Methods: Due to its porous and appropriate structure and for putting cells in 3D space, the kidney of BALB/c mouse was selected and decellulalized using physical and chemical methods. After decellularization, the scaffold was impregnated in chondroitin sulfate solution (CS) for 24 hr. Then, 60×105 human adipose-derived mesenchymal stem cells were seeded on the scaffold to assess their behavior on day 5, 10, 15, 20, and 25. Results: After 48 hr, DAPI staining approved completed decellularized kidney by 1% SDS (sodium dodecyl sulfate). Migration and establishment of a number of cells to the remaining area of the glomerulus was observed. In addition, cell accumulation on the scaffold surface as well as cells migration to the depth of kidney formed an epithelium-like structure. Up to the day 15, microscopic study of different days of seeding showed the gradual adhesion of large number of cells to the scaffold. Conclusion: Glycosaminoglycan could be a right option for impregnation. It is used for smartification and strengthening of natural scaffolds and induction of some behaviors in stem cells. PMID:26557968
NASA Technical Reports Server (NTRS)
Adler, Robert F.; Tao, Wei-Kuo; Simpson, Joanne; Prasad, N.; Yeh, H.-Y. M.
1990-01-01
The relationship between the rain rate and the brightness temperature (Tb) was investigated using a cloud model/microwave radiative transfer model combination to obtain the rain-rate/Tb relations for four different frequencies: 10, 19, 37, and 86 GHz. The results at 19, 37, and 86 GHz were found to be significantly affected by ice in the modeled convective system, while the results at 10 GHz showed very little effect. Nonprecipitating cloud water was found to affect Tb in two ways. First, at low rain rates, the presence of significant cloud water produced higher Tb values than in cases with little cloud water. The second effects occurs at 19, 37, and 86 GHz at higher rainrates associated with significant ice formation; the scattering by ice lowered the Tb.
Natural-convection heat transfer of a spherical lighting fixture
Ikeda, Takamasa; Fujii, Tetsu
1994-09-01
The surface temperatures of the inner lamp and the outer globe of a spherical lighting fixture, the surfaces of which are painted black, were measured. From the results, the average convective heat-transfer coefficients between the inner lamp and the outer globe and on the outer surface of the globe were obtained. These data are correlated with the aid of existing equations for two concentric spheres and the outer surface of a single sphere. The relationships between the maximum and mean temperatures on the lamp and the globe were also obtained. By the use of these equations, a method for the optimal thermal design of spherical lighting fixtures is proposed.
NASA Technical Reports Server (NTRS)
Seybert, C. D.; Evans, J. W.; Leslie, F.; Jones, W. K., Jr.
2000-01-01
Natural convection, driven by temperature-or concentration gradients or both, is an inherent phenomenon during solidification of materials on Earth. This convection has practical consequences (e.g effecting macrosegregation) but also renders difficult the scientific examination of diffusive/conductive phenomena during solidification. It is possible to halt, or even reverse, natural convection by exploiting the variation (with temperature, for example) of the susceptibility of a material. If the material is placed in a vertical magnetic field gradient, a buoyancy force of magnetic origin arises and, at a critical field gradient, can balance the normal buoyancy forces to halt convection. At higher field gradients the convection can be reversed. The effect has been demonstrated in experiments at Marshall Space Flight Center where flow was measured by PIV in MnCl2 solution in a superconducting magnet. In auxiliary experiments the field in the magnet and the properties of the solution were measured. Computations of the natural convection, its halting and reversal, using the commercial software FLUENT were in good agreement with the measurements.
NASA Astrophysics Data System (ADS)
Ballmer, M. D.; Conrad, C. P.; Harmon, N.; Bianco, T. A.; Smith, E. I.
2011-12-01
Volcanism far from plate boundaries is typically ascribed to mantle plumes. Related hotspot melting can indeed explain a wide range of observations at a fair number of age-progressive volcano lineaments. Various other sites of intraplate volcanism, however, such as in the western US and the south Pacific, require alternative explanations, for which candidates are asthenospheric flow and lithospheric cracking. While cracking mechanisms presume abundant metastable melt in the asthenosphere, upward flow in the asthenosphere can generate significant amounts of decompression melting to support intraplate volcanism. Shear-driven upwellings (SDU) arise from rheological heterogeneity in the asthenosphere, a configuration that can redirect vigorous horizontal shearing into localized upwellings. Thermally-driven instabilities (small-scale convection, or SSC) develop at the base of cool and relatively thick lithosphere. It has been shown that areas with abundant non-hotspot volcanism indeed statistically correlate with regions of large asthenospheric shear [Conrad et al., 2011, see URL below]. In three-dimensional thermochemical numerical simulations, we investigate the relevance of SDU and SSC for mantle melting and intraplate volcanism. Depending on ambient mantle viscosity, plate thickness, asthenospheric shear, as well as the geometry and viscosity contrast of rheological heterogeneity (i.e., low-viscosity pockets), either SDU or SSC dominates asthenospheric flow and decompression melting. SDU is promoted over SSC by high viscosities, young plates, high shear, as well as wide, shallow, and strongly contrasting heterogeneity. Once melting initiates, it is amplified through the positive feedback effects of melt buoyancy (buoyant melting instability, or BMI), and lubrication. However, degrees of melting are effectively capped by stiffening of the residue upon melt and water extraction. In some cases, the dominance of SDU over SSC or vice-versa is transient, since
Thermocapillary flow and natural convection in a melt column with an unknown melt/solid interface
NASA Technical Reports Server (NTRS)
Lan, C. W.; Kou, Sindo
1991-01-01
A vertical melt column set up between an upper heating rod and a lower sample rod, i.e., the so-called half-zone system, is a convenient experimental tool for studying convection in the melt in floating-zone crystal growth. In order to help understand the convection observed in the melt column, a computer model has been developed to describe steady state, axisymmetrical thermocapillary flow and natural convection in the melt. The governing equations and boundary conditions are expressed in general non-orthogonal curvilinear coordinates in order to accurately treat the unknown melt/solid interface as well as all other physical boundaries in the system. The effects of key dimensionless variables on the following items are discussed: (1) convection and temperature distribution in the melt; (2) the shape of the melt/solid interface; (3) the height of the melt column. These dimensionless variables are the Grashof, Marangoni and Prandtl numbers.
Transient natural convection inside rigid drops in a liquid-liquid direct-contact heat exchanger
Hutchins, J.F.
1988-01-01
Natural convection was simulated inside spherical container and drops. The transient Navier-Stokes and energy equations were solved by employing finite-difference techniques. Pseudosteady-state natural convection inside spheres was simulated. Pseudosteady state was maintained by keeping the driving force for natural convection constant. To obtain pseudosteady state conditions, the temperature at the inside surface of the sphere was steadily increased so that the temperature difference between the surface and the center remained constant. The results were compared to experimental data found in the literature. It was found that the Nusselt number (Pr > 0.7) for pseudosteady state correlated to the Raleigh number by the following relation: Nu = 1.19Ra{sup .2215}, 10{sup 5} < Ra < 10{sup 8}. The simulation results were compared to experimental data of two other researchers who measured drop-temperature profiles in direct-contact heat-exchange columns. The simulation results demonstrate good correlation to the experimental data.
Transient testing of the FFTF for decay-heat removal by natural convection
Beaver, T R; Johnson, H G; Stover, R L
1982-06-01
This paper reports on the series of transient tests performed in the FFTF as a major part of the pre-operations testing program. The structure of the transient test program was designed to verify the capability of the FFTF to safely remove decay heat by natural convection. The series culminated in a scram from full power to complete natural convection in the plant, simulating a loss of all electrical power. Test results and acceptance criteria related to the verification of safe decay heat removal are presented.
Analysis and measurements of interzonal natural convection heat transfer in buildings
Hill, D.; Kirkpatrick, A.; Burns, P.
1986-08-01
Natural convection heat transfer through doorways can be an important process by which thermal energy is transferred from one zone to another zone of a building. The topic of this paper is interzonal natural convection in a two zone and a three zone multilevel full scale building. Aperture velocity and temperature distributions are measured and the experimental interzonal mass flow rate and heat transfer are determined. A Bernoulli model is derived to predict the neutral heights, velocity profiles, and interzonal heat transfer. The measured and predicted interzonal flow rate and heat transfer are compared and found to be in good agreement.
Natural Convection Cooling of the Advanced Stirling Radioisotope Generator Engineering Unit
NASA Technical Reports Server (NTRS)
Lewandowski, Edward J.; Hill, Dennis
2011-01-01
After fueling and prior to launch, the Advanced Stirling Radioisotope Generator (ASRG) will be stored for a period of time then moved to the launch pad for integration with the space probe and mounting on the launch vehicle. During this time, which could be as long as 3 years, the ASRG will operate continuously with heat rejected from the housing and fins. Typically, the generator will be cooled by forced convection using fans. During some of the ground operations, maintaining forced convection may add significant complexity, so allowing natural convection may simplify operations. A test was conducted on the ASRG Engineering Unit (EU) to quantify temperatures and operating parameters with natural convection only and determine if the EU could be safely operated in such an environment. The results show that with natural convection cooling the ASRG EU Stirling convertor pressure vessel temperatures and other parameters had significant margins while the EU was operated for several days in this configuration. Additionally, an update is provided on ASRG EU testing at NASA Glenn Research Center, where the ASRG EU has operated for over 16,000 hr and underwent extensive testing.
Validation of PARET for the modeling of heat transfer under natural convection core cooling
Ibrahim, J.K.; Kassim, M.S.; Mohammed, F.
1995-12-31
The PARET code is a one-dimensional, coupled thermal-hydraulic and point-kinetics code, which was originally developed for the analysis of SPERT-I transients and later adapted for the analysis of transient behavior in research reactors. Due to its ease of transportability and relative simplicity of input preparation, it is widely used internationally and is particularly attractive for research reactors with limited computational facilities. The thermal-hydraulic modeling of the current version of PARET accounts for buoyancy forces in the core and external pressure gradients that may arise from density differences between the core inlet and outlet. This feature of PARET makes it a useful tool for the analysis of research reactors cooled by natural convection as well as those cooled by forced convection. Since PARET has been applied to the analysis of the International Atomic Energy Agency 10-MW benchmark cores for protected and unprotected transients and also for the analysis of SPERT-I transients, its forced convection heat-removal model is reliable. However, there has been little experience with the capability of PARET to model heat removal in cores cooled by natural convection. This paper reports the results of some experiments performed at the Malaysian PUSPATI reactor to compare PARET predictions for power increases under natural convection core cooling to measured data.
Natural convection in a horizontal cylinder with axial rotation.
Sánchez, Odalys; Mercader, Isabel; Batiste, Oriol; Alonso, Arantxa
2016-06-01
We study the problem of thermal convection in a laterally heated horizontal cylinder rotating about its axis. A cylinder of aspect ratio Γ=H/2R=2 containing a small Prandtl number fluid (σ=0.01) representative of molten metals and molten semiconductors at high temperature is considered. We focus on a slow rotation regime (Ω<8), where the effects of rotation and buoyancy forces are comparable. The Navier-Stokes and energy equations with the Boussinesq approximation are solved numerically to calculate the basic states, analyze their linear stability, and compute several secondary flows originated from the instabilities. Due to the confined cylindrical geometry-the presence of lateral walls and lids-all the flows are completely three dimensional, even the basic steady states. Results characterizing the basic states as the rotation rate increases are presented. As it occurred in the nonrotating case for higher values of the Prandtl number, two curves of steady states with the same symmetric character coexist for moderate values of the Rayleigh number. In the range of Ω considered, rotation has a stabilizing effect only for very small values. As the value of the rotation rate approaches Ω=3.5 and Ω=4.5, the scenario of bifurcations becomes more complex due to the existence in both cases of very close bifurcations of codimension 2, which in the latter case involve both curves of symmetric solutions. PMID:27415364
Natural convection in a horizontal cylinder with axial rotation
NASA Astrophysics Data System (ADS)
Sánchez, Odalys; Mercader, Isabel; Batiste, Oriol; Alonso, Arantxa
2016-06-01
We study the problem of thermal convection in a laterally heated horizontal cylinder rotating about its axis. A cylinder of aspect ratio Γ =H /2 R =2 containing a small Prandtl number fluid (σ =0.01 ) representative of molten metals and molten semiconductors at high temperature is considered. We focus on a slow rotation regime (Ω <8 ), where the effects of rotation and buoyancy forces are comparable. The Navier-Stokes and energy equations with the Boussinesq approximation are solved numerically to calculate the basic states, analyze their linear stability, and compute several secondary flows originated from the instabilities. Due to the confined cylindrical geometry—the presence of lateral walls and lids—all the flows are completely three dimensional, even the basic steady states. Results characterizing the basic states as the rotation rate increases are presented. As it occurred in the nonrotating case for higher values of the Prandtl number, two curves of steady states with the same symmetric character coexist for moderate values of the Rayleigh number. In the range of Ω considered, rotation has a stabilizing effect only for very small values. As the value of the rotation rate approaches Ω =3.5 and Ω =4.5 , the scenario of bifurcations becomes more complex due to the existence in both cases of very close bifurcations of codimension 2, which in the latter case involve both curves of symmetric solutions.
Weight and water loss in the neonate in natural and forced convection.
Thompson, M H; Stothers, J K; McLellan, N J
1984-01-01
We describe a simple method of determining weight loss and hence water loss of infants in incubators. Unlike previously reported methods, it does not interfere with the microenvironment surrounding the infant. Weight loss of 16 term and 32 preterm infants was measured in both forced and natural convection. No significant increase in water loss was observed in the term infants but in the preterm infants the mean loss in natural convection was 0.85 g/kg/hour compared with 1.26 g/kg/hour in forced convection: in the most extreme situation it was doubled. This water loss represents a substantial energy loss and suggestions to minimise it are discussed. Images Fig. 1 PMID:6497432
Pebay, Cécile; Sella, Catherine; Thouin, Laurent; Amatore, Christian
2013-12-17
Mass transport at infinite regular arrays of microband electrodes was investigated theoretically and experimentally in unstirred solutions. Even in the absence of forced hydrodynamics, natural convection limits the convection-free domain up to which diffusion layers may expand. Hence, several regimes of mass transport may take place according to the electrode size, gap between electrodes, time scale of the experiment, and amplitude of natural convection. They were identified through simulation by establishing zone diagrams that allowed all relative contributions to mass transport to be delineated. Dynamic and steady-state regimes were compared to those achieved at single microband electrodes. These results were validated experimentally by monitoring the chronoamperometric responses of arrays with different ratios of electrode width to gap distance and by mapping steady-state concentration profiles above their surface through scanning electrochemical microscopy. PMID:24283775
Johnson, Christian D.; Truex, Michael J.
2006-07-25
This document describes a suite of MNA/EA reaction modules that were developed for addressing complex chlorinated solvent reactions using RT3D. As an introduction, an overview of these MNA/EA reaction modules is presented, including discussions of similarities between reaction modules, the purpose of key reaction parameters, and important considerations for using the reaction modules. Subsequent sections provide the details of the reaction kinetics (conceptual model and equations), data input requirements, and example (batch reactor) results for each reaction module. This document does not discuss reaction module implementation or validation; such information will accompany the software in the form of release notes or a supplement to the RT3D manual.
Tsai, Tsung-Yuan; Lu, Tung-Wu; Chen, Chung-Ming; Kuo, Mei-Ying; Hsu, Horng-Chaung
2010-03-15
Purpose: Accurate measurement of the three-dimensional (3D) rigid body and surface kinematics of the natural human knee is essential for many clinical applications. Existing techniques are limited either in their accuracy or lack more realistic experimental evaluation of the measurement errors. The purposes of the study were to develop a volumetric model-based 2D to 3D registration method, called the weighted edge-matching score (WEMS) method, for measuring natural knee kinematics with single-plane fluoroscopy to determine experimentally the measurement errors and to compare its performance with that of pattern intensity (PI) and gradient difference (GD) methods. Methods: The WEMS method gives higher priority to matching of longer edges of the digitally reconstructed radiograph and fluoroscopic images. The measurement errors of the methods were evaluated based on a human cadaveric knee at 11 flexion positions. Results: The accuracy of the WEMS method was determined experimentally to be less than 0.77 mm for the in-plane translations, 3.06 mm for out-of-plane translation, and 1.13 deg. for all rotations, which is better than that of the PI and GD methods. Conclusions: A new volumetric model-based 2D to 3D registration method has been developed for measuring 3D in vivo kinematics of natural knee joints with single-plane fluoroscopy. With the equipment used in the current study, the accuracy of the WEMS method is considered acceptable for the measurement of the 3D kinematics of the natural knee in clinical applications.
The Fractional Step Method Applied to Simulations of Natural Convective Flows
NASA Technical Reports Server (NTRS)
Westra, Douglas G.; Heinrich, Juan C.; Saxon, Jeff (Technical Monitor)
2002-01-01
This paper describes research done to apply the Fractional Step Method to finite-element simulations of natural convective flows in pure liquids, permeable media, and in a directionally solidified metal alloy casting. The Fractional Step Method has been applied commonly to high Reynold's number flow simulations, but is less common for low Reynold's number flows, such as natural convection in liquids and in permeable media. The Fractional Step Method offers increased speed and reduced memory requirements by allowing non-coupled solution of the pressure and the velocity components. The Fractional Step Method has particular benefits for predicting flows in a directionally solidified alloy, since other methods presently employed are not very efficient. Previously, the most suitable method for predicting flows in a directionally solidified binary alloy was the penalty method. The penalty method requires direct matrix solvers, due to the penalty term. The Fractional Step Method allows iterative solution of the finite element stiffness matrices, thereby allowing more efficient solution of the matrices. The Fractional Step Method also lends itself to parallel processing, since the velocity component stiffness matrices can be built and solved independently of each other. The finite-element simulations of a directionally solidified casting are used to predict macrosegregation in directionally solidified castings. In particular, the finite-element simulations predict the existence of 'channels' within the processing mushy zone and subsequently 'freckles' within the fully processed solid, which are known to result from macrosegregation, or what is often referred to as thermo-solutal convection. These freckles cause material property non-uniformities in directionally solidified castings; therefore many of these castings are scrapped. The phenomenon of natural convection in an alloy under-going directional solidification, or thermo-solutal convection, will be explained. The
Levin, A.E. ); Montgomery, B.H. )
1990-01-01
The Thermal-Hydraulic Out of Reactor Safety (THORS) Program at Oak Ridge National Laboratory (ORNL) had as its objective the testing of simulated, electrically heated liquid metal reactor (LMR) fuel assemblies in an engineering-scale, sodium loop. Between 1971 and 1985, the THORS Program operated 11 simulated fuel bundles in conditions covering a wide range of normal and off-normal conditions. The last test series in the Program, THORS-SHRS Assembly 1, employed two parallel, 19-pin, full-length, simulated fuel assemblies of a design consistent with the large LMR (Large Scale Prototype Breeder -- LSPB) under development at that time. These bundles were installed in the THORS Facility, allowing single- and parallel-bundle testing in thermal-hydraulic conditions up to and including sodium boiling and dryout. As the name SHRS (Shutdown Heat Removal System) implies, a major objective of the program was testing under conditions expected during low-power reactor operation, including low-flow forced convection, natural convection, and forced-to-natural convection transition at various powers. The THORS-SHRS Assembly 1 experimental program was divided up into four phases. Phase 1 included preliminary and shakedown tests, including the collection of baseline steady-state thermal-hydraulic data. Phase 2 comprised natural convection testing. Forced convection testing was conducted in Phase 3. The final phase of testing included forced-to-natural convection transition tests. Phases 1, 2, and 3 have been discussed in previous papers. The fourth phase is described in this paper. 3 refs., 2 figs.
NASA Astrophysics Data System (ADS)
Missoum, Abdelkrim; Elmir, Mohamed; Bouanini, Mohamed; Belkacem, Abdellah; Draoui, Belkacem
2016-03-01
This study focuses on the numerical simulation of heat transfer by natural convection in a rectangular enclosure, filled with a liquid metal (low Prandtl number) partially heated from below with a sinusoidal temperature. The value of the study lies in its involvement in the crystal growth for the manufacture of semiconductors and electronics cooling. Indeed, the occurrence of convection during crystal growth can lead to in homogeneities that lead to striations and defects that affect the quality of the crystals obtained by the Bridgman techniques or Chochrawlski. Temperature of the oscillations, due to the instabilities of the convective flow in the liquid metal, also induces non-uniform cooling in the solidification front. Convection is then studied in order to reduce it. A modelling of the problem in two dimensions was conducted using Comsol computer code that is based on the finite element method, by varying the configuration of the control parameters, namely, the Rayleigh number, the nature of fluid (Prandtl number) and amplitude of temperature on heat transfer rate (Nusselt number) on convective structures that appear.
Numerical investigation of natural convection heat transfer in a three-dimensional annular enclosure
NASA Astrophysics Data System (ADS)
Yung, Chain-Nan; de Witt, Kenneth J.; Keith, Theo G., Jr.
Natural convective flow and heat transfer in a three-dimensional annular enclose have been investigated numerically. The analysis uses dimensionless equations of continuity, momentum, and energy in Cartesian coordinates, which are cast into a generalized curvilinear system and solved by using a prediction-correction algorithm. For short horizontal cylinders, the local heat transfer rate is found to decrease sharply near the end walls due to convective velocity suppression; the overall heat transfer rate is less than that predicted by a two-dimensional model. Heat transfer rates are presented as a function of the Rayleigh number and compared with the available experimental data.
Tracking naturally occurring indoor features in 2-D and 3-D with lidar range/amplitude data
Adams, M.D.; Kerstens, A.
1998-09-01
Sensor-data processing for the interpretation of a mobile robot`s indoor environment, and the manipulation of this data for reliable localization, are still some of the most important issues in robotics. This article presents algorithms that determine the true position of a mobile robot, based on real 2-D and 3-D optical range and intensity data. The authors start with the physics of the particular type of sensor used, so that the extraction of reliable and repeatable information (namely, edge coordinates) can be determined, taking into account the noise associated with each range sample and the possibility of optical multiple-path effects. Again, applying the physical model of the sensor, the estimated positions of the mobile robot and the uncertainty in these positions are determined. They demonstrate real experiments using 2-D and 3-D scan data taken in indoor environments. To update the robot`s position reliably, the authors address the problem of matching the information recorded in a scan to, first, an a priori map, and second, to information recorded in previous scans, eliminating the need for an a priori map.
NASA Astrophysics Data System (ADS)
Danaila, Ionut; Moglan, Raluca; Hecht, Frédéric; Le Masson, Stéphane
2014-10-01
We present a new numerical system using finite elements with mesh adaptivity for the simulation of solid-liquid phase change systems. In the liquid phase, the natural convection flow is simulated by solving the incompressible Navier-Stokes equations with Boussinesq approximation. A variable viscosity model allows the velocity to progressively vanish in the solid phase, through an intermediate mushy region. The phase change is modeled by introducing an implicit enthalpy source term in the heat equation. The final system of equations describing the liquid-solid system by a single domain approach is solved using a Newton iterative algorithm. The space discretization is based on a P2-P1 Taylor-Hood finite elements and mesh adaptivity by metric control is used to accurately track the solid-liquid interface or the density inversion interface for water flows. The numerical method is validated against classical benchmarks that progressively add strong non-linearities in the system of equations: natural convection of air, natural convection of water, melting of a phase-change material and water freezing. Very good agreement with experimental data is obtained for each test case, proving the capability of the method to deal with both melting and solidification problems with convection. The presented numerical method is easy to implement using FreeFem++ software using a syntax close to the mathematical formulation.
Campbell, A N
2015-07-14
When any exothermic reaction proceeds in an unstirred vessel, natural convection may develop. This flow can significantly alter the heat transfer from the reacting fluid to the environment and hence alter the balance between heat generation and heat loss, which determines whether or not the system will explode. Previous studies of the effects of natural convection on thermal explosion have considered reactors where the temperature of the wall of the reactor is held constant. This implies that there is infinitely fast heat transfer between the wall of the vessel and the surrounding environment. In reality, there will be heat transfer resistances associated with conduction through the wall of the reactor and from the wall to the environment. The existence of these additional heat transfer resistances may alter the rate of heat transfer from the hot region of the reactor to the environment and hence the stability of the reaction. This work presents an initial numerical study of thermal explosion in a spherical reactor under the influence of natural convection and external heat transfer, which neglects the effects of consumption of reactant. Simulations were performed to examine the changing behaviour of the system as the intensity of convection and the importance of external heat transfer were varied. It was shown that the temporal development of the maximum temperature in the reactor was qualitatively similar as the Rayleigh and Biot numbers were varied. Importantly, the maximum temperature in a stable system was shown to vary with Biot number. This has important consequences for the definitions used for thermal explosion in systems with significant reactant consumption. Additionally, regions of parameter space where explosions occurred were identified. It was shown that reducing the Biot number increases the likelihood of explosion and reduces the stabilising effect of natural convection. Finally, the results of the simulations were shown to compare favourably with
Multidimensional hydrodynamic convection in full amplitude RR Lyrae models
NASA Astrophysics Data System (ADS)
Deupree, R.; Geroux, C.
2016-05-01
Multidimensional (both 2D and 3D) hydrodynamic calculations have been performed to compute full amplitude RR Lyrae models. The multi- dimensional nature allows convection to be treated in a more realistic way than simple 1D formulations such as the local mixing length theory. We focus on some aspects of multidimensional calculations and on the model for treating convection.
Luo, Kang; Yi, Hong-Liang Tan, He-Ping
2014-05-15
Transitions and bifurcations of transient natural convection in a horizontal annulus with radiatively participating medium are numerically investigated using the coupled lattice Boltzmann and direct collocation meshless (LB-DCM) method. As a hybrid approach based on a common multi-scale Boltzmann-type model, the LB-DCM scheme is easy to implement and has an excellent flexibility in dealing with the irregular geometries. Separate particle distribution functions in the LBM are used to calculate the density field, the velocity field and the thermal field. In the radiatively participating medium, the contribution of thermal radiation to natural convection must be taken into account, and it is considered as a radiative term in the energy equation that is solved by the meshless method with moving least-squares (MLS) approximation. The occurrence of various instabilities and bifurcative phenomena is analyzed for different Rayleigh number Ra and Prandtl number Pr with and without radiation. Then, bifurcation diagrams and dual solutions are presented for relevant radiative parameters, such as convection-radiation parameter Rc and optical thickness τ. Numerical results show that the presence of volumetric radiation changes the static temperature gradient of the fluid, and generally results in an increase in the flow critical value. Besides, the existence and development of dual solutions of transient convection in the presence of radiation are greatly affected by radiative parameters. Finally, the advantage of LB-DCM combination is discussed, and the potential benefits of applying the LB-DCM method to multi-field coupling problems are demonstrated.
Critical heat flux in natural convection cooled TRIGA reactors with hexagonal bundle
Yang, J.; Avery, M.; De Angelis, M.; Anderson, M.; Corradini, M.; Feldman, E. E.; Dunn, F. E.; Matos, J. E.
2012-07-01
A three-rod bundle Critical Heat Flux (CHF) study at low flow, low pressure, and natural convection condition has been conducted, simulating TRIGA reactors with the hexagonally configured core. The test section is a custom-made trefoil shape tube with three identical fuel pin heater rods located symmetrically inside. The full scale fuel rod is electrically heated with a chopped-cosine axial power profile. CHF experiments were carried out with the following conditions: inlet water subcooling from 30 K to 95 K; pressure from 110 kPa to 230 kPa; mass flux up to 150 kg/m{sup 2}s. About 50 CHF data points were collected and compared with a few existing CHF correlations whose application ranges are close to the testing conditions. Some tests were performed with the forced convection to identify the potential difference between the CHF under the natural convection and forced convection. The relevance of the CHF to test parameters is investigated. (authors)
Natural convection in tunnels at Yucca Mountain and impact on drift seepage
Halecky, N.; Birkholzer, J.T.; Peterson, P.
2010-04-15
The decay heat from radioactive waste that is to be disposed in the once proposed geologic repository at Yucca Mountain (YM) will significantly influence the moisture conditions in the fractured rock near emplacement tunnels (drifts). Additionally, large-scale convective cells will form in the open-air drifts and will serve as an important mechanism for the transport of vaporized pore water from the fractured rock in the drift center to the drift end. Such convective processes would also impact drift seepage, as evaporation could reduce the build up of liquid water at the tunnel wall. Characterizing and understanding these liquid water and vapor transport processes is critical for evaluating the performance of the repository, in terms of water-induced canister corrosion and subsequent radionuclide containment. To study such processes, we previously developed and applied an enhanced version of TOUGH2 that solves for natural convection in the drift. We then used the results from this previous study as a time-dependent boundary condition in a high-resolution seepage model, allowing for a computationally efficient means for simulating these processes. The results from the seepage model show that cases with strong natural convection effects are expected to improve the performance of the repository, since smaller relative humidity values, with reduced local seepage, form a more desirable waste package environment.
NASA Astrophysics Data System (ADS)
Wang, Zheng; Zhang, Yeshun; Zhang, Jinxiang; Huang, Lei; Liu, Jia; Li, Yongkui; Zhang, Guozheng; Kundu, Subhas C.; Wang, Lin
2014-11-01
Sericin, a major component of silk, has a long history of being discarded as a waste during silk processing. The value of sericin for tissue engineering is underestimated and its potential application in regenerative medicine has just begun to be explored. Here we report the successful fabrication and characterization of a covalently-crosslinked 3D pure sericin hydrogel for delivery of cells and drugs. This hydrogel is injectable, permitting its implantation through minimally invasive approaches. Notably, this hydrogel is found to exhibit photoluminescence, enabling bioimaging and in vivo tracking. Moreover, this hydrogel system possesses excellent cell-adhesive capability, effectively promoting cell attachment, proliferation and long-term survival of various types of cells. Further, the sericin hydrogel releases bioactive reagents in a sustained manner. Additionally, this hydrogel demonstrates good elasticity, high porosity, and pH-dependent degradation dynamics, which are advantageous for this sericin hydrogel to serve as a delivery vehicle for cells and therapeutic drugs. With all these unique features, it is expected that this sericin hydrogel will have wide utility in the areas of tissue engineering and regenerative medicine.
Wang, Zheng; Zhang, Yeshun; Zhang, Jinxiang; Huang, Lei; Liu, Jia; Li, Yongkui; Zhang, Guozheng; Kundu, Subhas C.; Wang, Lin
2014-01-01
Sericin, a major component of silk, has a long history of being discarded as a waste during silk processing. The value of sericin for tissue engineering is underestimated and its potential application in regenerative medicine has just begun to be explored. Here we report the successful fabrication and characterization of a covalently-crosslinked 3D pure sericin hydrogel for delivery of cells and drugs. This hydrogel is injectable, permitting its implantation through minimally invasive approaches. Notably, this hydrogel is found to exhibit photoluminescence, enabling bioimaging and in vivo tracking. Moreover, this hydrogel system possesses excellent cell-adhesive capability, effectively promoting cell attachment, proliferation and long-term survival of various types of cells. Further, the sericin hydrogel releases bioactive reagents in a sustained manner. Additionally, this hydrogel demonstrates good elasticity, high porosity, and pH-dependent degradation dynamics, which are advantageous for this sericin hydrogel to serve as a delivery vehicle for cells and therapeutic drugs. With all these unique features, it is expected that this sericin hydrogel will have wide utility in the areas of tissue engineering and regenerative medicine. PMID:25412301
Wang, Zheng; Zhang, Yeshun; Zhang, Jinxiang; Huang, Lei; Liu, Jia; Li, Yongkui; Zhang, Guozheng; Kundu, Subhas C; Wang, Lin
2014-01-01
Sericin, a major component of silk, has a long history of being discarded as a waste during silk processing. The value of sericin for tissue engineering is underestimated and its potential application in regenerative medicine has just begun to be explored. Here we report the successful fabrication and characterization of a covalently-crosslinked 3D pure sericin hydrogel for delivery of cells and drugs. This hydrogel is injectable, permitting its implantation through minimally invasive approaches. Notably, this hydrogel is found to exhibit photoluminescence, enabling bioimaging and in vivo tracking. Moreover, this hydrogel system possesses excellent cell-adhesive capability, effectively promoting cell attachment, proliferation and long-term survival of various types of cells. Further, the sericin hydrogel releases bioactive reagents in a sustained manner. Additionally, this hydrogel demonstrates good elasticity, high porosity, and pH-dependent degradation dynamics, which are advantageous for this sericin hydrogel to serve as a delivery vehicle for cells and therapeutic drugs. With all these unique features, it is expected that this sericin hydrogel will have wide utility in the areas of tissue engineering and regenerative medicine. PMID:25412301
On the convectively unstable nature of optimal streaks in boundary layers
NASA Astrophysics Data System (ADS)
Brandt, Luca; Cossu, Carlo; Chomaz, Jean-Marc; Huerre, Patrick; Henningson, Dan S.
2003-06-01
The objective of the study is to determine the absolute/convective nature of the secondary instability experienced by finite-amplitude streaks in the flat-plate boundary layer. A family of parallel streaky base flows is defined by extracting velocity profiles from direct numerical simulations of nonlinearly saturated optimal streaks. The computed impulse response of the streaky base flows is then determined as a function of streak amplitude and streamwise station. Both the temporal and spatio-temporal instability properties are directly retrieved from the impulse response wave packet, without solving the dispersion relation or applying the pinching point criterion in the complex wavenumber plane. The instability of optimal streaks is found to be unambiguously convective for all streak amplitudes and streamwise stations. It is more convective than the Blasius boundary layer in the absence of streaks; the trailing edge-velocity of a Tollmien Schlichting wave packet in the Blasius boundary layer is around 35% of the free-stream velocity, while that of the wave packet riding on the streaky base flow is around 70%. This is because the streak instability is primarily induced by the spanwise shear and the associated Reynolds stress production term is located further away from the wall, in a larger velocity region, than for the Tollmien Schlichting instability. The streak impulse response consists of the sinuous mode of instability triggered by the spanwise wake-like profile, as confirmed by comparing the numerical results with the absolute/convective instability properties of the family of two-dimensional wakes introduced by Monkewitz (1988). The convective nature of the secondary streak instability implies that the type of bypass transition studied here involves streaks that behave as amplifiers of external noise.
Dudek, D.; Fletcher, T.H.
1987-02-01
When a heated solid sphere is introduced into an ambient fluid, a natural convective flow occurs which results in a drag force on the sphere. This study involves the numerical calculation of both the steady-state and the transient natural convective drag force around spheres at low Grashof numbers. Numerical techniques are taken from Geoola and Cornish. An empirical expression is suggested for the total drag coefficient for Grashof numbers ranging from 4 x 10/sup -4/ to 0.5 and Prandtl number = 0.72: log C/sub DT/ = 1.25 + 0.31 log Gr - 0.097(log Gr)/sup 2/. The dimensionless time required to reach 90% of the steady-state drag force can be approximated by the second-order polynomial: log t/sub 90%/ = 1.32 - log Gr - 0.11(Gr)/sup 2/.
Tzeng, P.Y.; Soong, C.Y.; Sheu, T.S.
1997-02-07
The present work is concerned with a numerical investigation of transient laminar natural convection and the associated flow-mode transition in a two-dimensional rectangular enclosure. Navier-Stokes/Boussinesq equations for fluid flow and energy balance are solved by using the SIMPLE-C algorithm. Air of Pr = 0.71 in a differentially heated enclose of length-to-height aspect ratio As = 4 and at Ra = 5,000 is chosen as the flow model to examine the influences of the inclination. Calculations of time accuracy are performed to investigate the transient procedure of the flow-mode transition with increasing or decreasing inclination. The present results reveal that, at some critical situations, natural convection in inclined enclosures is very sensitive to the change in tilt angle, and the associated heat transfer rates are closely related to the correspondent cellular flow patterns.
Scaling of the turbulent natural convection flow in a heated square cavity
NASA Astrophysics Data System (ADS)
Henkes, R. A. W. M.; Hoogendoorn, C. J.
1994-05-01
By numerically solving the Reynolds equations for air and water in a square cavity, with differentially heated vertical walls, at Rayleigh numbers up to 10(exp 20) the scalings of the turbulent natural convection flow are derived. Turbulence is modeled by the standard k-epsilon model and by the low-Reynolds-number k-epsilon models of Chien and of Jones and Launder. Both the scalings with respect to the Rayleigh number (based on the cavity size H) and with respect to the local height (y/H) are considered. The scalings are derived for the inner layer, outer layer, and core region. The Rayleigh number scalings are almost the same as the scalings for the natural convection boundary layer along a hot vertical plate. The scalings found are almost independent of the k-epsilon model used.
Tuncay, K.; Romer, S.; Ortoleva, P.; Hoak, T.; Sundberg, K.
1998-12-31
The power of the reaction, transport, mechanical (RTM) modeling approach is that it directly uses the laws of geochemistry and geophysics to extrapolate fracture and other characteristics from the borehole or surface to the reservoir interior. The objectives of this facet of the project were to refine and test the viability of the basin/reservoir forward modeling approach to address fractured reservoir in E and P problems. The study attempts to resolve the following issues: role of fracturing and timing on present day location and characteristics; clarifying the roles and interplay of flexure dynamics, changing rock rheological properties, fluid pressuring and tectonic/thermal histories on present day reservoir location and characteristics; and test the integrated RTM modeling/geological data approach on a carbonate reservoir. Sedimentary, thermal and tectonic data from Andector Field, West Texas, were used as input to the RTM basin/reservoir simulator to predict its preproduction state. The results were compared with data from producing reservoirs to test the RTM modeling approach. The effects of production on the state of the field are discussed in a companion report. The authors draw the following conclusions: RTM modeling is an important new tool in fractured reservoir E and P analysis; the strong coupling of RTM processes and the geometric and tensorial complexity of fluid flow and stresses require the type of fully coupled, 3-D RTM model for fracture analysis as pioneered in this project; flexure analysis cannot predict key aspects of fractured reservoir location and characteristics; fracture history over the lifetime of a basin is required to understand the timing of petroleum expulsion and migration and the retention properties of putative reservoirs.
Calculation of Post-Closure Natural Convection Heat and Mass Transfer in Yucca Mountain Drifts
S. Webb; M. Itamura
2004-03-16
Natural convection heat and mass transfer under post-closure conditions has been calculated for Yucca Mountain drifts using the computational fluid dynamics (CFD) code FLUENT. Calculations have been performed for 300, 1000, 3000, and 10,000 years after repository closure. Effective dispersion coefficients that can be used to calculate mass transfer in the drift have been evaluated as a function of time and boundary temperature tilt.
Natural convection within a vertical finite-length channel in free space
Lin, S.C.; Chang, K.P.; Hung, Y.H. )
1994-04-01
Natural convection within a vertical finite length channel in free space is studied in this article to remove assumptions that need to be made on velocity and temperature profiles at the channel entrance. For small channel aspect ratios and low Rayleigh numbers, significant deviations of the Nusselt number and temperature distributions exist due to the effects of vertical thermal diffusion and free space stratification in the channel. A new correlation was proposed on induced Reynolds number for vertical finite length channel. 8 refs.
CFD numerical simulation of air natural convection over a heated cylindrical surface
NASA Astrophysics Data System (ADS)
Flori, M.; Vîlceanu, L.
2015-06-01
In this study a CFD numerical simulation is used to describe the fluid flow and heat transfer in air surrounding a heated horizontal cylinder. The model is created in 2D space dimension involving a finite element solver of Navier-Stokes equations. As natural convection phenomenon is induced by a variable fluid density field with temperature rising, the Boussinesq approximation was coupled to the model.
Simulation of natural convection in a rectangular loop using finite elements
Pepper, D W; Hamm, L L; Kehoe, A B
1984-01-01
A two-dimensional finite-element analysis of natural convection in a rectangular loop is presented. A psi-omega formulation of the Boussinesque approximation to the Navier-Stokes equation is solved by the false transient technique. Streamlines and isotherms at Ra = 10/sup 4/ are shown for three different modes of heating. The results indicate that corner effects should be considered when modeling flow patterns in thermosyphons.
Delmas, A.A.; Wilkes, K.E.
1992-04-01
A two-dimensional code for solving equations of convective heat transfer in porous media is used to analyze heat transfer by conduction and convection in the attic insulation configuration. The particular cases treated correspond to loose-fill fiberglass insulation, which is characterized by high porosity and air permeability. The effects of natural convection on the thermal performance of the insulation are analyzed for various densities, permeabilities, and thicknesses of insulation. With convection increasing the total heat transfer through the insulation, the thermal resistance was found to decrease as the temperature difference across the insulating material increases. The predicted results for the thermal resistance are compared with data obtained in the large-scale climate simulator at the Roof Research Center using the attic test module, where the same phenomenon has already been observed. The way the wood joists within the insulation influence the start of convection is studied for differing thermophysical and dynamic properties of the insulating material. The presence of wood joists induces convection at a lower temperature difference.
NASA Astrophysics Data System (ADS)
Bower, S. M.; Saylor, J. R.
2009-11-01
Presented are the results from an experimental investigation of the effects of surface conditions at an air/water interface on transport phenomena within the context of natural convection-driven evaporation. Experiments were conducted using tanks of heated water under several different surface conditions: 1) contamination with an oleyl alcohol monolayer, 2) contamination with a stearic acid monolayer, and 3) ``clean'' or surfactant-free. These surface conditions create the following hydrodynamic boundary conditions: 1) constant elasticity, 2) no-slip, and 3) shear-free. The effect of these boundary conditions on evaporation and air-side natural convection heat transfer is presented via the power law relationships between the Sherwood and Rayleigh numbers (for evaporation) and the Nusselt and Rayleigh numbers (for natural convection heat transfer). Additionally, infrared imagery of the water surface was collected during these experiments, yielding qualitative information on the effect of these boundary conditions on the flow near the interface. Few studies exist in which the effects of surface conditions on interfacial heat and mass transfer are investigated, making this work particularly relevant.
Emergency cooling down of fast-neutron reactors by natural convection (a review)
NASA Astrophysics Data System (ADS)
Zhukov, A. V.; Sorokin, A. P.; Kuzina, Yu. A.
2013-05-01
Various methods for emergency cooling down of fast-neutron reactors by natural convection are discussed. The effectiveness of using natural convection for these purposes is demonstrated. The operating principles of different passive decay heat removal systems intended for cooling down a reactor are explained. Experimental investigations carried out in Russia for substantiating the removal of heat in cooling down fast-neutron reactors are described. These investigations include experimental works on studying thermal hydraulics in small-scale simulation facilities containing the characteristic components of a reactor (reactor core elements, above-core structure, immersed and intermediate heat exchangers, pumps, etc.). It is pointed out that a system that uses leaks of coolant between fuel assemblies holds promise for fast-neutron reactor cooldown purposes. Foreign investigations on this problem area are considered with making special emphasis on the RAMONA and NEPTUN water models. A conclusion is drawn about the possibility of using natural convection as the main method for passively removing heat in cooling down fast-neutron reactors, which is confirmed experimentally both in Russia and abroad.
Air cooling of a vented enclosure by combined conduction, natural convection and radiation
Yu, E.; Joshi, Y.K.
1996-12-31
A three-dimensional investigation of combined conduction, natural convection and radiation in vented enclosures is carried out. A discrete flush type heat source mounted on a vertical substrate is used to simulate an electronic component. A uniform volumetric generation rate is assumed within the heat source. Combined natural convection in the air, conduction in the heat source, the substrate and the enclosure walls, and surface radiation are solved for Rayleigh numbers at 2.6 {times} 10{sup 6} and 2.0 {times} 10{sup 7}. Radiation is incorporated based on the radiosity/irradiation approach. The resulting flow and temperature patterns are discussed, focusing on radiation and three-dimensional effects. The relative contributions of natural convection and radiation are investigated for different emissivities of internal surface of the substrate. Heat transfer rates from the substrate and other internal walls are presented to illustrate conjugate heat transfer due to combined modes. The numerical solutions are found in reasonably good agreement with the data.
Zhang, Liwei; Anderson, Nicole; Dilmore, Robert; Soeder, Daniel J; Bromhal, Grant
2014-09-16
Potential natural gas leakage into shallow, overlying formations and aquifers from Marcellus Shale gas drilling operations is a public concern. However, before natural gas could reach underground sources of drinking water (USDW), it must pass through several geologic formations. Tracer and pressure monitoring in formations overlying the Marcellus could help detect natural gas leakage at hydraulic fracturing sites before it reaches USDW. In this study, a numerical simulation code (TOUGH 2) was used to investigate the potential for detecting leaking natural gas in such an overlying geologic formation. The modeled zone was based on a gas field in Greene County, Pennsylvania, undergoing production activities. The model assumed, hypothetically, that methane (CH4), the primary component of natural gas, with some tracer, was leaking around an existing well between the Marcellus Shale and the shallower and lower-pressure Bradford Formation. The leaky well was located 170 m away from a monitoring well, in the Bradford Formation. A simulation study was performed to determine how quickly the tracer monitoring could detect a leak of a known size. Using some typical parameters for the Bradford Formation, model results showed that a detectable tracer volume fraction of 2.0 × 10(-15) would be noted at the monitoring well in 9.8 years. The most rapid detection of tracer for the leak rates simulated was 81 days, but this scenario required that the leakage release point was at the same depth as the perforation zone of the monitoring well and the zones above and below the perforation zone had low permeability, which created a preferred tracer migration pathway along the perforation zone. Sensitivity analysis indicated that the time needed to detect CH4 leakage at the monitoring well was very sensitive to changes in the thickness of the high-permeability zone, CH4 leaking rate, and production rate of the monitoring well. PMID:25144442
NASA Astrophysics Data System (ADS)
Gholizadeh Doonechaly, N.; Rahman, S. S.
2012-05-01
Simulation of naturally fractured reservoirs offers significant challenges due to the lack of a methodology that can utilize field data. To date several methods have been proposed by authors to characterize naturally fractured reservoirs. Among them is the unfolding/folding method which offers some degree of accuracy in estimating the probability of the existence of fractures in a reservoir. Also there are statistical approaches which integrate all levels of field data to simulate the fracture network. This approach, however, is dependent on the availability of data sources, such as seismic attributes, core descriptions, well logs, etc. which often make it difficult to obtain field wide. In this study a hybrid tectono-stochastic simulation is proposed to characterize a naturally fractured reservoir. A finite element based model is used to simulate the tectonic event of folding and unfolding of a geological structure. A nested neuro-stochastic technique is used to develop the inter-relationship between the data and at the same time it utilizes the sequential Gaussian approach to analyze field data along with fracture probability data. This approach has the ability to overcome commonly experienced discontinuity of the data in both horizontal and vertical directions. This hybrid technique is used to generate a discrete fracture network of a specific Australian gas reservoir, Palm Valley in the Northern Territory. Results of this study have significant benefit in accurately describing fluid flow simulation and well placement for maximal hydrocarbon recovery.
Aksenova, A.E.; Chudanov, V.V.; Strizhov, V.F.; Vabishchevich, P.N.
1995-09-01
Unsteady natural convection of a heat-generating fluid with phase transitions in the enclosures of a square section with isothermal rigid walls is investigated numerically for a wide range of dimensionless parameters. The quasisteady state solutions of conjugate heat and mass transfer problem are compared with available experimental results. Correlation relations for heat flux distributions at the domain boundaries depending on Rayleigh and Ostrogradskii numbers are obtained. It is shown that generally heat transfer is governed both by natural circulation and crust formation phenomena. Results of this paper may be used for analysis of experiments with prototypic core materials.
Finite element, stream function-vorticity solution of steady laminar natural convection
NASA Astrophysics Data System (ADS)
Stevens, W. N. R.
1982-12-01
Stream function-vorticity finite element solution of two-dimensional incompressible viscous flow and natural convection is considered. Steady state solutions of the natural convection problem have been obtained for a wide range of the two independent parameters. Use of boundary vorticity formulae or iterative satisfaction of the no-slip boundary condition is avoided by application of the finite element discretization and a displacement of the appropriate discrete equations. Solution is obtained by Newton-Raphson iteration of all equations simultaneously. The method then appears to give a steady solution whenever the flow is physically steady, but it does not give a steady solution when the flow is physically unsteady. In particular, no form of asymmetric differencing is required. The method offers a degree of economy over primitive variable formulations. Physical results are given for the square cavity convection problem. The paper also reports on earlier work in which the most commonly used boundary vorticity formula was found not to satisfy the no-slip condition, and in which segregated solution procedures were attempted with very minimal success.
Triplett, C.E.
1996-12-01
This thesis presents the results of an experimental investigation of natural convection heat transfer in a staggered array of heated cylinders, oriented horizontally within a rectangular enclosure. The main purpose of this research was to extend the knowledge of heat transfer within enclosed bundles of spent nuclear fuel rods sealed within a shipping or storage container. This research extends Canaan`s investigation of an aligned array of heated cylinders that thermally simulated a boiling water reactor (BWR) spent fuel assembly sealed within a shipping or storage cask. The results are presented in terms of piecewise Nusselt-Rayleigh number correlations of the form Nu = C(Ra){sup n}, where C and n are constants. Correlations are presented both for individual rods within the array and for the array as a whole. The correlations are based only on the convective component of the heat transfer. The radiative component was calculated with a finite-element code that used measured surface temperatures, rod array geometry, and measured surface emissivities as inputs. The correlation results are compared to Canaan`s aligned array results and to other studies of natural convection in horizontal tube arrays.
NASA Astrophysics Data System (ADS)
Pletinckx, D.
2011-09-01
The current 3D hype creates a lot of interest in 3D. People go to 3D movies, but are we ready to use 3D in our homes, in our offices, in our communication? Are we ready to deliver real 3D to a general public and use interactive 3D in a meaningful way to enjoy, learn, communicate? The CARARE project is realising this for the moment in the domain of monuments and archaeology, so that real 3D of archaeological sites and European monuments will be available to the general public by 2012. There are several aspects to this endeavour. First of all is the technical aspect of flawlessly delivering 3D content over all platforms and operating systems, without installing software. We have currently a working solution in PDF, but HTML5 will probably be the future. Secondly, there is still little knowledge on how to create 3D learning objects, 3D tourist information or 3D scholarly communication. We are still in a prototype phase when it comes to integrate 3D objects in physical or virtual museums. Nevertheless, Europeana has a tremendous potential as a multi-facetted virtual museum. Finally, 3D has a large potential to act as a hub of information, linking to related 2D imagery, texts, video, sound. We describe how to create such rich, explorable 3D objects that can be used intuitively by the generic Europeana user and what metadata is needed to support the semantic linking.
An analysis of natural convection film boiling from spheres using the spherical coordinate system
Tso, C.P.; Leong, K.C.; Tan, H.S.
1995-11-01
The problem of natural convection film boiling on a sphere was analyzed by solving the momentum and energy equations in spherical coordinates. These solutions were compared to the analytical model of Frederking and Clark based on the Cartesian coordinate system, empirical correlation of Frederking and Clark and recent experimental data of Tso et al. for boiling in various refrigerants and liquid nitrogen. For the average Nusselt number, good agreement with Frederking and Clark`s model was obtained. Results using spherical coordinates yield a limiting value of 2 for the average Nusselt number near a modified Rayleigh number of 1 which could not be extracted from Frederking and Clark`s model.
NASA Astrophysics Data System (ADS)
Raju, S. Suresh Kumar; Narahari, M.; Pendyala, Rajashekhar
2014-10-01
In this paper, a numerical solution of the unsteady two-dimensional natural convection along a vertical plate in the presence of Soret and chemical reaction effects is presented. The governing non-dimensional coupled non-linear partial differential equations have been evaluated by using an implicit finite-difference technique of Crank-Nicolson scheme. Numerical predictions for the velocity, concentration, local and average skin-friction and Sherwood number for distinct values of chemical reaction parameter and Soret number are plotted graphically. It is found that the fluid velocity and concentration decreases while increasing chemical reaction parameter whereas an increase in the Soret number increases the fluid velocity and concentration.
Effect of free surface shape on combined thermocapillary and natural convection
NASA Technical Reports Server (NTRS)
Kamotani, Yasuhiro; Platt, Jonathan
1992-01-01
Combined thermocapillary and natural convection in an open square cavity with differentially-heated side walls is studied numerically as well as experimentally. The test fluid is silicone oil with Prandtl number of 105. The shape of fluid-free surface is made either flat or curved to study its effect on the flow. A finite difference scheme to deal with a curved free surface is developed. The experimental results shown agree with the numerical results. With the curved-free surface, the flow and local heat transfer rate are reduced in the corner regions, and a sharp peak in heat transfer rate at the top edge of the cold wall disappears.
Passive decay heat removal by natural air convection after severe accidents
Erbacher, F.J.; Neitzel, H.J.; Cheng, X.
1995-09-01
The composite containment proposed by the Research Center Karlsruhe and the Technical University Karlsruhe is to cope with severe accidents. It pursues the goal to restrict the consequences of core meltdown accidents to the reactor plant. One essential of this new containment concept is its potential to remove the decay heat by natural air convection and thermal radiation in a passive way. To investigate the coolability of such a passive cooling system and the physical phenomena involved, experimental investigations are carried out at the PASCO test facility. Additionally, numerical calculations are performed by using different codes. A satisfying agreement between experimental data and numerical results is obtained.
NASA Astrophysics Data System (ADS)
Li, Chunggang; Tsubokura, Makoto; Complex Phenomena Unified Simulation Research Team
2014-11-01
The complete transition from laminar to turbulent natural convection in a long channel is investigated using compressible direct numerical simulation (DNS). Numerical methods of Roe scheme with precontioning and dual time stepping are used for addressing the flow field which is low speed but the density is variable. During the transient development, there are four stages which are laminar, unstable process, relaminarization and turbulence can be obviously identified. After reaching the quasi steady state, the laminar, transition and turbulence simultaneously coexist in the same flow field. Additionally, the comparisons of the statistics with the experimental data are also well consistent.
Anomalies of the natural convection of water near 3.98°C
NASA Astrophysics Data System (ADS)
Baturov, L. N.; Govor, I. N.
2016-02-01
Natural convection of water in a cylindrical cavity with an open surface at a temperature of about 3.98°C (temperature of the maximum water density) is accompanied by typical anomalies on time dependences of temperatures of water layers. In particular, stabilization of temperature T st is observed in the bottom region of the cavity and duration of such stabilization t st may reach several hours depending on the experimental conditions. The results for solutions of sodium chloride and ethanol at a relatively low rate of water cooling show that temperature T st coincides with temperature T max corresponding to the maximum density of solutions.
NASA Astrophysics Data System (ADS)
Haddad, Zoubida; Abid, Chérifa; Mohamad, A. A.; Rahli, O.; Bawazer, S.
2016-08-01
An experimental and numerical study was performed to investigate the effect of different formulas for nanofluid thermal conductivity and dynamic viscosity on natural convective heat transfer. It was found that the heat transfer across the enclosure using different models can be enhanced or deteriorated with respect to the base fluid. Also, it was found that the inconsistencies in the reported thermal conductivity and dynamic viscosity from different research groups are mainly due to the characterization of the nanofluid, including determination of colloidal stability and particle size, (i.e., aggregates size) within nanofluid.
Drying characteristic of barley under natural convection in a mixed-mode type solar grain dryer
Basunia, M.A.; Abe, T.
1999-07-01
Thin-layer solar drying characteristics of barley were determined at average natural air flow temperature ranging from 43.4 to 51.7 C and for relative humidities ranging from 16.5% to 37.5%. A mixed-mode type natural convection solar dryer was used for this experiment. The data of sample weight, and dry and wet bulb temperatures of the drying air were recorded continuously throughout the drying period for each test. The drying data were then fitted to the Page model. The model gave a good fit for the moisture content with an average standard error of 0.305% dry basis. The parameter N in Page's equation was assumed as a product-dependent constant which made it easy to compare the effects of independent variables on the natural convection solar drying rate without causing considerable error in predicting the drying rate for barley. A linear relationship was found between the parameter K, temperature T, and relative humidity R{sub H}.
Lee, Il S.; Yu, Yong H.; Son, Hyoung M.; Hwang, Jin S.; Suh, Kune Y.
2006-07-01
An experimental study is performed to investigate the natural convection heat transfer characteristics with subcooled coolant to create engineering database for basic applications in a lead alloy cooled reactor. Tests are performed in the ALTOS (Applied Liquid-metal Thermal Operation Study) apparatus as part of MITHOS (Metal Integrated Thermo Hydrodynamic Operation System). A relationship is determined between the Nusselt number Nu and the Rayleigh number Ra in the liquid metal rectangular pool. Results are compared with correlations and experimental data in the literature. Given the similar Ra condition, the present test results for Nu of the liquid metal pool with top subcooling are found to be similar to those predicted by the existing correlations or experiments. The current test results are utilized to develop natural convection heat transfer correlations applicable to low Prandtl number Pr fluids that are heated from below and cooled by the external coolant above. Results from this study are slated to be used in designing BORIS (Battery Optimized Reactor Integral System), a small lead cooled modular fast reactor for deployment at remote sites cycled with MOBIS (Modular Optimized Brayton Integral System) for electricity generation, tied with NAVIS (Naval Application Vessel Integral System) for ship propulsion, joined with THAIS (Thermochemical Hydrogen Acquisition Integral System) for hydrogen production, and coupled with DORIS (Desalination Optimized Reactor Integral System) for seawater desalination. Tests are performed with Wood's metal (Pb-Bi-Sn-Cd) filling a rectangular pool whose lower surface is heated and upper surface cooled by forced convection of water. The test section is 20 cm long, 11.3 cm high and 15 cm wide. The simulant has a melting temperature of 78 deg. C. The constant temperature and heat flux condition was realized for the bottom heating once the steady state had been met. The test parameters include the heated bottom surface temperature
Mayor, T S; Couto, S; Psikuta, A; Rossi, R M
2015-12-01
The ability of clothing to provide protection against external environments is critical for wearer's safety and thermal comfort. It is a function of several factors, such as external environmental conditions, clothing properties and activity level. These factors determine the characteristics of the different microclimates existing inside the clothing which, ultimately, have a key role in the transport processes occurring across clothing. As an effort to understand the effect of transport phenomena in clothing microclimates on the overall heat transport across clothing structures, a numerical approach was used to study the buoyancy-driven heat transfer across horizontal air layers trapped inside air impermeable clothing. The study included both the internal flow occurring inside the microclimate and the external flow occurring outside the clothing layer, in order to analyze the interdependency of these flows in the way heat is transported to/from the body. Two-dimensional simulations were conducted considering different values of microclimate thickness (8, 25 and 52 mm), external air temperature (10, 20 and 30 °C), external air velocity (0.5, 1 and 3 m s(-1)) and emissivity of the clothing inner surface (0.05 and 0.95), which implied Rayleigh numbers in the microclimate spanning 4 orders of magnitude (9 × 10(2)-3 × 10(5)). The convective heat transfer coefficients obtained along the clothing were found to strongly depend on the transport phenomena in the microclimate, in particular when natural convection is the most important transport mechanism. In such scenario, convective coefficients were found to vary in wavy-like manner, depending on the position of the flow vortices in the microclimate. These observations clearly differ from data in the literature for the case of air flow over flat-heated surfaces with constant temperature (which shows monotonic variations of the convective heat transfer coefficients, along the length of the surface). The flow
NASA Astrophysics Data System (ADS)
Mayor, T. S.; Couto, S.; Psikuta, A.; Rossi, R. M.
2015-12-01
The ability of clothing to provide protection against external environments is critical for wearer's safety and thermal comfort. It is a function of several factors, such as external environmental conditions, clothing properties and activity level. These factors determine the characteristics of the different microclimates existing inside the clothing which, ultimately, have a key role in the transport processes occurring across clothing. As an effort to understand the effect of transport phenomena in clothing microclimates on the overall heat transport across clothing structures, a numerical approach was used to study the buoyancy-driven heat transfer across horizontal air layers trapped inside air impermeable clothing. The study included both the internal flow occurring inside the microclimate and the external flow occurring outside the clothing layer, in order to analyze the interdependency of these flows in the way heat is transported to/from the body. Two-dimensional simulations were conducted considering different values of microclimate thickness (8, 25 and 52 mm), external air temperature (10, 20 and 30 °C), external air velocity (0.5, 1 and 3 m s-1) and emissivity of the clothing inner surface (0.05 and 0.95), which implied Rayleigh numbers in the microclimate spanning 4 orders of magnitude (9 × 102-3 × 105). The convective heat transfer coefficients obtained along the clothing were found to strongly depend on the transport phenomena in the microclimate, in particular when natural convection is the most important transport mechanism. In such scenario, convective coefficients were found to vary in wavy-like manner, depending on the position of the flow vortices in the microclimate. These observations clearly differ from data in the literature for the case of air flow over flat-heated surfaces with constant temperature (which shows monotonic variations of the convective heat transfer coefficients, along the length of the surface). The flow patterns and
Stanton, M M; Samitier, J; Sánchez, S
2015-08-01
Three-dimensional (3D) bioprinting has recently emerged as an extension of 3D material printing, by using biocompatible or cellular components to build structures in an additive, layer-by-layer methodology for encapsulation and culture of cells. These 3D systems allow for cell culture in a suspension for formation of highly organized tissue or controlled spatial orientation of cell environments. The in vitro 3D cellular environments simulate the complexity of an in vivo environment and natural extracellular matrices (ECM). This paper will focus on bioprinting utilizing hydrogels as 3D scaffolds. Hydrogels are advantageous for cell culture as they are highly permeable to cell culture media, nutrients, and waste products generated during metabolic cell processes. They have the ability to be fabricated in customized shapes with various material properties with dimensions at the micron scale. 3D hydrogels are a reliable method for biocompatible 3D printing and have applications in tissue engineering, drug screening, and organ on a chip models. PMID:26066320
Experimental study of natural convection enhancement using a Fe3O4-water based magnetic nanofluid.
Stoian, Floriana D; Holotescu, Sorin
2012-10-01
The effect of nanoparticles dispersed in a carrier fluid on the natural convection heat transfer is still raising controversies. While the reported experimental results show no improvement or even worsening of the heat transfer performance of nanofluids, the numerical simulations show an increase of the heat transfer coefficient, at least for certain ranges of Ra number. We report an experimental investigation regarding the natural convection heat transfer performance of a Fe3O4-water based nanofluid, in a cylindrical enclosure. The fluid was heated linearly from the bottom wall using an electric heater and cooled from the upper wall by a constant flow of water, such that a constant temperature difference between the upper and bottom walls was obtained at steady-state. The experiment was also carried out using water, in order to observe the effect of the addition of Fe3O4 nanoparticles on the heat transfer coefficient. Several regimes were tested, both for water and nanofluid. The experimental results showed that values obtained for the heat transfer coefficient for Fe3O4-water nanofluid were higher than those for water, at the same temperature difference. The present experimental results are also compared with our previous work and the reference literature. PMID:23421199
A new look at natural convection from isothermal vertical parallel plates
Li, H.H.; Chung, B.T.F.
1996-12-31
Natural convection between isothermal plates is solved numerically by applying the full Navier-Stokes equations. The elliptic formulation allows separating the effect of the Rayleigh number, Ra, and the aspect ratio, L/B. Calculations are made on a wide range of the Rayleigh number and the aspect ratio, and the Nusselt number is provided as a function of both Ra and B/L. The conventional correlations in the literature presenting the Nusselt number in terms of a single parameter, RaB/L, have been found inaccurate. At a small value of RaB/L, multiple values of Nusselt number are obtained for different combinations of Ra and B/L. Previous results are found to be the special cases of the present study. A minimum Rayleigh number is also obtained above which a fully-developed flow is possible. To simulate the natural convective flow, the ambient pressure is given at the exit while the pressure at the entrance is related to the ambient pressure by the Bernoulli equation. Velocities at the entrance and exit are also solved from the Navier-Stokes equations.
NASA Astrophysics Data System (ADS)
Ambrosini, Dario; Tanda, Giovanni
2006-01-01
In this work, natural convection heat transfer in vertical channels is experimentally investigated by applying different optical techniques, namely holographic interferometry and schlieren. Both these techniques are based on the temperature dependence of the air refractive index but they detect different optical quantities and their use involves different instrumentation and optical components. Optical methods, non-intrusive in nature, are particularly suitable for the visualization of flow and thermal fields as witnessed by their increasing use in a range of scientific and engineering disciplines; for this reason, the introduction of these experimental tools into a laboratory course can be of high value. Physics and engineering students can get familiarized with optical techniques, grasp the basics of thermal phenomena, usually elusive, which can be more easily understood if they are made visible, and begin to master digital image analysis, a key skill in laboratory activities. A didactic description of holographic interferometry and schlieren is provided and experimental results obtained for vertical, smooth and rib-roughened channels with asymmetrical heating are presented. A comparison between distributions of the local heat transfer coefficient (or its dimensionless counterpart, the Nusselt number) revealed good agreement between the results separately obtained by the two techniques, thus proving their suitability for investigating free convection heat transfer in channels.
Kang, S.; Ha, K. S.; Lee, S. W.; Park, S. D.; Kim, S. M.; Seo, H.; Kim, J. H.; Bang, I. C.
2012-07-01
The safety issues of the SFRs are important due to the fact that it uses sodium as a nuclear coolant, reacting vigorously with water and air. For that reason, there are efforts to seek for alternative candidates of liquid metal coolants having excellent heat transfer property and to adopt improved safety features to the SFR concepts. This study considers gallium as alternative liquid metal coolant applicable to safety features in terms of chemical activity issue of the sodium and aims to experimentally investigate the natural convection capability of gallium as a feasibility study for the development of gallium-based passive safety features in SFRs. In this paper, the design and construction of the liquid gallium natural convection loop were carried out. The experimental results of heat transfer coefficient of liquid gallium resulting in heat removal {approx}2.53 kW were compared with existing correlations and they were much lower than the correlations. To comparison of the experimental data with computer code analysis, gallium property code was developed for employing MARS-LMR (Korea version of RELAP) based on liquid gallium as working fluid. (authors)
3d-3d correspondence revisited
NASA Astrophysics Data System (ADS)
Chung, Hee-Joong; Dimofte, Tudor; Gukov, Sergei; Sułkowski, Piotr
2016-04-01
In fivebrane compactifications on 3-manifolds, we point out the importance of all flat connections in the proper definition of the effective 3d {N}=2 theory. The Lagrangians of some theories with the desired properties can be constructed with the help of homological knot invariants that categorify colored Jones polynomials. Higgsing the full 3d theories constructed this way recovers theories found previously by Dimofte-Gaiotto-Gukov. We also consider the cutting and gluing of 3-manifolds along smooth boundaries and the role played by all flat connections in this operation.
Natural and mixed convection in the cylindrical pool of TRIGA reactor
NASA Astrophysics Data System (ADS)
Henry, R.; Tiselj, I.; Matkovič, M.
2016-05-01
Temperature fields within the pool of the JSI TRIGA MARK II nuclear research reactor were measured to collect data for validation of the thermal hydraulics computational model of the reactor tank. In this context temperature of the coolant was measured simultaneously at sixty different positions within the pool during steady state operation and two transients. The obtained data revealed local peculiarities of the cooling water dynamics inside the pool and were used to estimate the coolant bulk velocity above the reactor core. Mixed natural and forced convection in the pool were simulated with a Computational Fluid Dynamics code. A relatively simple CFD model based on Unsteady RANS turbulence model was found to be sufficient for accurate prediction of the temperature fields in the pool during the reactor operation. Our results show that the simple geometry of the TRIGA pool reactor makes it a suitable candidate for a simple natural circulation benchmark in cylindrical geometry.
The program FANS-3D (finite analytic numerical simulation 3-dimensional) and its applications
NASA Technical Reports Server (NTRS)
Bravo, Ramiro H.; Chen, Ching-Jen
1992-01-01
In this study, the program named FANS-3D (Finite Analytic Numerical Simulation-3 Dimensional) is presented. FANS-3D was designed to solve problems of incompressible fluid flow and combined modes of heat transfer. It solves problems with conduction and convection modes of heat transfer in laminar flow, with provisions for radiation and turbulent flows. It can solve singular or conjugate modes of heat transfer. It also solves problems in natural convection, using the Boussinesq approximation. FANS-3D was designed to solve heat transfer problems inside one, two and three dimensional geometries that can be represented by orthogonal planes in a Cartesian coordinate system. It can solve internal and external flows using appropriate boundary conditions such as symmetric, periodic and user specified.
Plume generation in natural thermal convection at high Rayleigh and Prandtl numbers
NASA Astrophysics Data System (ADS)
Lithgow-Bertelloni, C.; Richards, M. A.; Conrad, C. P.; Griffiths, R. W.
2001-05-01
We study natural thermal convection of a fluid (corn syrup) with a large Prandtl number (103 107) and temperature-dependent viscosity. The experimental tank (1 × 1 × 0.3m) is heated from below with insulating top and side boundaries, so that the fluid experiences secular heating as experiments proceed. This setup allows a focused study of thermal plumes from the bottom boundary layer over a range of Rayleigh numbers relevant to convective plumes in the deep interior of the Earth's mantle. The effective value of Ra, based on the viscosity of the fluid at the interior temperature, varies from 105 at the beginning to almost 108 toward the end of the experiments. Thermals (plumes) from the lower boundary layer are trailed by continuous conduits with long residence times. Plumes dominate flow in the tank, although there is a weaker large-scale circulation induced by material cooling at the imperfectly insulating top and sidewalls. At large Ra convection is extremely time-dependent and exhibits episodic bursts of plumes, separated by periods of quiescence. This bursting behaviour probably results from the inability of the structure of the thermal boundary layer and its instabilities to keep pace with the rate of secular change in the value of Ra. The frequency of plumes increases and their size decreases with increasing Ra, and we characterize these changes via in situ thermocouple measurements, shadowgraph videos, and videos of liquid crystal films recorded during several experiments. A scaling analysis predicts observed changes in plume head and tail radii with increasing Ra. Since inertial effects are largely absent no transition to ‘hard’ thermal turbulence is observed, in contrast to a previous conclusion from numerical calculations at similar Rayleigh numbers. We suggest that bursting behaviour similar to that observed may occur in the Earth's mantle as it undergoes secular cooling on the billion-year time scale.
TACO3D. 3-D Finite Element Heat Transfer Code
Mason, W.E.
1992-03-04
TACO3D is a three-dimensional, finite-element program for heat transfer analysis. An extension of the two-dimensional TACO program, it can perform linear and nonlinear analyses and can be used to solve either transient or steady-state problems. The program accepts time-dependent or temperature-dependent material properties, and materials may be isotropic or orthotropic. A variety of time-dependent and temperature-dependent boundary conditions and loadings are available including temperature, flux, convection, and radiation boundary conditions and internal heat generation. Additional specialized features treat enclosure radiation, bulk nodes, and master/slave internal surface conditions (e.g., contact resistance). Data input via a free-field format is provided. A user subprogram feature allows for any type of functional representation of any independent variable. A profile (bandwidth) minimization option is available. The code is limited to implicit time integration for transient solutions. TACO3D has no general mesh generation capability. Rows of evenly-spaced nodes and rows of sequential elements may be generated, but the program relies on separate mesh generators for complex zoning. TACO3D does not have the ability to calculate view factors internally. Graphical representation of data in the form of time history and spatial plots is provided through links to the POSTACO and GRAPE postprocessor codes.
NASA Astrophysics Data System (ADS)
Meulien Ohlmann, Odile
2013-02-01
Today the industry offers a chain of 3D products. Learning to "read" and to "create in 3D" becomes an issue of education of primary importance. 25 years professional experience in France, the United States and Germany, Odile Meulien set up a personal method of initiation to 3D creation that entails the spatial/temporal experience of the holographic visual. She will present some different tools and techniques used for this learning, their advantages and disadvantages, programs and issues of educational policies, constraints and expectations related to the development of new techniques for 3D imaging. Although the creation of display holograms is very much reduced compared to the creation of the 90ies, the holographic concept is spreading in all scientific, social, and artistic activities of our present time. She will also raise many questions: What means 3D? Is it communication? Is it perception? How the seeing and none seeing is interferes? What else has to be taken in consideration to communicate in 3D? How to handle the non visible relations of moving objects with subjects? Does this transform our model of exchange with others? What kind of interaction this has with our everyday life? Then come more practical questions: How to learn creating 3D visualization, to learn 3D grammar, 3D language, 3D thinking? What for? At what level? In which matter? for whom?
Kumar, Akhil; Roy, Sudeep; Tripathi, Shubhandra; Sharma, Ashok
2016-01-01
Beta-site APP cleaving enzyme1 (BACE1) catalyzes the rate determining step in the generation of Aβ peptide and is widely considered as a potential therapeutic drug target for Alzheimer's disease (AD). Active site of BACE1 contains catalytic aspartic (Asp) dyad and flap. Asp dyad cleaves the substrate amyloid precursor protein with the help of flap. Currently, there are no marketed drugs available against BACE1 and existing inhibitors are mostly pseudopeptide or synthetic derivatives. There is a need to search for a potent inhibitor with natural scaffold interacting with flap and Asp dyad. This study screens the natural database InterBioScreen, followed by three-dimensional (3D) QSAR pharmacophore modeling, mapping, in silico ADME/T predictions to find the potential BACE1 inhibitors. Further, molecular dynamics of selected inhibitors were performed to observe the dynamic structure of protein after ligand binding. All conformations and the residues of binding region were stable but the flap adopted a closed conformation after binding with the ligand. Bond oligosaccharide interacted with the flap as well as catalytic dyad via hydrogen bond throughout the simulation. This led to stabilize the flap in closed conformation and restricted the entry of substrate. Carbohydrates have been earlier used in the treatment of AD because of their low toxicity, high efficiency, good biocompatibility, and easy permeability through the blood-brain barrier. Our finding will be helpful in identify the potential leads to design novel BACE1 inhibitors for AD therapy. PMID:25707809
Developing natural convection in a fluid layer with localized heating and large viscosity variation
NASA Astrophysics Data System (ADS)
Hickox, C. E.; Chu, Tze Yao
Numerical simulations and laboratory experiments are used to elucidate aspects of transient natural convection in a magma chamber. The magma chamber is modeled as a horizontal fluid layer confined within an enclosure of square planform and heated from below by a strip heater centered on the lower boundary of the enclosure. The width of the strip heater and the depth of the fluid layer are one-fourth of the layer width. Corn syrup is used as the working fluid in order to approximate the large viscosity variation with temperature and the large Prandtl number typical of magma. The quiescent, uniform, fluid layer is subjected to instantaneous heating from the strip heater producing a transient flow which is dominated by two counter-rotating convective cells. Experimentally determined characteristics of the developing flow are compared with numerical simulations carried out with a finite element computer program. The results of numerical simulations are in essential agreement with experimental data. Differences between the numerical simulations and experimental measurements are conjectured to result from non-ideal effects present in the experiment which are difficult to represent accurately in a numerical simulation.
Nature's Grand Experiment: Linkage between magnetospheric convection and the radiation belts
NASA Astrophysics Data System (ADS)
Rodger, Craig J.; Cresswell-Moorcock, Kathy; Clilverd, Mark A.
2016-01-01
The solar minimum of 2007-2010 was unusually deep and long lived. In the later stages of this period the electron fluxes in the radiation belts dropped to extremely low levels. The flux of relativistic electrons (>1 MeV) was significantly diminished and at times was below instrument thresholds both for spacecraft located in geostationary orbits and also those in low-Earth orbit. This period has been described as a natural "Grand Experiment" allowing us to test our understanding of basic radiation belt physics and in particular the acceleration mechanisms which lead to enhancements in outer belt relativistic electron fluxes. Here we test the hypothesis that processes which initiate repetitive substorm onsets drive magnetospheric convection, which in turn triggers enhancement in whistler mode chorus that accelerates radiation belt electrons to relativistic energies. Conversely, individual substorms would not be associated with radiation belt acceleration. Contrasting observations from multiple satellites of energetic and relativistic electrons with substorm event lists, as well as chorus measurements, show that the data are consistent with the hypothesis. We show that repetitive substorms are associated with enhancements in the flux of energetic and relativistic electrons and enhanced whistler mode wave intensities. The enhancement in chorus wave power starts slightly before the repetitive substorm epoch onset. During the 2009/2010 period the only relativistic electron flux enhancements that occurred were preceded by repeated substorm onsets, consistent with enhanced magnetospheric convection as a trigger.
Developing natural convection in a fluid layer with localized heating and large viscosity variation
Hickox, C.E.; Chu, Tze Yao.
1991-01-01
Numerical simulations and laboratory experiments are used to elucidate aspects of transient natural convection in a magma chamber. The magma chamber is modeled as a horizontal fluid layer confined within an enclosure of square planform and heated from below by a strip heater centered on the lower boundary of the enclosure. The width of the strip heater and the depth of the fluid layer are one-fourth of the layer width. Corn syrup is used as the working fluid in order to approximate the large viscosity variation with temperature and the large Prandtl number typical of magma. The quiescent, uniform, fluid layer is subjected to instantaneous heating from the strip heater producing a transient flow which is dominated by two counter-rotating convective cells. Experimentally determined characteristics of the developing flow are compared with numerical simulations carried out with a finite element computer program. The results of numerical simulations are in essential agreement with experimental data. Differences between the numerical simulations and experimental measurements are conjectured to result from non-ideal effects present in the experiment which are difficult to represent accurately in a numerical simulation.
NASA Technical Reports Server (NTRS)
Diaguila, Anthony J; Freche, John C
1951-01-01
Blade-to-coolant heat-transfer data and operating data were obtained with a natural-convection water-cooled turbine over range of turbine speeds and inlet-gas temperatures. The convective coefficients were correlated by the general relation for natural-convection heat transfer. The turbine data were displaced from a theoretical equation for natural convection heat transfer in the turbulent region and from natural-convection data obtained with vertical cylinders and plates; possible disruption of natural convection circulation within the blade coolant passages was thus indicated. Comparison of non dimensional temperature-ratio parameters for the blade leading edge, midchord, and trailing edge indicated that the blade cooling effectiveness is greatest at the midchord and least at the trailing edge.
NASA Astrophysics Data System (ADS)
Magee, N. B.; Boaggio, K.; Bancroft, L.; Bandamede, M.
2015-12-01
Recent work has highlighted micro-scale roughness on the surfaces of ice crystals grown and imaged in-situ within the chambers of environmental scanning electron microscopes (ESEM). These observations appear to align with theoretical and satellite observations that suggest a prevalence of rough ice in cirrus clouds. However, the atmospheric application of the lab observations are indeterminate because the observations have been based only on crystals grown on substrates and in pure-water vapor environments. In this work, we present details and results from the development of a transfer technique which allows natural and lab-grown ice and snow crystals to be captured, preserved, and transferred into the ESEM for 3D imaging. Ice crystals were gathered from 1) natural snow, 2) a balloon-borne cirrus particle capture device, and 3) lab-grown ice crystals from a diffusion chamber. Ice crystals were captured in a pre-conditioned small-volume (~1 cm3) cryo-containment cell. The cell was then sealed closed and transferred to a specially-designed cryogenic dewer (filled with liquid nitrogen or crushed dry ice) for transport to a new Hitachi Field Emission, Variable Pressure SEM (SU-5000). The cryo-cell was then removed from the dewer and quickly placed onto the pre-conditioned cryo transfer stage attached to the ESEM (Quorum 3010T). Quantitative 3D topographical digital elevation models of ice surfaces are reported from SEM for the first time, including a variety of objective measures of statistical surface roughness. The surfaces of the transported crystals clearly exhibit signatures of mesoscopic roughening that are similar to examples of roughness seen in ESEM-grown crystals. For most transported crystals, the habits and crystal edges are more intricate that those observed for ice grown directly on substrates within the ESEM chamber. Portions of some crystals do appear smooth even at magnification greater than 1000x, a rare observation in our ESEM-grown crystals. The
NASA Astrophysics Data System (ADS)
Oztop, Hakan F.; Abu-Nada, Eiyad; Varol, Yasin; Al-Salem, Khaled
2011-04-01
In this study, the problem of steady state natural convection in an enclosure filled with a nanofluid has been analyzed numerically by using heating and cooling by sinusoidal temperature profiles on one side. The governing partial differential equations, in terms of the dimensionless stream function-vorticity and temperature, are solved numerically using the finite volume method for various inclination angles 0∘≤ϕ≤90∘, different types of nanoparticles (TiO 2 and Al 2O 3) and fractions of nanoparticles 0≤φ≤0.1, whereas the range of the Rayleigh number Ra is 10 3-10 5. It is found that the addition of nanoparticles into water affects the fluid flow and temperature distribution especially for higher Rayleigh numbers. An enhancement in heat transfer rate was registered for the whole range of Rayleigh numbers. However, low Rayleigh numbers show more enhancement compared to high Rayleigh numbers.
A numerical solution of variable porosity effects on natural convection in a packed-sphere cavity
David, E.; Lauriat, G. ); Cheng, P. )
1991-05-01
The problem of natural convection in differentially heated vertical cavities filled with spherical particles saturated with Newtonian fluids is investigated numerically. The Brinkman-Darcy-Ergun equation is used as the momentum equation, and the wall effect on porosity variation is approximated by an exponential function. The effect of variable stagnant thermal conductivities is taken into consideration in the energy equation. The formulation of the problem shows that the flow and heat transfer characteristics depend on six dimensionless parameters, namely, the Rayleigh and Prandtl numbers of the fluid phase, the dimensionless particle diameter, the conductivity ratio of the two phases, the bulk porosity, and the aspect ratio of the cavity. The influences of these parameters on the heat transfer rate are thoroughly investigated. The predicted Nusselt numbers are compared with existing experimental results. It is found that the computed Nusselt numbers based on the present model compare the best with experimental data.
Natural convection mass transfer on a vertical steel structure submerged in a molten aluminum pool
Cheung, F.B.; Yang, B.C.; Shiah, S.W.; Cho, D.H.; Tan, M.J.
1995-02-01
The process of dissolution mass transport along a vertical steel structure submerged in a large molten aluminum pool is studied theoretically. A mathematical model is developed from the conservation laws and thermodynamic principles, taking full account of the density variation in the dissolution boundary layer due to concentration differences. Also accounted for are the influence of the solubility of the wall material on species transfer and the motion of the solid/liquid interface at the dissolution front. The governing equations are solved by a combined analytical-numerical technique to determine the characteristics of the dissolution boundary layer and the rate of natural convection mass transfer. Based upon the numerical results, a correlation for the average Sherwood number is obtained. It is found that the Sherwood number depends strongly on the saturated concentration of the substrate at the moving dissolution front but is almost independent of the freestream velocity.
Design and Scaling of the Natural Convection Shutdown Heat Removal Test Facility
Lisowski, Darius D.; Gerardi, Craig D.; Bremer, Nathan C.; Farmer, Mitchell T.
2014-01-01
The Natural convection Shutdown heat removal Test Facility (NSTF) at Argonne National Laboratory (ANL) reflects a 1/2 scale model of one conceptual design for passive safety in advanced reactors. The project was initiated in 2010 primarily to conduct ex-vessel, passive decay heat removal experiments in support of the Advanced Reactor Concepts (ARC), Small Modular Reactor (SMR), and Next Generation Nuclear Plant (NGNP) programs while also generating data for code validation purposes. The facility successfully demonstrated scoping objectives in late 2013, and is expected to begin testing by early 2014. The following paper summarizes some of the key design and scaling considerations used in construction of the experimental facility, along with an overview of the current instrumentation and data acquisition methods. Details of the distributed fiber optic temperature system will be presented, which introduces a level of data density suitable for CFD validation and is a first-of-its-kind for largescale thermal hydraulics facilities.
Natural convection flow in porous enclosure with localized heating from below with heat flux
NASA Astrophysics Data System (ADS)
Siddiki, Md. Noor-A.-Alam; Molla, Md. Mamun; Saha, Suvash C.
2016-07-01
Unsteady natural convection flow in a two dimensional fluid saturated porous enclosure with localized heating from below with heat flux, symmetrical cooling from the sides and the insulated top wall has been investigated numerically. The governing equations are the Darcy's law for the porous media and the energy equation for the temperature field has been considered. The non-dimensional Darcy's law in terms of the stream function is solved by finite difference method using the successive over-relaxation (SOR) scheme and the energy equation is solved by Alternative Direction Alternative (ADI) scheme. The uniform heat flux source is located centrally at the bottom wall. The numerical results are presented in terms of the streamlines and isotherms, as well as the local and average rate of heat transfer for the wide range of the Darcy's Rayleigh number and the length of the heat flux source at the bottom wall.
MHD natural convection flow along a vertical wavy surface with heat generation and pressure work
NASA Astrophysics Data System (ADS)
Alim, M. A.; Kabir, K. H.; Andallah, L. S.
2016-07-01
In this paper, the influence of pressure work on MHD natural convection flow of viscous incompressible fluid along a uniformly heated vertical wavy surface with heat generation has been investigated. The governing boundary layer equations are first transformed into a non-dimensional form using suitable set of dimensionless variables. The resulting nonlinear system of partial differential equations are mapped into the domain of a vertical flat plate and then solved numerically employing the implicit finite difference method, known as Keller-box scheme. The numerical results for the velocity profiles, temperature profiles, skin friction coefficient, the rate of heat transfers, the streamlines and the isotherms are shown graphically and skin friction coefficient and rate of heat transfer have been shown in tabular form for different values of the selective set of parameters consisting of pressure work parameter Ge, the magnetic parameter M, Prandtl number Pr, heat generation parameter Q and the amplitude of the wavy surface.
NASA Astrophysics Data System (ADS)
Liviu, Pascu; Adriana, Putan; Vasile, Putan; Alina, Lascutoni
2012-09-01
The similarity between steel ladles and hot water model regarding natural convection phenomena has been analyzed through examination of the numerical solutions of turbulent Navier-Stokes partial differential equations governing the phenomena in question. Key similarity criteria for non-isothermal physical modeling of steel ladles with hot-water models have been derived as Frm = Frp and (β∇T)m = (β∇T)p where the subscript m and p stand for the water model and the prototype steel ladle, respectively. Accordingly, appropriate conditions fulfilling the above criteria, such as model size, water temperature, time scale factor and the scale factor of boundary heat loss fluxes, have been proposed and discussed.
Weaver, J.A.; Viskanta, R. )
1992-01-01
An investigation of natural convection is presented to examine the influence of a horizontal temperature gradient and a concentration gradient occurring from the bottom to the cold wall in a cavity. As the solutal buoyancy force changes from augmenting to opposing the thermal buoyancy force, the fluid motion switches from unicellular to multicellular flow (fluid motion is up the cold wall and down the hot wall for the bottom counterrotating flow cell). Qualitatively, the agreement between predicted streamlines and smoke flow patterns is generally good. In contrast, agreement between measured and predicted temperature and concentration distributions ranges from fair to poor. Part of the discrepancy can be attributed to experimental error. However, there remains considerable discrepancy between data and predictions due to the idealizations of the mathematical model, which examines only first-order physical effects. An unsteady flow, variable thermophysical properties, conjugate effects, species interdiffusion, and radiation were not accounted for in the model. 31 refs.
Dunn, T.A.; McCallen, R.C.
2000-10-17
The Galerkin Finite Element Method was used to predict a natural convection flow in an enclosed cavity. The problem considered was a differentially heated, tall (8:1), rectangular cavity with a Rayleigh number of 3.4 x 10{sup 5} and Prandtl number of 0.71. The incompressible Navier-Stokes equations were solved using a Boussinesq approximation for the buoyancy force. The algorithm was developed for efficient use on massively parallel computer systems. Emphasis was on time-accurate simulations. It was found that the average temperature and velocity values can be captured with a relatively coarse grid, while the oscillation amplitude and period appear to be grid sensitive and require a refined computation.
Domanus, H.M.; Sha, W.T.
1981-01-01
The single-phase COMMIX (COMponent MIXing) computer code performs fully three-dimensional, transient, thermal-hydraulic analyses of liquid-sodium LMFBR components. It solves the conservation equations of mass, momentum, and energy as a boundary-value problem in space and as an initial-value problem in time. The concepts of volume porosity, surface permeability and distributed resistance, and heat source have been employed in quasi-continuum (rod-bundle) applications. Results from three transient simulations involving forced and natural convection are presented: (1) a sodium-filled horizontal pipe initially of uniform temperature undergoing an inlet velocity rundown transient, as well as an inlet temperature transient; (2) a 19-pin LMFBR rod bundle undergoing a velocity transient; and, (3) a simulation of a water test of a 1/10-scale outlet plenum undergoing both velocity and temperature transients.
NASA Astrophysics Data System (ADS)
Barakos, G.; Mitsoulis, E.; Assimacopoulos, D.
1994-04-01
Numerical simulations have been undertaken for the benchmark problem of natural convection flow in a square cavity. The control volume method is used to solve the conservation equations for laminar and turbulent flows for a series of Rayleigh numbers (Ra) reaching values up to 10(exp 10). The k-epsilon model has been used for turbulence modelling with and without logarithmic wall functions. Uniform and non-uniform (stretched) grids have been employed with increasing density to guarantee accurate solutions, especially near the walls for high Ra-values. ADI and SIP solvers are implemented to accelerate convergence. Excellent agreement is obtained with previous numerical solutions, while some discrepancies with others for high Ra-values may be due to a possibly different implementation of the wall functions. Comparisons with experimental data for heat transfer (Nusselt number) clearly demonstrates the limitations of the standard k-epsilon model with logarithmic wall functions, which gives significant overpredictions.
NASA Astrophysics Data System (ADS)
Slezak, Ondrej; Yasuhara, Ryo; Lucianetti, Antonio; Vojna, David; Mocek, Tomas
2015-06-01
Thermal birefringence-induced depolarization in terbium gallium garnet (TGG) ceramic rods has been numerically evaluated for the geometry and heating conditions in a previous experiment. In this model, the spatially resolved heat transfer coefficient corresponding to natural convection cooling and the offset of the beam from the rotational axis of the rod have been incorporated and the realistic beam profile used in the experiment has been considered. A resulting beam depolarization ratio of 4.3 × 10-4 has been calculated for an input power of 117 W. The results were found to be in good agreement with the measured values. Furthermore, a parametric study of the depolarization ratio for higher input powers has been performed leading to a depolarization ratio of 3.3 × 10-2 for 1 kW input power.
Natural convection heat transfer on two horizontal cylinders in liquid sodium
Hata, K.; Shiotsu, M.; Takeuchi, Y.
1995-09-01
Natural convection heat transfer on two horizontal 7.6 mm diameter test cylinders assembled with the ratio of the distance between each cylinder axis to the cylinder diameter, S/D, of 2 in liquid sodium was studied experimentally and theoretically. The heat transfer coefficients on the cylinder surface due to the same heat inputs ranging from 1.0 X 10{sup 7} to 1.0 x 10{sup 9} W/m{sup 3} were obtained experimentally for various setting angeles, {gamma}, between vertical direction and the plane including both of these cylinder axis over the range of zero to 90{degrees}. Theoretical equations for laminar natural convection heat transfer from the two horizontal cylinders were numerically solved for the same conditions as the experimental ones considering the temperature dependence of thermophysical properties concerned. The average Nusselt numbers, Nu, values on the Nu versus modified Rayleigh number, R{sub f}, graph. The experimental values of Nu for the upper cylinder are about 20% lower than those for the lower cylinder at {gamma} = 0{degrees} for the range of R{sub f} tested here. The value of Nu for the upper cylinder becomes higher and approaches that for the lower cylinder with the increase in {gamma} over range of 0 to 90{degrees}. The values of Nu for the lower cylinder at each {gamma} are almost in agreement with those for a single cylinder. The theoretical values of Nu on two cylinders except those for R{sub f}<4 at {gamma} = 0{degrees} are in agreement with the experimental data at each {gamma} with the deviations less than 15%. Correlations for Nu on the upper and lower cylinders were obtained as functions of S/D and {gamma} based n the theoretical solutions for the S/D ranged over 1.5 to 4.0.
NASA Astrophysics Data System (ADS)
Carr, Simon; Spencer, Kate; James, Tempest; Lucy, Diggens
2015-04-01
Saltmarshes are globally important environments which, though occupying < 4% of the Earth's surface, provide a range of ecosystem services. Yet, they are threatened by sea level rise, human population growth, urbanization and pollution resulting in degradation. To compensate for this habitat loss many coastal restoration projects have been implemented over the last few decades, largely driven by legislative requirements for improved biodiversity e.g. the EU Habitats Directive and Birds Directive. However, there is growing evidence that restored saltmarshes, recreated through the return to tidal inundation of previously drained and defended low-lying coastal land, do not have the same species composition even after 100 years and while environmental enhancement has been achieved, there may be consequences for ecosystem functioning This study presents the findings of a comparative analysis of detailed sediment structure and hydrological functioning of equivalent natural and de-embanked saltmarsh sediments at Orplands Farm, Essex, UK. 3D x-ray CT scanning of triplicate undisturbed sediment cores recovered in 2013 have been used to derive detailed volumetric reconstructions of macropore structure and networks, and to infer differences in bulk microporosity between natural and de-embanked saltmarshes. These volumes have been further visualised for qualitative analysis of the main sediment components, and extraction of key macropore space parameters for quantified analysis including total porosity and connectivity, as well as structure, organisation and efficiency (tortuosity) of macropore networks. Although total porosity was significantly greater within the de-embanked saltmarsh sediments, pore networks in these samples were less organised and more tortuous, and were also inferred to have significantly lower micro-porosity than those of the natural saltmarsh. These datasets are applied to explain significant differences in the hydraulic behaviour and functioning
ERIC Educational Resources Information Center
Hastings, S. K.
2002-01-01
Discusses 3 D imaging as it relates to digital representations in virtual library collections. Highlights include X-ray computed tomography (X-ray CT); the National Science Foundation (NSF) Digital Library Initiatives; output peripherals; image retrieval systems, including metadata; and applications of 3 D imaging for libraries and museums. (LRW)
ERIC Educational Resources Information Center
Saatadjian, Esteban; Lesage, Francois; Mota, Jose Paulo B.
2013-01-01
A project that involves the numerical simulation of transport phenomena is an excellent method to teach this subject to senior/graduate chemical engineering students. The subject presented here has been used in our senior/graduate course, it concerns the study of natural convection heat transfer between two concentric, horizontal, saturated porous…
Qi, Cong; He, Yurong; Yan, Shengnan; Tian, Fenglin; Hu, Yanwei
2013-01-01
Considering interaction forces (gravity and buoyancy force, drag force, interaction potential force, and Brownian force) between nanoparticles and a base fluid, a two-phase Lattice Boltzmann model for natural convection of nanofluid is developed in this work. It is applied to investigate the natural convection in a square enclosure (the left wall is kept at a high constant temperature (TH), and the top wall is kept at a low constant temperature (TC)) filled with Al2O3/H2O nanofluid. This model is validated by comparing numerical results with published results, and a satisfactory agreement is shown between them. The effects of different nanoparticle fractions and Rayleigh numbers on natural convection heat transfer of nanofluid are investigated. It is found that the average Nusselt number of the enclosure increases with increasing nanoparticle volume fraction and increases more rapidly at a high Rayleigh number. Also, the effects of forces on nanoparticle volume fraction distribution in the square enclosure are studied in this paper. It is found that the driving force of the temperature difference has the biggest effect on nanoparticle volume fraction distribution. In addition, the effects of interaction forces on flow and heat transfer are investigated. It is found that Brownian force, interaction potential force, and gravity-buoyancy force have positive effects on the enhancement of natural convective heat transfer, while drag force has a negative effect. PMID:23374509
NASA Technical Reports Server (NTRS)
Parker, E. N.
1974-01-01
It had been pointed out by Parker (1974) that the basic cause of the sunspot phenomenon is the enhanced heat transport in the magnetic field of the sunspot. The enhanced transport occurs through convective overstability which operates as a heat engine generating Alfven waves. The characteristics of the convective forces present are investigated along with questions concerning overstability and convectively driven Alfven waves. Relations regarding instability and convectively driven surface waves are discussed and attention is given to individual overstable Alfven modes. It is found that the form of an Alfven wave in the absence of convective forces is entirely arbitrary, so that waves with any arbitrary profile can be fitted into a vertical column of the field without disturbing the fluid outside. With the introduction of convective forces the situation changes so that the presence of lateral boundaries alters the form of the basic wave modes.