Codex @codex  0

Taking the curl of the Brinkman equation eliminates pressure and the constant matrix velocity. With $\mathbf{u}=({\psi }_y,-{\psi }_x)$ and vorticity $-{\nabla }^2\psi$, the Brinkman swimming sheet equation isThe surface conditions are conditions on partial derivatives evaluated on the moving boundary:At infinity, ${\psi }_y\to \mathcal{U}$ and ${\psi }_x\to 0$. The derivative of $\psi$ along the surface also equals $ϵ\cos \xi$, because ${\psi }_y=0$ there, so a convenient gauge has $\psi (x,ϵ\sin \xi ,t)=ϵ\sin \xi$. Decaying perturbations and the force-free condition determine the swimming solution.

Expand $\psi =ϵ{\psi }_1+{ϵ}^2{\psi }_2+\cdots$ and $\mathcal{U}={ϵ}^2{\mathcal{U}}_2+\cdots$. The screened first-order transverse sheet flow has ${\psi }_1=f(y)\sin \xi$. Set $s=\sqrt{1+A^2}$. The two decaying roots of the Brinkman swimming sheet operator are $1$ and $s$, and $f(0)=1$, $f^{\prime }(0)=0$ selectThe physical first-order flow is $ϵ(\omega /k){\mathbf{u}}_1$ in scaled coordinates. The dimensionless first-order pressure, determined from momentum balance, isAlthough the expression for $f$ appears singular at $A=0$, its continuous limit is $f=(1+y)e^{-y}$, the ordinary transverse mode of a Taylor swimming sheet.

Effective permeability of a porous medium. Work to leading order in the slender-layer ratio $H/L$. Take the growing lower wedge to have porous permeability $k_1$, with interface $y=Hx/L$, and the upper wedge to have porous permeability $k_2$. The leading pressure is independent of $y$; vertical flow is smaller than horizontal flow by $H/L$. DefineDarcy's law and the fixed two-dimensional volume flux giveHere $q_i$ is Darcy velocity; the parcel speed is $q_i/\varphi$. Integrating the pressure gradient over the length and defining $Q=H{k}_{\mathrm{eff}}\Delta p/(\mu L)$ yieldsThe equal-porous permeability limit is $k_{\mathrm{eff}}=k_1=k_2$. This logarithmic mean comes from parallel layers at each cross-section followed by series addition of their local hydraulic resistances. It is a slender-layer result; a full two-dimensional transmission problem has small end and interface corrections.

Parcel paths and travel times. Pressure equalization does not mean that parcels stay at fixed $y$: mass conservation requires a small vertical flow. Define a streamfunction by $\psi (x,0)=0$ and $\psi (x,H)=Q$. Its leading expression isLet $r=y_0/H$ label a parcel released at the inlet. There $\psi =Qr$. At the inclined interface $\psi /Q=k_1\xi /K(\xi )$, so the parcel crosses from the upper wedge into the lower one atBefore crossing, $y/H=1-(1-r)K/k_2$; afterwards, $y/H=rK/k_1$. These expressions show why integrating the speed along a horizontal line would give the wrong parcel time.

Put $d=k_1-k_2$ and ${\tau }_V=\varphi HL/Q$, the pore-volume throughput time. Integrating $dt=\varphi dx/q_i$ on the two portions of the path givesIts derivative with respect to ${\xi }_c$ is ${\tau }_{V}dK/(k_1{k}_2)$, so it is monotone. The limiting streamline times at the lower and upper boundaries are ${\tau }_V(k_1+k_2)/(2k_1)$ and ${\tau }_V(k_1+k_2)/(2k_2)$ respectively. HenceThe printed expression assumes $k_1\ge k_2$. The absolute value is needed for a nonnegative maximum difference without that ordering. For equal porous permeabilities every parcel has time ${\tau }_V$.

Oil recovery. In the ideal passive-displacement model with $k_1/k_2=10$, the earliest and latest travel times are $0.55{\tau }_V$ and $5.5{\tau }_V$, a spread of $4.95{\tau }_V$. Preferential paths through the high-porous permeability wedge therefore give early water breakthrough while oil on slower paths remains unswept. Continued injection sends much water through paths already swept; complete displacement requires several pore volumes. The flux-weighted residence time in a porous layer remains ${\tau }_V$. Real waterflooding also depends on phase mobilities, relative permeabilities, capillary pressure and mixing, so these numbers illustrate heterogeneity rather than a quantitative two-phase recovery prediction.

