Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear system”

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 1,297 records · Page 72
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

A scalable matrix-free spectral element approach for unsteady PDE constrained optimization using PETSc/TAO

In this work, we provide a new approach for the efficient matrix-free application of the transpose of the Jacobian for the spectral element method for the adjoint-based solution of partial differential equation (PDE) constrained optimization. This results in optimizations of nonlinear PDEs using explicit integrators where the integration of the adjoint problem is not more expensive than the forward simulation. Solving PDE constrained optimization problems entails combining expertise from multiple areas, including simulation, computation of derivatives, and optimization. The Portable, Extensible Toolkit for Scientific computation (PETSc) together with its companion package, the Toolkit for Advanced Optimization (TAO), is an integrated numerical software library that contains an algorithmic/software stack for solving linear systems, nonlinear systems, ordinary differential equations, differential algebraic equations, and large-scale optimization problems and, as such, is an ideal tool for performing PDE-constrained optimization. This paper describes an efficient approach in which the software stack provided by PETSc/TAO can be used for large-scale nonlinear time-dependent problems. Time integration can involve a range of high-order methods, both implicit and explicit. The PDE-constrained optimization algorithm used is gradient-based and seamlessly integrated with the simulation of the physical problem.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

A memory-distributed quasi-Newton solver for nonlinear programming problems with a small number of general constraints

We approach the problem of parallelizing state-of-the-art nonlinear programming optimization algorithms.Specifically, we focus on parallelizing quasi-Newton interior-point methods that use limited-memory secant Hessian approximations. Such interior-point methods are known to have better convergence properties and to be more effective on large-scale problems than gradient-based and derivative-free optimization algorithms. We target nonlinear and potentially nonconvex optimization problems with an arbitrary number of bound constraints and a small number of general equality and inequality constraints on the optimization variables. These problems occur for example in the form of optimal control, optimal design, and inverse problems governed by ordinary or partial differential equations, whenever they are expressed in a “reduced-space” optimization approach. We introduce and analyze the time and space complexity of a decomposition method for solving the quasi-Newton linear systems that leverages the fact that the quasi-Newton Hessian matrix has a small number of dense blocks that border a low-rank update of a diagonal matrix. This enables an efficient parallelization on memory-distributed computers of the iterations of the optimization algorithm, a state-of-the-art filter line-search interior-point algorithm by Wächter et. al. We illustrate the efficiency of the proposed method by solving structural topology optimization problems on up to 4608 cores on a parallel machine.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Vibration mitigation in interlocking metasurfaces

Interlocking metasurfaces (ILMs) are a new class of mechanical metamaterials with patterned surface features that enable joining of two adjacent bodies, serving as an alternative to welds, fasteners, or adhesives. While the strength of ILMs has been explored previously, it is also possible to envision ILMs that possess the ability to dampen vibration or tailor energy transmission across the interface. In this study, we explore the vibration transmission under frequency sweep and random excitation of three ILM designs. This is the first study to characterize the response of ILMs under vibration. We find that ILMs which are tightly coupled (e.g., small interference fit) behave as a linear system and depart only marginally from the behavior of a solid body. However, through strategic manipulation of the ILM topology, and in particular the fit clearances, we show that vibrations can be damped by over an order of magnitude. Design strategies to employ ILMs for vibration mitigation are developed using a first-order model.

