
Under the classical nonrelativistic consideration of the spacetime we
propose the model of the laws of gravity and Electrodynamics, invariant under
the galilean transformations and moreover, under every change of noninertial
cartesian coordinate system. Being in the frames of nonrelativistic model of
the spacetime, we adopt some general ideas of the General Theory of
Relativity, like the assumption of invariance of the most general physical laws
in every inertial and noninertial coordinate system and equivalence of
factious forces in noninertial coordinate systems and the force of gravity.
Moreover, in the frames of our model, we obtain that the laws of
Nonrelativistic Quantum Mechanics also invariant under the change of inertial
or noninertial cartesian coordinate system.

Motivated by the formula, due to Bourgain, Brezis and Mironescu,
\begin{equation*}
\lim_{\varepsilon\to 0^+} \int_\Omega\int_\Omega
\frac{u(x)u(y)^q}{xy^q}\,\rho_\varepsilon(xy)\,dx\,dy=K_{q,N}\\nabla
u\_{L^{q}}^q\,, \end{equation*} that characterizes the functions in $L^q$ that
belong to $W^{1,q}$ (for $q>1$) and $BV$ (for $q=1$), respectively, we study
what happens when one replaces the denominator in the expression above by
$xy$. It turns out that, for $q>1$ the corresponding functionals "see" only
the jumps of the $BV$ function. We further identify the function space relevant
to the study of these functionals, the space $BV^q$, as the Besov space
$B^{1/q}_{q,\infty}$. We show, among other things, that $BV^q(\Omega)$ contains
both the spaces $BV(\Omega)\cap L^\infty(\Omega)$ and $W^{1/q,q}(\Omega)$. We
also present applications to the study of singular perturbation problems of
AvilesGiga type.

Motivated by the study of non abelian Chern Simons vortices of non
topological type in Gauge Field Theory, we analyse the solvability of planar
Liouville systems of Toda type in presence of singular sources. We identify
necessary and sufficient conditions on the "flux" pair which ensure the radial
solvability of the system. Since the given system includes the (integrable) 2 X
2 Toda system as a particular case, thus we recover the existence result
available in this case. Our method relies on a blowup analysis, which even in
the radial setting, takes new turns compared with the single equation case.

We use the technique developed in [32][33] to construct the upper and the
lower bounds for classes of problems containing nonlocal terms, including
problems in micromagnetics and problems arising in the variational study of the
Method of Vanishing Viscosity for systems of conservation laws. We reduced
these problems to the problems considered in [32][33], with the appropriate
prescribed differential constraint.

In part II we constructed the lower bound, in the spirit of $\Gamma$
$\liminf$ for some general classes of singular perturbation problems, with or
without the prescribed differential constraint, taking the form
E_\e(v):=\int_\Omega \frac{1}{\e}F\Big(\e^n\nabla^n v,...,\e\nabla
v,v\Big)dx\quad\text{for}\;\; v:\Omega\subset\R^N\to\R^k\;\;\text{such
that}\;\; A\cdot\nabla v=0, where the function $F\geq 0$ and $A:\R^{k\times
N}\to\R^m$ is a prescribed linear operator (for example, $A:\equiv 0$,
$A\cdot\nabla v:=\text{curl}\, v$ and $A\cdot\nabla v=\text{div} v$).
Furthermore, we studied the cases where we can easy prove the coinciding of
this lower bound and the upper bound obtained in [33]. In particular we find
the formula for the $\Gamma$limit for the general class of anisotropic
problems without a differential constraint (i.e., in the case $A:\equiv 0$).

Using our results in [15], we provided existence theorems for the general
classes of nonlinear evolutions. Finally, we give examples of applications of
our results to parabolic, hyperbolic, Shr\"{o}dinger, NavierStokes and other
timedependent systems of equations.

We develop a variational technique for some wide classes of nonlinear
evolutions. The novelty here is that we derive the main information directly
from the corresponding EulerLagrange equations. In particular, we prove that
not only the minimizer of the appropriate energy functional but also any
critical point must be a solution of the corresponding evolutional system.

In Part I we construct the upper bound, in the spirit of $\Gamma$ $\limsup$,
achieved by multidimensional profiles, for some general classes of singular
perturbation problems, with or without the prescribed differential constraint,
taking the form $$E_\e(v):=\int_\Omega \frac{1}{\e}F\Big(\e^n\nabla^n
v,...,\e\nabla v,v\Big)dx\quad\text{for} v:\Omega\subset\R^N\to\R^k \text{such
that} A\cdot\nabla v=0,$$ where the function $F\geq 0$ and $A:\R^{k\times
N}\to\R^m$ is a prescribed linear operator (for example, $A:\equiv 0$,
$A\cdot\nabla v:=\text{curl}v$ and $A\cdot\nabla v=\text{div}\,v$) which
includes, in particular, the problems considered in [27]. This bound is in
general sharper then one obtained in [27].