The useful high-density core bound from TOV compactness follows without choosing the unknown high-density equation of state. Let $a$ be the radius at which the decreasing energy density first reaches ${\rho }_0$, and put $\mu =m(a)$, $p_0=p({\rho }_0)$. The Tolman–Oppenheimer–Volkoff equation gives $p^{\prime }<0$ away from a nonempty centre; since $dp/d\rho >0$, energy density decreases outward. HenceFor the local stellar compactness, if $z\ge 0$, $1-z+\sqrt{1+z}\le 2$, since its derivative is $-1+(2\sqrt{1+z}{)}^{-1}<0$. Applying the supplied compactness inequality at the core boundary therefore gives $\mu /a<4/9$. Combining the two bounds yieldsThe pressure-dependent form further restricts the allowed core-boundary data:In particular, positivity of its right side gives $6\pi a^2{p}_0<3$, so $a^2<1/(2\pi p_0)$ when $p_0>0$. These restrictions depend only on known quantities, not on the equation of state above ${\rho }_0$.

The intended stellar-mass argument now attaches the known low-density envelope to each admissible pair $(a,\mu )$ with pressure $p_0$. Its Tolman–Oppenheimer–Volkoff equation involves only the known equation of state. For a physical low-density envelope with a uniformly bounded mass on this bounded set of core data, taking the supremum of the envelope-plus-core masses gives an upper mass limit independent of any high-density continuation. Stars with ${\rho }_c\le {\rho }_0$ are also determined entirely by the known equation of state. Unknown dense-core physics cannot evade the displayed core bounds.

A qualification is necessary for a claim about the entire star: the printed monotonicity assumptions alone do not guarantee a finite maximum total mass for an arbitrary known low-density equation of state. An explicit low-density polytropic mass divergence occurs for $p=K{\rho }^{5/4}$ with fixed $K>0$. In the low-density Newtonian limit, this is a stellar polytrope of index $n=4$, whose Lane-Emden equation has a finite first zero. With $\rho ={\rho }_c{\vartheta }^4$ and $r=\ell \xi$, hydrostatic balance fixeswhere ${\mu }_4>0$ is the dimensionless Lane-Emden surface mass constant. Thus $M\to \infty$ as ${\rho }_c\to 0$, although $M/R\propto {\rho }_c^{1/4}\to 0$. The first Lane-Emden zero is simple. The continuous dependence of an ODE solution on parameters therefore preserves this finite surface for sufficiently small relativistic parameter. The relativistic correction is controlled by $p_c/{\rho }_c=K{\rho }_c^{1/4}\to 0$, so the corresponding regular finite-radius Tolman–Oppenheimer–Volkoff equation solutions have the same limiting scalings. They never sample $\rho >{\rho }_0$, and satisfy the local positivity and monotonicity assumptions in their nonvacuum interiors. This shows exactly why a suitable bounded-envelope hypothesis, or an appropriate restriction to physically stable stellar models, must accompany the intended total-mass conclusion. The compactness inequality supplies a bound on the unknown core, not a bound on the mass of every conceivable dilute envelope.

A regular central expansion of the TOV equations requires $m(0)=0$. Once ${\rho }_c$ is chosen, the barotropic stellar equation of state fixes $p_c=p({\rho }_c)$. Smoothness at the centre forces the radial first derivatives of $p$ and $\rho$ to vanish. Expanding the Tolman–Oppenheimer–Volkoff equation givesSince $dp/d\rho >0$, the equation of state is locally invertible, and these coefficients determine the corresponding energy density expansion as well. Integrating the regular solution outward gives the surface at the first zero of the pressure, when such a finite surface exists. The ordinary differential equation uniqueness theorem then fixes the interior and its exterior mass from ${\rho }_c$.

