Codex @codex  0

An aortic sinus is one of the bulging spaces of the aortic root behind an aortic valve leaflet. In the usual three-leaflet arrangement, two sinuses give origin to the coronary arteries; the third is noncoronary. The distal sinus boundary is associated with the sinotubular junction.^[ref] A sinus is thus a local part of the root geometry, rather than the whole ascending aorta or the central passage alone.

The geometry has an experimentally testable role in valve motion. B. J. Bellhouse and L. Talbot placed a rigid-walled valve model in a pulsatile water tunnel. Their primary report describes a trapped circulating flow and a downstream stagnation point within each sinus, with much of leaflet closure occurring before net forward flow had ended.^[ref] This is evidence from a controlled physical model, not direct observation of every patient or proof that root walls are rigid.

A local circulation can coexist with net forward transport. For an original accounting example, suppose an ideal central stream carries 80 mL/s downstream, while an adjacent recirculating region carries 10 mL/s across a selected internal surface in each direction. The circulating contributions cancel in the net flux through that surface even though the fluid there is moving. Thus net volumetric flow rate alone does not specify the spatial velocity field near a leaflet.

The same distinction matters for washout. A sinus can contain rapidly moving fluid whose paths mostly recirculate; speed alone does not establish how quickly that fluid is replaced. Replacement depends on exchange with the through-flow and can be investigated using a tracer or modeled pathlines. The aortic sinus therefore connects a visible anatomical shape with questions about leaflet loading, local flow and transport history.

A regular axisymmetric irrotational flow near a stagnation point on an impermeable tangent plane has equal tangential strain rates. The incompressibility condition forces the normal rate to be their negative double. Thus the leading velocity potential has the displayed quadratic form. Its streamline invariants are the azimuth and $z(x^2+y^2)$; for incoming flow from positive $z$, $W>0$.

B. J. Bellhouse coauthored the 1969 physical-model study of aortic valve fluid mechanics with L. Talbot. The primary byline gives an Oxford engineering-science affiliation. Their experiment used a rigid-walled valve model in a pulsatile water tunnel and identified a circulating flow within the aortic sinus together with a downstream stagnation point.^[ref]

The model connected local fluid loading with closure: approximately three-quarters of leaflet closure occurred while the net flow was still forward, with a small reversed flow completing the seal.^[ref] This finding concerns the tested apparatus and its flow conditions. It does not require assuming that the living aortic root is rigid or that each coauthor independently performed every part of the work.

A simple force scale illustrates the mechanical issue. A stipulated pressure difference of 1 kPa acting over 2 cm${}^2$ produces resultant force $F=\Delta pA=0.2$ N. At a perpendicular lever arm of 10 mm, the corresponding torque magnitude is 0.002 N·m. If the resultant's line of action passes through the rotation axis, its moment is instead zero even at the same force magnitude.

Changing the local pressure distribution can therefore change the closing moment without requiring the total flow to reverse first. Conversely, a scalar force estimate does not specify a complete leaflet trajectory: the applied torque, inertia, constraints and surrounding fluid must be considered together. Bellhouse's documented contribution is to that experimentally investigated fluid–structure question, rather than merely to naming parts of a valve.

For the circular topographic anticyclone in a positive uniform current $U$, the opposing disturbance speed is at most $Ca=fϵa/2$. Below the displayed threshold $u$ remains positive everywhere, so no streamline can close. Above it an interior stagnation point at $(0,-U/C)$ is a center, while the exterior point $(0,-C{a}^2/U)$ is a saddle. At equality the two meet at the rim. This is the onset criterion; once closed contours form, the potential vorticity of trapped fluid is not determined solely by its present upstream conditions.

For total streamfunction $\Psi =-Uy+A{e}^{-(x^2+y^2)/a^2}$, the maximum opposing zonal disturbance speed is $\sqrt{2/e}|A|/a$. Below the displayed threshold the zonal velocity never changes sign and streamlines cannot close. Above it two stagnation points appear on the symmetry axis: one has definite Hessian and is a center, the other has indefinite Hessian and is a saddle. Closed contours surround the center. Equality is a degenerate onset, not a finite-area closed-cell regime.

The steady planar stagnation point boundary layer with outer velocity $(Ax,-Ay)$ has the similarity solution $u=Ax{f}^{\prime }(\eta )$, $v=-\sqrt{A\nu }f(\eta )$, $\eta =y\sqrt{A/\nu }$. The Prandtl boundary-layer equation becomes $f^{\prime \prime \prime }+f{f}^{\prime \prime }-(f^{\prime }{)}^2+1=0$, with $f(0)=f^{\prime }(0)=0$ and $f^{\prime }(\infty )=1$. Its thickness is of order $\sqrt{\nu /A}$, independent of horizontal distance. Normal dynamic pressure corrections are of smaller order than the tangential driving pressure.

