## Rate of change of mass of a parameterized region

### Remember to vote up questions/answers you find interesting or helpful (requires 15 reputation points)

Let $R_t$ be a family of compact, simply connected regions in the plane defined by

`$R_t = \{x\in\mathbb{R}^2 : h(x) \leq t\}$`

for all $t$, where $h(x)$ is some nicely behaved smooth function. Suppose $f(x)$ is a probability density on $\mathbb{R}^2$ and define

$M(t) = \iint_{R_t} f(x) dA$

for all $t$. Is it true that

$\frac{dM}{dt}|_{t=t_0} = \int_{\partial R_{t_0}} f(x) /\|\nabla h\| ds$

where $ds$ denotes integration with respect to arc length? If not, what is the right expression for $\frac{dM}{dt}$? I assume this is some well-known first-year calculus-type problem but I can't find it stated in any context (it may very well be a common homework problem though I've not seen it).

-

2
This is not a good question for MO. See en.wikipedia.org/wiki/Coarea_formula – Anton Petrunin May 9 2012 at 0:51
This is a special case of the coarea formula. – Liviu Nicolaescu May 9 2012 at 17:58
Thanks to both of you, I wasn't familiar with that result! That clears this question up for me nicely. – Jennifer Gao May 11 2012 at 0:26

## 1 Answer

With $H$ the Heaviside function (characteristic function of $\mathbb R_+$), you have $$M(t)=\int_{\mathbb R^n} f(x) H(t-h(x)) dx$$ and thus, at least formally, $$\dot M(t)=\int_{\mathbb R^n} f(x) \delta_0(t-h(x)) dx= \int_{\mathbb R^n} f(x)\Vert\nabla h(x)\Vert ^{-1} \underbrace{\delta_0(t-h(x))\Vert\nabla h(x)\Vert dx}_{d\sigma} ,$$ where $d\sigma$ is the Euclidean surface measure on {$x, h(x) =t$}. Here, it is important to assume that $h$ is say $C^1$ such that $dh\not=0$ at $h=t$ for $t$ in neighborhood of some distinguished value $t_0$. As a result, $$\dot M(t)=\int_{\partial R(t)}f(x)\Vert\nabla h(x)\Vert ^{-1} d\sigma.$$ The previous computation can be justified by Green's formula: for $X$ a $C^1_c$ vector field on $\Omega$ ( an open subset of $\mathbb R^n$ with a $C^1$ boundary) $$\int_\Omega \text{div} X\ dx=\int_{\partial \Omega} X\cdot \nu\ d\sigma,$$ where $\nu$ is the exterior unit normal to $\partial \Omega$, and $d\sigma$ is the Euclidean surface measure on $\partial \Omega$. That formula can be proven by showing $$d\sigma=\lim_{\epsilon\rightarrow 0} \phi(\frac{\rho(x)}{\epsilon})\epsilon^{-1}\Vert\nabla \rho(x)\Vert\quad\text{(distribution sense),}$$ for $\phi\in C^\infty_c(\mathbb R),\int \phi(t) dt=1$, $\Omega=${$x, \rho(x)<0$}, $d\rho\not=0$ at $\rho=0$.

-
{"extraction_info": {"found_math": true, "script_math_tex": 0, "script_math_asciimath": 0, "math_annotations": 0, "math_alttext": 0, "mathml": 0, "mathjax_tag": 0, "mathjax_inline_tex": 34, "mathjax_display_tex": 5, "mathjax_asciimath": 0, "img_math": 0, "codecogs_latex": 0, "wp_latex": 0, "mimetex.cgi": 0, "/images/math/codecogs": 0, "mathtex.cgi": 0}, "config": {"markdown_headings": true, "markdown_code": true, "boilerplate_config": {"ratio_threshold": 0.18, "absolute_threshold": 10, "end_threshold": 5, "enable": true}, "remove_buttons": true, "remove_image_figures": true, "remove_link_clusters": true, "table_config": {"min_rows": 2, "min_cols": 3, "format": "github"}, "remove_chinese": true, "remove_edit_buttons": true, "extract_latex": true}, "english_confidence": 0.9266827702522278, "perplexity_flag": "head"}