
We initiate the study of the finiteness condition
$\int_{\Omega}u(x)^{\beta}\,dx\leq C(\Omega,\beta)<+\infty$ where
$\Omega\subseteq{\mathbb{R}}^n$ is an open set and $u$ is the solution of the
Saint Venant problem $\Delta u=1$ in $\Omega$, $u=0$ on $\partial\Omega$. The
central issue which we address is that of determining the range of values of
the parameter $\beta>0$ for which the aforementioned condition holds under
various hypotheses on the smoothness of $\Omega$ and demands on the nature of
the constant $C(\Omega,\beta)$. Classes of domains for which our analysis
applies include bounded piecewise $C^1$ domains in ${\mathbb{R}}^n$, $n\geq 2$,
with conical singularities (in particular polygonal domains in the plane),
polyhedra in ${\mathbb{R}}^3$, and bounded domains which are locally of class
$C^2$ and which have (finitely many) outwardly pointing cusps. For example, we
show that if $u_N$ is the solution of the Saint Venant problem in the regular
polygon $\Omega_N$ with $N$ sides circumscribed by the unit disc in the plane,
then for each $\beta\in(0,1)$ the following asymptotic formula holds: %
{eqnarray*} \int_{\Omega_N}u_N(x)^{\beta}\,dx=\frac{4^\beta\pi}{1\beta}
+{\mathcal{O}}(N^{\beta1})\quad{as}\,\,N\to\infty. {eqnarray*} % One of the
original motivations for addressing the aforementioned issues was the study of
sublevel set estimates for functions $v$ satisfying $v(0)=0$, $\nabla v(0)=0$
and $\Delta v\geq c>0$.

We consider two types of multilinear pseudodifferential operators. First, we
prove the boundedness of multilinear pseudodifferential operators with symbols
which are only measurable in the spatial variables in weighted Lebesgue spaces.
These results generalise earlier work of the present authors concerning linear
pseudopseudodifferential operators. Secondly, we investigate the boundedness
of bilinear pseudodifferential operators with symbols in the H\"ormander
$S^{m}_{\rho, \delta}$ classes. These results are new in the case $\rho < 1$,
that is, outwith the scope of multilinear Calder\'onZygmund theory.

We prove the global $L^2 \times L^2 \to L^1$ boundedness of bilinear Fourier
integral operators with amplitudes in $S^0_{1,0} (n,2)$. To achieve this, we
require that the phase function can be written as $(x,\xi,\eta) \mapsto
\phase_1(x,\xi) + \phase_2(x,\eta)$ where each $\phase_j$ belongs to the class
$\Phi^2$ and satisfies the strong nondegeneracy condition. This result extends
that of R. Coifman and Y. Meyer regarding pseudodifferential operators to the
case of Fourier integral operators.