an online method for interpolating linear parametric reduced order models a two step online method is proposed for interpolating projection based linear parametric reduced order models roms in order to construct a new rom for a new set of parameter values the first step of this method transforms each precomputed rom into a consistent set of generalized coordinates the second step interpolates the associated linear operators on their appropriate matrix manifold real time performance is achieved by precomputing inner products between the reduced order bases underlying the precomputed roms the proposed method is illustrated by applications in mechanical and aeronautical engineering in particular its robustness is demonstrated by its ability to handle the case where the sampled parameter set values exhibit a mode veering phenomenon

asymptotic expansion of solutions to the drift diffusion equation with large initial data abstract we consider the large time behavior of the solution to the initial value problem for the nernst planck type drift diffusion equation in whole spaces in the l p framework the global existence and the decay of the solution were shown moreover the second order asymptotic expansion of the solution as t was derived we also deduce the higher order asymptotic expansion of the solution especially we discuss the contrast between the odd dimensional case and the even dimensional case

on the existence of local strong solutions to chemotaxis shallow water system with large data and vacuum abstract in this paper motivated by the chemotaxis navier stokes system arising from mathematical biology 43 a modified shallow water type chemotactic model is derived for large initial data allowing vacuum the local existence of strong solutions together with the blow up criterion is established

chemotaxis with logistic source very weak global solutions and their boundedness properties abstract we consider the chemotaxis system u t u u v g u x t 0 0 v v u x t 0 in a smooth bounded domain r n where 0 and g generalizes the logistic function g u a u b u with 1 a 0 and b 0 a concept of very weak solutions is introduced and global existence of such solutions for any nonnegative initial data u 0 l 1 is proved under the assumption that 2 1 n moreover boundedness properties of the constructed solutions are studied inter alia it is shown that if b is sufficiently large and u 0 l has small norm in l for some n 2 then the solution is globally bounded finally in the case that additionally n 2 holds a bounded set in l can be found which eventually attracts very weak solutions emanating from arbitrary l 1 initial data the paper closes with numerical experiments that illustrate some of the theoretically established results

computer model calibration using high dimensional output this work focuses on combining observations from field experiments with detailed computer simulations of a physical process to carry out statistical inference of particular interest here is determining uncertainty in resulting predictions this typically involves calibration of parameters in the computer simulator as well as accounting for inadequate physics in the simulator the problem is complicated by the fact that simulation code is sufficiently demanding that only a limited number of simulations can be carried out we consider applications in characterizing material properties for which the field data and the simulator output are highly multivariate for example the experimental data and simulation output may be an image or may describe the shape of a physical object we make use of the basic framework of kennedy and o hagan however the size and multivariate nature of the data lead to computational challenges in implementing the framework to overcome these challenges we make use of basis repre

nonlinear aggregation diffusion equations radial symmetry and long time asymptotics we analyze under which conditions equilibration between two competing effects repulsion modeled by nonlinear diffusion and attraction modeled by nonlocal interaction occurs this balance leads to continuous compactly supported radially decreasing equilibrium configurations for all masses all stationary states with suitable regularity are shown to be radially symmetric by means of continuous steiner symmetrization techniques calculus of variations tools allow us to show the existence of global minimizers among these equilibria finally in the particular case of newtonian interaction in two dimensions they lead to uniqueness of equilibria for any given mass up to translation and to the convergence of solutions of the associated nonlinear aggregation diffusion equations towards this unique equilibrium profile up to translations as t rightarrow infty

global boundedness in a fully parabolic quasilinear chemotaxis system with singular sensitivity abstract in this paper we study the global boundedness of solutions to the quasilinear fully parabolic chemotaxis system u t d u u s u v v t v v u where bounded domain r n n 2 subject to the non flux boundary conditions the diffusivity fulfills d u a 0 u 1 with a 0 0 and 0 while the density signal governed sensitivity satisfies 0 s u b 0 u 1 and 0 v v k for b 0 0 and k r it is shown that the solution is globally bounded provided 1 and k 1 this result demonstrates the effect of signal dependent sensitivity on the blow up prevention

definability in the lattice of equational theories of commutative semigroups in this paper we study first order definability in the lattice of equational theories of commutative semigroups in a series of papers j jezek solving problems posed by a tarski and r mckenzie has proved in particular that each equational theory is first order definable in the lattice of equational theories of a given type up to automorphism and that such lattices have no automorphisms besides the obvious syntactically defined ones with exceptions for special unary types he has proved also that the most important classes of theories of a given type are so definable in a later paper jezek and mckenzie have almost proved the same facts for the lattice of equational theories of semigroups there were good reasons to believe that the same can be proved for the lattice of equational theories of commutative semigroups in this paper however we show that the case of commutative semigroups is different

singularity formation in chemotaxis a conjecture of nagai consider the initial boundary value problem for the system s ut uxx uvx x vt u av on an interval 0 1 for t 0 where a 0 with ux 0 t ux 1 t 0 suppose mu v0 are positive constants the corresponding spatially homogeneous global solution u t mu v t mu a v0 mu a exp at is stable in the sense that if mu v0 are positive constants the corresponding spatially homogeneous solution will be uniformly close to u cdot v cdot we consider in sequence space an approximate system s which is related to s in the following sense the chemotactic term uvx x is replaced by the inverse fourier transform of the finite part of the convolution integral for the fourier transform of uvx x here the finite part of the convolution on the line at a point x of two functions f g is defined as int 0 x f y g y x dy we prove the following if mu a then in every neighborhood of mu v0 there are spatially nonconstant initial data for which the solution of proble

