On the natural convection boundary layer adjacent to an inclined flat plate subject to ramp heating

An investigation of the natural convection boundary layer adjacent to an inclined semi-infinite plate subject to a temperature boundary condition which follows a ramp function up until some specified time and then remains constant is reported. The development of the flow from start-up to a steadystate has been described based on scaling analyses and verified by numerical simulations. Attention in this study has been given to fluids having a Prandtl number Pr less than unity. The boundary layer flow depends on the comparison of the time at which the ramp heating is completed and the time at which the boundary layer completes its growth. If the ramp time is long compared with the steady state time, the layer reaches a quasi steady mode in which the growth of the layer is governed solely by the thermal balance between convection and conduction. On the other hand, if the ramp is completed before the layer becomes steady; the subsequent growth is governed by the balance between buoyancy and inertia, as for the case of instantaneous heating.

Scaling analysis of natural convection boundary layers : Application to attic space

The natural convection thermal boundary layer adjacent to an inclined flat plate and inclined walls of an attic space subject to instantaneous and ramp heating and cooling is investigated. A scaling analysis has been performed to describe the flow behaviour and heat transfer. Major scales quantifying the flow velocity, flow development time, heat transfer and the thermal and viscous boundary layer thicknesses at different stages of the flow development are established. Scaling relations of heating-up and cooling-down times and heat transfer rates have also been reported for the case of attic space. The scaling relations have been verified by numerical simulations over a wide range of parameters. Further, a periodic temperature boundary condition is also considered to show the flow features in the attic space over diurnal cycles.

An analytical framework for quantifying aquifer response time scales associated with transient boundary conditions

A major challenge in studying coupled groundwater and surface-water interactions arises from the considerable difference in the response time scales of groundwater and surface-water systems affected by external forcings. Although coupled models representing the interaction of groundwater and surface-water systems have been studied for over a century, most have focused on groundwater quantity or quality issues rather than response time. In this study, we present an analytical framework, based on the concept of mean action time (MAT), to estimate the time scale required for groundwater systems to respond to changes in surface-water conditions. MAT can be used to estimate the transient response time scale by analyzing the governing mathematical model. This framework does not require any form of transient solution (either numerical or analytical) to the governing equation, yet it provides a closed form mathematical relationship for the response time as a function of the aquifer geometry, boundary conditions, and flow parameters. Our analysis indicates that aquifer systems have three fundamental time scales: (i) a time scale that depends on the intrinsic properties of the aquifer; (ii) a time scale that depends on the intrinsic properties of the boundary condition, and; (iii) a time scale that depends on the properties of the entire system. We discuss two practical scenarios where MAT estimates provide useful insights and we test the MAT predictions using new laboratory-scale experimental data sets.

The effect of closed boundary conditions on a stationary dynamo

One of two boundary conditions generally assumed in solutions of the dynamo equation is related to the disappearance of the azimuthal field at the boundary. Parker (1984) points out that for the realization of this condition the field must escape freely through the surface. Escape requires that the field be detached from the gas in which it is embedded. In the case of the sun, this can be accomplished only through reconnection in the tenuous gas above the visible surface. Parker concludes that the observed magnetic activity on the solar surface permits at most three percent of the emerging flux to escape. He arrives at the conclusion that, instead of B(phi) = 0, the partial derivative of B(phi) to r is equal to zero. The present investigation is concerned with the effect of changing the boundary condition according to Parker's conclusion. Implications for the solar convection zone are discussed.

On compressible boundary layer on a flat plate with uniform suction

A quartic profile in terms of the normal distance from the wall has been taken and coefficients are evaluated by satisfying one more boundary condition on the wall than the usual one. By doing so, the limitations about the Reynolds number of the quartic profile adopted by Lew (1949) has been removed. The Kármán (1921) Momentum Integral Equation has been used to evaluate the various characteristics of the flow. A comparative study of Lew's quartic profile and exponential profile together with the quartic profile of the present paper has been undertaken and the graphs for the various characteristics of the flow for a number of Mach numbers and suction coefficients have been drawn. At the end, certain conclusions of general nature about the velocity profiles have been recorded.

Exact analysis of Burgers equation on semiline with flux condition at the origin

The initial boundary value problem for the Burgers equation in the domain x greater-or-equal, slanted 0, t > 0 with flux boundary condition at x = 0 has been solved exactly. The behaviour of the solution as t tends to infinity is studied and the “asymptotic profile at infinity” is obtained. In addition, the uniqueness of the solution of the initial boundary value problem is proved and its inviscid limit as var epsilon → 0 is obtained.

Continuum theory of edge states of topological insulators: variational principle and boundary conditions

