TY - GEN A1 - Wolf, Thomas T1 - A comparison of four approaches to the calculation of conservation laws N2 - The paper compares computational aspects of four approaches to compute conservation laws of single differential equations or systems of them, ODEs and PDEs. The only restriction, required by two of the four corresponding computer algebra programs, is that each DE has to be solvable for a leading derivative. Extra constraints may be given. Examples of new conservation laws include non-polynomial expressions, an explicit variable dependence and conservation laws involving arbitrary functions. Examples involve the following equations: Ito, Liouville, Burgers, Kadomtsev-Petviashvili, Karney-Sen-Chu-Verheest, Boussinesq, Tzetzeica, Benney. T3 - ZIB-Report - SC-99-01 Y1 - 1999 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-3896 ER - TY - GEN A1 - Sokolov, Vladimir V. A1 - Wolf, Thomas T1 - A symmetry test for quasilinear coupled systems N2 - It is well known that the following class of systems of evolution equations \begin{eqnarray} \label{nsgen} \cases{ u_{t}=u_{xx}+F(u,v,u_x,v_x),\cr v_{t}=-v_{xx}+G(u,v,u_x,v_x),\cr} \end{eqnarray} is very rich in integrable cases. The complete classification problem is very difficult. Here we consider only the most interesting (from our opinion) subclass of systems (1). Namely, we consider equations linear in all derivatives of the form \begin{eqnarray} \label{kvazgen} \cases{ u_t = u_{xx} + A_{1}(u,v) u_x + A_{2}(u,v) v_x + A_{0}(u,v)\cr v_t = - v_{xx} + B_{1}(u,v) v_x + B_{2}(u,v) u_x + B_{0}(u,v). \cr} \end{eqnarray} without any restrictions on the functions $A_{i}(u,v), B_{i}(u,v)$. T3 - ZIB-Report - SC-99-02 Y1 - 1999 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-3901 ER - TY - GEN A1 - Wolf, Thomas A1 - Brandt, Andreas A1 - Mohammadzadeh, Majid T1 - Computer algebra algorithms and routines for the computation of conservation laws and fixing of gauge in differential expressions N2 - Three different approaches for the determination of conservation laws of differential equations are presented. For three corresponding REDUCE computer algebra programs CONLAW1/2/3 the necessary subroutines are discribed. One of them simplifies general solutions of overdetermined PDE systems so that all remaining free functions and constants correspond to independent conservation laws. It determines redundant functions and constants in differential expressions and is equally useful for the determination of symmetries or the fixing of gauge freedom in differential expressions. T3 - ZIB-Report - SC-98-05 Y1 - 1998 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-3482 ER - TY - GEN A1 - Wolf, Thomas A1 - Efimovskaya, Olga T1 - On integrability of the Kontsevich non-abelian ODE system N2 - We consider systems of ODEs with the right hand side being Laurent polynomials in several non-commutative unknowns. In particular, these unknowns could be matrices of arbitrary size. An important example of such a system was proposed by M. Kontsevich. We prove the integrability of the Kontsevich system by finding a Lax pair, corresponding first integrals and commuting flows. We also provide a pre-Hamiltonian operator which maps gradients of integrals for the Kontsevich system to symmetries. T3 - ZIB-Report - 11-24 KW - integrability, Lax pairs, noncommutative ODE, Laurent ODE Y1 - 2011 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-13113 ER - TY - GEN A1 - Wolf, Thomas A1 - Schrüfer, Eberhard A1 - Webster, Kenneth T1 - Solving large linear algebraic systems in the context of integrable non-abelian Laurent ODEs N2 - The paper reports on a computer algebra program {\sc LSSS} (Linear Selective Systems Solver) for solving linear algebraic systems with rational coefficients. The program is especially efficient for very large sparse systems that have a solution in which many variables take the value zero. The program is applied to the symmetry investigation of a non-abelian Laurent ODE introduced recently by M.\ Kontsevich. The computed symmetries confirmed that a Lax pair found for this system earlier generates all first integrals of degree at least up to 14. T3 - ZIB-Report - 12-32 KW - linear algebraic systems, integrability, Lax pairs, noncommutative ODE, Laurent ODE Y1 - 2012 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-16102 SN - 1438-0064 ER - TY - GEN A1 - Vierhaus, Ingmar A1 - Fügenschuh, Armin A1 - Gottwald, Robert Lion A1 - Grösser, Stefan N. T1 - Modern Nonlinear Optimization Techniques for an Optimal Control of System Dynamics Models N2 - We study System Dynamics models with several free parameters that can be altered by the user. We assume that the user's goal is to achieve a certain dynamic behavior of the model by varying these parameters. In order to the find best possible combination of parameter settings, several automatic parameter tuning methods are described in the literature and readily available within existing System Dynamic software packages. We give a survey on the available techniques in the market and describe their theoretical background. Some of these methods are already six decades old, and meanwhile newer and more powerful optimization methods have emerged in the mathematical literature. One major obstacle for their direct use are tabled data in System Dynamics models, which are usually interpreted as piecewise linear functions. However, modern optimization methods usually require smooth functions which are twice continuously differentiable. We overcome this problem by a smooth spline interpolation of the tabled data. We use a test set of three complex System Dynamic models from the literature, describe their individual transition into optimization problems, and demonstrate the applicability of modern optimization algorithms to these System Dynamics Optimization problems. T3 - ZIB-Report - 14-08 KW - System Dynamics KW - Optimal Control KW - Nonlinear Optimization KW - Spline Interpolation Y1 - 2014 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-48159 SN - 1438-0064 ER - TY - GEN A1 - Fügenschuh, Armin A1 - Vierhaus, Ingmar T1 - A Global Approach to the Optimal Control of System Dynamics Models N2 - The System Dynamics (SD) methodology is a framework for modeling and simulating the dynamic behavior of socioeconomic systems. Characteristic for the description of such systems is the occurrence of feedback loops together with stocks and flows. The mathematical equations that describe the system are usually ordinary differential equations and nonlinear algebraic constraints. Therefore seemingly simple systems can show a nonintuitive, unpredictable behavior over time. Controlling a dynamical system means to specify potential interventions from outside that should keep the system on the desired track, and to define an evaluation schema to compare different controls among each other, so that a "best" control can be defined in a meaningful way. The central question is how to compute such globally optimal control for a given SD model, that allows the transition of the system into a desired state with minimum effort. We propose a mixed-integer nonlinear programming (MINLP) reformulation of the System Dynamics Optimization (SDO) problem. MINLP problems can be solved by linear programming based branch-and-bound approach. We demonstrate that standard MINLP solvers are not able to solve SDO problem. To overcome this obstacle, we introduce a special-tailored bound propagation method. We apply our new method to a predator-prey model with additional hunting activity as control, and to a mini-world model with the consumption level as control. Numerical results for these test cases are presented. T3 - ZIB-Report - 13-28 KW - System Dynamics KW - Global Optimal Control KW - Mixed-Integer Nonlinear Optimization KW - Bounds Strengthening Y1 - 2013 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-18600 SN - 1438-0064 ER - TY - GEN A1 - Moualeu-Ngangue, Dany Pascal A1 - Röblitz, Susanna A1 - Ehrig, Rainald A1 - Deuflhard, Peter T1 - Parameter Identification for a Tuberculosis Model in Cameroon N2 - A deterministic model of tuberculosis in sub-Saharan Africa in general and Cameroon in particular including lack of access to the treatment and weak diagnose capacity is designed and analyzed with respect to its transmission dynamics. The model includes both frequency- and density-dependent transmissions. It is shown that the model is mathematically well-posed and epidemiologically reasonable. Solutions are non-negative and bounded whenever the initial values are non-negative. A sensitivity analysis of model parameters is performed and most sensitive parameters of the model are identified using a state-of-the-art Gauss-Newton Method. In particular, parameters representing the proportion of individuals having access to medical facilities have a large impact on the dynamics of the disease. It has been shown that an increase of these parameter values over the time can significantly reduce the disease burden in the population within the next 15 years. T3 - ZIB-Report - 13-72 KW - Tuberculosis, parameter identification, sensitivity analysis, Gauss-Newton method, Nonlinear dynamical systems. Y1 - 2013 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-43087 SN - 1438-0064 ER - TY - GEN A1 - Moualeu-Ngangue, Dany Pascal A1 - Weiser, Martin A1 - Ehrig, Rainald A1 - Deuflhard, Peter T1 - Optimal control for a tuberculosis model with undetected cases in Cameroon N2 - This paper considers the optimal control of tuberculosis through education, diagnosis campaign and chemoprophylaxis of latently infected. A mathematical model which includes important components such as undiagnosed infectious, diagnosed infectious, latently infected and lost-sight infectious is formulated. The model combines a frequency dependent and a density dependent force of infection for TB transmission. Through optimal control theory and numerical simulations, a cost-effective balance of two different intervention methods is obtained. Seeking to minimize the amount of money the government spends when tuberculosis remain endemic in the Cameroonian population, Pontryagin's maximum principle is used to characterize the optimal control. The optimality system is derived and solved numerically using the forward-backward sweep method (FBSM). Results provide a framework for designing cost-effective strategies for diseases with multiple intervention methods. It comes out that combining chemoprophylaxis and education, the burden of TB can be reduced by 80 % in 10 years T3 - ZIB-Report - 13-73 KW - Tuberculosis KW - Optimal control KW - Nonlinear dynamical systems Y1 - 2013 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-43142 SN - 1438-0064 ER - TY - GEN A1 - Berg, Mascha A1 - Plöntzke, Julia A1 - Leonhard-Marek, Sabine A1 - Müller, Kerstin-Elisabeth A1 - Röblitz, Susanna T1 - A dynamic model to simulate potassium balance in dairy cows. N2 - High performing dairy cows require a particular composition of nutritional ingredients, adapted to their individual requirements and depending on their production status. The optimal dimensioning of minerals in the diet, one of them being potassium, is indispensable for the prevention of imbalances. The potassium balance in cows is the result of potassium intake, distribution in the organism, and excretion, it is closely related with the glucose and electrolyte metabolism. In this paper, we present a dynamical model for the potassium balance in lactating and non-lactating dairy cows based on ordinary differential equations. Parameter values are obtained from clinical trial data and from the literature. To verify the consistency of the model, we present simulation outcomes for three different scenarios: potassium balance in (i) non-lactating cows with varying feed intake, (ii) non-lactating cows with varying potassium fraction in the diet, and (iii) lactating cows with varying milk production levels. The results give insights into the short and long term potassium metabolism, providing an important step towards the understanding of the potassium network, the design of prophylactic feed additives, and possible treatment strategies. T3 - ZIB-Report - 17-47 KW - cow, potassium, model, ode Y1 - 2017 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-64756 SN - 1438-0064 ER - TY - THES A1 - Moualeu-Ngangue, Dany Pascal T1 - A Mathematical Tuberculosis Model in Cameroon N2 - This thesis firstly presents a nonlinear extended deterministic model for the transmission dynamics of tuberculosis, based on realistic assumptions and data collected from the WHO. This model enables a comprehensive qualitative analysis of various aspects in the outbreak and control of tuberculosis in Sub-Saharan Africa countries and successfully reproduces the epidemiology of tuberculosis in Cameroon for the past (from 1994-2010). Some particular properties of the model and its solution have been presented using the comparison theorem applied to the theory of differential equations. The existence and the stability of a disease free equilibrium has been discussed using the Perron-Frobenius theorem and Metzler stable matrices. Furthermore, we computed the basic reproduction number, i.e. the number of cases that one case generates on average over the course of its infectious period. Rigorous qualitative analysis of the model reveals that, in contrast to the model without reinfections, the full model with reinfection exhibits the phenomenon of backward bifurcation, where a stable disease-free equilibrium coexists with a stable endemic equilibrium when a certain threshold quantity, known as the basic reproduction ratio (R0), is less than unity. The global stability of the disease-free equilibrium has been discussed using the concepts of Lyapunov stability and bifurcation theory. With the help of a sensitivity analysis using data of Cameroon, we identified the relevant parameters which play a key role for the transmission and the control of the disease. This was possible applying sophisticated numerical methods (POEM) developed at ZIB. Using advanced approaches for optimal control considering the costs for chemoprophylaxis, treatment and educational campaigns should provide a framework for designing realistic cost effective strategies with different intervention methods. The forward-backward sweep method has been used to solve the numerical optimal control problem. The numerical result of the optimal control problem reveals that combined effort in education and chemoprophylaxis may lead to a reduction of 80\% in the number of infected people in 10 years. The mathematical and numerical approaches developed in this thesis could be similarly applied in many other Sub-Saharan countries where TB is a public health problem. Y1 - 2013 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:188-fudissthesis000000095327-7 UR - http://www.diss.fu-berlin.de/diss/receive/FUDISS_thesis_000000095327 ER - TY - GEN A1 - Anco, Stephen A1 - Ali, Sajid A1 - Wolf, Thomas T1 - Symmetry analysis and exact solutions of semilinear heat flow in multi-dimensions N2 - A symmetry group method is used to obtain exact solutions for a semilinear radial heat equation in $n>1$ dimensions with a general power nonlinearity. The method involves an ansatz technique to solve an equivalent first-order PDE system of similarity variables given by group foliations of this heat equation, using its admitted group of scaling symmetries. This technique yields explicit similarity solutions as well as other explicit solutions of a more general (non-similarity) form having interesting analytical behavior connected with blow up and dispersion. In contrast, standard similarity reduction of this heat equation gives a semilinear ODE that cannot be explicitly solved by familiar integration techniques such as point symmetry reduction or integrating factors. T3 - ZIB-Report - 11-26 KW - semilinear heat equation, similarity reduction, exact solutions, group foliation, symmetry Y1 - 2011 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-13213 ER - TY - GEN A1 - Anco, Stephen A1 - Ali, Sajid A1 - Wolf, Thomas T1 - Exact solutions of nonlinear partial differential equations by the method of group foliation reduction N2 - A novel symmetry method for finding exact solutions to nonlinear PDEs is illustrated by applying it to a semilinear reaction-diffusion equation in multi-dimensions. The method is based on group foliation reduction and employs a separation ansatz to solve an equivalent first-order group foliation system whose independent and dependent variables respectively consist of the invariants and differential invariants of a given one-dimensional group of point symmetries for the reaction-diffusion equation. With this method, solutions of the reaction-diffusion equation are obtained in an explicit form, including group-invariant similarity solutions and travelling-wave solutions, as well as dynamically interesting solutions that are not invariant under any of the point symmetries admitted by this equation. T3 - ZIB-Report - 11-25 KW - semilinear heat equation, similarity reduction, exact solutions, group foliation, symmetry Y1 - 2011 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-13171 ER -