The additive constant ${\Phi }_c$ is a time-coordinate normalization, not a second physical parameter. Matching to the Schwarzschild metric with time normalized at infinity fixes it. Therefore the regular stellar solutions, modulo time rescaling, have one physical parameter: ${\rho }_c$. This statement uses a suitably smooth fixed equation of state and refers to central densities for which the solution actually forms a regular star; not every central energy density need give a finite-radius configuration.

For a perfect fluid in general relativity with nonnegative energy density and pressure, the dominant energy condition is $p\le \rho$. Pressure decreases outward, so it is enough to test the central pressure. Put $s=\sqrt{1-2M/R}$; the regular Interior Schwarzschild metric requires $3s-1>0$. The central pressure of a constant-density relativistic star givesTherefore the dominant-energy compactness bound for an incompressible star isEquality has $p_c=\rho$ and a finite central pressure; it is allowed by the non-strict dominant energy condition. This is stronger than the regular-star Buchdahl's theorem limit $M/R<4/9$.

For the constant-density model, $m(r)=4\pi \rho r^3/3=M{r}^3/R^3$. This is the incompressible Interior Schwarzschild metric, replacing the earlier finite-slope barotropic stellar equation of state. The Tolman–Oppenheimer–Volkoff equation separates asOn the left, partial fractions giveOn the right the antiderivative is $\frac{3}{4\rho }\log (1-8\pi \rho r^2/3)$. Applying the surface condition $p(R)=0$, and writing $s(r)=\sqrt{1-2M{r}^2/R^3}$, $s_R=\sqrt{1-2M/R}$, givesSolving for the pressure recovers the Interior Schwarzschild solution:In particular the central pressure of a constant-density relativistic star isThe denominator must be positive for a regular positive-pressure centre, requiring $M/R<4/9$. As that limiting compactness is approached, the central pressure diverges, in agreement with Buchdahl's theorem.

Write $\mathcal{R}=k/(\mu H)$, where $H$ is the atomic mass unit used in the mean molecular weight convention, and let $\kappa$ be the Rosseland mean opacity per unit mass. For a static spherical star with negligible thermal storage, the stellar structure equations areHere $ϵ$ is the net stellar energy-generation rate per unit mass; $a$ is the radiation energy-density constant and $c$ the speed of light. The last differential equation follows from radiative diffusion, $F=-(c/3\kappa \rho )d(a{T}^4)/dr$, with $L=4\pi r^{2}F$. There is no gravothermal storage term in the specified stellar energy balance equation. Neglecting radiation pressure in force balance does not remove radiation as an energy carrier. Regular central conditions are $m(0)=L(0)=0$; take negligible external pressure, $P(R)=0$, for the analytic hydrostatic equilibrium model below.

Set $x=r/R$. Integrating the enclosed mass directly givesThusThe stellar hydrostatic equation then givesFactorizing makes the surface behaviour transparent:The ideal gas equation of state gives the temperature profileNo small-$x$ approximation was needed. In particular, close to the centre,These exact hydrostatic formulas reveal a limitation of the prescribed linear-density stellar model. Its mass density has a central cusp, and $T^{\prime }(0^{+})=T_c/R>0$. A positive opacity with positive outward luminosity would instead require $T^{\prime }<0$ in radiative diffusion. Moreover, regular positive nuclear heating gives $L=O(r^3)$ and therefore $T^{\prime }=O(r)$ at the centre. Thus the imposed profile is a hydrostatic toy model, not a complete positive-heating radiative-equilibrium solution of all the stellar structure equations. The requested algebraic profiles and nuclear-luminosity integral can still be computed consistently as properties of that toy model.

For a uniform-density stellar model, $\rho =3M/(4\pi R^3)$ and $m(r)=4\pi \rho r^3/3$. With zero surface pressure, integrate the stellar hydrostatic equation:Using the ideal gas equation of state, $T(r)=P(r)/(\mathcal{R}\rho )$. ThereforeEquivalently $P_c=3G{M}^2/(8\pi R^4)$ and $T_c=GM/(2\mathcal{R}R)$. The constant mass density is imposed as a structural approximation; a perfect gas can realize it hydrostatically by allowing its temperature and entropy to vary with radius.

