Stationary solutions of liquid two-layer thin film models

We investigate stationary solutions of a thin-film model for liquid two-layer flows in an energetic formulation that is motivated by its gradient flow structure. The goal is to achieve a rigorous understanding of the contact-angle conditions for such two-layer systems. We pursue this by investigating a corresponding energy that favors the upper liquid to dewet from the lower liquid substrate, leaving behind a layer of thickness $h_*$. After proving existence of stationary solutions for the resulting system of thin-film equations we focus on the limit $h_*\to 0$ via matched asymptotic analysis. This yields a corresponding sharp-interface model and a matched asymptotic solution that includes logarithmic switch-back terms. We compare this with results obtained using $\Gamma$-convergence, where we establish existence and uniqueness of energetic minimizers in that limit

Assessment of the potential impacts of plant traits across environments by combining global sensitivity analysis and dynamic modeling in wheat

A crop can be viewed as a complex system with outputs (e.g. yield) that are affected by inputs of genetic, physiology, pedo-climatic and management information. Application of numerical methods for model exploration assist in evaluating the major most influential inputs, providing the simulation model is a credible description of the biological system. A sensitivity analysis was used to assess the simulated impact on yield of a suite of traits involved in major processes of crop growth and development, and to evaluate how the simulated value of such traits varies across environments and in relation to other traits (which can be interpreted as a virtual change in genetic background). The study focused on wheat in Australia, with an emphasis on adaptation to low rainfall conditions. A large set of traits (90) was evaluated in a wide target population of environments (4 sites x 125 years), management practices (3 sowing dates x 2 N fertilization) and $CO_2$ (2 levels). The Morris sensitivity analysis method was used to sample the parameter space and reduce computational requirements, while maintaining a realistic representation of the targeted trait x environment x management landscape ($\sim$ 82 million individual simulations in total). The patterns of parameter x environment x management interactions were investigated for the most influential parameters, considering a potential genetic range of +/- 20% compared to a reference. Main (i.e. linear) and interaction (i.e. non-linear and interaction) sensitivity indices calculated for most of APSIM-Wheat parameters allowed the identifcation of 42 parameters substantially impacting yield in most target environments. Among these, a subset of parameters related to phenology, resource acquisition, resource use efficiency and biomass allocation were identified as potential candidates for crop (and model) improvement.Comment: 22 pages, 8 figures. This work has been submitted to PLoS On

Numerical convergence for semilinear parabolic equations

We present a convergence result for finite element discretisations of semilinear parabolic equations, in which the evaluation of the nonlinearity requires some high order of regularity of the solution. For example a coefficient might depend on derivatives or pointevaluation of the solution. We do not rely on high regularity of the exact solution itself and as a payoff we can not deduce convergence rates. As an example the convergence result is applied to a nonlinear Fokker--Planck type battery model

Evaluating and predicting agricultural management effects under tile drainage using modified APSIM

An accurate and management sensitive simulation model for tile-drained Midwestern soils is needed to optimize the use of agricultural management practices (e.g., winter cover crops) to reduce nitrate leaching without adversely affecting corn yield. Our objectives were to enhance the Agricultural Production Systems Simulator (APSIM) for tile drainage, test the modified model for several management scenarios, and then predict nitrate leaching with and without winter wheat cover crop. Twelve years of data (1990–2001) from northeast Iowa were used for model testing. Management scenarios included continuous corn and corn–soybean rotations with single or split N applications. For 38 of 44 observations, yearly drain flow was simulated within 50 mm of observed for low drainage (\u3c 100 mm) or within 30% of observed for high drain flow. Corn yield was simulated within 1500 kg/ha for 12 of 24 observations. For 30 of 45 observations yearly nitrate-N loss in tile drains was simulated within 10 kg N/ha for low nitrate-N loss (\u3c 20 kg N/ha) or within 30% of observed for high nitrate-N loss. Several of the poor yield and nitrate-N loss predictions appear related to poor N-uptake simulations. The model accurately predicted greater corn yield under split application (140–190 kg N/ha) compared to single 110 kg N/ha application and higher drainage and nitrate-N loss under continuous corn compared to corn/soybean rotations. A winter wheat cover crop was predicted to reduce nitrate-N loss 38% (341 vs. 537 kg N/ha with and without cover) under 41-years of corn-soybean rotations and 150 kg N/ha applied to corn. These results suggest that the modified APSIM model is a promising tool to help estimate the relative effect of alternative management practices under fluctuating high water tables

Hysteresis in the context of hydrogen storage and lithium-ion batteries

