Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “homotopy”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

HomPINNs: Homotopy physics-informed neural networks for learning multiple solutions of nonlinear elliptic differential equations

Physics-informed neural networks (PINNs) based machine learning is an emerging framework for solving nonlinear differential equations. However, due to the implicit regularity of neural network structure, PINNs can only find the flattest solution in most cases by minimizing the loss functions. In this paper, we combine PINNs with the homotopy continuation method, a classical numerical method to compute isolated roots of polynomial systems, and propose a new deep learning framework, named homotopy physics-informed neural networks (HomPINNs), for solving multiple solutions of nonlinear elliptic differential equations. The implementation of an HomPINN is a homotopy process that is composed of the training of a fully connected neural network, named the starting neural network, and training processes of several PINNs with different tracking parameters. The starting neural network is to approximate a starting function constructed by the trivial solutions, while other PINNs are to minimize the loss functions defined by boundary condition and homotopy functions, varying with different tracking parameters. These training processes are regraded as different steps of a homotopy process, and a PINN is initialized by the well-trained neural network of the previous step, while the first starting neural network is initialized using the default initialization method. Finally, several numerical examples are presented to show the efficiency of our proposed HomPINNs, including reaction-diffusion equations with a heart-shaped domain.

97 MATHEMATICS AND COMPUTING↗

Calculating Probability Densities with Homotopy and Applications to Particle Filters

In this work, we explore a homotopy sampling procedure and its generalization, loosely based on importance sampling, known as annealed importance sampling. The procedure makes use of a known probability distribution to find, via homotopy, the unknown normalization of a target distribution, as well as samples of the target distribution. We propose a reformulation of the method that leads to a rejection sampling alternative. Estimates of the error as a function of homotopy stages and sample averages are derived for the algorithmic version of the method. These estimates are useful in making computational efficiency decisions on how the calculation should proceed, given a computer architecture. Consideration is given to how the procedure can be adapted to Bayesian stationary and nonstationary estimation problems. The connection between homotopy sampling and thermodynamic integration is made. Emphasis is placed on the non-stationary problems, and in particular, on a sequential estimation technique known as particle filtering. It is shown that a modification of the particle filter framework to include the homotopy process can improve the computational robustness of particle filters.

97 MATHEMATICS AND COMPUTING↗

A modified homotopy method for solving the principal-agent bilevel programming problem

In this paper, the principal-agent bilevel programming problem with integral operator is considered, in which the upper-level object is the agent that maximizes its expected utility with respect to an agreed compensation contract. The constraints are the principal’s participation and the agent’s incentive compatibility. The latter is a lower-level optimization problem with respect to its private action. To solve an equivalent single-level nonconvex programming problem with integral operator, a modified homotopy method for solving the Karush–Kuhn–Tucker system is proposed. This method requires only an interior point and, not necessarily, a feasible initial approximation for the constraint shifting set. Global convergence is proven under some mild conditions. Numerical experiments were performed by our homotopy method as well as by fmincon in Matlab, LOQO and MINOS. The experiments showed that: designing a piecewise linear contract is much better than designing a piecewise constant contract and only needs to solve a much lower-dimensional optimization problem and hence needs much less computation time; the optimal value of the principal-agent model with designing piecewise linear contract tends to a limitation, while the discrete segments gradually increase; and finally, the proposed modified homotopy method is feasible and effective.

97 MATHEMATICS AND COMPUTING↗

Neutron star under homotopy perturbation method

Highlights: • We study isotropic neutron star under Einstein’s General Relativity. • We get mass function by solving TOV equation via homotopy perturbation method. • Using causality and stability conditions we find allowed ranges of parameters. • We provide stable, physically valid and non-singular solution set. - Abstract: We obtain a mass function solving the Tolman–Oppenheimer–Volkoff (TOV) equation for isotropic and spherically symmetric system via homotopy perturbation method (HPM). Using the mass function we construct a stellar model which can be determined from the equation of state (EOS) parameter (ω) and a model parameter (n). With the help of Einstein field equations we develop three solutions which can describe different properties and the core–crust structure of neutron star (NS). Solution I is valid for NS having the inner and outer radius near the surface of the star. The star is physical up to the inner radius whereas negative density occurs and the energy conditions are violated in the upper region from the inner to outer radius. Solution II represents NS with high gravitational redshift as well as compactification factor. All the features of NS can be given by solution III which involves only the EOS parameter. Our model predicts maximum mass for a NS with the central density 5.5×10{sup 15} g/cm{sup 3} and surface redshift 0.69 is to be 2.01 M{sub ⊙} for the EOS parameter ω=0.73. The predicted range for the surface redshift is 0.57<Z{sub s}<1.95 for the allowed ranges 8.4<n<10.9 and 1/3<ω<1 in the presented NS model.