aggregation vs global diffusive behavior in the higher dimensional keller segel model abstract we consider the classical parabolic parabolic keller segel system u t u u v x t 0 v t v v u x t 0 under homogeneous neumann boundary conditions in a smooth bounded domain r n it is proved that in space dimension n 3 for each q n 2 and p n one can find e 0 0 such that if the initial data u 0 v 0 satisfy u 0 l q e and v 0 l p e then the solution is global in time and bounded and asymptotically behaves like the solution of a discoupled system of linear parabolic equations in particular u v approaches the steady state m m as t where m is the total mass m u 0 of the population moreover we shall show that if is a ball then for arbitrary prescribed m 0 there exist unbounded solutions emanating from initial data u 0 v 0 having total mass u 0 m

galerkin proper orthogonal decomposition methods for a general equation in fluid dynamics error estimates for galerkin proper orthogonal decomposition pod methods for nonlinear parabolic systems arising in fluid dynamics are proved for the time integration the backward euler scheme is considered the asymptotic estimates involve the singular values of the pod snapshot set and the grid structure of the time discretization as well as the snapshot locations

existence of global bounded classical solution to a quasilinear attraction repulsion chemotaxis system with logistic source abstract this paper concerns the global existence of bounded classical solution to a quasilinear attraction repulsion chemotaxis system with logistic source u t d u u u v u w u u 2 in v t v u v in w t w u w in under homogeneous neumann boundary condition with positive parameters and d s c 0 s m 1 for s 0 d 0 0 we prove that in dimension three there exists a unique global bounded classical solution provided that m 6 5 and 0 with an additional assumption d s c 0 s m 1 1 for s 0 we prove that for any m 1 6 5 there exists a constant 0 0 m such that for all 0 the above problem admits a unique global bounded classical solution

analysis for a system of coupled reaction diffusion parabolic equations arising in biology we study a system where the equations are coupled by transport and reaction terms in the general case we prove local existence and uniqueness of the solution and in a particular case for a priori fixed time t global existence on 0 t

explicit lower bound of blow up time for an attraction repulsion chemotaxis system abstract in this paper we study classical solutions to the zero flux attraction repulsion chemotaxis system u t u u v u w in 0 t 0 v u v in 0 t 0 w u w in 0 t where is a smooth and bounded domain of r 2 t is the blow up time and are positive real numbers from the literature it is known that under a proper interplay between the above parameters and suitable assumptions on the initial data u x 0 u 0 c 0 system has a unique classical solution which becomes unbounded as t t the main result of this investigation is to provide an explicit lower bound for t estimated in terms of u 0 2 d x and attained by means of well established techniques based on ordinary differential inequalities

a new look at proper orthogonal decomposition we investigate some basic properties of the proper orthogonal decomposition pod method as it is applied to data compression and model reduction of finite dimensional nonlinear systems first we provide an analysis of the errors involved in solving a nonlinear ode initial value problem using a pod reduced order model then we study the effects of small perturbations in the ensemble of data from which the pod reduced order model is constructed on the reduced order model we explain why in some applications this sensitivity is a concern while in others it is not we also provide an analysis of computational complexity of solving an ode initial value problem and study the computational savings obtained by using a pod reduced order model we provide several examples to illustrate our theoretical results

reduced basis approximation for nonlinear parametrized evolution equations based on empirical operator interpolation we present a new approach to treating nonlinear operators in reduced basis approximations of parametrized evolution equations our approach is based on empirical interpolation of nonlinear differential operators and their frechet derivatives efficient offline online decomposition is obtained for discrete operators that allow an efficient evaluation for a certain set of interpolation functionals an a posteriori error estimate for the resulting reduced basis method is derived and analyzed numerically we introduce a new algorithm the podei greedy algorithm which constructs the reduced basis spaces for the empirical interpolation and for the numerical scheme in a synchronized way the approach is applied to nonlinear parabolic and hyperbolic equations based on explicit or implicit finite volume discretizations we show that the resulting reduced scheme is able to capture the evolution of both smooth and discontinuous solutions in case of symmetries of the problem the approach realizes an automatic and intuitive space compression or even space dimensionality reduction we perform empirical investigations of the error convergence and run times in all cases we obtain a good run time acceleration

global existence versus blow up in a high dimensional keller segel equation with degenerate diffusion and nonlocal aggregation abstract in this paper we study the degenerate keller segel equation with nonlocal aggregation u t u m u b u in r d r where m 1 d 3 and b u 2 u with 2 d by analyzing the interaction between the degenerate diffusion and the nonlocal aggregation we determine the conditions for initial data under which weak solutions globally exist or blow up in finite time with m 1 d d d furthermore a sharper criterion is given for global existence and finite time blow up of weak solutions with m in the subrange 2 d 2 d d d 1 d d

dynamics in a parabolic elliptic two species population competition model with cross diffusion for one species abstract this work deals with a simplified parabolic elliptic population competition model with cross diffusion for one species we first claim that this model has at least a non trivial global smooth solution and then we detect explicit parameter ranges under which the coexistence steady state and the two semi trivial steady states are globally and asymptotically stable

