Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Algebraic models”

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,459 records · Page 81
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Compositional modeling in porous media using constant volume flash and flux computation without the need for phase identification

The paper deals with the numerical solution of a compositional model describing compressible two-phase flow of a mixture composed of several components in porous media with species transfer between the phases. The mathematical model is formulated by means of the extended Darcy's laws for all phases, components continuity equations, constitutive relations, and appropriate initial and boundary conditions. The splitting of components among the phases is described using a new formulation of the local thermodynamic equilibrium which uses volume, temperature, and moles as specification variables. The problem is solved numerically using a combination of the mixed-hybrid finite element method for the total flux discretization and the finite volume method for the discretization of transport equations. A new approach to numerical flux approximation is proposed, which does not require the phase identification and determination of correspondence between the phases on adjacent elements. The time discretization is carried out by the backward Euler method. The resulting large system of nonlinear algebraic equations is solved by the Newton–Raphson iterative method. We provide eight examples of different complexity to show reliability and robustness of our approach.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Chern-Simons-matter dualities with SO and USp gauge groups

In the last few years several dualities were found between the low-energy behaviors of Chern-Simons-matter theories with unitary gauge groups coupled to scalars, and similar theories coupled to fermions. In this paper we generalize those dualities to orthogonal and symplectic gauge groups. In particular, we conjecture dualities between SO(N) k Chern-Simons theories coupled to N f real scalars in the fundamental representation, and SO(k)- N+N f /2 coupled to N f real (Majorana) fermions in the fundamental. For N f = 0 these are just level-rank dualities of pure Chern-Simons theories, whose precise form we clarify. They lead us to propose new gapped boundary states of topological insulators and superconductors. As a result, for k = 1 we get an interesting low-energy duality between N f free Majorana fermions and an SO( N ) 1 Chern-Simons theory coupled to N f scalar fields (with N f ≤ N -2).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Fermion geometry and the renormalization of the Standard Model Effective Field Theory

The geometry of field space governs on-shell scattering amplitudes. We formulate a geometric description of effective field theories which extends previous results for scalars and gauge fields to fermions. The field-space geometry reorganizes and simplifies the computation of quantum loop corrections. Using this geometric framework, we calculate the fermion loop contributions to the renormalization group equations for bosonic operators in the Standard Model Effective Field Theory up to mass dimension eight.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2018

A multigrid waveform relaxation method for solving the poroelasticity equations

In this work, a multigrid waveform relaxation method is proposed for solving a collocated finite difference discretization of the linear Biot’s model. This gives rise to the first space–time multigrid solver for poroelasticity equations in the literature. The waveform relaxation iteration is based on a point-wise Vanka smoother that couples the pressure variable at a grid-point with the displacements around it. A semi-algebraic mode analysis is proposed to theoretically analyze the convergence of the multigrid waveform relaxation algorithm. This analysis is novel since it combines the semi-algebraic analysis, suitable for parabolic problems, with the non-standard analysis for overlapping smoothers. The practical utility of the method is illustrated through several numerical experiments in one and two dimensions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

A robust hierarchical solver for ill-conditioned systems with applications to ice sheet modeling

A hierarchical solver is proposed for solving sparse ill-conditioned linear systems in parallel. The solver is based on a modification of the LoRaSp method, but employs a deferred-compression technique, which provably reduces the approximation error and significantly improves efficiency. Moreover, the deferred-compression technique introduces minimal overhead and does not affect parallelism. As a result, the new solver achieves linear computational complexity under mild assumptions and excellent parallel scalability. To demonstrate the performance of the new solver, we focus on applying it to solve sparse linear systems arising from ice sheet modeling. The strong anisotropic phenomena associated with the thin structure of ice sheets creates serious challenges for existing solvers. To address the anisotropy, we additionally developed a customized partitioning scheme for the solver, which captures the strong-coupling direction accurately. In general, the partitioning can be computed algebraically with existing software packages, and thus the new solver is generalizable for solving other sparse linear systems. Our results show that ice sheet problems of about 300 million degrees of freedom have been solved in just a few minutes using a thousand processors.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Enabling Extended-Term Simulation of Power Systems with High PV Penetration. Final Report

This is the final Technical Report for DOE-SETO Project Award # DE-EE0036461. The goal of this project is to advance the understanding of the grid impact of high penetration of photovoltaic (PV) generation by developing novel numerical methods to solve the differential algebraic equations (DAEs) that define power systems. This will overcome the limitations of current software packages – namely that they only consider fast dynamics over brief time periods. The work presented in this final project report covers results over the entire period of the project. This includes results on model development, code development for the PST repository, datasets in the PST repository, algorithm development and results from variable time-step simulations, development and results from multirate simulations, and sensitivity analysis of key parameter in variable time-step methods. In addition, this report discusses project outreach activities to stakeholders, and a summary of project products. Also covered in this final report is the writing of two conference papers (one of which has already been accepted) and a journal paper. In addition, the updating of two inverter models (both grid forming and grid following) to be compatible with the latest version of PST software is discussed.

14 SOLAR ENERGY↗
NASA NTRS1974

