A combined integrating and differentiating matrix formulation for boundary value problems on rectangular domains

Integrating and differentiating matrices allow the numerical integration and differential of functions whose values are known at points of a discrete grid. Previous derivations of these matrices were restricted to one dimensional grids or to rectangular grids with uniform spacing in at least one direction. Integrating and differentiating matrices were developed for grids with nonuniform spacing in both directions. The use of these matrices as operators to reformulate boundary value problems on rectangular domains as matrix problems for a finite dimensional solution vector is considered. The method requires nonuniform grids which include near boundary points. An eigenvalue problem for the transverse vibrations of a simply supported rectangular plate is solved to illustrate the method

Two-dimensional integrating matrices on rectangular grids

The use of integrating matrices in solving differential equations associated with rotating beam configurations is examined. In vibration problems, by expressing the equations of motion of the beam in matrix notation, utilizing the integrating matrix as an operator, and applying the boundary conditions, the spatial dependence is removed from the governing partial differential equations and the resulting ordinary differential equations can be cast into standard eigenvalue form. Integrating matrices are derived based on two dimensional rectangular grids with arbitrary grid spacings allowed in one direction. The derivation of higher dimensional integrating matrices is the initial step in the generalization of the integrating matrix methodology to vibration and stability problems involving plates and shells

Stability of the laminar boundary layer in a streamwise corner

The stability of viscous, incompressible flow along a streamwise corner, often called the corner boundary layer problem is examined. The semi-infinite boundary value problem satisfied by small amplitude disturbances in the "bending boundary layer' region is obtained. The mean secondary flow induced by the corner exhibits a flow reversal in this region. Uniformly valid "first approximations' to solutions of the governing differential equations are derived. Uniformity at infinity is achieved by a suitable choice of the large parameter and use of an approximate Langer variable. Approximations to solutions of balanced type have a phase shift across the critical layer which is associated with instabilities in the case of two dimensional boundary layer profiles

Integrating matrix formulations for vibrations of rotating beams including the effects of concentrated masses

By expressing partial differential equations of motion in matrix notation, utilizing the integrating matrix as a spatial operator, and applying the boundary conditions, the resulting ordinary differential equations can be cast into standard eigenvalue form upon assumption of the usual time dependence. As originally developed, the technique was limited to beams having continuous mass and stiffness properties along their lengths. Integrating matrix methods are extended to treat the differential equations governing the flap, lag, or axial vibrations of rotating beams having concentrated masses. Inclusion of concentrated masses is shown to lead to the same kind of standard eigenvalue problem as before, but with slightly modified matrices

Higher modes of the Orr-Sommerfeld problem for boundary layer flows

The discrete spectrum of the Orr-Sommerfeld problem of hydrodynamic stability for boundary layer flows in semi-infinite regions is examined. Related questions concerning the continuous spectrum are also addressed. Emphasis is placed on the stability problem for the Blasius boundary layer profile. A general theoretical result is given which proves that the discrete spectrum of the Orr-Sommerfeld problem for boundary layer profiles (U(y), 0,0) has only a finite number of discrete modes when U(y) has derivatives of all orders. Details are given of a highly accurate numerical technique based on collocation with splines for the calculation of stability characteristics. The technique includes replacement of 'outer' boundary conditions by asymptotic forms based on the proper large parameter in the stability problem. Implementation of the asymptotic boundary conditions is such that there is no need to make apriori distinctions between subcases of the discrete spectrum or between the discrete and continuous spectrums. Typical calculations for the usual Blasius problem are presented

Existence and non-uniqueness of similarity solutions of a boundary layer problem

A Blasius boundary value problem with inhomogeneous lower boundary conditions f(0) = 0 and f'(0) = - lambda with lambda strictly positive was considered. The Crocco variable formulation of this problem has a key term which changes sign in the interval of interest. It is shown that solutions of the boundary value problem do not exist for values of lambda larger than a positive critical value lambda. The existence of solutions is proven for 0 lambda lambda by considering an equivalent initial value problem. It is found however that for 0 lambda lambda, solutions of the boundary value problem are nonunique. Physically, this nonuniqueness is related to multiple values of the skin friction