parameter and state model reduction for large scale statistical inverse problems a greedy algorithm for the construction of a reduced model with reduction in both parameter and state is developed for an efficient solution of statistical inverse problems governed by partial differential equations with distributed parameters large scale models are too costly to evaluate repeatedly as is required in the statistical setting furthermore these models often have high dimensional parametric input spaces which compounds the difficulty of effectively exploring the uncertainty space we simultaneously address both challenges by constructing a projection based reduced model that accepts low dimensional parameter inputs and whose model evaluations are inexpensive the associated parameter and state bases are obtained through a greedy procedure that targets the governing equations model outputs and prior information the methodology and results are presented for groundwater inverse problems in one and two dimensions

unification in commutative idempotent monoids abstract embedding sets as a datastructure into resolution based deduction requires a unification algorithm respecting associativity commutativity and idempotence of the constructor function union we shall present such an algorithm and prove its correctness and termination although this algorithm is very simple the enormous proliferation of unifiers restricts its value to theoretical considerations for practical purposes the more special unitary theory of boolean expressions seem to be better suited

application of an adams type inequality to a two chemical substances chemotaxis system abstract this paper deals with positive solutions of the fully parabolic system u t u u v in 0 1 v t v v w in 0 2 w t w w u in 0 under homogeneous neumann boundary conditions or mixed boundary conditions no flux and dirichlet conditions in a smooth bounded domain r n n 4 with positive parameters 1 2 0 and nonnegative smooth initial data u 0 v 0 w 0 in the lower dimensional case n 3 it is proved that for all reasonable initial data solutions of the system exist globally in time and remain bounded in the case n 4 it is shown that in the radially symmetric setting solutions to the neumann boundary value problem of the system exist globally in time and remain bounded if u 0 l 1 8 2 as to the mixed boundary value problem we will establish global existence and boundedness of solutions if u 0 l 1 8 2 without radial symmetry the key ingredients are a lyapunov functional and an adams type inequality a lyapunov functional of the above problems will be constructed and the constant 8 2 is deduced from the critical constant in the adams type inequality this result is regarded as a generalization of the well known 8 problem in the keller segel system to higher dimensions

on explosions of solutions to a system of partial differential equations modelling chemotaxis a system of partial differential equations modelling chemotactic aggregation is analysed keller segel model conditions on the system of paramaters are given implying global existence of smooth solutions in two space dimensions and radially symmetric situations explosion of the bacteria concentration in finite time is shown for a class of initial values

dynamics of a bio reactor model with chemotaxis long time dynamics of solutions to a strongly coupled system of parabolic equations modeling the competition in bio reactors with chemotaxis will be studied in particular we show that the dynamical system possesses a global attractor and that it is strongly uniformly persistent if the trivial steady state is unstable using a result of smith and waltman on perturbation of global attractors we also show that the positive steady state is unique and globally attracting

blow up and decay criteria for a model of chemotaxis for a system of equations introduced by jager and luckhaus 1992 6 as a model of chemotaxis the questions of blow up and global existence criteria are investigated specifically for a convex region a lower bound for the blow up time is derived if the solution blows up and explicit criteria to ensure non blow up are determined

global existence of solutions to a parabolic elliptic chemotaxis system with critical degenerate di ffusion abstract this paper is devoted to the analysis of non negative solutions for a degenerate parabolic elliptic patlak keller segel system with critical nonlinear diffusion in a bounded domain with homogeneous neumann boundary conditions our aim is to prove the existence of a global weak solution under a smallness condition on the mass of the initial data there by completing previous results on nite blow up for large masses under some higher regularity condition on solutions the uniqueness of solutions is proved by using a classical duality technique

stabilization of projection based reduced order models summary a rigorous method for stabilizing projection based linear reduced order models without significantly affecting their accuracy is proposed unlike alternative approaches this method is computationally efficient it requires primarily the solution of a small scale convex optimization problem furthermore it is nonintrusive in the sense that it operates directly on readily available reduced order operators these can be precomputed using any data compression technique including balanced truncation balanced proper orthogonal decomposition proper orthogonal decomposition or moment matching the proposed method is illustrated with three applications the stabilization of the reduction of the computational fluid dynamics based model of a linearized unsteady supersonic flow the reduction of a computational structural dynamics system and the stabilization of the reduction of a coupled computational fluid dynamics computational structural dynamics model of a linearized aeroelastic system in the transonic flow regime copyright 2012 john wiley sons ltd

a method for interpolating on manifolds structural dynamics reduced order models a rigorous method for interpolating a set of parameterized linear structural dynamics reduced order models roms is presented by design this method does not operate on the underlying set of parameterized full order models hence it is amenable to an online real time implementation it is based on mapping appropriately the rom data onto a tangent space to the manifold of symmetric positive definite matrices interpolating the mapped data in this space and mapping back the result to the aforementioned manifold algorithms for computing the forward and backward mappings are offered for the case where the roms are derived from a general galerkin projection method and the case where they are constructed from modal reduction the proposed interpolation method is illustrated with applications ranging from the fast dynamic characterization of a parameterized structural model to the fast evaluation of its response to a given input in all cases good accuracy is demonstrated at real time processing speeds copyright 2009 john wiley sons ltd