A Computing Method for Sound Propagation Through a Nonuniform Jet Stream

Understanding the principles of jet noise propagation is an essential ingredient of systematic noise reduction research. High speed computer methods offer a unique potential for dealing with complex real life physical systems whereas analytical solutions are restricted to sophisticated idealized models. The classical formulation of sound propagation through a jet flow was found to be inadequate for computer solutions and a more suitable approach was needed. Previous investigations selected the phase and amplitude of the acoustic pressure as dependent variables requiring the solution of a system of nonlinear algebraic equations. The nonlinearities complicated both the analysis and the computation. A reformulation of the convective wave equation in terms of a new set of dependent variables is developed with a special emphasis on its suitability for numerical solutions on fast computers. The technique is very attractive because the resulting equations are linear in nonwaving variables. The computer solution to such a linear system of algebraic equations may be obtained by well-defined and direct means which are conservative of computer time and storage space. Typical examples are illustrated and computational results are compared with available numerical and experimental data.

Padula, S. L.↗
NASA NTRS2000

TOA Lightning Location Retrieval on Spherical and Oblate Spheroidal Earth Geometries

A simple linear algebraic solution is introduced for retrieving the location and time-of-occurrence of lightning ground strikes on a spherical Earth from a network of four or more time-of-arrival (TOA) sensors. Since the solution accounts for Earth curvature, it represents an extension to earlier planar model results described by Koshak et al. A test of the retrieval method is provided using computer-simulated data sets. The method is easy to comprehend and completely avoids reference to the mathematics of spherical hyperbolas such as discussed by Lewis. A quasi-analytic extension to the spherical Earth solution is provided for an oblate spheroidal Earth geometry, and the importance/relevance of oblate effects are discussed. Future application of these methods in support of the North American National Lightning Detection Network (NALDN) described by Cummins et al. is desirable, but additional theoretical investigations are required to incorporate magnetic bearing information into the present solution process.

Koshak, W. J.↗
DOE OSTI2014

Magnetization stability analysis of the Stoner-Wohlfarth model under a spin-polarized current with a tilted polarization

The stationary-state solutions of magnetization dynamics under a spin-polarized current that was polarized in an arbitrary direction were investigated by solving the Landau-Lifshitz-Gilbert-Slonczewski equation for a single-domain magnet. Taking into consideration the uniaxial magnetic anisotropy, the equilibrium directions of the magnetization vectors were analytically obtained by solving an algebraic cubic equation. It was found that one to three pairs of magnetization equilibrium states existed, depending on the current intensity and the direction of the spin polarization. By numerically analyzing the stabilities of these equilibrium states, the threshold switching current for the reversing the magnetic vector was obtained under different current polarization configurations, which may be useful for use in future spintronics devices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1992

Geometry modeling and multi-block grid generation for turbomachinery configurations

An interactive 3D grid generation code, Turbomachinery Interactive Grid genERation (TIGER), was developed for general turbomachinery configurations. TIGER features the automatic generation of multi-block structured grids around multiple blade rows for either internal, external, or internal-external turbomachinery flow fields. Utilization of the Bezier's curves achieves a smooth grid and better orthogonality. TIGER generates the algebraic grid automatically based on geometric information provided by its built-in pseudo-AI algorithm. However, due to the large variation of turbomachinery configurations, this initial grid may not always be as good as desired. TIGER therefore provides graphical user interactions during the process which allow the user to design, modify, as well as manipulate the grid, including the capability of elliptic surface grid generation.

Shih, Ming H.↗
DOE OSTI2018

Two-Level Least Squares Methods in Krylov Subspaces

Two-level least squares acceleration approaches are applied to the Chebyshev acceleration method and the restarted conjugate residual method in solving systems of linear algebraic equations with sparse unsymmetric coefficient matrices arising from finite volume or finite element approximations of boundary-value problems on irregular grids. Application of the proposed idea to other iterative restarted processes also is considered. The efficiency of the algorithms suggested is investigated numerically on a set of model Dirichlet problems for the convection-diffusion equation.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1988

Calculation of symmetric and asymmetric vortex seperation on cones and tangent ogives based on discrete vortex models

An inviscid discrete vortex model, with newly derived expressions for the tangential velocity imposed at the separation points, is used to investigate the symmetric and asymmetric vortex separation on cones and tangent ogives. The circumferential locations of separation are taken from experimental data. Based on a slender body theory, the resulting simultaneous nonlinear algebraic equations in a cross-flow plane are solved with Broyden's modified Newton-Raphson method. Total force coefficients are obtained through momentum principle with new expressions for nonconical flow. It is shown through the method of function deflation that multiple solutions exist at large enough angles of attack, even with symmetric separation points. These additional solutions are asymmetric in vortex separation and produce side force coefficients which agree well with data for cones and tangent ogives.

Chin, S.↗
NASA NTRS1977

Frequency domain compensation of a DYNGEN turbofan engine model