The processes of reversible storage of hydrogen in a metal by loading and unloading and of charging and discharging of lithium-ion batteries have many things in common. The both processes are accompanied by a phase transition and loading and unloading run along different paths, so that hysteretic behavior is observed. For hydrogen storage we consider a fine powder of magnesium (Mg) particles and lithium storage is studied for iron phosphate (FePO$_4$) particles forming the cathode of a lithium-ion battery. The mathematical models that are established in \cite{DGJ08} and \cite{DGH09a}, describe phase transitions and hysteresis exclusively in a single particle and on that basis they can predict the observed hysteretic plots with almost horizontal plateaus. Interestingly the models predict that the coexistence of a 2-phase system in an individual particle disappears, if its size is below a critical value. However, measurements reveal that this is qualitatively not reflected by the mentioned hysteretic plots of loading and unloading. In other words: The behavior of a storage system consisting of many particles is qualitatively independent of the fact whether the individual particles itself develop a 2-phase system or if they remain in a single phase state. This apparent paradoxical observation will be resolved in this article. It will be shown that if each of the individual particles homogeneously distributes the supplied matter, nevertheless the many particle ensemble exhibits phase transition and hysteresis, because one of the two phases is realized in some part of the particles while the remaining part is in the other phase

The behavior of a many particle cathode in a lithium-ion battery

We study the almost reversible storage process of charging and discharging of lithium-ion batteries. That process is accompanied by a phase transition and charging and discharging run along different paths, so that hysteretic behavior is observed. We are interested in the storage problem of the cathode of a lithium-ion battery consisting of a system of many iron phosphate (FePO4) particles. There are mathematical models, see [DGJ08], [DGGHJ09] and [DG09], that describe phase transitions and hysteresis exclusively in a single storage particle and they can describe the observed hysteretic voltage-charge plots with almost horizontal plateaus. Interestingly the models predict that the coexistence of a 2-phase system in an individual particle disappears, if its size is below a critical value. The disappearance of the phase transition in the single particle model implies the disappearance of the hysteresis. However, in the experiment hysteretic behavior survives. In other words: The behavior of a storage system consisting of many particles is qualitatively independent of the fact whether the individual particles itself develop a 2-phase system or if they remain in a single phase state. This apparent paradoxical observation will be resolved in this article by a many particle model. It will be shown that if each of the individual particles is in a homogeneous state, nevertheless the many particle ensemble exhibits phase transition and hysteresis ..

Report of the Dark Energy Task Force

Dark energy appears to be the dominant component of the physical Universe, yet there is no persuasive theoretical explanation for its existence or magnitude. The acceleration of the Universe is, along with dark matter, the observed phenomenon that most directly demonstrates that our theories of fundamental particles and gravity are either incorrect or incomplete. Most experts believe that nothing short of a revolution in our understanding of fundamental physics will be required to achieve a full understanding of the cosmic acceleration. For these reasons, the nature of dark energy ranks among the very most compelling of all outstanding problems in physical science. These circumstances demand an ambitious observational program to determine the dark energy properties as well as possible

Gradient flow perspective of thin-film bilayer flows

We study gradient flow formulations of thin-film bilayer flows with triple-junctions between liquid/liquid/air. First we highlight the gradient structure in the Stokes free-boundary flow and identify its solutions with the well known PDE with boundary conditions. Next we propose a similar gradient formulation for the corresponding thin-film model and formally identify solutions with those of the corresponding free-boundary problem. A robust numerical algorithm for the thin-film gradient flow structure is then provided. Using this algorithm we compare the sharp triple-junction model with precursor models. For their stationary solutions a rigorous connection is established using Gamma-convergence. For time-dependent solutions the comparison of numerical solutions shows a good agreement for small and moderate times. Finally we study spreading in the zero-contact angle case, where we compare numerical solutions with asymptotically exact source-type solutions

Gradient flow perspective of thin-film bilayer flows

We study gradient flow formulations of thin-film bilayer flows with triple-junctions between liquid/liquid/air. First we highlight the gradient structure in the Stokes free-boundary flow and identify its solutions with the well-known PDE with boundary conditions. Next we propose a similar gradient formulation for the corresponding thin-film model and formally identify solutions with those of the corresponding free-boundary problem. A robust numerical algorithm for the thin-film gradient flow structure is then provided. Using this algorithm we compare the sharp triple-junction model with precursor models. For their stationary solutions a rigorous connection is established using [Gamma]-convergence. For time-dependentsolutions the comparison of numerical solutions shows a good agreement for small and moderate times. Finally we study spreading in the zero-contact angle case, where we compare numerical solutions with asymptotically exact source-type solutions

Blow-up versus boundedness in a nonlocal and nonlinear Fokker-Planck equation

Literaturverz