model reduction for large scale systems with high dimensional parametric input space a model constrained adaptive sampling methodology is proposed for the reduction of large scale systems with high dimensional parametric input spaces our model reduction method uses a reduced basis approach which requires the computation of high fidelity solutions at a number of sample points throughout the parametric input space a key challenge that must be addressed in the optimization control and probabilistic settings is the need for the reduced models to capture variation over this parametric input space which for many applications will be of high dimension we pose the task of determining appropriate sample points as a pde constrained optimization problem which is implemented using an efficient adaptive algorithm that scales well to systems with a large number of parameters the methodology is demonstrated using examples with parametric input spaces of dimension 11 and 21 which describe thermal analysis and design of a heat conduction fin and compared with statistically based sampling methods for these examples the model constrained adaptive sampling leads to reduced models that for a given basis size have error several orders of magnitude smaller than that obtained using the other methods

nonlinear component analysis as a kernel eigenvalue problem a new method for performing a nonlinear form of principal component analysis is proposed by the use of integral operator kernel functions one can efficiently compute principal components in high dimensional feature spaces related to input space by some nonlinear map for instance the space of all possible five pixel products in 16 16 images we give the derivation of the method and present experimental results on polynomial feature extraction for pattern recognition

reduced dimensional gaussian process emulators of parametrized partial differential equations based on isomap in this paper isomap and kernel isomap are used to dramatically reduce the dimensionality of the output space to efficiently construct a gaussian process emulator of parametrized partial differential equations the output space consists of spatial or spatio temporal fields that are functions of multiple input variables for such problems standard multi output gaussian process emulation strategies are computationally impractical and or make restrictive assumptions regarding the correlation structure the method we develop can be applied without modification to any problem involving vector valued targets and vector valued inputs it also extends a method based on linear dimensionality reduction to response surfaces that cannot be described accurately by a linear subspace of the high dimensional output space comparisons to the linear method are made through examples that clearly demonstrate the advantages of nonlinear dimensionality reduction

boundedness of solutions to a quasilinear chemotaxis haptotaxis model we study global solutions of a class of chemotaxis haptotaxis systems generalizing the prototype u t u 1 m 1 u u u 1 q 1 v u u 1 p 1 w h u w 0 v v u w t v w in a bounded domain r n n 1 with smooth boundary h u w u 1 u r 1 w with parameters m 1 r 1 and positive constants p q it is shown that either max q 1 p 2 p m max m 2 n r or max q 1 p 2 p m r and b 0 is large enough then for any sufficiently smooth initial data there exists a classical solution which is global in time and bounded the results of this paper improve the results of tao and winkler 2014 46 51

global boundedness in a quasilinear chemotaxis system with general density signal governed sensitivity abstract in this paper we study the global boundedness of solutions to the quasilinear parabolic chemotaxis system u t d u u s u v 0 v v u subject to homogeneous neumann boundary conditions and the initial data u 0 in a bounded and smooth domain r n n 2 where the diffusivity d u is supposed to satisfy d u a 0 u 1 with a 0 0 and r while the density signal governed sensitivity fulfills 0 s u b 0 u 1 and 0 v v k for b 0 0 and k r it is shown that the solution is globally bounded if 1 2 n k 2 n with n 3 and k 1 or 1 for k 1 this implies that the large k benefits the global boundedness of solutions due to the weaker chemotactic migration of the signal dependent sensitivity at high signal concentrations moreover when arrives at the critical value we establish the global boundedness of solutions for the coefficient properly small it should be emphasized that the smallness of under k 1 is positively related to the total cellular mass u 0 d x which is attributed to the stronger singularity of v at v 0 for k 1 and the fact that v can be estimated from below by a multiple of u 0 d x in addition distinctive phenomena concerning this model are observed by comparison with the known results

repulsion effects on boundedness in the higher dimensional fully parabolic attraction repulsion chemotaxis system abstract this paper deals with an attraction repulsion chemotaxis system u t d u u u v u w x t 0 1 v t v u v x t 0 2 w t w u w x t 0 under homogeneous neumann boundary conditions in a smooth bounded domain r n n 2 where parameters i i 1 2 and are positive and diffusion coefficient d u c 2 0 satisfies d u 0 for u 0 d u d u m 1 with d 0 and m 1 for all u 0 it is proved that the corresponding initial boundary value problem possesses a unique global bounded classical solution for m 2 2 n in particular in the case 1 2 and the solution is globally bounded if m 2 2 n n 2 n 2 n 2 therefore due to the inhibition of repulsion to the attraction the range of m 2 2 n of boundedness is enlarged and the results of 21 is thus extended to the higher dimensional attraction repulsion chemotaxis system with nonlinear diffusion

global asymptotic stability of constant equilibria in a fully parabolic chemotaxis system with strong logistic dampening abstract this paper deals with nonnegative solutions of the neumann initial boundary value problem for the parabolic chemotaxis system u t u u v u u 2 x t 0 v t v v u x t 0 in bounded convex domains r n n 1 with smooth boundary it is shown that if the ratio is sufficiently large then the unique nontrivial spatially homogeneous equilibrium given by u v 1 is globally asymptotically stable in the sense that for any choice of suitably regular nonnegative initial data u 0 v 0 such that u 0 0 the above problem possesses a uniquely determined global classical solution u v with u v t 0 u 0 v 0 which satisfies u t 1 l 0 and v t 1 l 0 as t

self similar solutions to a parabolic system modeling chemotaxis we study the forward self similar solutions to a parabolic system modeling chemotaxisut u u v vt v u in the whole space r2 where is a positive constant using the liouville type result and the method of moving planes it is proved that self similar solutions u v must be radially symmetric about the origin then the structure of the set of self similar solutions is investigated as a consequence it is shown that there exists a threshold in r2u for the existence of self similar solutions in particular for 0 1 2 there exists a self similar solution u v if and only if r2u 8

