An Optimal Order Interior Penalty Discontinuous Galerkin Discretization of the Compressible Navier-Stokes Equations

In this article we propose a new symmetric version of the interior penalty discontinuous Galerkin finite element method for the numerical approximation of the compressible Navier-Stokes equations. Here, particular emphasis is devoted to the construction of an optimal numerical method for the evaluation of certain target functionals of practical interest, such as the lift and drag coefficients of a body immersed in a viscous fluid. With this in mind, the key ingredients in the construction of the method include: (i) An adjoint consistent imposition of the boundary conditions; (ii) An adjoint consistent reformulation of the underlying target functional of practical interest; (iii) Design of appropriate interior-penalty stabilization terms. Numerical experiments presented within this article clearly indicate the optimality of the proposed method when the error is measured in terms of both the L_2-norm, as well as for certain target functionals. Computational comparisons with other discontinuous Galerkin schemes proposed in the literature, including the second scheme of Bassi & Rebay, cf. [11], the standard SIPG method outlined in [25], and an NIPG variant of the new scheme will be undertaken.

Estimation of whole body muscle, adipose/fat mass, validation in health and during weight loss. Development of prediction equations using MRI as reference method

Background: Body composition is affected by diseases, and affects responses to medical treatments, dosage of medicines, etc., while an abnormal body composition contributes to the causation of many chronic diseases. While we have reliable biochemical tests for certain nutritional parameters of body composition, such as iron or iodine status, and we have harnessed nuclear physics to estimate the body’s content of trace elements, the very basic quantification of body fat content and muscle mass remains highly problematic. Both body fat and muscle mass are vitally important, as they have opposing influences on chronic disease, but they have seldom been estimated as part of population health surveillance. Instead, most national surveys have merely reported BMI and waist, or sometimes the waist/hip ratio; these indices are convenient but do not have any specific biological meaning. Anthropometry offers a practical and inexpensive method for muscle and fat estimation in clinical and epidemiological settings; however, its use is imperfect due to many limitations, such as a shortage of reference data, misuse of terminology, unclear assumptions, and the absence of properly validated anthropometric equations. To date, anthropometric methods are not sensitive enough to detect muscle and fat loss. Aims: The aim of this thesis is to estimate Adipose/fat and muscle mass in health disease and during weight loss through; 1. evaluating and critiquing the literature, to identify the best-published prediction equations for adipose/fat and muscle mass estimation; 2. to derive and validate adipose tissue and muscle mass prediction equations; and 3.to evaluate the prediction equations along with anthropometric indices and the best equations retrieved from the literature in health, metabolic illness and during weight loss. Methods: a Systematic review using Cochrane Review method was used for reviewing muscle mass estimation papers that used MRI as the reference method. Fat mass estimation papers were critically reviewed. Mixed ethnic, age and body mass data that underwent whole body magnetic resonance imaging to quantify adipose tissue and muscle mass (dependent variable) and anthropometry (independent variable) were used in the derivation/validation analysis. Multiple regression and Bland-Altman plot were applied to evaluate the prediction equations. To determine how well the equations identify metabolic illness, English and Scottish health surveys were studied. Statistical analysis using multiple regression and binary logistic regression were applied to assess model fit and associations. Also, populations were divided into quintiles and relative risk was analysed. Finally, the prediction equations were evaluated by applying them to a pilot study of 10 subjects who underwent whole-body MRI, anthropometric measurements and muscle strength before and after weight loss to determine how well the equations identify adipose/fat mass and muscle mass change. Results: The estimation of fat mass has serious problems. Despite advances in technology and science, prediction equations for the estimation of fat mass depend on limited historical reference data and remain dependent upon assumptions that have not yet been properly validated for different population groups. Muscle mass does not have the same conceptual problems; however, its measurement is still problematic and reference data are scarce. The derivation and validation analysis in this thesis was satisfactory, compared to prediction equations in the literature they were similar or even better. Applying the prediction equations in metabolic illness and during weight loss presented an understanding on how well the equations identify metabolic illness showing significant associations with diabetes, hypertension, HbA1c and blood pressure. And moderate to high correlations with MRI-measured adipose tissue and muscle mass before and after weight loss. Conclusion: Adipose tissue mass and to an extent muscle mass can now be estimated for many purposes as population or groups means. However, these equations must not be used for assessing fatness and categorising individuals. Further exploration in different populations and health surveys would be valuable.

