An approximate solution method for boundary layer flow of a power law fluid over a flat plate

The work in this paper deals with the development of momentum and thermal boundary layers when a power law fluid flows over a flat plate. At the plate we impose either constant temperature, constant flux or a Newton cooling condition. The problem is analysed using similarity solutions, integral momentum and energy equations and an approximation technique which is a form of the Heat Balance Integral Method. The fluid properties are assumed to be independent of temperature, hence the momentum equation uncouples from the thermal problem. We first derive the similarity equations for the velocity and present exact solutions for the case where the power law index n = 2. The similarity solutions are used to validate the new approximation method. This new technique is then applied to the thermal boundary layer, where a similarity solution can only be obtained for the case n = 1.

Bifurcations and chaos in single-roll natural convection with low Prandtl number

Convective flows of a small Prandtl number fluid contained in a two-dimensional cavity subject to a lateral thermal gradient are numerically studied by using different techniques. The aspect ratio (length to height) is kept at around 2. This value is found optimal to make the flow most unstable while keeping the basic single-roll structure. Two cases of thermal boundary conditions on the horizontal plates are considered: perfectly conducting and adiabatic. For increasing Rayleigh numbers we find a transition from steady flow to periodic oscillations through a supercritical Hopf bifurcation that maintains the centrosymmetry of the basic circulation. For a Rayleigh number of about ten times that of the Hopf bifurcation the system initiates a complex scenario of bifurcations. In the conductive case these include a quasiperiodic route to chaos. In the adiabatic one the dynamics is dominated by the interaction of two Neimark-Sacker bifurcations of the basic periodic solutions, leading to the stable coexistence of three incommensurate frequencies, and finally to chaos. In all cases, the complex time-dependent behavior does not break the basic, single-roll structure.

Improving the accuracy of heat balance integral methods applied to thermal problems with time dependent boundary conditions

In this paper the two main drawbacks of the heat balance integral methods are examined. Firstly we investigate the choice of approximating function. For a standard polynomial form it is shown that combining the Heat Balance and Refined Integral methods to determine the power of the highest order term will either lead to the same, or more often, greatly improved accuracy on standard methods. Secondly we examine thermal problems with a time-dependent boundary condition. In doing so we develop a logarithmic approximating function. This new function allows us to model moving peaks in the temperature profile, a feature that previous heat balance methods cannot capture. If the boundary temperature varies so that at some time t & 0 it equals the far-field temperature, then standard methods predict that the temperature is everywhere at this constant value. The new method predicts the correct behaviour. It is also shown that this function provides even more accurate results, when coupled with the new CIM, than the polynomial profile. Analysis primarily focuses on a specified constant boundary temperature and is then extended to constant flux, Newton cooling and time dependent boundary conditions.

Heat transfer between nanoparticles: Thermal conductance for near-field interactions

We analyze the heat transfer between two nanoparticles separated by a distance lying in the near-field domain in which energy interchange is due to the Coulomb interactions. The thermal conductance is computed by assuming that the particles have charge distributions characterized by fluctuating multipole moments in equilibrium with heat baths at two different temperatures. This quantity follows from the fluctuation-dissipation theorem for the fluctuations of the multipolar moments. We compare the behavior of the conductance as a function of the distance between the particles with the result obtained by means of molecular dynamics simulations. The formalism proposed enables us to provide a comprehensive explanation of the marked growth of the conductance when decreasing the distance between the nanoparticles.

Generalized backward doubly stochastic differential equations and SPDEs with nonlinear Neumann boundary conditions

In this paper, a new class of generalized backward doubly stochastic differential equations is investigated. This class involves an integral with respect to an adapted continuous increasing process. A probabilistic representation for viscosity solutions of semi-linear stochastic partial differential equations with a Neumann boundary condition is given.

Thermal tests on the LISA pathfinder engineering model

Report for the scientific sojourn carried out at Albert Einstein Institut in Germany, from April to July 2006.

Non-hyperbolic boundary equilibrium bifurcations in planar Filippov systems: a case study approach

Boundary equilibrium bifurcations in piecewise smooth discontinuous systems are characterized by the collision of an equilibrium point with the discontinuity surface. Generically, these bifurcations are of codimension one, but there are scenarios where the phenomenon can be of higher codimension. Here, the possible collision of a non-hyperbolic equilibrium with the boundary in a two-parameter framework and the nonlinear phenomena associated with such collision are considered. By dealing with planar discontinuous (Filippov) systems, some of such phenomena are pointed out through specific representative cases. A methodology for obtaining the corresponding bi-parametric bifurcation sets is developed.

