Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numeric invariants”

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

Numeric invariants from multidimensional persistence

Topological data analysis is the study of data using techniques from algebraic topology. Often, one begins with a finite set of points representing data and a “filter” function which assigns a real number to each datum. Using both the data and the filter function, one can construct a filtered complex for further analysis. For example, applying the homology functor to the filtered complex produces an algebraic object known as a “one-dimensional persistence module”, which can often be interpreted as a finite set of intervals representing various geometric features in the data. If one runs the above process incorporating multiple filter functions simultaneously, one instead obtains a multidimensional persistence module. Unfortunately, these are much more difficult to interpret. In this article, we analyze the space of multidimensional persistence modules from the perspective of algebraic geometry. First we build a moduli space of a certain subclass of easily analyzed multidimensional persistence modules, which we construct specifically to capture much of the information which can be gained by using multidimensional persistence instead of one-dimensional persistence. Fruthermore, we argue that the global sections of this space provide interesting numeric invariants when evaluated against our subclass of multidimensional persistence modules. Finally, we extend these global sections to the space of all multidimensional persistence modules and discuss how the resulting numeric invariants might be used to study data. This paper extends the results of Adcock et al. (Homol Homotopy Appl 18(1), 381–402, 2016) by constructing numeric invariants from the computation of a multidimensional persistence module as given by Carlsson et al. (J Comput Geom 1(1), 72–100, 2010).

58 GEOSCIENCES↗

Design of double- and multi-bend achromat lattices with large dynamic aperture and approximate invariants

A numerical method to design nonlinear double- and multi-bend achromat (DBA and MBA) lattices with approximate invariants of motion is investigated. The search for such nonlinear lattices is motivated by Fermilab’s Integrable Optics Test Accelerator (IOTA), whose design is based on an integrable Hamiltonian system with two invariants of motion. While it may not be possible to design an achromatic lattice for a dedicated synchrotron light source storage ring with one or more exact invariants of motion, it is possible to tune the sextupoles and octupoles in existing DBA and MBA lattices to produce approximate invariants. In our procedure, the lattice is tuned while minimizing the turn-by-turn fluctuations of the Courant-Snyder actions J x and J y at several distinct amplitudes, while simultaneously minimizing diffusion of the on-energy betatron tunes. The resulting lattices share some important features with integrable ones, such as a large dynamic aperture, trajectories confined to invariant tori, robustness to resonances and errors, and a large amplitude-dependent tune-spread. Compared to the nominal NSLS-II lattice, the single- and multi-bunch instability thresholds are increased and the bunch-by-bunch feedback gain can be reduced.

43 PARTICLE ACCELERATORS↗

On an origin of numerical diffusion: Violation of invariance under space-time inversion

The invariant properties of the convection equation du/dt + adu/dx = 0 (where d is the partial differential operator) with respect to spatial reflection, time reversal, and space-time inversion are studied. Generally, a finite-difference analog of this equation may possess some or none of these properties. It is shown that, under certain conditions, the von Neumann amplification factor of an analog satisfies a special relation for each invariant property this analog possesses. Particularly, an analog is neutrally stable and thus free of numerical diffusion if it possesses the invariant property related to space-time inversion. It is also explained why generally (1) an upwind scheme possesses neither the invariant property related to spatial reflection nor that related to space-time inversion, and (2) an explicit scheme possesses neither the invariant property related to time reversal nor that related to space-time inversion. Extension to the viscous case and a remarkable connection between the current work and a new numerical framework for solving conservation laws are also discussed.

Chang, Sin-Chung↗

Numerical solution of the azimuthal-invariant thin-layer Navier-Stokes equations

The paper reports a numerical procedure developed for a two-dimensional azimuthal (or planar) invariant form of the thin-layer Navier-Stokes equations. Generalization of the governing equations is described, and the equations are solved with an implicit approximate factorization finite-difference scheme. Inviscid and viscous results are presented for both external and internal flows and for spinning and nonspinning bodies.

Nietubicz, C. J.↗

Thermal Decomposition of an Impure (Roxbury) Siderite: Relevance to the Presence of Chemically Pure Magnetite Crystals in ALH84001 Carbonate Disks

The question of the origin of nanophase magnetite in Martian meteorite ALH84001 has been widely debated for nearly a decade. Golden et al. have reported producing nearly chemically pure magnetite from thermal decomposition of chemically impure siderite [(Fe, Mg, Mn)CO3]. This claim is significant for three reasons: first, it has been argued that chemically pure magnetite present in the carbonate disks in Martian meteorite ALH84001 could have formed by the thermal decomposition of the impure carbonate matrix in which they are embedded; second, the chemical purity of magnetite has been previously used to identify biogenic magnetite; and, third, previous studies of thermal decomposition of impure (Mg,Ca,Mn)-siderites, which have been investigated under a wide variety of conditions by numerous researchers, invariably yields a mixed metal oxide phase as the product and not chemically pure magnetite. The explanation for this observation is that these siderites all possess the same crystallographic structure (Calcite; R3c) so solid solutions between these carbonates are readily formed and can be viewed on an atomic scale as two chemically different but structurally similar lattices.