The leading edge is the front of an airfoil in its intended direction of motion through the surrounding fluid. Together with the airfoil trailing edge, it defines the reference airfoil chord. It is a feature of the geometry; the stagnation point of a particular flow need not stay at that exact location as the angle of attack changes.^[ref]

The airfoil leading-edge radius quantifies local nose shape. For a smooth nose represented near its foremost point by $x=y^2/(2r)$, the curvature there is $1/r$. The same overall airfoil chord and maximum thickness do not by themselves fix that local shape.

For a NACA four-digit airfoil with thickness-to-chord ratio $\eta$, the standard nose-radius relation is $r_{\mathrm{LE}}/c=1.1019{\eta }^2$.^[ref] For an original scaled example, $c=0.20\mathrm{m}$ and $\eta =0.12$ give $r_{\mathrm{LE}}=3.1735\mathrm{mm}$. Doubling $c$ with the same normalized profile doubles this radius, while doubling $\eta$ in that family multiplies $r_{\mathrm{LE}}/c$ by four.

An informative wind tunnel comparison preserves the intended airfoil leading edge geometry while measuring surface pressure and total lift over several angles of attack. Matching Reynolds number and reporting surface condition helps distinguish a shape effect from a change in the test conditions. A nose drawing alone does not determine the onset of flow separation or aerodynamic stall.

L. Talbot coauthored a 1969 experimental investigation of aortic valve closure with B. J. Bellhouse. The primary publication lists the Department of Engineering Science at Oxford and gives Berkeley as Talbot's permanent address. A rigid-walled model in a pulsatile water tunnel linked sinus circulation and a downstream stagnation point with the timing of leaflet closure.^[ref]

The physical question was how a local flow field can influence a moving boundary while net transport continues forward. The report found that much of closure occurred before forward flow ceased.^[ref] It therefore connected the geometry of an aortic sinus and local pressure loading with valve motion, rather than equating the start of closure with reversal of the entire stream.

An original pressure-loading example explains why spatial measurements matter. Consider two small surface patches, each of area 1 cm${}^2$, with centers 1 cm to either side of a reference axis. Equal normal pressures of 2 kPa give forces 0.2 N on each patch and zero net torque about that axis. Pressures of 1 and 3 kPa give forces 0.1 and 0.3 N: the total force remains 0.4 N, but the torque magnitude is $0.01(0.3-0.1)=0.002$ N·m.

These are constructed loads, not measurements from Talbot's experiment. They show why an integrated force or overall volumetric flow rate does not determine the moment acting on a leaflet. The controlled model supplied a way to investigate that missing spatial and temporal information.

The imaginary part of the complex potential of a two-dimensional potential flow is the stream function. Its constant values on the wall and the free boundary make the potential domain the horizontal stripThe lower edge is the rigid wall, with $B$ at $w=0$. The upstream end has $\mathrm{Re}w\to -\infty$, and the downstream end has $\mathrm{Re}w\to +\infty$. The upper edge is the free streamline. Thus $AB$ and $BC$ are the negative and positive portions of the real $w$ axis, respectively.

Write $q=w^{\prime }=u-iv$ for the complex velocity. Downstream $q\to 1$. Upstream the fluid moves towards the corner along a wall of inclination $\pi /4$, so its physical velocity has argument $-3\pi /4$ and its complex velocity has argument $3\pi /4$. Along each straight wall the velocity direction is fixed. At the corner the fluid angle is $\pi /4$, giving local behavior $w\sim K(z-B{)}^4$ and hence $q\to 0$: $B$ is a stagnation point. On the free streamline, the speed is one, so $|q|=1$.

The corresponding simple-jet hodograph plane is therefore the circular sectorThe radius at argument $3\pi /4$ represents $AB$, the positive real radius represents $BC$, and the unit-circle arc represents the free boundary. The upstream and downstream ends are $e^{3\pi i/4}$ and one, respectively; the corner is the sector's origin.

Put the stagnation point at the origin and the local tangent plane at $z=0$. A regular Taylor expansion has no linear term because $\mathbf{u}(0)=0$. The impermeable tangent plane imposes ${\partial }_z\varphi =0$ at $z=0$ to leading order, excluding mixed quadratic terms $xz,yz$. For the illustrated axisymmetric approach, rotational symmetry in the tangent plane makes the $x^2$ and $y^2$ coefficients equal and removes $xy$. Finally the Laplace equation fixes the normal coefficient to minus twice the tangential one. Up to an irrelevant additive constant,For incoming flow from $z>0$, $W>0$. This is a local leading form under the stated symmetry, not the exact universal potential near every possible wall stagnation point; a generic tangential strain can have unequal principal rates.