We develop a continuum theory to model low energy excitations of a generic four-band time reversal invariant electronic system with boundaries. We propose a variational energy functional for the wavefunctions which allows us to derive natural boundary conditions valid for such systems. Our formulation is particularly suited for developing a continuum theory of the protected edge/surface excitations of topological insulators both in two and three dimensions. By a detailed comparison of our analytical formulation with tight binding calculations of ribbons of topological insulators modelled by the Bernevig-Hughes-Zhang (BHZ) Hamiltonian, we show that the continuum theory with a natural boundary condition provides an appropriate description of the low energy physics.

Boundary perturbations and the Helmholtz equation in three dimensions

We propose an analytic perturbative scheme in the spirit of Lord Rayleigh's work for determining the eigenvalues of the Helmholtz equation in three dimensions inside an arbitrary boundary where the eigenfunction satisfies either the Dirichlet boundary condition or the Neumann boundary condition. Although numerous works are available in the literature for arbitrary boundaries in two dimensions, to the best of our knowledge the formulation in three dimensions is proposed for the first time. In this novel prescription, we have expanded the arbitrary boundary in terms of spherical harmonics about an equivalent sphere and obtained perturbative closed-form solutions at each order for the problem in terms of corrections to the equivalent spherical boundary for both the boundary conditions. This formulation is in parallel with the standard time-independent Rayleigh-Schrodinger perturbation theory. The efficacy of the method is tested by comparing the perturbative values against the numerically calculated eigenvalues for spheroidal, superegg and superquadric shaped boundaries. It is shown that this perturbation works quite well even for wide departure from spherical shape and for higher excited states too. We believe this formulation would find applications in the field of quantum dots and acoustical cavities.

Euler Solution Using Adaptive Cartesian Grid With A Gridless Boundary Treatment

A quadtree-based adaptive Cartesian grid generator and flow solver were developed. The grid adaptation based on pressure or density gradient was performed and a gridless method based on the least-square fashion was used to treat the wall surface boundary condition, which is generally difficult to be handled for the common Cartesian grid. First, to validate the technique of grid adaptation, the benchmarks over a forward-facing step and double Mach reflection were computed. Second, the flows over the NACA 0012 airfoil and a two-element airfoil were calculated to validate the developed gridless method. The computational results indicate the developed method is reasonable for complex flows.

Direct numerical simulation of turbulent flows over superhydrophobic surfaces with gas pockets using linearized boundary conditions

Superhydrophobic surfaces are shown to be effective for surface drag reduction under laminar regime by both experiments and simulations (see for example, Ou and Rothstein, Phys. Fluids 17:103606, 2005). However, such drag reduction for fully developed turbulent flow maintaining the Cassie-Baxter state remains an open problem due to high shear rates and flow unsteadiness of turbulent boundary layer. Our work aims to develop an understanding of mechanisms leading to interface breaking and loss of gas pockets due to interactions with turbulent boundary layers. We take advantage of direct numerical simulation of turbulence with slip and no-slip patterned boundary conditions mimicking the superhydrophobic surface. In addition, we capture the dynamics of gas-water interface, by deriving a proper linearized boundary condition taking into account the surface tension of the interface and kinematic matching of interface deformation and normal velocity conditions on the wall. We will show results from our simulations predicting the dynamical behavior of gas pocket interfaces over a wide range of dimensionless surface tensions.

Boundary conditions at the solid-liquid surface over the multiscale channel size from nanometer to micron

The boundary condition at the solid surface is one of the important problems for the microfluidics. In this paper we study the effects of the channel sizes on the boundary conditions (BC), using the hybrid computation scheme adjoining the molecular dynamics (MD) simulations and the continuum fluid mechanics. We could reproduce the three types of boundary conditions (slip, no-slip and locking) over the multiscale channel sizes. The slip lengths are found to be mainly dependent on the interfacial parameters with the fixed apparent shear rate. The channel size has little effects on the slip lengths if the size is above a critical value within a couple of tens of molecular diameters. We explore the liquid particle distributions nearest the solid walls and found that the slip boundary condition always corresponds to the uniform liquid particle distributions parallel to the solid walls, while the no-slip or locking boundary conditions correspond to the ordered liquid structures close to the solid walls. The slip, no-slip and locking interfacial parameters yield the positive, zero and negative slip lengths respectively. The three types of boundary conditions existing in "microscale" still occur in "macroscale". However, the slip lengths weakly dependent on the channel sizes yield the real shear rates and the slip velocity relative to the solid wall traveling speed approaching those with the no-slip boundary condition when the channel size is larger than thousands of liquid molecular diameters for all of the three types of interfacial parameters, leading to the quasi-no-slip boundary conditions.

An augmented method for free boundary problems with moving contact lines

