Adaptive Regularization for Nonconvex Optimization Using Inexact Function Values and Randomly Perturbed Derivatives

A regularization algorithm allowing random noise in derivatives and inexact function values is proposed for computing approximate local critical points of any order for smooth unconstrained optimization problems. For an objective function with Lipschitz continuous $p$-th derivative and given an arbitrary optimality order $q \leq p$, it is shown that this algorithm will, in expectation, compute such a point in at most $O\left(\left(\min_{j\in\{1,\ldots,q\}}\epsilon_j\right)^{-\frac{p+1}{p-q+1}}\right)$ inexact evaluations of $f$ and its derivatives whenever $q\in\{1,2\}$, where $\epsilon_j$ is the tolerance for $j$th order accuracy. This bound becomes at most $O\left(\left(\min_{j\in\{1,\ldots,q\}}\epsilon_j\right)^{-\frac{q(p+1)}{p}}\right)$ inexact evaluations if $q>2$ and all derivatives are Lipschitz continuous. Moreover these bounds are sharp in the order of the accuracy tolerances. An extension to convexly constrained problems is also outlined.Comment: 22 page

Adaptive Regularization Algorithms with Inexact Evaluations for Nonconvex Optimization

A regularization algorithm using inexact function values and inexact derivatives is proposed and its evaluation complexity analyzed. This algorithm is applicable to unconstrained problems and to problems with inexpensive constraints (that is constraints whose evaluation and enforcement has negligible cost) under the assumption that the derivative of highest degree is $\beta$-H\"{o}lder continuous. It features a very flexible adaptive mechanism for determining the inexactness which is allowed, at each iteration, when computing objective function values and derivatives. The complexity analysis covers arbitrary optimality order and arbitrary degree of available approximate derivatives. It extends results of Cartis, Gould and Toint (2018) on the evaluation complexity to the inexact case: if a $q$th order minimizer is sought using approximations to the first $p$ derivatives, it is proved that a suitable approximate minimizer within $\epsilon$ is computed by the proposed algorithm in at most $O(\epsilon^{-\frac{p+\beta}{p-q+\beta}})$ iterations and at most $O(|\log(\epsilon)|\epsilon^{-\frac{p+\beta}{p-q+\beta}})$ approximate evaluations. An algorithmic variant, although more rigid in practice, can be proved to find such an approximate minimizer in $O(|\log(\epsilon)|+\epsilon^{-\frac{p+\beta}{p-q+\beta}})$ evaluations.While the proposed framework remains so far conceptual for high degrees and orders, it is shown to yield simple and computationally realistic inexact methods when specialized to the unconstrained and bound-constrained first- and second-order cases. The deterministic complexity results are finally extended to the stochastic context, yielding adaptive sample-size rules for subsampling methods typical of machine learning.Comment: 32 page

Dynamic energy release rate in couple-stress elasticity

This paper is concerned with energy release rate for dynamic steady state crack problems in elastic materials with microstructures. A Mode III semi-infinite crack subject to loading applied on the crack surfaces is considered. The micropolar behaviour of the material is described by the theory of couple-stress elasticity developed by Koiter. A general expression for the dynamic J-integral including both traslational and micro-rotational inertial contributions is derived, and the conservation of this integral on a path surrounding the crack tip is demonstrated

Remarks on the energy release rate for an antiplane moving crack in couple stress elasticity

This paper is concerned with the steady-state propagation of an antiplane semi-infinite crack in couple stress elastic materials. A distributed loading applied at the crack faces and moving with the same velocity of the crack tip is considered, and the influence of the loading profile variations and microstructural effects on the dynamic energy release rate is investigated. The behaviour of both energy release rate and maximum total shear stress when the crack tip speed approaches the critical speed (either that of the shear waves or that of the localised surface waves) is studied. The limit case corresponding to vanishing characteristic scale lengths is addressed both numerically and analytically by means of a comparison with classical elasticity results.Comment: 37 pages, 13 figure

The strange case of negative reflection