a new approach toward stabilization in a two species chemotaxis model with logistic source abstract this paper aims at providing an alternative approach to study global dynamic properties for a two species chemotaxis model with the main novelty being that both populations mutually compete with the other on account of the lotka volterra dynamics more precisely we consider the following neumann initial boundary value problem u t d 1 u 1 u w 1 u 1 u a 1 v x t 0 v t d 2 v 2 v w 2 v 1 a 2 u v x t 0 0 d 3 w w u v x t 0 in a bounded domain r n n 1 with smooth boundary where d 1 d 2 d 3 1 2 1 2 a 1 a 2 are positive constants when a 1 0 1 and a 2 0 1 it is shown that under some explicit largeness assumptions on the logistic growth coefficients 1 and 2 the corresponding neumann initial boundary value problem possesses a unique global bounded solution which moreover approaches a unique positive homogeneous steady state u v w of above system in the large time limit the respective decay rate of this convergence is shown to be exponential when a 1 1 and a 2 0 1 if 2 is suitable large for all sufficiently regular nonnegative initial data u 0 and v 0 with u 0 0 and v 0 0 the globally bounded solution of above system will stabilize toward 0 1 1 as t in algebraic

interpolation method for adapting reduced order models and application to aeroelasticity reduced order models are usually thought of as computationally inexpensive mathematical representations that offer the potential for near real time analysis although most reduced order models can operate in near real time their construction can be computationally expensive as it requires accumulating a large number of system responses to input excitations furthermore reduced order models usually lack robustness with respect to parameter changes and therefore must often be rebuilt for each parameter variation together these two issues underline the need for a fast and robust method for adapting precomputed reduced order models to new sets of physical or modeling parameters to this effect this paper presents an interpolation method based on the grassmann manifold and its tangent space at a point that is applicable to structural aerodynamic aeroelastic and many other reduced order models based on projection schemes this method is illustrated here with the adaptation of computational fluid dynamics based aeroelastic reduced order models of complete fighter configurations to new values of the freestream mach number good correlations with results obtained from direct reduced order model reconstruction high fidelity nonlinear and linear simulations are reported thereby highlighting the potential of the proposed reduced order model adaptation method for near real time aeroelastic predictions using precomputed reduced order model databases

boundedness in a quasilinear parabolic parabolic keller segel system with subcritical sensitivity abstract we consider the quasilinear parabolic parabolic keller segel system u t d u u s u v x t 0 v t v v u x t 0 under homogeneous neumann boundary conditions in a convex smooth bounded domain r n with n 1 it is proved that if s u d u c u with 2 n and some constant c 0 for all u 1 then the classical solutions to the above system are uniformly in time bounded provided that d u satisfies some technical conditions such as algebraic upper and lower growth resp decay estimates as u this boundedness result is optimal according to a recent result by the second author winkler 2010 27 which says that if s u d u c u for u 1 with c 0 and some 2 n n 2 then for each mass m 0 there exist blow up solutions with mass u 0 m in addition this paper also proves a general boundedness result for quasilinear non uniformly parabolic equations by modifying the iterative technique of moser alikakos alikakos 1979 1

complexity of makanin s algorithm the exponent of periodicity is an important factor in estimates of complexity of word unification algorithms we prove that the exponent of periodicity of a minimal solution of a word equation is of order 2 1 07d where d is the length of the equation we also give a lower bound 2 0 29d so our upper bound is almost optimal and exponentially better than the original bound 6d 22d4 2 consequently our result implies an exponential improvement of known upper bounds on complexity of word unification algorithms

model reduction for compressible flows using pod and galerkin projection abstract we present a framework for applying the method of proper orthogonal decomposition pod and galerkin projection to compressible fluids for incompressible flows only the kinematic variables are important and the techniques are well known in a compressible flow both the kinematic and thermodynamic variables are dynamically important and must be included in the configuration space we introduce an energy based inner product which may be used to obtain pod modes for this configuration space we then obtain an approximate version of the navier stokes equations valid for cold flows at moderate mach number and project these equations onto a pod basis the resulting equations of motion are quadratic and are much simpler than projections of the full compressible navier stokes equations

the basin of attraction of the steady states for a chemotaxis model in mathbb r 2 with critical mass we consider the cauchy problem for a parabolic elliptic system in mathbb r 2 which is amathematical model of chemotaxis and also amodel of self attracting particles in the critical mass case we determine the basin of attraction of the steady states for the cauchy problem through a lyapunov functional

the diffusion limit of transport equations ii chemotaxis equations in this paper we use the diffusion limit expansion of transport equations developed earlier t hillen and h g othmer siam j appl math 61 2000 pp 751 775 to study the limiting equation under a variety of external biases imposed on the motion when applied to chemotaxis or chemokinesis these biases produce modification of the turning rate the movement speed or the preferred direction of movement depending on the strength of the bias it leads to anisotropic diffusion to a drift term in the flux or to both in the parabolic limit we show that the classical chemotaxis equation which we call the patlak keller segel alt pksa equation arises only when the bias is sufficiently small using this general framework we derive phenomenological models for chemotaxis of flagellated bacteria of slime molds and of myxobacteria we also show that certain results derived earlier for one dimensional motion can easily be generalized to two or three dimensional motion as well