Bolmin, O. [Carnegie Mellon Univ., Pittsburgh, PA ↗
DOE OSTI2022

Characterization of nonlinear ultrasonic waves behavior while interacting with poor interlayer bonds in large-scale additive manufactured materials

Over the past decades, researchers have developed several nonlinear ultrasonic techniques for quality control of materials commonly used in different applications. Owing to the superior sensitivity of nonlinear ultrasound waves to small defects such as micro-cracks, their applicability in different nondestructive testing (NDT) problems has been investigated in numerous studies. These studies utilize frequency domain analysis to detect the generation of higher harmonics because of the formation of defects in the inspected medium. Frequency domain analysis based on the Fourier transform is a significant approach used in linear systems; however, it may not perform adequately on nonlinear systems. Hence, studies on nonlinear dynamics and physics consider analyzing systems' behavior in the phase-space domain. In contrast to the frequency domain analysis, which can result in information loss, analysis in the phase-space domain retains all the information regarding a system's states. Here, we investigate the nonlinearities induced by poor interlayer bonding in polymer-based additive manufactured parts in the phase-space domain. It is convenient to characterize the nonlinearity in the phase-space domain because it provides a geometrical representation of a system's states. Two types of low quality interlayer bond are considered. The first type is simulated artificially while the second type is manufactured by reducing the bond quality during the printing process. The analysis verified that the received ultrasonic signals exhibit classical nonlinear behavior in the phase-space domain while interacting with simulated poor interlayer bonds. In addition, the results showed that the behavior of ultrasonic waves is amplitude-dependent and evolves into models that have not been previously reported. Furthermore, Largest Lyapunov Exponent (LLE) is used to quantify the behavior of nonlinear ultrasonic waves while interacting with poor interlayer bonds. Using LLE, it was observed that the divergence rate of the phase-space trajectories depends on the amplitude of the excitation. This observation quantitatively proves that nonlinear behavior of ultrasound while interacting with poor interlayer bonds can be amplitude-dependent. The results of both simulated and inherent poor interlayer bond cases showed that LLE can be used as a reliable quantitative damage-sensitive feature to detect and potentially characterize weak bonds, which are difficult to detect using conventional approaches. Additionally, the reported results in the phase-space domain provide a basis for proposing a new mathematical model for ultrasonic waves interacting with poor interlayer bonds.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Incomplete Sparse Approximate Inverses for Parallel Preconditioning

In this study, we propose a new preconditioning method that can be seen as a generalization of block-Jacobi methods, or as a simplification of the sparse approximate inverse (SAI) preconditioners. The “Incomplete Sparse Approximate Inverses” (ISAI) is in particular efficient in the solution of sparse triangular linear systems of equations. Those arise, for example, in the context of incomplete factorization preconditioning. ISAI preconditioners can be generated via an algorithm providing fine-grained parallelism, which makes them attractive for hardware with a high concurrency level. Finally, in a study covering a large number of matrices, we identify the ISAI preconditioner as an attractive alternative to exact triangular solves in the context of incomplete factorization preconditioning.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Scalable line and plane relaxation in a parallel structured multigrid solver

The efficient solution of sparse, linear systems that arise through the discretization of partial differential equations remains a key challenge for a range of high performance scientific simulations. One approach for reducing data movement and improving performance is by exposing and exploiting structure in a problem through the use of robust structured multilevel solvers. By choosing coarsening that preserves the structure of the problem, these methods maintain efficient structured computation and communication throughout the multigrid hierarchy. However, when coarsening is not permitted to be dependent on the operator, anisotropy must be addressed by the smoother — producing error compatible for coarse-grid correction with structured coarsening. Here, the components required in a scalable parallel structured solver are described with a focus on memory and communication efficiency of robust smoothers. While the implementation of communication and memory reduction techniques in smoothers integrated in a complete 3D solver present a significant engineering challenge, a novel approach is proposed that addresses these challenges systematically through a change to the solver’s execution model. Enabled by user-level threading paired with a set of data and communication abstractions, this approach permits seamless aggregation of communication in plane smoothers — directly reusing code for a 2D distributed multilevel cycle. Results show an effective reduction in communication costs for coarse-grid problems, and result in a speedup of 8.7x in smoothing routines shown in Fig. 12 using this approach. This produces a significant improvement to strong scalability while maintaining favorable weak scaling behavior. Finally, a parallel scaling study using a series of refined meshes is included that demonstrates the effectiveness of this approach in an application of interest.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

A Comparison of Linear Solvers for Resolving Flow in Three-Dimensional Discrete Fracture Networks

We compare various methods for resolving steady flow within three-dimensional discrete fracture networks, including direct methods, Krylov subspace methods with and without preconditioning, and multi-grid methods. We compared the performance of the methods based on compute times and scaling of the solution as a function of the number of grid nodes and log-variance of the hydraulic aperture. The methods are applied to three test cases: (a) variable density of networks with a truncated power-law distribution of fracture lengths, (b) a fixed network composed of monodisperse fracture sizes but varied permeability/aperture heterogeneity, (c) and a network based on field site in Nevada, US. We chose these cases to allow us to study the impact of the mesh size and flow properties, as well as to demonstrate our conclusions on a large-scale, realistic problem (more than 40 million mesh nodes). A direct solution using Cholesky factorization outperformed other methods for every example but was closely followed in performance by some algebraic multigrid (AMG) preconditioned Krylov subspace methods. Among the Krylov methods, conjugate gradients (CG) with an AMG preconditioner performs the best. Generally, Cholesky factorization is recommended, but CG with an AMG preconditioner may be suitable for very large problems beyond 40 million nodes where the entire linear system cannot reside in memory.

58 GEOSCIENCES↗
DOE OSTI2024

Addressing quantum’s “fine print” with efficient state preparation and information extraction for quantum algorithms and geologic fracture networks

Abstract Quantum algorithms provide an exponential speedup for solving certain classes of linear systems, including those that model geologic fracture flow. However, this revolutionary gain in efficiency does not come without difficulty. Quantum algorithms require that problems satisfy not only algorithm-specific constraints, but also application-specific ones. Otherwise, the quantum advantage carefully attained through algorithmic ingenuity can be entirely negated. Previous work addressing quantum algorithms for geologic fracture flow has illustrated core algorithmic approaches while incrementally removing assumptions. This work addresses two further requirements for solving geologic fracture flow systems with quantum algorithms: efficient system state preparation and efficient information extraction. Our approach to addressing each is consistent with an overall exponential speed-up.

58 GEOSCIENCES↗
DOE OSTI2019

Inertial extended-Lagrangian scheme for solving charge equilibration models

The inertial extended Lagrangian/self-consistent field scheme (iEL–SCF) has been adopted for solving charge equilibration in LAMMPS as part of the reactive force field ReaxFF, which due to the charge conservation constraint requires solving two sets of linear system of equations for the new charges at each molecular dynamics time-step. Therefore, the extended Lagrangian for charge equilibration is comprised of two auxiliary variables for the intermediate charges which serve as an initial guess for the real charges. We demonstrate that the iEL–SCF is able to reduce the number of SCF cycles by 50–80% of the original conjugate gradient self-consistent field solver as tested across diverse systems including water, ferric hydroxide, nitramine RDX, and hexanitrostilbene.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

High-Fidelity Large-Signal Order Reduction Approach for Composite Load Model

With the increasing penetration of electronic loads and distributed energy resources, conventional load models cannot capture their dynamics. Therefore, a new comprehensive composite load model is developed by Western Electricity Coordinating Council (WECC). However, this model is a complex high-order non-linear system with multi-time-scale property, which poses challenges on power system studies with heavy computational burden. In order to reduce the model complexity, the authors firstly develop a large-signal order reduction (LSOR) method using singular perturbation theory. In this method, the fast dynamics are integrated into the slow ones to preserve transient characteristics of the former. Then, accuracy assessment conditions are proposed and embedded into the LSOR to improve and guarantee the accuracy of reduced-order model. Finally, the reduced-order WECC composite load model is derived by using the proposed algorithm. Overall, simulation results show that the reduced-order large-signal model significantly alleviates the computational burden while maintaining similar dynamic responses as the original composite load model.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2012

Dirac matrices for Chern-Simons gravity

A genuine gauge theory for the Poincare, de Sitter or anti-de Sitter algebras can be constructed in (2n- 1)-dimensional spacetime by means of the Chern-Simons form, yielding a gravitational theory that differs from General Relativity but shares many of its properties, such as second order field equations for the metric. The particular form of the Lagrangian is determined by a rank n, symmetric tensor invariant under the relevant algebra. In practice, the calculation of this invariant tensor can be reduced to the computation of the trace of the symmetrized product of n Dirac Gamma matrices {Gamma}{sub ab} in 2n-dimensional spacetime. While straightforward in principle, this calculation can become extremely cumbersome in practice. For large enough n, existing computer algebra packages take an inordinate long time to produce the answer or plainly fail having used up all available memory. In this talk we show that the general formula for the trace of the symmetrized product of 2n Gamma matrices {Gamma}{sub ab} can be written as a certain sum over the integer partitions s of n, with every term being multiplied by a numerical cofficient {alpha}{sub s}. We then give a general algorithm that computes the {alpha}-coefficients as the solution of a linear system of equations generated by evaluating the general formula for different sets of tensors B{sup ab} with random numerical entries. A recurrence relation between different coefficients is shown to hold and is used in a second, 'minimal' algorithm to greatly speed up the computations. Runtime of the minimal algorithm stays below 1 min on a typical desktop computer for up to n = 25, which easily covers all foreseeable applications of the trace formula.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2013

Note: Spectral motional Stark effect diagnostic for measurement of magnetic fields below 0.3 T

The paper reports on development of the spectral motional Stark effect (MSE) diagnostic in the midplane of the gas dynamic trap (GDT) linear system for magnetic confinement of anisotropic hot-ion plasma. The axially symmetric GDT vacuum magnetic field has a minimum value in the midplane, which varies from 0.2 to 0.35 T in different regimes of operation. Buildup of 15 keV ion population generates a diamagnetic reduction of magnetic field in the plasma core of up to 30% in the maximum density region, as measured by the existing eight-line MSE diagnostic. Commissioning of the midplane MSE provided first direct measurements of diamagnetic modifications in the minimum magnetic field GDT section, a necessary complement to the understanding of equilibrium and self-organization of high-β plasmas in GDT. Making use of the stable short-pulse diagnostic beam and calibration of the apparent spectral width of beam emission lines allow for the measurement of the plasma magnetic field of 0.29 ± 0.007 T with the integration time of 200 μs.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2014

A generalized tool for accurate time-domain separation of excited modes in spin-torque oscillators

We propose and develop an advanced signal processing technique that, combined with micromagnetic simulations, is able to deeply describe the non-stationary behavior of spin-torque oscillators, both in terms of time domain and spatial distribution of the magnetization dynamics. The Hilbert-Huang Transform is used for the identification of the time traces of each oscillation in a multimode excitation and enhanced with masking signals and the Ensemble Empirical Mode Decomposition. We emphasize that the technique developed here is general and can be used for any physical non-linear system in the presence of multimode dynamical excitation or intermittence.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

Cluster synchronization in networks of neurons with chemical synapses

In this work, we study the cluster synchronization of chemically coupled and generally formulated networks which are allowed to be nonidentical. The sufficient condition for the existence of stably synchronous clusters is derived. Specifically, we only need to check the stability of the origins of m decoupled linear systems. Here, m is the number of subpopulations. Examples of nonidentical networks such as Hindmarsh-Rose (HR) neurons with various choices of parameters in different subpopulations, or HR neurons in one subpopulation and FitzHugh-Nagumo neurons in the other subpopulation are provided. Explicit threshold for the coupling strength that guarantees the stably cluster synchronization can be obtained.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Resonances and oscillatory behavior near multi-species plasma equilibria

We consider dynamic multi-species plasma equilibria whose variables depend on a single spatial coordinate and linear perturbations of these. The linearized system may be reduced to a second-order one satisfied by the respective fluid streamfunctions. For the two-species case, the electron mass is a parameter small enough for a WKB asymptotic analysis to be justified. It turns out that the points where either the ion or electron equilibrium velocity equals the ratio between the temporal and transversal frequencies of the perturbation are turning or singular points of the system, connecting exponentially increasing or decreasing solutions to oscillatory ones. The crucial role of singular points in the balance between the different contributions to the electron kinetic energy is explored.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Characterization and calibration of 8-channel E-band heterodyne radiometer system for SST-1 tokamak

An 8-channel E-band heterodyne radiometer system (74–86 GHz) is designed, characterized, and calibrated to measure the radial electron temperature profile by measuring Electron Cyclotron Emission spectrum at SST-1 Tokamak. The developed radiometer has a noise equivalent temperature of 1 eV and sensitivity of 5 × 10{sup 9} V/W. In order to precisely measure the absolute value of electron temperature, a calibration measurement of the radiometer system is performed using hot-cold Dicke switch method, which confirms the system linearity.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2015

Stochastic linear multistep methods for the simulation of chemical kinetics

In this paper, we introduce the Stochastic Adams-Bashforth (SAB) and Stochastic Adams-Moulton (SAM) methods as an extension of the τ-leaping framework to past information. Using the Θ-trapezoidal τ-leap method of weak order two as a starting procedure, we show that the k-step SAB method with k ≥ 3 is order three in the mean and correlation, while a predictor-corrector implementation of the SAM method is weak order three in the mean but only order one in the correlation. These convergence results have been derived analytically for linear problems and successfully tested numerically for both linear and non-linear systems. A series of additional examples have been implemented in order to demonstrate the efficacy of this approach.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2015

Isotropic octupolar second harmonic generation response in LaBGeO{sub 5} glass-ceramic with spherulitic precipitation

A spherulitic crystallization of the crystalline phase LaBGeO{sub 5} is generated in the 25La{sub 2}O{sub 3}-25B{sub 2}O{sub 3}-50 GeO{sub 2} glass system. Linear and nonlinear optical properties of lanthanum borogermanate glass-ceramic have been investigated at both macroscopic and microscopic scales. Polarized μ-Raman analysis has evidenced a radial distribution of the crystallites along the c-axis inside spherulites, whereas polarized μ-Second Harmonic Generation (SHG) analysis revealed intensity maxima perpendicularly to the c-axis crystallites orientation. Polarized SHG mapping of a spherulite indicate that no dipolar response along the c-axis oriented crystallites occurs despite the individual dipolar symmetry C{sub 3} of the crystallites. At a larger mm scale, the isotropic scattering of spherulites recorded from macroscopic SHG experiment in the forward direction is consistent with an average coherent octupolar response per spherulite. These SHG analyses at different scale are both in accordance with radial antiferroelectric orientation along the c-axis of crystallites inside each spherulite.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