Search Results
A Bayesian approach to parameter identification in gas networks
The inverse problem of identifying the friction coefficient in an isothermal semilinear Euler system is considered. Adopting a Bayesian approach, the goal is to identify the distribution of the quantity of interest based on a finite number of noisy measurements of the pressure at the boundaries of the domain. First well-posedness of the underlying non-linear PDE system is shown using semigroup theory, and then Lipschitz continuity of the solution operator with respect to the friction coefficient is established. Based on the Lipschitz property, well-posedness of the resulting Bayesian inverse problem for the identification of the friction coefficient is inferred. Numerical tests for scalar and distributed parameters are performed to validate the theoretical results.
A Redlich-Kister type free energy model for Li-intercalation compounds with variable lattice occupation numbers
One of the central quantities of a lithium ion intercalation compound is the open circuit potential, the voltage a battery material delivers in thermodynamic equilibrium. This voltage is related to the chemical potential of lithium in the insertion material and in general a non-linear function of the mole fraction of intercalated lithium. Experimental data shows further that it is specific for various materials. The open circuit voltage is a central ingredient for mathematical models of whole battery cells, which are used to investigate and simulate the charge and discharge behavior and to interpret experimental data on non-equilibrium processes. However, since no overall predictive theoretical method presently exists for the open circuit voltage, it is commonly fitted to experimental data. Simple polynomial fitting approaches are widely used, but they lack any thermodynamic interpretation. More recently systematically and thermodynamically motivated approaches are used to model the open circuit potential. We provide here an explicit free energy density which accounts for variable occupation numbers of Li on the intercalation lattice as well as RedlichKister-type enthalpy contributions. The derived chemical potential is validated by experimental data of Liy(Ni1/3Mn1/3Co1/3)O2 and we show that only two parameters are sufficient to obtain an overall agreement of the non-linear open circuit potential within the experimental error.
Modeling and simulations of beam stabilization in edge-emitting broad area semiconductor devices
A 2+1 dimensional PDE traveling wave model describing spatial-lateral dynamics of edge-emitting broad area semiconductor devices is considered. A numerical scheme based on a split-step Fourier method is presented and implemented on a parallel compute cluster. Simulations of the model equations are used for optimizing of existing devices with respect to the emitted beam quality, as well as for creating and testing of novel device design concepts
Bistability and hysteresis in an optically injected two-section semiconductor laser
The effect of coherent single frequency injection on two-section semiconductor lasers is studied numerically using a model based on a set of delay differential equations. The existence of bistability between different CW and non-stationary regimes of operation is demonstrated in the case of sufficiently large linewidth enhancement factors.
Impact of slippage on the morphology and stability of a dewetting rim
In this study lubrication theory is used to describe the stability and morphology of the rim that forms as a thin polymer film dewets from a hydrophobized silicon wafer. Thin film equations are derived from the governing hydrodynamic equations for the polymer to enable the systematic mathematical and numerical analysis of the properties of the solutions for different regimes of slippage and for a range of time scales. Dewetting rates and the cross sectional profiles of the evolving rims are derived for these models and compared to experimental results. Experiments also show that the rim is typically unstable in the spanwise direction and develops thicker and thinner parts that may grow into ``fingers''. Linear stability analysis as well as nonlinear numerical solutions are presented to investigate shape and growth rate of the rim instability. It is demonstrated that the difference in morphology and the rate at which the instability develops can be directly attributed to the magnitude of slippage. Finally, a derivation is given for the dominant wavelength of the bulges along the unstable rim.
Error control for the approximation of Allen-Cahn and Cahn-Hilliard equations with a logarithmic potential
A fully computable upper bound for the finite element approximation error of Allen-Cahn and Cahn-Hilliard equations with logarithmic potentials is derived. Numerical experiments show that for the sharp interface limit this bound is robust past topological changes. Modifications of the abstract results to derive quasi-optimal error estimates in different norms for lowest order finite element methods are discussed and lead to weaker conditions on the residuals under which the conditional error estimates hold.
Moment asymptotics for branching random walks in random environment
We consider the long-time behaviour of a branching random walk in random environment on the lattice Zd. The migration of particles proceeds according to simple random walk in continuous time, while the medium is given as a random potential of spatially dependent killing/branching rates. The main objects of our interest are the annealed moments m_np , i.e., the p-th moments over the medium of the n-th moment over the migration and killing/branching, of the local and global population sizes. For n = 1, this is well-understood citeGM98, as m_1 is closely connected with the parabolic Anderson model. For some special distributions, citeA00 extended this to ngeq2, but only as to the first term of the asymptotics, using (a recursive version of) a Feynman-Kac formula for m_n. In this work we derive also the second term of the asymptotics, for a much larger class of distributions. In particular, we show that m_n^p m_1^np are asymptotically equal, up to an error e^o(t). The cornerstone of our method is a direct Feynman-Kac-type formula for mn, which we establish using the spine techniques developed in citeHR1.1
Anisotropic finite element mesh adaptation via higher dimensional embedding
In this paper we provide a novel anisotropic mesh adaptation technique for adaptive finite element analysis. It is based on the concept of higher dimensional embedding, which was exploited in [1-4] to obtain an anisotropic curvature adapted mesh that fits a complex surface in ℝ3. In the context of adaptive finite element simulation, the solution (which is an unknown function ƒ: Ω ⊂; ℝd → ℝ) is sought by iteratively modifying a finite element mesh according to a mesh sizing field described via a (discrete) metric tensor field that is typically obtained through an error estimator. We proposed to use a higher dimensional embedding, Φf(x) := (x1, …, xd, s f (x1, …, xd), s ∇ f (x1, …, xd))t, instead of the mesh sizing field for the mesh adaption. This embedding contains both informations of the function ƒ itself and its gradient. An isotropic mesh in this embedded space will correspond to an anisotropic mesh in the actual space, where the mesh elements are stretched and aligned according to the features of the function ƒ. To better capture the anisotropy and gradation of the mesh, it is necessary to balance the contribution of the components in this embedding. We have properly adjusted Φf(x) for adaptive finite element analysis. To better understand and validate the proposed mesh adaptation strategy, we first provide a series of experimental tests for piecewise linear interpolation of known functions. We then applied this approach in an adaptive finite element solution of partial differential equations. Both tests are performed on two-dimensional domains in which adaptive triangular meshes are generated. We compared these results with the ones obtained by the software BAMG - a metric-based adaptive mesh generator. The errors measured in the L2 norm are comparable. Moreover, our meshes captured the anisotropy more accurately than the meshes of BAMG.
On the evolution of the empirical measure for the hard-sphere dynamics
We prove that the evolution of marginals associated to the empirical measure of a finite system of hard spheres is driven by the BBGKY hierarchical expansion. The usual hierarchy of equations for L1 measures is obtained as a corollary. We discuss the ambiguities arising in the corresponding notion of microscopic series solution to the Boltzmann-Enskog equation.
From nonlinear to linear elasticity in a coupled rate-dependent/independent system for brittle delamination
We revisit the weak, energetic-type existence results obtained in [RT15] for a system for rateindependent, brittle delamination between two visco-elastic, physically nonlinear bulk materials and explain how to rigorously extend such results to the case of visco-elastic, linearly elastic bulk materials. Our approximation result is essentially based on deducing the MOSCO-convergence of the functionals involved in the energetic formulation of the system. We apply this approximation result in two different situations at small strains: Firstly, to pass from a nonlinearly elastic to a linearly elastic, brittle model on the time-continuous level, and secondly, to pass from a time-discrete to a time-continuous model using an adhesive contact approximation of the brittle model, in combination with a vanishing, super-quadratic regularization of the bulk energy. The latter approach is beneficial if the model also accounts for the evolution of temperature.