79 ASTRONOMY AND ASTROPHYSICS↗

Tool employing homotopy-based approaches in finding the controlling unstable equilibrium point in the electric power grid

The analysis tool employs a computer-implemented algorithm that uses homotopy-based approaches to map the solution from the exit point to the controlling unstable equilibrium point (UEP). The computational time is reduced significantly by using an approximate exit point rather than computing an accurate exit point as it is required in finding the controlling UEPs using traditional transient stability direct methods. In addition, this method eliminates the necessity of computing the minimum gradient point (MGP) which is a key element in using Newton methods. These properties provide an advantage to homotopy-based approaches over traditional iterative methods in terms of both speed of computation and reliability of finding solutions.

97 MATHEMATICS AND COMPUTING↗

Approximation of the Frequency-Amplitude Curve Using the Homotopy Analysis Method

The Frequency-Amplitude (F-A) curve on power system oscillation under a large disturbance characterizes how a natural oscillation mode transitions to nonlinear oscillations with growing amplitudes and decaying frequencies. The existing formulation of the F-A curve is derived by solving elliptical integrals on oscillation of a single-machine-infinite-bus equivalent about the targeted oscillation mode. The formula is in a form of infinite series and needs to sum a large number of terms for satisfactory accuracy. This paper introduces an explicit, approximate expression obtained from the Homotopy Analysis Method on the F-A curve. The proposed F-A curve expression is derived from an SMIB system and verified on the IEEE 3-machine 9-bus system to show how the oscillation frequency of a dominant mode varies with oscillation amplitude under large disturbances.

frequency-amplitude (F-A) curve↗

Newton-Raphson AC Power Flow Convergence Based on Deep Learning Initialization and Homotopy Continuation

Power flow forms the basis of many power system studies. With the increased penetration of renewable energy, grid planners tend to perform multiple power flow simulations under various operating conditions and not just selected snapshots at peak or light load conditions. Getting a converged AC power flow (ACPF) case remains a significant challenge for grid planners especially in large power grid networks. This paper proposes a two-stage approach to improve Newton-Raphson ACPF convergence and was applied to a 6102 bus Electric Reliability Council of Texas (ERCOT) system. The first stage utilizes a deep learning-based initializer with data re-training. Here a deep neural network (DNN) initializer is developed to provide better initial voltage magnitude and angle guesses to aid in power flow convergence. This is because Newton-Raphson ACPF is quite sensitive to the initial conditions and bad initialization could lead to divergence. The DNN initializer includes a data re-training framework that improves the initializer's performance when faced with limited training data. The DNN initializer successfully solved 3,285 cases out of 3,899 non-converging dispatch and performed better than random forest and DC power flow initialization methods. ACPF cases not solved in this first stage are then passed through a hot-starting algorithm based on homotopy continuation with switched shunt control. The hot-starting algorithm successfully converged 416 cases out of the remaining 614 non-converging ACPF dispatch. In conclusion, the combined two-stage approach achieved a 94.9% success rate, by converging a total of 3,701 cases out of the initial 3,899 unsolved cases.

Deep learning↗

Power System Transient Stability Analysis through a Homotopy Analysis Method

As an important function of energy management systems (EMSs), online contingency analysis plays an important role in providing power system security warnings of instability. At present, N-1 contingency analysis still relies on time-consuming numerical integration. To save computational cost, the paper proposes a quasi-analytical method to evaluate transient stability through time domain periodic solutions’ frequency sensitivities against initial values. First, dynamic systems described in classical models are modified into damping free systems whose solutions are either periodic or expanded (non-convergent). Second, because the sensitivities experience sharp changes when periodic solutions vanish and turn into expanded solutions, transient stability is assessed using the sensitivity. Third, homotopy analysis is introduced to extract frequency information and evaluate the sensitivities only from initial values so that time consuming numerical integration is avoided. Finally, a simple case is presented to demonstrate application of the proposed method, and simulation results show that the proposed method is promising.

Wang, Shaobu↗

Some aspects of the five-dimensional Lovelock black hole spacetime: Strong homotopy retract, perihelion precession and quasistationary levels