Numerical simulation of combined heat transfer phenomena (conduction, convection, and radiation). Application to he optimisation of thermal systems

Thermal systems interchanging heat and mass by conduction, convection, radiation (solar and thermal ) occur in many engineering applications like energy storage by solar collectors, window glazing in buildings, refrigeration of plastic moulds, air handling units etc. Often these thermal systems are composed of various elements for example a building with wall, windows, rooms, etc. It would be of particular interest to have a modular thermal system which is formed by connecting different modules for the elements, flexibility to use and change models for individual elements, add or remove elements without changing the entire code. A numerical approach to handle the heat transfer and fluid flow in such systems helps in saving the full scale experiment time, cost and also aids optimisation of parameters of the system. In subsequent sections are presented a short summary of the work done until now on the orientation of the thesis in the field of numerical methods for heat transfer and fluid flow applications, the work in process and the future work.

What is the Jacobian of a Riemann surface with boundary?

We define the Jacobian of a Riemann surface with analytically parametrized boundary components. These Jacobians belong to a moduli space of "open abelian varieties" which satisfies gluing axioms similar to those of Riemann surfaces, and therefore allows a notion of "conformal field theory" to be defined on this space. We further prove that chiral conformal field theories corresponding to even lattices factor through this moduli space of open abelian varieties.

Near-linear dynamics in KdV with periodic boundary conditions

Near linear evolution in Korteweg de Vries (KdV) equation with periodic boundary conditions is established under the assumption of high frequency initial data. This result is obtained by the method of normal form reduction.

Optimal exponent heat balance and refined integral methods applied to Stefan problems

When using a polynomial approximating function the most contentious aspect of the Heat Balance Integral Method is the choice of power of the highest order term. In this paper we employ a method recently developed for thermal problems, where the exponent is determined during the solution process, to analyse Stefan problems. This is achieved by minimising an error function. The solution requires no knowledge of an exact solution and generally produces significantly better results than all previous HBI models. The method is illustrated by first applying it to standard thermal problems. A Stefan problem with an analytical solution is then discussed and results compared to the approximate solution. An ablation problem is also analysed and results compared against a numerical solution. In both examples the agreement is excellent. A Stefan problem where the boundary temperature increases exponentially is analysed. This highlights the difficulties that can be encountered with a time dependent boundary condition. Finally, melting with a time-dependent flux is briefly analysed without applying analytical or numerical results to assess the accuracy.

A one-dimensional Keller-Segel equation with a drift issued from the boundary

We investigate in this note the dynamics of a one-dimensional Keller-Segel type model on the half-line. On the contrary to the classical configuration, the chemical production term is located on the boundary. We prove, under suitable assumptions, the following dichotomy which is reminiscent of the two-dimensional Keller-Segel system. Solutions are global if the mass is below the critical mass, they blow-up in finite time above the critical mass, and they converge to some equilibrium at the critical mass. Entropy techniques are presented which aim at providing quantitative convergence results for the subcritical case. This note is completed with a brief introduction to a more realistic model (still one-dimensional).

Mesophyll conductance to CO2 and leaf morphological characteristics under drought stress during Quercus ilex L resprouting

• Quercus ilex L., the dominant species in Mediterranean forests and one with a great capacity for resprouting after disturbances, is threatened by the expected increase in fire frequency and drought associated with climate change. • The aim of this study was to determine the contribution of photosynthesis limitants, especially mesophyll conductance (gmes ) during this species’ resprouting and under summer drought. • Resprouts showed 5.3-fold increased gmes and 3.8-fold increased stomatal conductance (gs) atmidday with respect to leaves of undisturbed individuals. With increased drought, structural changes (decreased density and increased thickness) in resprouts contributed to the observed higher photosynthesis and increased gmes. However, gmes only partially depended on leaf structure, and was also under physiological control. Resprouts also showed lower non-stomatal limitations (around 50% higher carboxylation velocity (Vc,max) and capacity for ribulose-1,5-bisphosphate regeneration (Jmax)). A significant contribution of gmes to leaf carbon isotope discrimination values was observed. • gmes exhibits a dominant role in photosynthesis limitation in Q. ilex and is regulated by factors other than morphology. During resprouting after disturbances, greater capacity to withstand drought, as evidenced by higher gmes , gs and lower non-stomatal limitants, enables increased photosynthesis and rapid growth.