In this paper, we show the phenomenon of negative reflection occurring in a mechanical phononic structure, consisting of a grating of fixed inclusions embedded in a linear elastic matrix. The negative reflection is not due to the introduction of a subwavelength metastructure or materials with negative mechanical properties. Numerical analyses for out-of-plane shear waves demonstrate that there exist frequencies at which most of the incident energy is reflected at negative angles. The effect is symmetric with respect to a line that is not parallel to the normal direction to the grating structure. Simulations at different angles of incidence and computations of the energy fluxes show that negative reflection is achievable in a wide range of loading conditions

Integral identities for a semi-infinite interfacial crack in anisotropic elastic bimaterials

The focus of the article is on the analysis of a semi-infinite crack at the interface between two dissimilar anisotropic elastic materials, loaded by a general asymmetrical system of forces acting on the crack faces. Recently derived symmetric and skew-symmetric weight function matrices are introduced for both plane strain and antiplane shear cracks, and used together with the fundamental reciprocal identity (Betti formula) in order to formulate the elastic fracture problem in terms of singular integral equations relating the applied loading and the resulting crack opening. The proposed compact formulation can be used to solve many problems in linear elastic fracture mechanics (for example various classic crack problems in homogeneous and heterogeneous anisotropic media, as piezoceramics or composite materials). This formulation is also fundamental in many multifield theories, where the elastic problem is coupled with other concurrent physical phenomena.Comment: 29 pages, 4 figure

On fracture criteria for dynamic crack propagation in elastic materials with couple stresses

The focus of the article is on fracture criteria for dynamic crack propagation in elastic materials with microstructures. Steady-state propagation of a Mode III semi-infinite crack subject to loading applied on the crack surfaces is considered. The micropolar behavior of the material is described by the theory of couple-stress elasticity developed by Koiter. This constitutive model includes the characteristic lengths in bending and torsion, and thus it is able to account for the underlying microstructures of the material. Both translational and micro-rotational inertial terms are included in the balance equations, and the behavior of the solution near to the crack tip is investigated by means of an asymptotic analysis. The asymptotic fields are used to evaluate the dynamic J-integral for a couple-stress material, and the energy release rate is derived by the corresponding conservation law. The propagation stability is studied according to the energy-based Griffith criterion and the obtained results are compared to those derived by the application of the maximum total shear stress criterion.Comment: 31 pages, 6 figure

Electro-osmotic flows inside triangular microchannels

This work presents a numerical investigation of both pure electro-osmotic and combined electro-osmotic/pressure-driven flows inside triangular microchannels. A finite element analysis has been adopted to solve the governing equations for the electric potential and the velocity field, accounting for a finite thickness of the electric double layer. The influence of non-dimensional parameters such as the aspect ratio of the cross-section, the electrokinetic diameter and the ratio of the pressure force to the electric force on the flow behavior has been investigated. Numerical results point out that the velocity field is significantly influenced by the aspect ratio of the cross section and the electrokinetic diameter. More specifically, the aspect ratio plays an important role in determining the maximum volumetric flow rate, while the electrokinetic diameter is crucial to establishing the range of pressures that may be sustained by the electro-osmotic flow. Numerical results are also compared with two correlations available in the literature which enable to assess the volumetric flow rate and the pressure head for microchannels featuring a rectangular, a trapezoidal or an elliptical cross-section

Experimental Investigation on Latent Thermal Energy Storages (LTESs) Based on Pure and Copper-Foam-Loaded PCMs

In this work, a commercial paraffin PCM (RT35) characterized by a change range of the solid-liquid phase transition temperature Ts−l=29–36 °C and the low thermal conductivity λSL=0.2 W/m K is experimentally tested by submitting it to thermal charging/discharging cycles. The paraffin is contained in a case with a rectangular base and heated from the top due to electrical resistance. The aim of this research is to show the benefits that a 95% porous copper metal foam (pore density PD=20PPI) can bring to a PCM-based thermal storage system by simply loading it, due to the consequent increase in the effective thermal conductivity of the medium (λLOAD=7.03 W/m K). The experimental results highlight the positive effects of the copper foam presence, such as the heat conduction improvement throughout the system, and a significant reduction in time for the complete melting of the PCM. In addition, the experimental data highlight that in the copper-foam-loaded PCM the maximum temperature reached during the heating process is lower than 20K with respect to the test with pure PCM, imposing the same heat flux on the top (P=3.5 W/m2)