Adjoint optimization of natural convection problems: differentially heated cavity

Optimization of natural convection-driven flows may provide significant improvements to the performance of cooling devices, but a theoretical investigation of such flows has been rarely done. The present paper illustrates an efficient gradient-based optimization method for analyzing such systems. We consider numerically the natural convection-driven flow in a differentially heated cavity with three Prandtl numbers (\(Pr=0.15{-}7\)) at super-critical conditions. All results and implementations were done with the spectral element code Nek5000. The flow is analyzed using linear direct and adjoint computations about a nonlinear base flow, extracting in particular optimal initial conditions using power iteration and the solution of the full adjoint direct eigenproblem. The cost function for both temperature and velocity is based on the kinetic energy and the concept of entransy, which yields a quadratic functional. Results are presented as a function of Prandtl number, time horizons and weights between kinetic energy and entransy. In particular, it is shown that the maximum transient growth is achieved at time horizons on the order of 5 time units for all cases, whereas for larger time horizons the adjoint mode is recovered as optimal initial condition. For smaller time horizons, the influence of the weights leads either to a concentric temperature distribution or to an initial condition pattern that opposes the mean shear and grows according to the Orr mechanism. For specific cases, it could also been shown that the computation of optimal initial conditions leads to a degenerate problem, with a potential loss of symmetry. In these situations, it turns out that any initial condition lying in a specific span of the eigenfunctions will yield exactly the same transient amplification. As a consequence, the power iteration converges very slowly and fails to extract all possible optimal initial conditions. According to the authors’ knowledge, this behavior is illustrated here for the first time.     

Completely positive semidefinite rank

Mathematical Programming - An $$ntimes n$$ matrix X is called completely positive semidefinite (cpsd) if there exist $$dtimes d$$ Hermitian positive semidefinite matrices $${P_i}_{i=1}^n$$ (for...     

The effect of boundary conditions on guided wave propagation in two-dimensional models of healing bone

Guided wave propagation has recently drawn significant interest in the ultrasonic characterization of bone. In this work, we present a two-dimensional computational study of ultrasound propagation in healing bones aiming at monitoring the fracture healing process. In particular, we address the effect of fluid loading boundary conditions on the characteristics of guided wave propagation, using both time and time-frequency (t-f) signal analysis techniques, for three study cases. In the first case, the bone was assumed immersed in blood which occupied the semi-infinite spaces of the upper and lower surfaces of the plate. In the second case, the bone model was assumed to have the upper surface loaded by a 2mm thick layer of blood and the lower surface loaded by a semi-infinite fluid with properties close to those of bone marrow. The third case, involves a three-layer model in which the upper surface of the plate was again loaded by a layer of blood, whereas the lower surface was loaded by a 2mm layer of a fluid which simulated bone marrow. The callus tissue was modeled as an inhomogeneous material and fracture healing was simulated as a three-stage process. The results clearly indicate that the application of realistic boundary conditions has a significant effect on the dispersion of guided waves when compared to simplified models in which the bone's surfaces are assumed free.     

Thermal smearing and screening in a strong magnetic field for Dirac materials in comparison with the two dimensional electron liquid

We compute and compare the effects due to a uniform perpendicular magnetic field as well as temperature on the static polarization functions for monolayer graphene (MLG), associated with the Dirac point, with that for the two-dimensional electron liquid (2DEL) with the use of comprehensive numerical calculations. The relevance of our study to the Friedel oscillations for the screening of the potential for a dilute distribution of impurities is reported too. Our results show substantial differences due to screening for the 2DEL and MLG which have not been given adequate attention previously.     

Nonequilibrium thermodynamics and the optimal path to turbulence in shear flows

We determine the initial condition on the laminar-turbulent boundary closest to the laminar state using nonlinear optimization for plane Couette flow. Resorting to the general evolution criterion of nonequilibrium systems we optimize the route to the statistically steady turbulent state, i.e., the state characterized by the largest entropy production. This is the first time information from the fully turbulent state is included in the optimization procedure. We demonstrate that the optimal initial condition is localized in space for realistic flow domains.     

Congruences between the coefficients of the Tate curve via formal groups

Let be the Tate curve with canonical differential, . If the characteristic is , then the Hasse invariant, , of the pair should equal one. If , then calculation of leads to a nontrivial separable relation between the coefficients and . If or , Thakur related and via elementary methods and an identity of Ramanujan. Here, we treat uniformly all characteristics via explicit calculation of the formal group law of . Our analysis was motivated by the study of the invariant which is an infinite Witt vector generalizing the Hasse invariant.

     

The inverse conductivity problem via the calculus of functions of bounded variation

In this work, a novel approach for the solution of the inverse conductivity problem from one and multiple boundary measurements has been developed on the basis of the implication of the framework of $BV$ functions. The space of the functions of bounded variation is recommended here as the most appropriate functional space hosting the conductivity profile under reconstruction. For the numerical investigation of the inversion of the inclusion problem, we propose and implement a suitable minimization scheme of an enriched—constructed herein—functional, by exploiting the inner structure of $BV$ space. Finally, we validate and illustrate our theoretical results with numerical experiments.     

A Sharp Lp Inequality for Dyadic A1 Weights in Rn

The exact best possible range of p is determined such that anydyadic A₁ weight w on Rⁿ satisfies a reverse Hölder inequalityfor p, which depends on the dimension n and the correspondingA₁ constant of w. The proof is based on an effective linearizationof the dyadic maximal operator applied to dyadic step functions.2000 Mathematics Subject Classification 42B25.     

The Bellman functions of dyadic-like maximal operators and related inequalities

For each p>1 we precisely evaluate the main Bellman functions associated with the dyadic maximal operator on and the dyadic Carleson imbedding theorem. Actually, we do that in the more general setting of tree-like maximal operators. These provide refinements of the sharp L^p inequalities for those operators. For this we introduce an effective linearization for such maximal operators on an adequate set of functions.     

Higher-order aggregate networks in the analysis of temporal networks: path structures and centralities

Despite recent advances in the study of temporal networks, the analysis of time-stampednetwork data is still a fundamental challenge. In particular, recent studies have shownthat correlations in the ordering of links crucially alter causaltopologies of temporal networks, thus invalidating analyses based on static,time-aggregated representations of time-stamped data. These findings not only highlight animportant dimension of complexity in temporal networks, but also call for newnetwork-analytic methods suitable to analyze complex systems with time-varying topologies.Addressing this open challenge, here we introduce a novel framework for the study ofpath-based centralities in temporal networks. Studying betweenness,closeness and reach centrality, we first show than an application of these measures totime-aggregated, static representations of temporal networks yields misleading resultsabout the actual importance of nodes. To overcome this problem, we define path-basedcentralities in higher-order aggregate networks, a recently proposedgeneralization of the commonly used static representation of time-stamped data. Using dataon six empirical temporal networks, we show that the resulting higher-order measuresbetter capture the true, temporal centralities of nodes. Our resultsdemonstrate that higher-order aggregate networks constitute a powerful abstraction, withbroad perspectives for the design of new, computationally efficient data mining techniquesfor time-stamped relational data.