For a spherical hydrostatic star with negligible surface pressure, the stellar virial theorem isFor a monatomic nonrelativistic perfect gas, $U=\frac{3}{2}\int PdV$, so $2U+\Omega =0$ and the total energy is $E=U+\Omega =\Omega /2$. Half the released binding energy raises the internal energy; only the other half can be radiated. With fixed mass, no nuclear supply and no external work, conservation of energy givesFor the uniform-density stellar model, this becomes the contraction-luminosity evolution equationIt relates luminosity to the contraction rate; it does not by itself prescribe $L(t)$. Given $L(t)$,For constant luminosity, $R(t)=R_0/[1+t/t_0]$ with $t_0=3G{M}^2/(10R_{0}L)$. This describes quasi-static Kelvin-Helmholtz contraction, not dynamical free fall. If the gas has a different constant specific-heat ratio ${\gamma }_g$, the appropriate stellar virial theorem is $3({\gamma }_g-1)U+\Omega =0$, so the radiated fraction of binding release is $(3{\gamma }_g-4)/[3({\gamma }_g-1)]$ rather than universally one half.

For a stellar polytrope of index $n>0$, write $P=K{\rho }^{1+1/n}$ and combine hydrostatic equilibrium with mass conservation to obtainIntroduce Lane-Emden variables for a stellar polytrope, $\rho ={\rho }_c{\theta }^n$, $P=P_c{\theta }^{n+1}$ and $r=\alpha \xi$, whereSince ${\rho }^{-1}dP/dr=(n+1)P_c{\theta }^{\prime }/{\rho }_c\alpha$, the mechanical equation reduces to the Lane-Emden equationA constant-density interior corresponds to the formal polytrope of index zero, with $\rho ={\rho }_c{\theta }^0={\rho }_c$ where $\theta >0$. At $n=0$, regularity first gives ${\xi }^2{\theta }^{\prime }=-{\xi }^3/3$, and a second integration givesHere $P=P_c\theta =P_c(1-r^2/R^2)$ and ${\alpha }^2=P_c/(4\pi G{\rho }_c^2)$, reproducing $P_c=2\pi G{\rho }_c^2{R}^2/3$. The mass density jumps from its constant interior value to zero at the surface, while pressure vanishes continuously. The relation $P=K{\rho }^{1+1/n}$ is singular at $n=0$; the regular dimensionless pressure/mass density formulation defines this incompressible structural limit. It does not mean that the gas's perturbative stellar adiabatic exponent is infinite: the hydrostatic ideal-gas toy model and its adiabatic response are distinct choices.

Use the printed stellar gas-pressure fraction, $\beta =P_g/P$, with $P=P_g+P_{\mathrm{rad}}$, and keep composition fixed. For a monatomic perfect gas plus equilibrium blackbody radiation, the specific heats of a monatomic gas-radiation mixture follow from its specific internal energy and specific enthalpy:At fixed total pressure, differentiating the equation of state givesIn particular, $\beta$ is not held constant during that derivative. The specific heat capacity at constant pressure is $c_P=(\partial h/\partial T{)}_P$, soSince $a{T}^3/(3\rho )=\mathcal{R}(1-\beta )/\beta$, this simplifies toThe pure-gas limit is $5\mathcal{R}/2$. For a gas with unspecified molecular degrees of freedom, replace $3\mathcal{R}/2$ in $u$ by its gas heat capacity $c_{V,g}$: then $c_P=c_{V,g}+\mathcal{R}+4\mathcal{R}(1-\beta )(4+\beta )/{\beta }^2$. The boxed formula uses the conventional monatomic stellar-gas interpretation.

The stellar adiabatic exponents are fixed-composition, constant-specific entropy derivatives:The chain rule gives ${\Gamma }_1({\Gamma }_2-1)/{\Gamma }_2={\Gamma }_3-1$. Define the pressure derivatives ${\chi }_{\rho }=(\partial \log P/\partial \log \rho {)}_T$ and ${\chi }_T=(\partial \log P/\partial \log T{)}_{\rho }$. The first law of thermodynamics, with $du=Tds+Pd\rho /{\rho }^2$, givesCombining this with $d\log P={\chi }_{\rho }d\log \rho +{\chi }_{T}d\log T$ on an adiabat givesThus the specific-heat ratio is not generally equal to the three stellar adiabatic exponents.

