38,603 research outputs found
Elimination for Systems of Algebraic Differential Equations
We develop new upper bounds for several effective differential elimination techniques for systems of algebraic ordinary and partial differential equations. Differential elimination, also known as decoupling, is the process of eliminating a fixed subset of unknown functions from a system of differential equations in order to obtain differential algebraic consequences of the original system that do not depend on that fixed subset of unknowns. A special case of differential elimination, which we study extensively, is the question of consistency, that is, if the given system of differential equations has a solution. We first look solely at the ``algebraic data of the system of differential equations through the theory of differential kernels to provide a new upper bound for proving the consistency of the system. We then prove a new upper bound for the effective differential Nullstellensatz, which determines a sufficient number of times to differentiate the original system in order to prove its inconsistency. Finally, we study the Rosenfeld-Gröbner algorithm, which approaches differential elimination by decomposing the given system of differential equations into simpler systems. We analyze the complexity of the Rosenfeld-Gröbner algorithm by computing an upper bound for the orders of the derivatives in all intermediate steps and in the output of the algorithm
The GNAT method for nonlinear model reduction: effective implementation and application to computational fluid dynamics and turbulent flows
The Gauss--Newton with approximated tensors (GNAT) method is a nonlinear
model reduction method that operates on fully discretized computational models.
It achieves dimension reduction by a Petrov--Galerkin projection associated
with residual minimization; it delivers computational efficency by a
hyper-reduction procedure based on the `gappy POD' technique. Originally
presented in Ref. [1], where it was applied to implicit nonlinear
structural-dynamics models, this method is further developed here and applied
to the solution of a benchmark turbulent viscous flow problem. To begin, this
paper develops global state-space error bounds that justify the method's design
and highlight its advantages in terms of minimizing components of these error
bounds. Next, the paper introduces a `sample mesh' concept that enables a
distributed, computationally efficient implementation of the GNAT method in
finite-volume-based computational-fluid-dynamics (CFD) codes. The suitability
of GNAT for parameterized problems is highlighted with the solution of an
academic problem featuring moving discontinuities. Finally, the capability of
this method to reduce by orders of magnitude the core-hours required for
large-scale CFD computations, while preserving accuracy, is demonstrated with
the simulation of turbulent flow over the Ahmed body. For an instance of this
benchmark problem with over 17 million degrees of freedom, GNAT outperforms
several other nonlinear model-reduction methods, reduces the required
computational resources by more than two orders of magnitude, and delivers a
solution that differs by less than 1% from its high-dimensional counterpart
Two-branes with variable tension model and the effective Newtonian constant
It is shown that, in the two brane time variation model framework, if the
hidden brane tension varies according to the phenomenological Eotvos law, the
visible brane tension behavior is such that its time derivative is negative in
the past and positive after a specific time of cosmological evolution. This
behavior is interpreted in terms of an useful mechanical system analog and its
relation with the variation of the Newtonian (effective) gravitational
`constant' is explored.Comment: 15 pages, no figure, accepted for publication in Physical Review
Non-linear Dynamics in QED_3 and Non-trivial Infrared Structure
In this work we consider a coupled system of Schwinger-Dyson equations for
self-energy and vertex functions in QED_3. Using the concept of a
semi-amputated vertex function, we manage to decouple the vertex equation and
transform it in the infrared into a non-linear differential equation of
Emden-Fowler type. Its solution suggests the following picture: in the absence
of infrared cut-offs there is only a trivial infrared fixed-point structure in
the theory. However, the presence of masses, for either fermions or photons,
changes the situation drastically, leading to a mass-dependent non-trivial
infrared fixed point. In this picture a dynamical mass for the fermions is
found to be generated consistently. The non-linearity of the equations gives
rise to highly non-trivial constraints among the mass and effective (`running')
gauge coupling, which impose lower and upper bounds on the latter for dynamical
mass generation to occur. Possible implications of this to the theory of
high-temperature superconductivity are briefly discussed.Comment: 29 pages LATEX, 7 eps figures incorporated, uses axodraw style.
Discussion on the massless case (section 2) modified; no effect on
conclusions, typos correcte
- …