Highlights: • Some topological aspects are explored. • The perihelion precession is calculated. • The general Heun functions are applied to obtain analytical solutions. • The quasibound states and the radial wave eigenfunctions are obtained. In this work we explore some mathematical physics aspects of the spherically symmetric Lovelock black hole in high dimensions. Intended for this aim, we thoroughly consider the metric corresponding to the five-dimensional Lovelock black hole spacetime. We construct the strong retractions by the geodesic equations on the background under consideration. As a result, from the topological point of view, we construct the theory of strong homotopy retract, which will allow us, in principle, to better understand some of its suitable applications on astrophysics and cosmology, in particular, in the analysis of the spacetime singularities. We find the solutions of the equation of motion for both radial and angular coordinates, and then we describe the outer (“exterior”) and lower (“interior”) apparent horizons. Indeed, the outer apparent horizon is the last surface from which the light waves could still escape from the black hole. Thus, it is meaningful to analyze some physical phenomena related to quantum particles propagating outside the exterior apparent horizon, in particular, we discuss the quasistationary levels of scalar fields and their radial wave functions, which are given in terms of the general Heun functions. We also calculate the perihelion precession in this background.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Classical Optimal Control for Energy Minimization Based On Diffeomorphic Modulation under Observable-Response-Preserving Homotopy

In this work, we introduce the so-called “Classical Optimal Control Optimization” (COCO) method for global energy minimization based on the implementation of the diffeomorphic modulation under observable-response-preserving homotopy (DMORPH) gradient algorithm. A probe particle with time-dependent mass $m(t;β)$ and dipole $μ(r,t;β)$ is evolved classically on the potential energy surface $V(r)$ coupled to an electric field $E(t;β)$, as described by the time-dependent density of states represented on a grid, or otherwise as a linear combination of Gaussians generated by the k-means clustering algorithm. Control parameters $β$ defining $m(t;β), μ(r,t;β),$ and $E(t;β)$ are optimized by following the gradients of the energy with respect to $β,$ adapting them to steer the particle toward the global minimum energy configuration. We find that the resulting COCO algorithm is capable of resolving near-degenerate states separated by large energy barriers and successfully locates the global minima of golf potentials on flat and rugged surfaces, previously explored for testing quantum annealing methodologies and the quantum optimal control optimization (QuOCO) method. Preliminary results show successful energy minimization of multidimensional Lennard-Jones clusters. Beyond the analysis of energy minimization in the specific model systems investigated, we anticipate COCO should be valuable for solving minimization problems in general, including optimization of parameters in applications to machine learning and molecular structure determination.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Communication: Newton homotopies for sampling stationary points of potential energy landscapes

One of the most challenging and frequently arising problems in many areas of science is to find solutions of a system of multivariate nonlinear equations. There are several numerical methods that can find many (or all if the system is small enough) solutions but they all exhibit characteristic problems. Moreover, traditional methods can break down if the system contains singular solutions. Here, we propose an efficient implementation of Newton homotopies, which can sample a large number of the stationary points of complicated many-body potentials. We demonstrate how the procedure works by applying it to the nearest-neighbor ϕ{sup 4} model and atomic clusters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Calculation of the neutron diffusion equation by using Homotopy Perturbation Method

The distribution of the neutrons in a nuclear fuel element in the nuclear reactor core can be calculated by the neutron diffusion theory. It is the basic and the simplest approximation for the neutron flux function in the reactor core. In this study, the neutron flux function is obtained by the Homotopy Perturbation Method (HPM) that is a new and convenient method in recent years. One-group time-independent neutron diffusion equation is examined for the most solved geometrical reactor core of spherical, cubic and cylindrical shapes, in the frame of the HPM. It is observed that the HPM produces excellent results consistent with the existing literature.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Real homotopy theory and supersymmetric quantum mechanics

In the context of observing string backgrounds, much work has been devoted to the question of how similar a general quantum field theory (specifically, a two-dimensional superconformal field theory) is to a sigma model. Put differently, one would like to know how well or poorly one can understand the physics of string backgrounds in terms of concepts of classical geometry. Much attention has also been given of late to the question of how geometry can be encoded in a microscopic physical description that makes no explicit reference to space and time. We revisit the first question and review both well-known and less well-known results about geometry and sigma models from the perspective of dimensional reduction to supersymmetric quantum mechanics. The consequences of arising as the dimensional reduction of a d -dimensional theory for the resulting quantum mechanics are explored. In this context, we reinterpret the minimal models of rational (more precisely, complex) homotopy theory as certain supersymmetric Fock spaces, with unusual actions of the supercharges. The data of the Massey products appear naturally as supersymmetric vacuum states that are entangled between different degrees of freedom. This connection between entanglement and geometry is, as far as we know, not well known to physicists. In addition, we take note of an intriguing numerical coincidence in the context of string compactification on hyper-Kähler eight-manifolds.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Classification of fragile topology enabled by matrix homotopy