unboundedness for solutions to a degenerate drift diffusion equation with the l1 supercritical and the energy subcritical exponent abstract we consider large time behavior of weak solutions to a degenerate drift diffusion system related to keller segel system with the l 1 supercritical and the energy subcritical cases under relaxed weight condition it is known that the large time behavior of solutions is classified by the invariant norms of initial data for the l 1 critical case ogawa wakui proved that the negative entropy condition induces the unboundedness of corresponding solutions with the initial data decaying slowly at spacial infinity here the result is a continuous analogy of the l 1 critical case analogous results have been obtained in the theory of nonlinear schrodinger equations

on a parabolic elliptic system with gradient dependent chemotactic coefficient abstract we consider a second order pdes system of parabolic elliptic type with chemotactic terms the system describes the evolution of a biological species u moving towards a higher concentration of a chemical stimuli v in a bounded and open domain of r n in the system considered the chemotaxis sensitivity depends on the gradient of v i e the chemotaxis term has the following expression d i v u v p 2 v where is a positive constant and p satisfies p 1 if n 1 and p 1 n n 1 if n 2 we obtain uniform bounds in l and the existence of global in time solutions for the one dimensional case we prove the existence of infinitely many non constant steady states for p 1 2 for any positive and a given positive mass

on a constrained 2 d navier stokes equation the planar navier stokes equation exhibits in absence of external forces a trivial asymptotics in time nevertheless the appearence of coherent structures suggests non trivial intermediate asymptotics which should be explained in terms of the equation itself motivated by the separation of the different time scales observed in the dynamics of the navier stokes equation we study the well posedness and asymptotic behaviour of a constrained equation which neglects the variation of the energy and moment of inertia

measure valued solutions and the phenomenon of blow down in logarithmic diffusion abstract the existence of solutions of elliptic and parabolic equations with data a measure has always been quite important for the general theory a prominent example being the fundamental solutions of the linear theory in nonlinear equations the existence of such solutions may find special obstacles that can be either essential or otherwise they may lead to more general concepts of solution we give a particular review of results in the field of nonlinear diffusion as a new contribution we study in detail the case of logarithmic diffusion associated with ricci flow in the plane where we can prove existence of measure valued solutions the surprising thing is that these solutions become classical after a finite time in that general setting the standard concept of weak solution is not adequate but we can solve the initial value problem for the logarithmic diffusion equation in the plane with bounded nonnegative measures as initial data in a suitable class of measure solutions we prove that the problem is well posed the phenomenon of blow down in finite time is precisely described initial point masses diffuse into the medium and eventually disappear after a finite time t i m i 4

boundedness vs blow up in a chemotaxis system abstract we determine the critical blow up exponent for a keller segel type chemotaxis model where the chemotactic sensitivity equals some nonlinear function of the particle density assuming some growth conditions for the chemotactic sensitivity function we establish an a priori estimate for the solution of the problem considered and conclude the global existence and boundedness of the solution furthermore we prove the existence of solutions that become unbounded in finite or infinite time in that situation where this a priori estimate fails

stabilization of projection based reduced order models for linear time invariant systems via optimization based eigenvalue reassignment abstract a new approach for stabilizing unstable reduced order models roms for linear time invariant lti systems through an a posteriori post processing step applied to the algebraic rom system is developed the key idea is to modify the unstable eigenvalues of the rom system by moving these eigenvalues into the stable half of the complex plane it is demonstrated that this modification to the rom system eigenvalues can be accomplished using full state feedback a k a pole placement algorithms from control theory this approach ensures that the modified rom is stable provided the system s unstable poles are controllable and observable however the accuracy of the stabilized rom is not guaranteed to remedy this difficulty and guarantee an accurate stabilized rom a constrained nonlinear least squares optimization problem for the stabilized rom eigenvalues in which the error in the rom output is minimized is formulated this optimization problem is small and therefore computationally inexpensive to solve performance of the proposed algorithms is evaluated on two test cases for which roms constructed via the proper orthogonal decomposition pod galerkin method suffer from instabilities

permutability class of a semigroup given a permutation sn a semigroup s is called permutable if x1x2 xn x 1 x 2 x n for all n tuples x1 x2 xn s various permutability conditions were considered for special classes of permutations we define the permutability class of a semigroup s to be the infinite sequence g2 g3 g4 of the invariance groups gn gn s of the operation x1x2 xn on s the aim of this paper is to give a description of such sequences which generalizes a number of old results on permutability conditions in semigroups

infinite time aggregation for the critical patlak keller segel model in ℝ2 we analyze the two dimensional parabolic elliptic patlak keller segel model in the whole euclidean space r2 under the hypotheses of integrable initial data with finite second moment and entropy we first show local in time existence for any mass of free energy solutions namely weak solutions with some free energy estimates we also prove that the solution exists as long as the entropy is controlled from above the main result of the paper is to show the global existence of free energy solutions with initial data as before for the critical mass 8 actually we prove that solutions blow up as a delta dirac at the center of mass when t keeping constant their second moment at any time furthermore all moments larger than 2 blow up as t if initially bounded