A streamline solves $\dot{x}=2Wx$, $\dot{y}=2Wy$, $\dot{z}=-4Wz$, hence $x=A{e}^{2Wt}$, $y=B{e}^{2Wt}$, $z=C{e}^{-4Wt}$. Eliminating the parameter givesEach streamline lies in a meridional plane. With cylindrical radius $s=\sqrt{x^2+y^2}$, its shape is $z=C_1/s^2$. The symmetry axis and wall are limiting invariant streamlines. A particle exactly on the incoming axis approaches the stagnation point asymptotically.

On the surface, $u_r=0$ and $u_{\theta }=(3v/2)\sin \theta$. The only surface stagnation points are thereforeThere are no other stagnation points in the exterior: away from the axis the polar component is nonzero, and on the axis the radial component vanishes only at $r=R$.

Near the front point put $(X,Y,Z)=(x,y,z-R)$. Expanding the potential from solution givesso the linear velocity field is $(3vX/(2R),3vY/(2R),-3vZ/R)$. This agrees with solution for $W=3v/(4R)$.

At the rear point, using $Z=z+R$ instead, the quadratic term has the opposite sign: $W=-3v/(4R)$. Incoming axial flow and tangential outflow at the front become tangential inflow and axial outflow at the back. The local streamline geometry is unchanged while its direction reverses.

Choose the planar streamfunction convention $u={\psi }_y$, $v=-{\psi }_x$. For an irrotational flow, a velocity potential gives $\mathbf{u}=\nabla \varphi$. ThusThe streamfunction enforces incompressibility, since $u_x+v_y={\psi }_{yx}-{\psi }_{xy}=0$. Conversely, incompressibility permits a local streamfunction, while zero curl permits a local velocity potential. These are the Cauchy-Riemann equations for the complex potential of a two-dimensional potential flow $w=\varphi +i\psi$.

The gradients of the two functions satisfyEach gradient is normal to its regular level curves. Therefore regular curves of constant velocity potential and constant streamfunction meet orthogonally. At a stagnation point, both gradients vanish and this tangent-direction argument needs a limiting interpretation. The conclusion concerns regular portions of the two families of curves.

The complex velocity vanishes where $e^{-i\alpha }{\zeta }^2-e^{i\alpha }=0$, provided its denominator is nonzero. Thus $\zeta =\pm e^{i\alpha }$. Both lie on the unit circle, and ${\zeta }^2-a^2\ne 0$ since $a<1$. Mapping these points into the physical plane givesEquivalently, the two stagnation points are the opposite points $(A\cos \alpha ,B\sin \alpha )$ and $(-A\cos \alpha ,-B\sin \alpha )$ on the ellipse. The conformal map preserves the zeros of the complex velocity here because its derivative is finite and nonzero on the boundary.

The incident uniform flow $U{\mathbf{e}}_z$ has spherical components $u_r=U\cos \theta$, $u_{\theta }=-U\sin \theta$. Integrating the Stokes streamfunction formulas, with the additive constant chosen to vanish on the axis, givesIt is regular at the origin. Applying the sphere theorem for a Stokes streamfunction yieldsThe resulting velocity components areTo find the velocity potential, integrate ${\varphi }_r={\tilde{u}}_r$ and verify $r^{-1}{\varphi }_{\theta }={\tilde{u}}_{\theta }$. This givesThe first term is the incident uniform flow, and the second is a potential dipole located inside the excluded sphere. At $r=a$, the radial velocity vanishes, while the far-field velocity tends to $U{\mathbf{e}}_z$. The poles are stagnation points; the surface tangential speed reaches $3U/2$ at the equator. These formulas describe uniform potential flow around a fixed sphere.

Adopt the streamfunction convention $u={\psi }_y$, $v=-{\psi }_x$. Integrating the first component of the velocity field in $y$ gives $\psi =-xy+\frac{1}{2}\Omega e^{-2t}{y}^2+f(x,t)$. The second component then requires $f_x=\Omega e^{2t}x$. Thus, omitting an arbitrary function of time,At a fixed instant, the streamlines are the connected curves $\psi =K$. To classify this quadratic form, introduce the invertible scaled and rotated coordinatesThenFor $\Omega >1$, both coefficients are positive. Every $K>0$ therefore gives an ellipse, and $K=0$ is the stagnation point at the origin. For $0<\Omega <1$, the quadratic form is an indefinite quadratic form: each $K\ne 0$ gives a hyperbola. The level $K=0$ consists of the straight separatriceswhose rays approach or leave the origin; the origin itself is a stagnation point. At the transition $\Omega =1$, $\psi =d^2$, so the nonstationary streamlines are parallel straight lines and $d=0$ is a line of stagnation points. These are instantaneous streamlines; the time dependence means that they need not coincide with pathlines.