McKay, D.S.↗

Parallel-In-Time Magnus Integrators

Magnus integrators are a subset of geometric integration methods for the numerical solution of ordinary differential equations that conserve certain invariants in the numerical solution. Our report explores temporal parallelism of Magnus integrators, particularly in the context of nonlinear problems. The method combines the concurrent computation of matrix commutators and exponentials within a time step with a pipelined iteration applied to multiple time steps in parallel. The accuracy and efficiency of time parallel Magnus methods up to order six are highlighted through numerical examples and demonstrate that critical parallel speedup is possible compared to serial methods.

97 MATHEMATICS AND COMPUTING↗

A Numerical Method for Computing the State Transition Matrix Using Poincare Integral Invariants

The Poincare integral invariants describe the volumes of sets in Hamiltonian phase space. We use these invariants to derive a new numerical procedure for obtaining the state transition matrix (STM), which can be applied to both conservative and nonconservative systems. The method is analogous to a finite difference approximation of the STM, where perturbed states are numerically propagated along with the reference trajectory. We discuss the mathematical similarities between this new STM and existing methods, show numerical results for orbital motion and uncertainty propagation, and discuss new insights afforded by the Hamiltonian properties of phase flow.

state transition matrix↗

A consistent and conservative model and its scheme for N -phase- M -component incompressible flows

Here, we propose a consistent and conservative model for multiphase and multicomponent incompressible flows, where there can be arbitrary numbers of phases and components. Each phase has a background fluid called the pure phase, each pair of phases is immiscible, and components are dissolvable in some specific phases. The model is developed based on the multiphase Phase-Field model including the contact angle boundary condition, the diffuse domain approach, and the analyses on the proposed consistency conditions for multiphase and multicomponent flows. The model conserves the mass of individual pure phases, the amount of each component in its dissolvable region, and thus the mass of the fluid mixture, and the momentum of the flow. It ensures that no fictitious phases or components can be generated and that the summation of the volume fractions from the Phase-Field model is unity everywhere so that there is no local void or overfilling. It satisfies a physical energy law and it is Galilean invariant. A corresponding numerical scheme is developed for the proposed model, whose formal accuracy is 2nd-order in both time and space. It is shown to be consistent and conservative and its solution is demonstrated to preserve the Galilean invariance and energy law. Numerical tests indicate that the proposed model and scheme are effective and robust to study various challenging multiphase and multicomponent flows.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

IR-safe and UV-safe integrands in the EFTofLSS with exact time dependence

Because large-scale structure surveys may very well be the next leading sources of cosmological information, it is important to have a precise understanding of the cosmological observables; for this reason, the Effective Field Theory of Large-Scale Structure (EFTofLSS) was developed. So far, most results in the EFTofLSS have used the so-called Einstein-de Sitter approximation, an approximation of the time dependence which is known to be accurate to better than one percent. However, in order to reach even higher accuracy, the full time dependence must be used. The computation with exact time dependence is sensitive to both infrared (IR) and ultraviolet (UV) effects in the loop integrands, and while these effects must cancel because of diffeomorphism invariance, they make numerical computation much less efficient. We provide a formulation of the one-loop, equal-time exact-time-dependence power spectrum of density perturbations which is manifestly free of these spurious IR and UV divergences at the level of the integrand. We extend our results to the total matter mode with clustering quintessence, show that IR and UV divergences cancel, and provide the associated IR- and UV-safe integrand. This also establishes that the consistency conditions are satisfied in this system. In conclusion, we then use our one-loop result to do an improved precision comparison of the two-loop dark-matter power spectrum with the Dark Sky N-body simulation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Stochastic robustness of linear time-invariant control systems

A simple numerical procedure for estimating the stochastic robustness of a linear, time-invariant system is described. Monte Carlo evaluation of the system's eigenvalues allows the probability of instability and the related stochastic root locus to be estimated. This definition of robustness is an alternative to existing deterministic definitions that address both structured and unstructured parameter variations directly. This analysis approach treats not only Gaussian parameter uncertainties but also non-Gaussian cases, including uncertain-but-bounded variations. Trivial extensions of the procedure permit alternate discriminants to be considered. Thus the probabilities that stipulated degrees of instability will be exceeded or that closed-loop roots will leave desirable regions can also be estimated. Results are particularly amenable to graphical presentation.