Fractional differential equations with dependence on the Caputo-Katugampola derivative

In this paper we present a new type of fractional operator, the Caputo–Katugampola derivative. The Caputo and the Caputo–Hadamard fractional derivatives are special cases of this new operator. An existence and uniqueness theorem for a fractional Cauchy type problem, with dependence on the Caputo–Katugampola derivative, is proven. A decomposition formula for the Caputo–Katugampola derivative is obtained. This formula allows us to provide a simple numerical procedure to solve the fractional differential equation.

Modeling some real phenomena by fractional differential equations

This paper deals with fractional differential equations, with dependence on a Caputo fractional derivative of real order. The goal is to show, based on concrete examples and experimental data from several experiments, that fractional differential equations may model more efficiently certain problems than ordinary differential equations. A numerical optimization approach based on least squares approximation is used to determine the order of the fractional operator that better describes real data, as well as other related parameters.

Solutions of Tikhonov functional equations and applications to multiplication operators on Szegö spaces

We consider a natural representation of solutions for Tikhonov functional equations. This will be done by applying the theory of reproducing kernels to the approximate solutions of general bounded linear operator equations (when defined from reproducing kernel Hilbert spaces into general Hilbert spaces), by using the Hilbert-Schmidt property and tensor product of Hilbert spaces. As a concrete case, we shall consider generalized fractional functions formed by the quotient of Bergman functions by Szegö functions considered from the multiplication operators on the Szegö spaces.

Finite time extinction for nonlinear fractional evolution equations and related properties

The finite time extinction phenomenon (the solution reaches an equilibrium after a finite time) is peculiar to certain nonlinear problems whose solutions exhibit an asymptotic behavior entirely different from the typical behavior of solutions associated to linear problems. The main goal of this work is twofold. Firstly, we extend some of the results known in the literature to the case in which the ordinary time derivative is considered jointly with a fractional time differentiation. Secondly, we consider the limit case when only the fractional derivative remains. The latter is the most extraordinary case, since we prove that the finite time extinction phenomenon still appears, even with a non-smooth profile near the extinction time. Some concrete examples of quasi-linear partial differential operators are proposed. Our results can also be applied in the framework of suitable nonlinear Volterra integro-differential equations.

Finite time extinction for nonlinear fractional evolution equations and related properties.

The finite time extinction phenomenon (the solution reaches an equilibrium after a finite time) is peculiar to certain nonlinear problems whose solutions exhibit an asymptotic behavior entirely different from the typical behavior of solutions associated to linear problems. The main goal of this work is twofold. Firstly, we extend some of the results known in the literature to the case in which the ordinary time derivative is considered jointly with a fractional time differentiation. Secondly, we consider the limit case when only the fractional derivative remains. The latter is the most extraordinary case, since we prove that the finite time extinction phenomenon still appears, even with a non-smooth profile near the extinction time. Some concrete examples of quasi-linear partial differential operators are proposed. Our results can also be applied in the framework of suitable nonlinear Volterra integro-differential equations.

The Becker-Döring equations with exponentially size-dependent rate coefficients

This paper is concerned with an analysis of the Becker-Döring equations which lie at the heart of a number of descriptions of non-equilibrium phase transitions and related complex dynamical processes. The Becker-Döring theory describes growth and fragmentation in terms of stepwise addition or removal of single particles to or from clusters of similar particles and has been applied to a wide range of problems of physicochemical and biological interest within recent years. Here we consider the case where the aggregation and fragmentation rates depend exponentially on cluster size. These choices of rate coefficients at least qualitatively correspond to physically realistic molecular clustering scenarios such as occur in, for example, simulations of simple fluids. New similarity solutions for the constant monomer Becker-Döring system are identified, and shown to be generic in the case of aggregation and fragmentation rates that depend exponentially on cluster size.

Coagulation equations with mass loss