nonlinear model reduction via discrete empirical interpolation a dimension reduction method called discrete empirical interpolation is proposed and shown to dramatically reduce the computational complexity of the popular proper orthogonal decomposition pod method for constructing reduced order models for time dependent and or parametrized nonlinear partial differential equations pdes in the presence of a general nonlinearity the standard pod galerkin technique reduces dimension in the sense that far fewer variables are present but the complexity of evaluating the nonlinear term remains that of the original problem the original empirical interpolation method eim is a modification of pod that reduces the complexity of evaluating the nonlinear term of the reduced model to a cost proportional to the number of reduced variables obtained by pod we propose a discrete empirical interpolation method deim a variant that is suitable for reducing the dimension of systems of ordinary differential equations odes of a certain type as presented here it is applicable to odes arising from finite difference discretization of time dependent pdes and or parametrically dependent steady state problems however the approach extends to arbitrary systems of nonlinear odes with minor modification our contribution is a greatly simplified description of the eim in a finite dimensional setting that possesses an error bound on the quality of approximation an application of deim to a finite difference discretization of the one dimensional fitzhugh nagumo equations is shown to reduce the dimension from 1024 to order 5 variables with negligible error over a long time integration that fully captures nonlinear limit cycle behavior we also demonstrate applicability in higher spatial dimensions with similar state space dimension reduction and accuracy results

unification in primal algebras their powers and their varieties this paper examines the unification problem in the class of primal algebras and the varieties they generate an algebra is called primal if every function on its carrier can be expressed just in terms of the basic operations of the algebra the two element boolean algebra is the simplest nontrivial example every truth function can be realized in terms of the basic connectives for example negation and conjunction it is shown that unification in primal algebras is unitary that is if an equation has a solution it has a single most general one two unification algorithms based on equation solving techniques for boolean algebras due to boole and lo wenheim are studied in detail applications include certain finite post algebras and matrix rings over finite fields the former are algebraic models for many valued logics the latter cover in particular modular arithmetic then unification is extended from primal algebras to their direct powers which leads to unitary unification algorithms covering finite post algebras finite semisimple artinian rings and finite semisimple nonabelian groups finally the fact that the variety generated by a primal algebra coincides with the class of its subdirect powers is used this yields unitary unification algorithms for the equational theories of post algebras and p rings

localized discrete empirical interpolation method this paper presents a new approach to construct more efficient reduced order models for nonlinear partial differential equations with proper orthogonal decomposition and the discrete empirical interpolation method deim whereas deim projects the nonlinear term onto one global subspace our localized discrete empirical interpolation method ldeim computes several local subspaces each tailored to a particular region of characteristic system behavior then depending on the current state of the system ldeim selects an appropriate local subspace for the approximation of the nonlinear term in this way the dimensions of the local deim subspaces and thus the computational costs remain low even though the system might exhibit a wide range of behaviors as it passes through different regimes ldeim uses machine learning methods in the offline computational phase to discover these regions via clustering local deim approximations are then computed for each cluster in the online computational phase machine

goal oriented model constrained optimization for reduction of large scale systems optimization oriented reduced order models should target a particular output functional span an applicable range of dynamic and parametric inputs and respect the underlying governing equations of the system to achieve this goal we present an approach for determining a projection basis that uses a goal oriented model constrained optimization framework the mathematical framework permits consideration of general dynamical systems with general parametric variations and is applicable to both linear and nonlinear systems results for a simple linear model problem of the two dimensional heat equation demonstrate the ability of the goal oriented approach to target a particular output functional of interest application of the methodology to a more challenging example of a subsonic blade row governed by the unsteady euler flow equations shows a significant advantage of the new method over the proper orthogonal decomposition

local null controllability for a chemotaxis system of parabolic elliptic type abstract we consider the controllability of a chemotaxis system of parabolic elliptic type by linearizing the nonlinear system into two separate linear equations we can bypass the obstacle caused by the nonlinear drift term and establish local null controllability of the original nonlinear system this approach is different from the usual method for dealing with coupled parabolic systems

global existence of solutions for a strongly coupled quasilinear parabolic system we are interested in existence uniqueness and a priori estimates of non negative solutions to the following problem begin gathered u t delta left left 1 alpha v right u right uf u v hfill v t delta v vg u v quad in 0 infty times omega hfill u v 0 quad in 0 infty times partial omega hfill u 0 u 0 v 0 v 0 quad in omega hfill end gathered 1 0

stability of constant steady states and existence of unbounded solutions in time to a reaction diffusion equation modelling chemotaxis abstract we consider the initial boundary value problem to a reaction diffusion equation modelling chemotaxis in one space dimension under neumann boundary conditions we study the stability of constant steady states and show that there exists an unbounded solution in time induced by the instability of constant steady states

reduced order modelling of parameter dependent linear and nonlinear dynamic partial differential equation models in this paper we develop reduced order models for dynamic parameter dependent linear and nonlinear partial differential equations using proper orthogonal decomposition pod the main challenges are to accurately and efficiently approximate the pod bases for new parameter values and in the case of nonlinear problems to efficiently handle the nonlinear terms we use a bayesian nonlinear regression approach to learn the snapshots of the solutions and the nonlinearities for new parameter values computational efficiency is ensured by using manifold learning to perform the emulation in a low dimensional space the accuracy of the method is demonstrated on a linear and a nonlinear example with comparisons with a global basis approach