Stengel, Robert F.↗

Technical notes and correspondence: Stochastic robustness of linear time-invariant control systems

A simple numerical procedure for estimating the stochastic robustness of a linear time-invariant system is described. Monte Carlo evaluations of the system's eigenvalues allows the probability of instability and the related stochastic root locus to be estimated. This analysis approach treats not only Gaussian parameter uncertainties but non-Gaussian cases, including uncertain-but-bounded variation. Confidence intervals for the scalar probability of instability address computational issues inherent in Monte Carlo simulation. Trivial extensions of the procedure admit consideration of alternate discriminants; thus, the probabilities that stipulated degrees of instability will be exceeded or that closed-loop roots will leave desirable regions can also be estimated. Results are particularly amenable to graphical presentation.

Stengel, Robert F.↗

Narrow operator models of stellarator equilibria in Fourier Zernike basis

Numerical computation of the ideal magnetohydrodynamic (MHD) equilibrium magnetic field is at the base of stellarator optimisation and provides the starting point for solving more sophisticated partial differential equations like transport or turbulence models. Conventional approaches solve for a single stationary point of the ideal MHD equations, which is fully defined by three invariants and the numerical scheme employed by the solver. We present the first numerical approach that can solve for a continuous distribution of equilibria with fixed boundary and rotational transform, varying only the pressure invariant. This approach minimises the force residual by optimising parameters of multilayer perceptrons that map from a scalar pressure multiplier to the Fourier Zernike basis as implemented in the modern stellarator equilibrium solver DESC.

fusion plasma↗

Small-scale flow topologies in decaying isotropic turbulence laden with finite-size droplets

The topology of the fine-scale motions in decaying isotropic turbulence laden with droplets of super-Kolmogorov size is investigated using results from direct numerical simulations. The invariants of the velocity-gradient, rate-of-strain, and rate-of-rotation tensors are computed in the carrier phase. The joint probability density functions of the invariants are calculated and conditioned on different distances from the droplet surface. The results show that outside the viscous region near the interface, the flow topologies favor stable focus/stretching and unstable node/saddle/saddle structures, which is in agreement with those found in canonical homogeneous isotropic turbulence. Inside the viscous layer at the droplet surface, the flow topologies shift from a preference for high-enstrophy/low-dissipation motions to favoring low-enstrophy/high-dissipation. At the droplet surface, there is a strong tendency for boundary-layer-like and vortex-sheet flow topologies in which the strain and rotation rates are positively correlated. Furthermore, an interesting observation is that the shapes of the invariant distributions at the droplet surface are remarkably similar to those reported in the viscous sublayer of turbulent wall-bounded flows. Also, the results show that the smallest hydrodynamic length scale of the carrier fluid turbulence is located at the droplet interface and that this length scale is two to three times smaller than that of the surrounding bulk flow. From a computational viewpoint, this suggests a more stringent spatial resolution requirement for the direct numerical simulation of finite-size droplets in isotropic turbulence than its single-phase counterpart.

42 ENGINEERING↗

Approximating the linear quadratic optimal control law for hereditary systems with delays in the control

The fundamental control synthesis issue of establishing a priori convergence rates of approximation schemes for feedback controllers for a class of distributed parameter systems is addressed within the context of hereditary systems. Specifically, a factorization approach is presented for deriving approximations to the optimal feedback gains for the linear regulator-quadratic cost problem associated with time-varying functional differential equations with control delays. The approach is based on a discretization of the state penalty which leads to a simple structure for the feedback control law. General properties of the Volterra factors of Hilbert-Schmidt operators are then used to obtain convergence results for the controls, trajectories and feedback kernels. Two algorithms are derived from the basic approximation scheme, including a fast algorithm, in the time-invariant case. A numerical example is also considered.

Milman, Mark H.↗

Approximating the linear quadratic optimal control law for hereditary systems with delays in the control

The fundamental control synthesis issue of establishing a priori convergence rates of approximation schemes for feedback controllers for a class of distributed parameter systems is addressed within the context of hereditary schemes. Specifically, a factorization approach is presented for deriving approximations to the optimal feedback gains for the linear regulator-quadratic cost problem associated with time-varying functional differential equations with control delays. The approach is based on a discretization of the state penalty which leads to a simple structure for the feedback control law. General properties of the Volterra factors of Hilbert-Schmidt operators are then used to obtain convergence results for the controls, trajectories and feedback kernels. Two algorithms are derived from the basic approximation scheme, including a fast algorithm, in the time-invariant case. A numerical example is also considered.

Milman, Mark H.↗