Following Rosenbrock's ideas regarding the advantages of dominance in linear multivariable control systems, a new graphical technique is used for the design of compensators that achieve dominance. The technique is illustrated with an application to the problem of designing compensators for a linear turbofan-engine model. The resulting design is put into perspective by examining it in the light of two other multivariable frequency-domain methods. One, MacFarlane's method of characteristic loci, is used to realize a final design for stability and low interaction. The other is a direct technique based upon the algebraic expansion of the determinant of the return difference in terms of it's elements. Results from simulations carried out on the NASA DYNGEN software are included.

Schafer, R. M.↗
NASA NTRS1995

Elliptic surface grid generation on minimal and parametrized surfaces

An elliptic grid generation method, which generates boundary conforming grids in a two dimensional physical space, is presented. The method is based on the composition of an algebraic and elliptic transformation. The composite mapping obeys the Poisson grid generation system with control functions specified by the algebraic transformation. It is shown that the grid generation on a minimal surface in a three dimensional space is equivalent to the grid generation in a two dimensional domain in physical space. A second elliptic grid generation method, which generates boundary conforming grids on smooth surfaces, is presented. Concerning surface modeling, it is shown that bicubic Hermit interpolation is an excellent method to generate a smooth surface crossing a discrete set of control points.

Spekreijse, S. P.↗
NASA NTRS2009

A Module Language for Typing by Contracts

Assume-guarantee reasoning is a popular and expressive paradigm for modular and compositional specification of programs. It is becoming a fundamental concept in some computer-aided design tools for embedded system design. In this paper, we elaborate foundations for contract-based embedded system design by proposing a general-purpose module language based on a Boolean algebra allowing to define contracts. In this framework, contracts are used to negotiate the correctness of assumptions made on the definition of a component at the point where it is used and provides guarantees to its environment. We illustrate this presentation with the specification of a simplified 4-stroke engine model.

Glouche, Yann↗
NASA NTRS2016

Introduction to COFFE: The Next-Generation HPCMP CREATE-AV CFD Solver

HPCMP CREATE-AV Conservative Field Finite Element (COFFE) is a modular, extensible, robust numerical solver for the Navier-Stokes equations that invokes modularity and extensibility from its first principles. COFFE implores a flexible, class-based hierarchy that provides a modular approach consisting of discretization, physics, parallelization, and linear algebra components. These components are developed with modern software engineering principles to ensure ease of uptake from a user's or developer's perspective. The Streamwise Upwind/Petrov-Galerkin (SU/PG) method is utilized to discretize the compressible Reynolds-Averaged Navier-Stokes (RANS) equations tightly coupled with a variety of turbulence models. The mathematics and the philosophy of the methodology that makes up COFFE are presented.

Glasby, Ryan S.↗
DOE OSTI2021

Characterizing Topological Order with Matrix Product Operators

One of the most striking features of gapped quantum phases that exhibit topological order is the presence of long-range entanglement that cannot be detected by any local order parameter. The formalism of projected entangled-pair states is a natural framework for the parameterization of gapped ground state wavefunctions which allows one to characterize topological order in terms of the virtual symmetries of the local tensors that encode the wavefunction. In their most general form, these symmetries are represented by matrix product operators acting on the virtual level, which leads to a set of algebraic rules characterizing states with topological quantum order. This construction generalizes the concepts of $\mathsf G$- and twisted injectivity; the corresponding matrix product operators encode all topological features of the theory and provide a complete picture of the ground state manifold on the torus. We show how the string-net models of Levin and Wen fit within this formalism and in doing so provide a particularly intuitive interpretation of the pentagon equation for F-symbols as the pulling of matrix product operators through the string-net tensor network. Our approach paves the way to finding novel topological phases beyond string nets and elucidates the description of topological phases in terms of entanglement Hamiltonians and edge theories.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

General N=2 supersymmetric quantum mechanical model: Supervariable approach to its off-shell nilpotent symmetries

Using the supersymmetric (SUSY) invariant restrictions on the (anti-)chiral supervariables, we derive the off-shell nilpotent symmetries of the general one (0+1)-dimensional N=2 SUSY quantum mechanical (QM) model which is considered on a (1, 2)-dimensional supermanifold (parametrized by a bosonic variable t and a pair of Grassmannian variables θ and θ-bar with θ{sup 2}=(θ-bar){sup 2}=0,θ(θ-bar)+(θ-bar)θ=0). We provide the geometrical meanings to the two SUSY transformations of our present theory which are valid for any arbitrary type of superpotential. We express the conserved charges and Lagrangian of the theory in terms of the supervariables (that are obtained after the application of SUSY invariant restrictions) and provide the geometrical interpretation for the nilpotency property and SUSY invariance of the Lagrangian for the general N=2 SUSY quantum theory. We also comment on the mathematical interpretation of the above symmetry transformations. - Highlights: • A novel method has been proposed for the derivation of N=2 SUSY transformations. • General N=2 SUSY quantum mechanical (QM) model with a general superpotential, is considered. • The above SUSY QM model is generalized onto a (1, 2)-dimensional supermanifold. • SUSY invariant restrictions are imposed on the (anti-)chiral supervariables. • Geometrical meaning of the nilpotency property is provided.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