We derive and solve models for coagulation with mass loss arising, for example, from industrial processes in which growing inclusions are lost from the melt by colliding with the wall of the vessel. We consider a variety of loss laws and a variety of coagulation kernels, deriving exact results where possible, and more generally reducing the equations to similarity solutions valid in the large-time limit. One notable result is the effect that mass removal has on gelation: for small loss rates, gelation is delayed, whilst above a critical threshold, gelation is completely prevented. Finally, by forming an exact explicit solution for a more general initial cluster size distribution function, we illustrate how numerical results from earlier work can be interpreted in the light of the theory presented herein.

The Becker-Döring equations with monomer input, competition and inhibition

We investigate the Becker-Döring model of nucleation with three generalisations; an input of monomer, an input of inhibitor and finally, we allow the monomers to form two morphologies of cluster. We assume size-independent aggregation and fragmentation rates. Initially we consider the problem of constant monomer input and determine the steady-state solution approached in the large-time limit, and the manner in which it is approached. Secondly, in addition to a constant input of monomer we allow a constant input of inhibitor, which prevents clusters growing any larger and this removes them from the kinetics of the process; the inhibitor is consumed in the action of poisoning a cluster. We determine a critical ratio of poison to monomer input below which the cluster concentrations tend to a non-zero steady-state solution and the poison concentration tends to a finite value. Above the critical input ratio, the concentrations of all cluster sizes tend to zero and the poison concentration grows without limit. In both cases the solution in the large-time limit is determined. Finally we consider a model where monomers form two morphologies, but the inhibitor only acts on one morphology. Four cases are identified, depending on the relative poison to monomer input rates and the relative thermodynamic stability. In each case we determine the final cluster distribution and poison concentration. We find that poisoning the less stable cluster type can have a significant impact on the structure of the more stable cluster distribution; a counter-intuitive result. All results are shown to agree with numerical simulation.

Instabilities in threshold-diffusion equations with delay

The introduction of delays into ordinary or partial differential equation models is well known to facilitate the production of rich dynamics ranging from periodic solutions through to spatio-temporal chaos. In this paper we consider a class of scalar partial differential equations with a delayed threshold nonlinearity which admits exact solutions for equilibria, periodic orbits and travelling waves. Importantly we show how the spectra of periodic and travelling wave solutions can be determined in terms of the zeros of a complex analytic function. Using this as a computational tool to determine stability we show that delays can have very different effects on threshold systems with negative as opposed to positive feedback. Direct numerical simulations are used to confirm our bifurcation analysis, and to probe some of the rich behaviour possible for mixed feedback.

Elmar Schwinger: Literarische Erziehung und Gymnasium. Zur Entwicklung des bayerischen Gymnasiums in der Ära Niethammer/Thiersch. Mit einem Geleitwort von A. Reble. (Würzburger Arbeiten zur Erziehungswissenschaft. Bd. 18.) Bad Heilbrunn: Klinkhardt 1988 [Rezension]

Rezension von: Elmar Schwinger: Literarische Erziehung und Gymnasium. Zur Entwicklung des bayerischen Gymnasiums in der Ära Niethammer/Thiersch. Mit einem Geleitwort von A. Reble. (Würzburger Arbeiten zur Erziehungswissenschaft. Bd. 18.) Bad Heilbrunn: Klinkhardt 1988. 452 S.

"The X in balance", a serious game to learn how to solve mathematical equations

(Prefácio) This dissertation is submitted for the degree of Masters (Engenharia Informática) at University of Évora. Under the supervision of Professor Francisco Manuel Gonçalves Coelho, i have selected to work on game design. With the specific period of time and resources, an attempt has been made to make a serious educational game. While writing this thesis, the objective was to describe a math game for solving mathematical equations. Injecting learning factor in a game, is a main concern of this project. The document is about the description of ‘X in Balance’ game. This game provides a platform for school aged students to solve the equations by playing game. It also gives a unique dimension of putting fun and math in a same platform. The document describes full detail on the project. The first chapter gives an introduction about the problem faced by students in doing maths and the learning behavior of a game. It also points out the opportunities that this game might brings and the motivation behind doing this work. It describes the game concept and its genre too. Besides, the second chapter tells state of an art of serous educational game. It defines the concept of serious game and its types. Furthermore, it justifies the flexibility of serious games to adapt all learning styles. The impact of serious games on learning is also mentioned. It also includes the related work of other researchers.