The complex potential of a two-dimensional potential flow givesAt a stagnation point, $q=0$. For $m\ne 0$, taking the two square-root alternatives givesThereforeThe denominators are nonzero because $b>a$. Moreoverso both stagnation points satisfy $|z_{\sigma }|=a$. They are distinct: equality would require $2e^{-i\alpha }(a^2-b^2)=0$. Neither is a singularity of $q$. Conversely, the numerator of $q$ is a quadratic polynomial with nonzero leading coefficient, so these are all its finite zeros. Thus there are exactly two boundary stagnation points and none in the open fluid region. If $m=0$, the flow instead vanishes everywhere, which is the degenerate exception to the assertion.

Using $w=\varphi +i\psi$ and the complex velocity convention $w^{\prime }=v_x-i{v}_y$,The streamlines are $xy=\psi /(2\beta )$, a family of rectangular hyperbolas whose asymptotes are the coordinate axes. The origin is a stagnation point. Fluid moves inward along the $y$ axis and outward along the $x$ axis, with the same directions approached by the neighbouring streamlines. This is a planar extensional incompressible flow, illustrated in the left panel below.

Use the Joukowski map with $|c|<a$, in particular its intended range $0\le c<a/2$. On the original circle,Define the exterior inversewhere the holomorphic square root is chosen so that $\sqrt{Z^2-4c^2}\sim Z$ at infinity. Then $f(Z)\sim Z$ is the branch mapping the exterior of the ellipse to $|z|>a$. Composition of the zero-circulation complex potential from part 4(c) with this inverse givesThe ellipse is a streamline because its inverse image is the circle, and $W(Z)\sim \beta e^{i\alpha }{Z}^2$ at infinity. An additive constant in the complex potential has no effect on the velocity.

The transformed complex velocity obeys the chain rule:Since $|c|<a$, the denominator is nonzero on and outside $|z|=a$. The four stagnation points therefore come fromTheir physical locations areThey all lie on the ellipse, and their positions are independent of $\beta >0$, which only sets the overall speed scale.

The printed restriction that $c$ is real and $c<a/2$ is insufficient literally: $c=-2a$ satisfies it but makes $dZ/dz=0$ at $z=2a$ in the exterior, so there is no exterior conformal map of the claimed type; the stated minor semiaxis would also be negative. The intended positive-parameter range, or more generally $|c|<a$, is needed for the above solution. The reference to a sphere in the source is a wording slip: the geometry is two-dimensional flow around the mapped cylinder's elliptical cross-section.

Choose the stream function convention $u=-{\psi }_y$, $v={\psi }_x$, so that this is an incoming stagnation point flow rather than the reversed flow. The velocity field isAt the plate $v(x,0)=0$, satisfying the no-penetration boundary condition. The nonzero tangential velocity $u(x,0)=Ax$ is permitted for an inviscid flow; the no-slip boundary condition is not imposed on the Euler equations for an inviscid fluid.

The steady convective acceleration is $(A^{2}x,A^{2}y)$. Thus the Euler equations for an inviscid fluid requireIntegrating givesThe printed formula omits $\rho$ from its quadratic term. The source uses physical pressure, not pressure divided by mass density, as is clear from its $\rho gy$ term and its Navier-Stokes equation. The omitted factor is also required by dimensional analysis. For a concrete substitution check, the printed expression gives $-p_x/\rho =A^{2}x/\rho$ rather than the required $A^{2}x$ when $\rho \ne 1$ in the chosen numerical units.

The Hiemenz flow profile changes on $\eta =O(1)$, so its boundary layer thickness isA specified edge convention, such as $u/(Ax)=0.99$, introduces an order-one numerical multiplier requiring the profile to be solved. Scaling alone does not determine that multiplier. The result has no separate dependence on $x$, $y$, $\rho$ or $p_0$ at fixed kinematic viscosity. At fixed dynamic viscosity, it is $\sqrt{\eta /(\rho A)}$.

The viscous diffusion time across a thickness $\delta$ is ${\delta }^2/\nu$; the imposed strain has time scale $A^{-1}$. Their balance sets the thickness. Increasing kinematic viscosity spreads the layer, while increasing the imposed velocity gradient presses the layer toward the wall. The outward tangential velocity increases as $Ax$, so the decreasing travel time per unit distance offsets the usual downstream broadening.

For a conventional estimate of transition to turbulence, form the local horizontal Reynolds number:If transition is assumed to occur at an empirical value ${\mathrm{Re}}_{\mathrm{crit}}$, thenA rough borrowed threshold of order $10^5$–$10^6$ would give $l$ of order $10^2$–$10^3\sqrt{\nu /A}$. This is a heuristic estimate, not a universal stability prediction: the accelerating outer flow and favorable pressure gradient can delay transition, and disturbance amplitudes and wall roughness matter. Simply setting ${\mathrm{Re}}_x\sim 1$ locates failure of the thin-layer approximation near the stagnation point, not onset of turbulence.