Engineering the thermal emission of macroporous silicon

Todos los cuerpos emiten luz espontaneamente al ser calentados. El espectro de radiacion es una funcion de la temperatura y el material. Sin embargo, la mayoria de los materiales irradia, en general, en una banda espectral amplia. Algunas matereiales, por el contrario, son capaces de concentrar la radiacion termica en una banda espectral mucho mas estrecha. Estos materiales se conocen como emisores selectivos y su uso tiene un profundo impacto en la eficiencia de sistemas sistemas tales como iluminacion y conversion de energia termofotovoltaica. De los emisores selectivos se espera que sean capaces de operar a altas temperaturas y que emitan en una banda espectral muy concisa. Uno de los metodos mas prometedores para controlar y disenar el espectro de emision termico es la utilizacion de cristales fotonicos. Los cristales fotonicos son estructuras periodicas artificiales capaces de controlar y confinar la luz de formas sin precedentes. Sin embargo, la produccion de dichas estructuras con grandes superficies y capaces de soportar altas temperaturas sigue siendo una dificil tarea. Este trabajo esta dedicada al estudio de las propiedades de emision termica de estructuras 3D de silicio macroporoso en el rango espectral mid-IR (2-30 m). En particular, este trabajo se enfoca en reducir la elevada emisividad del silicio cristalino. Las muestras estudiadas en este trabajo tienen una periodicidad de 4 m, lo que limitan los resultados obtenidos a la banda del infrarrojo medio, aunque estructuras mucho mas pequenas son tecnologicamente realizables con el metodo de fabricacion utilizado. Hemos demostrado que el silicio macroporoso 3D puede inhibir completamente la emision termica en su superficie. Mas aun, esta banda se puede ajustar en un amplio margen mediante pequenos cambios durante la formacion de los macroporos. Tambien hemos demostrado que tanto el ancho como la frecuencia de la banda de inhibicion se puede doblar mediante la aplicacion de tecnicas de postprocesado adecuadas. Finalmente hemos mostrado que es posible crear bandas de baja emisividad arbitrariamente anchas mediante estructuras macroporosas aperiodicas.

Seasonal changes of whole root system conductance by a drought-tolerant grape root system

The role of root systems in drought tolerance is a subject of very limited information compared with above-ground responses. Adjustments to the ability of roots to supply water relative to shoot transpiration demand is proposed as a major means for woody perennial plants to tolerate drought, and is often expressed as changes in the ratios of leaf to root area (AL:AR). Seasonal root proliferation in a directed manner could increase the water supply function of roots independent of total root area (AR) and represents a mechanism whereby water supply to demand could be increased. To address this issue, seasonal root proliferation, stomatal conductance (gs) and whole root system hydraulic conductance (kr) were investigated for a drought-tolerant grape root system (Vitis berlandieri×V. rupestris cv. 1103P) and a non-drought-tolerant root system (Vitis riparia×V. rupestris cv. 101-14Mgt), upon which had been grafted the same drought-sensitive clone of Vitis vinifera cv. Merlot. Leaf water potentials (ψL) for Merlot grafted onto the 1103P root system (–0.91±0.02 MPa) were +0.15 MPa higher than Merlot on 101-14Mgt (–1.06±0.03 MPa) during spring, but dropped by approximately –0.4 MPa from spring to autumn, and were significantly lower by –0.15 MPa (–1.43±0.02 MPa) than for Merlot on 101-14Mgt (at –1.28±0.02 MPa). Surprisingly, gs of Merlot on the drought-tolerant root system (1103P) was less down-regulated and canopies maintained evaporative fluxes ranging from 35–20 mmol vine−1 s−1 during the diurnal peak from spring to autumn, respectively, three times greater than those measured for Merlot on the drought-sensitive rootstock 101-14Mgt. The drought-tolerant root system grew more roots at depth during the warm summer dry period, and the whole root system conductance (kr) increased from 0.004 to 0.009 kg MPa−1 s−1 during that same time period. The changes in kr could not be explained by xylem anatomy or conductivity changes of individual root segments. Thus, the manner in which drought tolerance was conveyed to the drought-sensitive clone appeared to arise from deep root proliferation during the hottest and driest part of the season, rather than through changes in xylem structure, xylem density or stomatal regulation. This information can be useful to growers on a site-specific basis in selecting rootstocks for grape clonal material (scions) grafted to them.