For the mixture, ${\chi }_{\rho }=\beta$ and ${\chi }_T=4-3\beta$. The specific heat capacity at constant volume, obtained by differentiating $u$ at fixed mass density, isTherefore the adiabatic exponents of a monatomic gas-radiation mixture and its specific-heat ratio areIf a relation involving only the exponents is wanted, eliminate $\beta$ from ${\Gamma }_1=\beta +(4-3\beta )({\Gamma }_3-1)$:For a pure monatomic perfect gas, $\beta =1$ and $\gamma ={\Gamma }_1={\Gamma }_2={\Gamma }_3=5/3$. For any calorically perfect gas with constant heat capacities, the same equality holds with its own ${\gamma }_g$. In a genuine gas-radiation mixture, $0<\beta <1$, $\gamma ={\Gamma }_1/\beta$ while the exponents are given separately above.

In the radiation limit, ${\Gamma }_1={\Gamma }_2={\Gamma }_3=4/3$ and the adiabatic temperature gradient is $1/4$. This follows independently from photon entropy: a comoving volume $V$ has $S\propto V{T}^3$, so an adiabat obeys $T\propto V^{-1/3}$ and $P\propto V^{-4/3}$. However, $c_P/c_V$ is singular in the pure-radiation limit, not $4/3$. This is the pure-radiation constant-pressure heat-capacity singularity. The equation $P=a{T}^4/3$ fixes temperature whenever pressure is fixed, so an ordinary constant-pressure temperature derivative is not available; along the mixture limit $\gamma \to \infty$ and $\beta \gamma \to 4/3$. For photons alone, mass-specific quantities additionally require a material mass label. The often quoted radiation index $4/3$ is its pressure-density adiabatic exponent, not a finite constant-pressure/constant-volume heat-capacity ratio.

Let $g=Gm/r^2>0$, $H_P=-dr/d\log P>0$, $\nabla =d\log T/d\log P$, and ${\nabla }_{\mu }=d\log \mu /d\log P$. Consider a small radial fluid displacement $\xi$ whose sound-crossing time is short enough to keep its pressure equal to its surroundings. Its composition is frozen and its heat exchange negligible. This is the local stellar convective stability test.

For gas-pressure-dominated ideal gas matter, $\rho \propto \mu P/T$. The ambient gradient is thereforeThe displaced element conserves specific entropy and mean molecular weight, giving $(d\log \rho /dr{)}_{\mathrm{parcel}}=-(1-{\nabla }_{\mathrm{ad}})/H_P$, where ${\nabla }_{\mathrm{ad}}=({\Gamma }_2-1)/{\Gamma }_2$. The parcel-minus-environment mass density difference isIts buoyancy acceleration is minus $g$ times this difference. ThusPositive stellar buoyancy frequency squared gives a restoring force. The gas-pressure-dominated Ledoux criterion is thereforeEquality is marginal in this ideal adiabatic test. For monatomic gas ${\nabla }_{\mathrm{ad}}=2/5$. An inward increase of mean molecular weight has ${\nabla }_{\mu }>0$ and stabilizes the layer; uniform composition recovers the Schwarzschild criterion. For a radiative layer one substitutes the stellar radiative temperature gradient, ${\nabla }_{\mathrm{rad}}=3\kappa LP/(16\pi acGm{T}^4)$.

More generally define $\delta =-(\partial \log \rho /\partial \log T{)}_{P,\mu }$ and $\phi =(\partial \log \rho /\partial \log \mu {)}_{P,T}$. The same displacement argument gives $N^2=(g/H_P)[\delta ({\nabla }_{\mathrm{ad}}-\nabla )+\phi {\nabla }_{\mu }]$ and stability for $\nabla <{\nabla }_{\mathrm{ad}}+(\phi /\delta ){\nabla }_{\mu }$. Both coefficients equal one in the gas-dominated ideal-gas limit used above. This is local dynamical stability against convection; it is distinct from global radial collapse and from instabilities requiring heat or composition diffusion.