POSITIVE FILTERED P_N METHOD FOR LINEAR TRANSPORT EQUATIONS AND THE ASSOCIATED OPTIMIZATION ALGORITHM

We propose a positive, accurate moment closure for linear kinetic transport equations based on a filtered spherical harmonic (FP_N) expansion in the angular variable. The FP_N moment equations are accurate approximations to linear kinetic equations, but they are known to suffer from the occurrence of unphysical, negative particle concentrations. The new positive filtered P_N (FP_N+) closure is developed to address this issue. The FP_N+ closure approximates the kinetic distribution by a spherical harmonic expansion that is non-negative on a finite, predetermined set of quadrature points. With an appropriate numerical PDE solver, the FP_N+ closure generates particle concentrations that are guaranteed to be non-negative. Under an additional, mild regularity assumption, we prove that as the moment order tends to infinity, the FP_N+ approximation converges, in the L2 sense, at the same rate as the FP_N approximation; numerical tests suggest that this assumption may not be necessary. By numerical experiments on the challenging line source benchmark problem, we confirm that the FP_N+ method indeed produces accurate and non-negative solutions. To apply the FP_N+ closure on problems at large temporal-spatial scales, we develop a positive asymptotic preserving (AP) numerical PDE solver. We prove that the propose AP scheme maintains stability and accuracy with standard mesh sizes at large temporal-spatial scales, while, for generic numerical schemes, excessive refinements on temporal-spatial meshes are required. We also show that the proposed scheme preserves positivity of the particle concentration, under some time step restriction. Numerical results confirm that the proposed AP scheme is capable for solving linear transport equations at large temporal-spatial scales, for which a generic scheme could fail. Constrained optimization problems are involved in the formulation of the FP_N+ closure to enforce non-negativity of the FP_N+ approximation on the set of quadrature points. These optimization problems can be written as strictly convex quadratic programs (CQPs) with a large number of inequality constraints. To efficiently solve the CQPs, we propose a constraint-reduced variant of a Mehrotra-predictor-corrector algorithm, with a novel constraint selection rule. We prove that, under appropriate assumptions, the proposed optimization algorithm converges globally to the solution at a locally q-quadratic rate. We test the algorithm on randomly generated problems, and the numerical results indicate that the combination of the proposed algorithm and the constraint selection rule outperforms other compared constraint-reduced algorithms, especially for problems with many more inequality constraints than variables.

THE STOCHASTIC NAVIER STOKES EQUATIONS FOR HEAT CONDUCTING, COMPRESSIBLE FLUIDS

This dissertation is devoted to the equations of motion governing the evolution of a fluid or gas at the macroscopic scale. The classical model is a PDE description known as the Navier-Stokes equations. The behavior of solutions is notoriously complex, leading many in the scientific community to describe fluid mechanics using a statistical language. In the physics literature, this is often done in an ad-hoc manner with limited precision about the sense in which the randomness enters the evolution equation. The stochastic PDE community has begun proposing precise models, where a random perturbation appears explicitly in the evolution equation. Although this has been an active area of study in recent years, the existing literature is almost entirely devoted to incompressible fluids. The purpose of this thesis is to take a step forward in addressing this statistical perspective in the setting of compressible fluids. In particular, we study the well posedness for the corresponding system of Stochastic Navier Stokes equations, satisfied by the density, velocity, and temperature. The evolution of the momentum involves a random forcing which is Brownian in time and colored in space. We allow for multiplicative noise, meaning that spatial correlations may depend locally on the fluid variables. Our main result is a proof of global existence of weak martingale solutions to the Cauchy problem set within a bounded domain, emanating from large initial datum. The proof involves a mix of deterministic and stochastic analysis tools. Fundamentally, the approach is based on weak compactness techniques from the deterministic theory combined with martingale methods. Four layers of approximate stochastic PDE's are built and analyzed. A careful study of the probability laws of our approximating sequences is required. We prove appropriate tightness results and appeal to a recent generalization of the Skorohod theorem. This ultimately allows us to deduce analogues of the weak compactness tools of Lions and Feireisl, appropriately interpreted in the stochastic setting.