On similarity solutions of a boundary layer problem with an upstream moving wall

The problem of a boundary layer on a flat plate which has a constant velocity opposite in direction to that of the uniform mainstream is examined. It was previously shown that the solution of this boundary value problem is crucially dependent on the parameter which is the ratio of the velocity of the plate to the velocity of the free stream. In particular, it was proved that a solution exists only if this parameter does not exceed a certain critical value, and numerical evidence was adduced to show that this solution is nonunique. Using Crocco formulation the present work proves this nonuniqueness. Also considered are the analyticity of solutions and the derivation of upper bounds on the critical value of wall velocity parameter

Self-Similar Blowup Solutions to the 2-Component Camassa-Holm Equations

In this article, we study the self-similar solutions of the 2-component Camassa-Holm equations% \begin{equation} \left\{ \begin{array} [c]{c}% \rho_{t}+u\rho_{x}+\rho u_{x}=0 m_{t}+2u_{x}m+um_{x}+\sigma\rho\rho_{x}=0 \end{array} \right. \end{equation} with \begin{equation} m=u-\alpha^{2}u_{xx}. \end{equation} By the separation method, we can obtain a class of blowup or global solutions for $\sigma=1$ or $-1$. In particular, for the integrable system with $\sigma=1$, we have the global solutions:% \begin{equation} \left\{ \begin{array} [c]{c}% \rho(t,x)=\left\{ \begin{array} [c]{c}% \frac{f\left( \eta\right) }{a(3t)^{1/3}},\text{ for }\eta^{2}<\frac {\alpha^{2}}{\xi} 0,\text{ for }\eta^{2}\geq\frac{\alpha^{2}}{\xi}% \end{array} \right. ,u(t,x)=\frac{\overset{\cdot}{a}(3t)}{a(3t)}x \overset{\cdot\cdot}{a}(s)-\frac{\xi}{3a(s)^{1/3}}=0,\text{ }a(0)=a_{0}% >0,\text{ }\overset{\cdot}{a}(0)=a_{1} f(\eta)=\xi\sqrt{-\frac{1}{\xi}\eta^{2}+\left( \frac{\alpha}{\xi}\right) ^{2}}% \end{array} \right. \end{equation} where $\eta=\frac{x}{a(s)^{1/3}}$ with $s=3t;$ $\xi>0$ and $\alpha\geq0$ are arbitrary constants.\newline Our analytical solutions could provide concrete examples for testing the validation and stabilities of numerical methods for the systems.Comment: 5 more figures can be found in the corresponding journal paper (J. Math. Phys. 51, 093524 (2010) ). Key Words: 2-Component Camassa-Holm Equations, Shallow Water System, Analytical Solutions, Blowup, Global, Self-Similar, Separation Method, Construction of Solutions, Moving Boundar

Observables in the Decays of B to Two Vector Mesons

In general there are nine observables in the decay of a B meson to two vector mesons defined in terms of polarization correlations of these mesons. Only six of these can be detected via the subsequent decay angular distributions because of parity conservation in those decays. The remaining three require the measurement of the spin polarization of one of the decay products.Comment: 12 pages, no figur

Learning Two-input Linear and Nonlinear Analog Functions with a Simple Chemical System

The current biochemical information processing systems behave in a predetermined manner because all features are defined during the design phase. To make such unconventional computing systems reusable and programmable for biomedical applications, adaptation, learning, and self-modification baaed on external stimuli would be highly desirable. However, so far, it haa been too challenging to implement these in real or simulated chemistries. In this paper we extend the chemical perceptron, a model previously proposed by the authors, to function as an analog instead of a binary system. The new analog asymmetric signal perceptron learns through feedback and supports MichaelisMenten kinetics. The results show that our perceptron is able to learn linear and nonlinear (quadratic) functions of two inputs. To the best of our knowledge, it is the first simulated chemical system capable of doing so. The small number of species and reactions allows for a mapping to an actual wet implementation using DNA-strand displacement or deoxyribozymes. Our results are an important step toward actual biochemical systems that can learn and adapt