a survey of projection based model reduction methods for parametric dynamical systems numerical simulation of large scale dynamical systems plays a fundamental role in studying a wide range of complex physical phenomena however the inherent large scale nature of the models often leads to unmanageable demands on computational resources model reduction aims to reduce this computational burden by generating reduced models that are faster and cheaper to simulate yet accurately represent the original large scale system behavior model reduction of linear nonparametric dynamical systems has reached a considerable level of maturity as reflected by several survey papers and books however parametric model reduction has emerged only more recently as an important and vibrant research area with several recent advances making a survey paper timely thus this paper aims to provide a resource that draws together recent contributions in different communities to survey the state of the art in parametric model reduction methods parametric model reduction targets the broad class of problems for wh

global solutions in a fully parabolic chemotaxis system with singular sensitivity and nonlinear signal production a fully parabolic chemotaxis system with singular sensitivity and nonlinear signal production is considered under suitable assumption on the sensitivity value and the signal production term it is shown that this system admits a global classical solution the proof is based on a series of a priori estimate techniques a fully parabolic chemotaxis system with singular sensitivity and nonlinear signal production is considered under suitable assumption on the sensitivity value and the signal production term it is shown that this system admits a global classical solution the proof is based on a series of a priori estimate techniques

eventual smoothness and stabilization of large data solutions in a three dimensional chemotaxis system with consumption of chemoattractant abstract this paper deals with positive solutions of u t u u v x t 0 v t v u v x t 0 under homogeneous neumann boundary conditions in bounded convex domains r 3 with smooth boundary it is shown that for arbitrarily large initial data this problem admits at least one global weak solution for which there exists t 0 such that u v is bounded and smooth in t moreover it is asserted that such solutions approach spatially constant equilibria in the large time limit

singularity patterns in a chemotaxis model the authors study a chemotactic model under certain assumptions and obtain the existence of a class of solutions which blow up at the center of an open disc in finite time such a finite time blow up of solutions implies chemotactic collapse namely concentration of species to form sporae the model studied is the limiting case of a basic chemotactic model when diffusion of the chemical approaches infinity which has the form ut u u v 0 v u 1 on r2 where is an open disc with no flux homogeneous neumann boundary conditions the initial conditions are continuous functions u x 0 u0 x 0 v x 0 v0 x 0 for x under these conditions the authors prove there exists a radially symmetric solution u r t which blows up at r 0 t t a specific description of such a solution is presented the authors also discuss the strong similarity between the chemotactic model they study and the classical stefan problem

unification in commutative semigroups abstract unification is one of the basic concepts of automated theorem proving it concerns such questions as finding solutions of finite sets of equations determining if every solution comes from a most general solution and if so determining how many most general solutions are needed to generate all solutions these solutions given in terms of substitutions are called more formally unifiers the unification type of a variery equational class of algebras is defined according to the cardinality or existence of minimal complete sets of most general unifiers of particular interest from a computational point of view are varieties of groups and semigroups so far the problem has been considered mainly for particular varieties in this paper we determine unification types for all varieties of commutative semigroups in particular we prove that for commutative semigroups the unification problem is solvable in the very strong sense that there is an algorithm which for any two finite sets 1and 2of semigroup equations produces the minimal complete set of the most general unifiers of 1over the variety of commutative semigroups generated by 2 it seems that this is the first so general decidability result in the area

energy type estimates and global solvability in a two dimensional chemotaxis haptotaxis model with remodeling of non diffusible attractant this paper deals with the coupled chemotaxis haptotaxis model ut u u v u w u 1 u w x t 0 0 v u v x t 0 wt vw w 1 w u x t 0 which was initially proposed by chaplain and lolas 2006 10 to describe the interactions between cancer cells the matrix degrading enzyme and the host tissue in a process of cancer cell invasion of tissue extracellular matrix here r2 is a bounded domain with smooth boundary and and are positive parameters as compared to previous mathematical studies the novelty here consists of allowing for positive values of reflecting processes of self remodeling of the extracellular matrix under zero flux boundary conditions it is shown that for any sufficiently smooth initial data u0 w0 satisfying first order compatibility conditions the model admits a unique global smooth solution a crucial ingredient in the proof is an energy like inequality which given t 0 yields boundedness of u t in llogl this serves as a starting point for a bootstrap argument used to derive higher regularity estimates sufficient for global extensibility of solutions

manifold learning for the emulation of spatial fields from computational models repeated evaluations of expensive computer models in applications such as design optimization and uncertainty quantification can be computationally infeasible for partial differential equation pde models the outputs of interest are often spatial fields leading to high dimensional output spaces although emulators can be used to find faithful and computationally inexpensive approximations of computer models there are few methods for handling high dimensional output spaces for gaussian process gp emulation approximations of the correlation structure and or dimensionality reduction are necessary linear dimensionality reduction will fail when the output space is not well approximated by a linear subspace of the ambient space in which it lies manifold learning can overcome the limitations of linear methods if an accurate inverse map is available in this paper we use kernel pca and diffusion maps to construct gp emulators for very high dimensional output spaces arising from pde model simulations for diffusion maps we develop a new inverse map approximation several examples are presented to demonstrate the accuracy of our approach

preventing blow up in a chemotaxis model abstract this paper deals with a chemotactic model of keller segel type the main feature of the keller segel model is the possibility of blow up of solutions in a finite time to eliminate the possibility of blow up a modified version of keller segel model is introduced the blow up control relies on the presence of a pressure function which increases faster than a logarithm for high enough cells densities for such a pressure function the solutions cannot blow up in a finite time some other conditions are introduced to ensure the global boundedness