Flat bands in twisted materials have attracted considerable attention due to the emergence of correlated phases that can be associated with the non-Wannier-representable nature of its single-particle states. Specifically, these bands can exhibit a class of topology that can be nullified by the addition of trivial bands, termed fragile topology, which has required an expansion of prior classification schemes. However, existing approaches for predicting fragile topology rely on momentum-space methods, e.g., Wilson loops, presenting a fundamental challenge for using fragile topology as a predictor of correlated phases in aperiodic systems, such as incommensurate twist angles in moiré materials. Here, we develop a ℤ 2 energy-resolved topological marker for classifying fragile phases using a system’s position-space description, enabling the direct classification of finite, disordered, and aperiodic materials. By translating the physical symmetries protecting the system’s fragile topological phase into matrix symmetries of the system’s Hamiltonian and position operators, we use matrix homotopy to construct our topological marker while simultaneously yielding a quantitative measure of topological robustness. We demonstrate our framework’s effectiveness in both a low-energy tight-binding model and a continuum photonic crystal model of 𝐶 2 ⁢𝒯-symmetric systems, and find that fragile topology can both persist under strong disorder and even exhibit disorder-induced reentrant phase transitions. Our photonic crystal results also demonstrate the robustness of fragile topology, and the applicability of our approach, to heterostructures lacking a bulk spectral gap. Overall, our framework serves as an efficient tool for elucidating fragile topology, offering guidance for the prediction and discovery of correlated phases in both crystalline and aperiodic materials.

Lee, Ki Young [Sandia National Laboratories (SNL-N↗

Homotopy characterization of non-Hermitian Hamiltonians

We revisit the problem of classifying topological band structures in non-Hermitian systems. Recently, a solution has been proposed, which is based on redefining the notion of energy band gap in two different ways, leading to the so-called “point-gap” and “line-gap” schemes. However, simple Hamiltonians without band degeneracies can be constructed which correspond to neither of the two schemes. Here, we resolve this shortcoming of the existing classifications by developing the most general topological characterization of non-Hermitian bands for systems without a symmetry. Our approach, which is based on homotopy theory, makes no particular assumptions on the band gap, and predicts significant extensions to the previous classification frameworks. In particular, we show that the one-dimensional invariant generalizes from $\mathbb{Z}$ winding number to the non-Abelian braid group, and that depending on the braid group invariants, the two-dimensional invariants can be cyclic groups $\mathbb{Z}_n$ (rather than $\mathbb{Z}$ Chern number). Finally, we interpret these results in terms of a correspondence with gapless systems, and we illustrate them in terms of analogies with other problems in band topology, namely, the fragile topological invariants in Hermitian systems and the topological defects and textures of nematic liquids.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Algebraic structures in closed superstring field theory, homotopy transfer, and effective actions

A consistent action for heterotic and type II superstring field theory was recently proposed by Sen. We give an algebraic formulation of this action in terms of certain twisted L ∞ algebra. We further show that Sen’s Wilsonian effective superstring field action can be obtained using homotopy transfer and the effective theory also possesses the algebraic structure of a twisted L ∞ algebra. Published by the American Physical Society 2024

Singh, Ranveer Kumar↗

A Homotopy-Based Method for Robust Computation of Controlling Unstable Equilibrium Points

This paper proposes a homotopy-based method for efficiently computing and precisely determining controlling unstable equilibrium points (CUEPs) in transient stability analysis and screening using direct methods. The proposed method is intended to overcome potential numerical convergence problems associated with computing CUEPs. These convergence problems are mostly due to irregularities in the region of convergence of CUEPs which make it difficult to determine suitable initial points. The most commonly used initial points in computing CUEPs are minimum gradient points (MGPs) on the stability boundary which in turn rely on the determination of exit points (EPs-points at which sustained-fault trajectories exit system stability boundaries). However, determination of MGPs and EPs are computationally involved and sophisticated numerical methods have been introduced such as stability-boundary-following procedure and shadowing method to correct the solution trajectory. The proposed method does not depend on MGPs and also does not require accurate EPs. The proposed method maps the solution from EPs to the CUEPs. The proposed method is applied on several known systems including the NE 39 bus system and the reduced WECC system and the results are provided. The results are compared with the traditional methods such as Boundary of stability region based CUEP method (BCU method).

42 ENGINEERING↗

Homotopy Solver

This software implements parallel versions of an interior-point solver, based on the publicly available ipopt solver. Here we have full control over the linear solver and our algorithm is fully parallel thus enabling scalability to large-scale optimization problems. This package also has a parallel implementation of a homotopy solver developed under the scalable methods for contact LDRD project 23-ERD-017. This solver is an mfem-based implementation of algorithm described in ``A filter trust-region Newton continuation method for nonlinear complementarity problems''. Cosmin G. Petra, Nai-Yuan Chiang, Jingyi Wang, Tucker Hartland, and Michael Puso (submitted), LLNL-JRNL-869761.

Hartland, Tucker [Lawrence Livermore National Labo↗