An augmented immersed interface method (IIM) is proposed for simulating one-phase moving contact line problems in which a liquid drop spreads or recoils on a solid substrate. While the present two-dimensional mathematical model is a free boundary problem, in our new numerical method, the fluid domain enclosed by the free boundary is embedded into a rectangular one so that the problem can be solved by a regular Cartesian grid method. We introduce an augmented variable along the free boundary so that the stress balancing boundary condition is satisfied. A hybrid time discretization is used in the projection method for better stability. The resultant Helmholtz/Poisson equations with interfaces then are solved by the IIM in an efficient way. Several numerical tests including an accuracy check, and the spreading and recoiling processes of a liquid drop are presented in detail. (C) 2010 Elsevier Ltd. All rights reserved.

Spectral functions of half-filled one-dimensional Hubbard rings with varying boundary conditions

We study the effect of varying the boundary condition on: the spectral function of a finite one-dimensional Hubbard chain, which we compute using direct (Lanczos) diagonalization of the Hamiltonian. By direct comparison with the two-body response functions and with the exact solution of the Bethe ansatz equations, we can identify both spinon and holon features in the spectra. At half-filling the spectra have the well-known structure of a low-energy holon band and its shadow-which spans the whole Brillouin zone-and a spinon band present for momenta less than the Fermi momentum. Features related to the twisted boundary condition are cusps in the spinon band. We show that the spectral building principle, adapted to account for both the finite system size and the twisted boundary condition, describes the spectra well in terms of single spinon and holon excitations. We argue that these finite-size effects are a signature of spin-charge separation and that their study should help establish the existence and nature of spin-charge separation in finite-size systems.

Boundary value problems for the elliptic sine-Gordon equation in a semi-strip

We study boundary value problems posed in a semistrip for the elliptic sine-Gordon equation, which is the paradigm of an elliptic integrable PDE in two variables. We use the method introduced by one of the authors, which provides a substantial generalization of the inverse scattering transform and can be used for the analysis of boundary as opposed to initial-value problems. We first express the solution in terms of a 2 by 2 matrix Riemann-Hilbert problem whose \jump matrix" depends on both the Dirichlet and the Neumann boundary values. For a well posed problem one of these boundary values is an unknown function. This unknown function is characterised in terms of the so-called global relation, but in general this characterisation is nonlinear. We then concentrate on the case that the prescribed boundary conditions are zero along the unbounded sides of a semistrip and constant along the bounded side. This corresponds to a case of the so-called linearisable boundary conditions, however a major difficulty for this problem is the existence of non-integrable singularities of the function q_y at the two corners of the semistrip; these singularities are generated by the discontinuities of the boundary condition at these corners. Motivated by the recent solution of the analogous problem for the modified Helmholtz equation, we introduce an appropriate regularisation which overcomes this difficulty. Furthermore, by mapping the basic Riemann-Hilbert problem to an equivalent modified Riemann-Hilbert problem, we show that the solution can be expressed in terms of a 2 by 2 matrix Riemann-Hilbert problem whose jump matrix depends explicitly on the width of the semistrip L, on the constant value d of the solution along the bounded side, and on the residues at the given poles of a certain spectral function denoted by h. The determination of the function h remains open.

Application of real and functional analysis to solve boundary value problems

This thesis is about using appropriate tools in functional analysis arid classical analysis to tackle the problem of existence and uniqueness of nonlinear partial differential equations. There being no unified strategy to deal with these equations, one approaches each equation with an appropriate method, depending on the characteristics of the equation. The correct setting of the problem in appropriate function spaces is the first important part on the road to the solution. Here, we choose the setting of Sobolev spaces. The second essential part is to choose the correct tool for each equation. In the first part of this thesis (Chapters 3 and 4) we consider a variety of nonlinear hyperbolic partial differential equations with mixed boundary and initial conditions. The methods of compactness and monotonicity are used to prove existence and uniqueness of the solution (Chapter 3). Finding a priori estimates is the main task in this analysis. For some types of nonlinearity, these estimates cannot be easily obtained, arid so these two methods cannot be applied directly. In this case, we first linearise the equation, using linear recurrence (Chapter 4). In the second part of the thesis (Chapter 5), by using an appropriate tool in functional analysis (the Sobolev Imbedding Theorem), we are able to improve previous results on a posteriori error estimates for the finite element method of lines applied to nonlinear parabolic equations. These estimates are crucial in the design of adaptive algorithms for the method, and previous analysis relies on, what we show to be, unnecessary assumptions which limit the application of the algorithms. Our analysis does not require these assumptions. In the last part of the thesis (Chapter 6), staying with the theme of choosing the most suitable tools, we show that using classical analysis in a proper way is in some cases sufficient to obtain considerable results. We study in this chapter nonexistence of positive solutions to Laplace's equation with nonlinear Neumann boundary condition. This problem arises when one wants to study the blow-up at finite time of the solution of the corresponding parabolic problem, which models the heating of a substance by radiation. We generalise known results which were obtained by using more abstract methods.