Rayleigh-Benard Convection with a Radial Ramp in Plate Separation

Pattern formation in Rayleigh-Benard convection in a large-aspect-ratio cylinder with a radial ramp in the plate separation is studied analytically and numerically by performing numerical simulations of the Boussinesq equations. A horizontal mean flow and a vertical large scale counterflow are quantified and used to understand the pattern wavenumber. Our results suggest that the mean flow, generated by amplitude gradients, plays an important role in the roll compression observed as the control parameter is increased. Near threshold the mean flow has a quadrupole dependence with a single vortex in each quadrant while away from threshold the mean flow exhibits an octupole dependence with a counter-rotating pair of vortices in each quadrant. This is confirmed analytically using the amplitude equation and Cross-Newell mean flow equation. By performing numerical experiments the large scale counterflow is also found to aid in the roll compression away from threshold but to a much lesser degree. Our results yield an understanding of the pattern wavenumbers observed in experiment away from threshold and suggest that near threshold the mean flow and large scale counterflow are not responsible for the observed shift to smaller than critical wavenumbers.Comment: 10 pages, 13 figure

Navigation

n/

Finite Size Scaling of Domain Chaos

Numerical studies of the domain chaos state in a model of rotating Rayleigh-Benard convection suggest that finite size effects may account for the discrepancy between experimentally measured values of the correlation length and the predicted divergence near onset

Extensive chaos in Rayleigh-Bénard convection

Using large-scale numerical calculations we explore spatiotemporal chaos in Rayleigh-Bénard convection for experimentally relevant conditions. We calculate the spectrum of Lyapunov exponents and the Lyapunov dimension describing the chaotic dynamics of the convective fluid layer at constant thermal driving over a range of finite system sizes. Our results reveal that the dynamics of fluid convection is truly chaotic for experimental conditions as illustrated by a positive leading-order Lyapunov exponent. We also find the chaos to be extensive over the range of finite-sized systems investigated as indicated by a linear scaling between the Lyapunov dimension of the chaotic attractor and the system size

Traveling waves in rotating Rayleigh-Bénard convection: Analysis of modes and mean flow

Numerical simulations of the Boussinesq equations with rotation for realistic no-slip boundary conditions and a finite annular domain are presented. These simulations reproduce traveling waves observed experimentally. Traveling waves are studied near threshhold by using the complex Ginzburg-Landau equation (CGLE): a mode analysis enables the CGLE coefficients to be determined. The CGLE coefficients are compared with previous experimental and theoretical results. Mean flows are also computed and found to be more significant as the Prandtl number decreases (from sigma=6.4 to sigma=1). In addition, the mean flow around the outer radius of the annulus appears to be correlated with the mean flow around the inner radius

Rayleigh-Benard Convection in Large-Aspect-Ratio Domains

The coarsening and wavenumber selection of striped states growing from random initial conditions are studied in a non-relaxational, spatially extended, and far-from-equilibrium system by performing large-scale numerical simulations of Rayleigh-B\'{e}nard convection in a large-aspect-ratio cylindrical domain with experimentally realistic boundaries. We find evidence that various measures of the coarsening dynamics scale in time with different power-law exponents, indicating that multiple length scales are required in describing the time dependent pattern evolution. The translational correlation length scales with time as $t^{0.12}$, the orientational correlation length scales as $t^{0.54}$, and the density of defects scale as $t^{-0.45}$. The final pattern evolves toward the wavenumber where isolated dislocations become motionless, suggesting a possible wavenumber selection mechanism for large-aspect-ratio convection.Comment: 5 pages, 6 figure

Power-Law Behavior of Power Spectra in Low Prandtl Number Rayleigh-Benard Convection

The origin of the power-law decay measured in the power spectra of low Prandtl number Rayleigh-Benard convection near the onset of chaos is addressed using long time numerical simulations of the three-dimensional Boussinesq equations in cylindrical domains. The power-law is found to arise from quasi-discontinuous changes in the slope of the time series of the heat transport associated with the nucleation of dislocation pairs and roll pinch-off events. For larger frequencies, the power spectra decay exponentially as expected for time continuous deterministic dynamics.Comment: (10 pages, 6 figures

Biplane Fluoroscopy for Hindfoot Motion Analysis during Gait: A Model-based Evaluation

The purpose of this study was to quantify the accuracy and precision of a biplane fluoroscopy system for model-based tracking of in vivo hindfoot motion during over-ground gait. Gait was simulated by manually manipulating a cadaver foot specimen through a biplane fluoroscopy system attached to a walkway. Three 1.6-mm diameter steel beads were implanted into the specimen to provide marker-based tracking measurements for comparison to model-based tracking. A CT scan was acquired to define a gold standard of implanted bead positions and to create 3D models for model-based tracking. Static and dynamic trials manipulating the specimen through the capture volume were performed. Marker-based tracking error was calculated relative to the gold standard implanted bead positions. The bias, precision, and root-mean-squared (RMS) error of model-based tracking was calculated relative to the marker-based measurements. The overall RMS error of the model-based tracking method averaged 0.43 ± 0.22 mm and 0.66 ± 0.43° for static and 0.59 ± 0.10 mm and 0.71 ± 0.12° for dynamic trials. The model-based tracking approach represents a non-invasive technique for accurately measuring dynamic hindfoot joint motion during in vivo, weight bearing conditions. The model-based tracking method is recommended for application on the basis of the study results

Enhanced tracer transport by the spiral defect chaos state of a convecting fluid

To understand how spatiotemporal chaos may modify material transport, we use direct numerical simulations of the three-dimensional Boussinesq equations and of an advection-diffusion equation to study the transport of a passive tracer by the spiral defect chaos state of a convecting fluid. The simulations show that the transport is diffusive and is enhanced by the spatiotemporal chaos. The enhancement in tracer diffusivity follows two regimes. For large Peclet numbers (that is, small molecular diffusivities of the tracer), we find that the enhancement is proportional to the Peclet number. For small Peclet numbers, the enhancement is proportional to the square root of the Peclet number. We explain the presence of these two regimes in terms of how the local transport depends on the local wave numbers of the convection rolls. For large Peclet numbers, we further find that defects cause the tracer diffusivity to be enhanced locally in the direction orthogonal to the local wave vector but suppressed in the direction of the local wave vector.Comment: 11 pages, 12 figure

Multilevel mesh partitioning for optimising domain shape

Multilevel algorithms are a successful class of optimisation techniques which address the mesh partitioning problem for mapping meshes onto parallel computers. They usually combine a graph contraction algorithm together with a local optimisation method which refines the partition at each graph level. To date these algorithms have been used almost exclusively to minimise the cut-edge weight in the graph with the aim of minimising the parallel communication overhead. However it has been shown that for certain classes of problem, the convergence of the underlying solution algorithm is strongly influenced by the shape or aspect ratio of the subdomains. In this paper therefore, we modify the multilevel algorithms in order to optimise a cost function based on aspect ratio. Several variants of the algorithms are tested and shown to provide excellent results