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,549 records · Page 86
All sourcesNASA NTRSDOE OSTI
NASA NTRS1984

Nonlinear evolution of magnetopause tearing modes

Since the magnetosheath plasma is highly turbulent, reconnection at the dayside magnetopause is likely to be temporally unsteady. The tearing mode can be viewed as a model for the unsteady development of a reconnecting magnetic topology. Magnetopause tearing occurs in the guide-field limit and has a wave packet spatial structure in the east-west direction. This paper solves for the nonlinear evolution of a single wavelength guide-field tearing mode including the effects of finite transit time on the Landau resonant electrons. Short wavelength modes evolve algebraically in time with perturbation amplitudes proportional to t-squared. Long wavelength modes are fully nonlinear, and the amplitude grows linearly in time.

Coroniti, F. V.↗
DOE OSTI2025

Renormalization-group equations of the LEFT at two loops: dimension-six baryon-number-violating operators

We present the second part of a systematic calculation of the two-loop anomalous dimensions for the low-energy effective field theory below the electroweak scale (LEFT): the baryon-number-violating sector at dimension six in the power counting. We obtain the results in two different schemes: in the algebraically consistent ’t Hooft-Veltman scheme for γ 5 , corrected for evanescent as well as chiral-symmetry-breaking effects through finite renormalizations; and in naive dimensional regularization, which in the considered sector of the theory does not lead to any ill-defined γ 5 -odd traces. Our results are of interest for a reanalysis of the constraints on physics beyond the Standard Model from proton-decay searches within an EFT framework at next-to-leading-logarithmic accuracy.

Baryon/Lepton Number Violation↗
DOE OSTI2021

Electromagnetic Transient Simulation Algorithms for Evaluation of Large-Scale Extreme Fast Charging Systems (Distribution Grid Models)

The distribution and transmission grids are observing an increased penetration of power electronics in loads and generations. For example, there is increasing interest in integrating in extreme fast charging (XFC) systems for fast charging of electrical vehicles. As these systems are integrated, developing high-fidelity electromagnetic transient model of XFC systems in distribution grids and evaluating their interactions with the power grid would be of significant interest. This model will be utilized for design of XFC systems, to identify upgrades in distribution and/or transmission grids, for planning purposes by transmission planners or operators or owners, among others. It can also be utilized in operations for improved reliable performance of the grid and/or XFC station. The challenge with simulating these models is the high computational complexity introduced by the large number of states present in the system and the time-step needed to simulate the system. In this paper, advanced simulations algorithms are applied to reduce the computational complexity of simulating large-scale XFC systems. The algorithms include numerical stiffness-based segregation, time constant-based segregation, clustering and aggregation on differential algebraic equations (DAEs), and multi-order integration approaches. While the first three algorithms split the matrix that needs to be inverted from a large matrix to much smaller matrices, the final algorithm reduces the computational burden of applying higher-order integration approaches in the complete system. The comparison made in the previous sentence is with respect to use of homogeneous integration approaches used in conventional electromagnetic transient simulators like power systems computer aided design (PSCAD). The approaches mentioned here have resulted in speed-up of 36x in the simulation of a single distribution system with 15 XFCs.

Debnath, Suman↗
NASA NTRS1985

Calculation of plasma properties in ion sources

Simple algebraic equations are developed that can be used to calculate the average primary electron density, the average primary to Maxwellian density ratio, and the average Maxwellian electron temperature in the discharge chamber of a cusped magnetic field ion source. These properties may be calculated as functions of the propellant gas flow rate and the propellant utilization efficiency. Each of the plasma properties calculated using the model agrees well with the corresponding experimental data for both argon and xenon gases over a wide range of propellant flow rates and utilization efficiencies.

Brophy, J. R.↗
NASA NTRS1990

Time-optical spinup maneuvers of flexible spacecraft

Attitude controllers for spacecraft have been based on the assumption that the bodies being controlled are rigid. Future spacecraft, however, may be quite flexible. Many applications require spinning up/down these vehicles. In this work the minimum time control of these maneuvers is considered. The time-optimal control is shown to possess an important symmetry property. Taking advantage of this property, the necessary and sufficient conditions for optimality are transformed into a system of nonlinear algebraic equations in the control switching times during one half of the maneuver, the maneuver time, and the costates at the mid-maneuver time. These equations can be solved using a homotopy approach. Control spillover measures are introduced and upper bounds on these measures are obtained. For a special case these upper bounds can be expressed in closed form for an infinite dimensional evaluation model. Rotational stiffening effects are ignored in the optimal control analysis. Based on a heuristic argument a simple condition is given which justifies the omission of these nonlinear effects. This condition is validated by numerical simulation.

Singh, G.↗
NASA NTRS1982

A generalization of algebraic surface drawing

An implicit surface mathematical description of three-dimensional space is defined in terms of all points which satisfy some equation F(x, y, z) equals 0. This form is ideal for space-shaded picture drawing, where the coordinates are substituted for x and y and the equation is solved for z. A new algorithm is presented which is applicable to functional forms other than those of first- and second-order polynomial functions, such as the summation of several Gaussian density distributions. The algorithm was created in order to model electron density maps of molecular structures, but is shown to be capable of generating shapes of esthetic interest.

Blinn, J. F.↗
DOE OSTI2018

Divisible Rigid Groups. II. Stability, Saturation, and Elementary Submodels

A group G is said to be rigid if it contains a normal series G = G{sub 1} > G{sub 2} > . . . > G{sub m} > G{sub m+1} = 1, whose quotients G{sub i}/G{sub i+1} are Abelian and, treated as right ℤ[G/G{sub i}]- modules, are torsion-free. A rigid group G is divisible if elements of the quotient G{sub i}/G{sub i+1} are divisible by nonzero elements of the ring ℤ[G/G{sub i}]. Every rigid group is embedded in a divisible one. Previously, it was stated that the theory 𝔗{sub m} of divisible m-rigid groups is complete. Here, it is proved that this theory is ω-stable. Furthermore, we describe saturated models, study elementary submodels of an arbitrary model, and find a representation for a countable saturated model in the form of a limit group in the Fraïssé system of all finitely generated m-rigid groups. Also, it is proved that the theory 𝔗{sub m} admits quantifier elimination down to a Boolean combination of ∀∃-formulas.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Repartitioning the Hamiltonian in many-body second-order Brillouin–Wigner perturbation theory: Uncovering new size-consistent models

Second-order Møller-Plesset perturbation theory is well-known as a computationally inexpensive approach to the electron correlation problem that is size-consistent with a size-consistent reference but fails to be regular. On the other hand, the less well-known many-body version of Brillouin-Wigner perturbation theory has the reverse properties: it is regular but fails to be size-consistent when used with the standard MP partitioning. Consequently, its widespread use remains limited. In this work, we analyze the ways in which it is possible to use alternative non-MP partitions of the Hamiltonian to yield variants of BW2 that are size-consistent as well as regular. We show that there is a vast space of such BW2 theories and also show that it is possible to define a repartitioned BW2 theory from the ground state density alone, which regenerates the exact correlation energy. We also provide a general recipe for deriving regular, size-consistent, and size-extensive partitions from physically meaningful components, and we apply the result to small model systems. The scope of these results appears to further set the stage for a revival of BW2 in quantum chemistry.

Ab initio perturbation↗
DOE OSTI2024

Qutrit circuits and algebraic relations: A pathway to efficient spin-1 Hamiltonian simulation

Quantum information processing has witnessed significant advancements through the application of qubit-based techniques within universal gate sets. Recently, exploration beyond the qubit paradigm to d -dimensional quantum units or qudits has opened new avenues for improving computational efficiency. This paper delves into the qudit-based approach, particularly addressing the challenges presented in the high-fidelity implementation of qudit-based circuits due to increased complexity. As an innovative approach towards enhancing qudit circuit fidelity, we explore algebraic relations, such as the Yang-Baxter-like turnover equation, which may enable circuit compression and optimization. The paper introduces the turnover relation for the three-qutrit time propagator and its potential use in reducing circuit depth. We further investigate whether this relation can be generalized for higher-dimensional quantum circuits, including a focused study on the one-dimensional spin-1 Heisenberg model. Our paper outlines both rigorous and numerically efficient approaches to potentially achieve this generalization, providing a foundation for further explorations in the field of qudit-based quantum computing.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

DEGAS: Dynamic Exascale Global Address Space Programming Environments

The Dynamic, Exascale Global Address Space programming environment (DEGAS) project will develop the next generation of programming models and runtime systems to meet the challenges of Exascale computing. The Berkeley part of the project concentrated on communication-optimal code generation to optimize speed and energy efficiency by reducing data movement. Our work developed communication lower bounds, and/or communication avoiding algorithms (that either meet the lower bound, or do much less communication than their conventional counterparts) for a variety of algorithms, including linear algebra, machine learning and genomics. The Berkeley part of the project concentrated on communication-optimal code generation to optimize speed and energy efficiency by reducing data movement. Our work developed communication lower bounds, and/or communication avoiding algorithms (that either meet the lower bound, or do much less communication than their conventional counterparts) for a variety of algorithms, including linear algebra, machine learning and genomics.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Qutrit Circuits and Algebraic Relations: A Pathway to Efficient Spin-1 Hamiltonian Simulation

Quantum information processing has witnessed significant advancements through the application of qubit- based techniques within universal gate sets. Recently, exploration beyond the qubit paradigm to d-dimensional quantum units or qudits has opened new avenues for improving computational efficiency. This paper delves into the qudit-based approach, particularly addressing the challenges presented in the high-fidelity implementation of qudit-based circuits due to increased complexity. As an innovative approach towards enhancing qudit circuit fidelity, we explore algebraic relations, such as the Yang-Baxter-like turnover equation, which may enable circuit compression and optimization. The paper introduces the turnover relation for the three-qutrit time propagator and its potential use in reducing circuit depth. We further investigate whether this relation can be generalized for higher-dimensional quantum circuits, including a focused study on the one-dimensional spin-1 Heisenberg model. Our paper outlines both rigorous and numerically efficient approaches to potentially achieve this generalization, providing a foundation for further explorations in the field of qudit-based quantum computing.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2018

Deploy threading in Nalu solver stack

The goal of the ExaWind project is to enable predictive simulations of wind farms composed of many MW-scale turbines situated in complex terrain. Predictive simulations will require computational fluid dynamics (CFD) simulations for which the mesh resolves the geometry of the turbines, and captures the rotation and large deflections of blades. Whereas such simulations for a single turbine are arguably petascale class, multi-turbine wind farm simulations will require exascale-class resources. The primary code in the ExaWind project is Nalu, which is an unstructured-grid solver for the acousticallyincompressible Navier-Stokes equations, and mass continuity is maintained through pressure projection. The model consists of the mass-continuity Poisson-type equation for pressure and a momentum equation for the velocity. For such modeling approaches, simulation times are dominated by linear-system setup and solution for the continuity and momentum systems. For the ExaWind challenge problem, the moving meshes greatly affect overall solver costs as re-initialization of matrices and re-computation of preconditioners is required at every time step In this Milestone, we examine the effect of threading on the solver stack performance against flat-MPI results obtained from previous milestones using Haswell performance data full-turbine simulations. Whereas the momentum equations are solved only with the Trilinos solvers, we investigate two algebraic-multigrid preconditioners for the continuity equations: Trilinos/Muelu and HYPRE/BoomerAMG. These two packages embody smoothed-aggregation and classical Ruge-Stiiben AMG methods, respectively. In our FY18 Q2 report, we described our efforts to improve setup and solve of the continuity equations under flat-MPI parallelism. While significant improvement was demonstrated in the solve phase, setup times remained larger than expected. Starting with the optimized settings described in the Q2 report, we explore here simulation performance where OpenMP threading is employed in the solver stack. For Trilinos, threading is acheived through the Kokkos abstraction where, whereas HYPRE/BoomerAMG employs straight OpenMP. We examined results for our mid-resolution baseline turbine simulation configuration (229M DOF). Simulations on 2048 Haswell cores explored the effect of decreasing the number of MPI ranks while increasing the number of threads. Both HYPRE and Trilinos exhibited similar overal solution times, and both showed dramatic increases in simulation time in the shift from MPI ranks to OpenMP threads. This increase is attributed to the large amount of work per MPI rank starting at the single-thread configuration. Decreasing MPI ranks, while increasing threads, may be increasing simulation time due to thread synchronization and start-up overhead contributing to the latency and serial time in the model. These result showed that an MPI+OpenMP parallel decomposition will be more effective as the amount per MPI rank computation per MPI rank decreases and the communication latency increases. This idea was demonstrated in a strong scaling study of our low-resolution baseline model (29M DOF) with the Trilinos-HYPRE configuration. While MPI-only results showed scaling improvement out to about 1536 cores, engaging threading carried scaling improvements out to 4128 cores — roughly 7000 DOF per core. This is an important result as improved strong scaling is needed for simulations to be executed over sufficiently long simulated durations (i.e., for many timesteps). In addition to threading work described above, the team examined solver-performance improvements by exploring communication-overhead in the HYPRE-GMRES implementation through a communicationoptimal- GMRE algorithm (CO-GMRES), and offloading compute-intensive solver actions to GPUs. To those ends, a HYPRE mini-app was allow us to easily test different solver approaches and HYPRE parameter settings without running the entire Nalu code. With GPU acceleration on the Summitdev supercomputer, a 20x speedup was achieved for the overall preconditioner and solver execution time for the mini-app. A study on Haswell processors showed that CO-GMRES provides benefits as one increases MPI ranks.

17 WIND ENERGY↗
DOE OSTI2022

Slip velocity boundary conditions for the lattice Boltzmann modeling of microchannel flows

Abstract Slip flows in ducts are important in numerous engineering applications, most notably in microchannel flows. Compared to the standard no‐slip Dirichlet condition, the case of slip formulates as a Robin‐type condition for the fluid tangential velocity. Such an increase in mathematical complexity is accompanied by a more challenging numerical transcription. The present work concerns with this topic, addressing the modeling of the slip velocity boundary condition in the lattice Boltzmann method (LBM) applied to steady slow viscous flows inside ducts of nontrivial shapes. As novelty, we extend the newly revised local second‐order boundary (LSOB) Dirichlet fluid flow method [ Philos. Trans. R. Soc. A 378, 20190404 (2020)] to implement the slip velocity condition within the two‐relaxation‐time (TRT) framework. The LSOB follows an in‐node philosophy where its operation principle seeks to explicitly reconstruct the unknown boundary populations in the form of a third‐order accurate Chapman–Enskog expansion, where the wall slip condition is built‐in as a normal Taylor‐type condition. The key point of this approach is that the required first‐ and second‐order momentum derivatives, rather than computed through nonlocal finite difference approximations, are locally determined through a simple local linear algebra procedure, whose formulation is particularly aided by the TRT symmetry argument. To express the obtained derivatives, two approaches are considered, called and , which operate with node and wall variables, respectively. These two formulations are developed to prescribe the physical slip condition over plane and curved walls, including the corners. Their consistency and accuracy characteristics are examined against alternative linkwise strategies to impose the wall slip velocity, such as the kinetic‐based diffusive bounce‐back scheme, the central linear interpolation slip scheme, and the multireflection slip scheme. The several slip schemes are tested over different 3D microchannel configurations, with walls not conforming with the LBM uniform mesh. Numerical tests confirm the advanced accuracy characteristics of the proposed LSOB slip boundary scheme, revealing the added challenge of the wall slip modeling, and that parabolic accuracy is a necessary requirement to reach second‐order accuracy within this problem class.

Silva, Goncalo↗
DOE OSTI2023

Nonperturbative constraints from symmetry and chirality on Majorana zero modes and defect quantum numbers in (2+1) dimensions

In (1+1)D topological phases, unpaired Majorana zero modes (MZMs) can arise only if the internal symmetry group G f of the ground state splits as G f = G b × $Z^f_2$, where $Z^f_2$ is generated by fermion parity, (–1) F . In contrast, (2+1)-dimensional [(2+1)D] topological superconductors (TSC) can host unpaired MZMs at defects even when G f is not of the form G b × $Z^f_2$. In this paper we study how G f together with the chiral central charge c – strongly constrain the existence of unpaired MZMs and the quantum numbers of symmetry defects. Our results utilize a recent algebraic characterization of (2+1)D invertible fermionic topological states, which provides a non-perturbative approach based on topological quantum field theory, beyond free fermions. We study physically relevant groups such as U(1) f $\rtimes$ H, SU (2) f × H, U(2) f $\rtimes$ H, generic Abelian groups, as well as more general compact Lie groups, antiunitary symmetries and crystalline symmetries. Further, we present an algebraic formula for the fermionic crystalline equivalence principle, which gives an equivalence between states with crystalline and internal symmetries. In light of our theory, we discuss several previously proposed realizations of unpaired MZMs in TSC materials such as Sr 2 RuO 4 , transition metal dichalcogenides and iron superconductors, in which crystalline symmetries are often important; in some cases we present additional predictions for the properties of these models.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2013

Three-cluster dynamics within an ab initio framework

In this study, we introduce a fully antisymmetrized treatment of three-cluster dynamics within the ab initio framework of the no-core shell model/resonating-group method. Energy-independent nonlocal interactions among the three nuclear fragments are obtained from realistic nucleon-nucleon interactions and consistent ab initio many-body wave functions of the clusters. The three-cluster Schrödinger equation is solved with bound-state boundary conditions by means of the hyperspherical-harmonic method on a Lagrange mesh. We discuss the formalism in detail and give algebraic expressions for systems of two single nucleons plus a nucleus. Using a soft similarity-renormalization-group evolved chiral nucleon-nucleon potential, we apply the method to a 4 He+n+n description of 6 He and compare the results to experiment and to a six-body diagonalization of the Hamiltonian performed within the harmonic-oscillator expansions of the no-core shell model. Differences between the two calculations provide a measure of core ( 4 He) polarization effects.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

A unified framework of stabilized finite element method for solving the Boltzmann transport equation

This paper presents a unified framework of stabilized finite element method for solving the Boltzmann transport equation. Unlike the traditional Petrov-Galerkin finite element method which modifies the test function to construct the stabilization term, we derive the stabilization methods from the standard Galerkin weak form with Sub-grid scale model. The basic idea of this method is to decompose the unknowns into its numerical solution and residual, with an approximation for the residual and embeds it in the Galerkin weak form to yield a stabilized variational formula. Different approximations of the residual lead to different stabilization methods, all the frequently used stabilized methods, including the Streamline Upwinding Petrov-Galerkin (SUPG) method, Galerkin/Least-Square (GLS) method, and Algebraic Sub-Grid Scale (ASGS) method can be obtained from this framework. The similarities and differences of the different approximations are compared in this paper. The numerical results show that the behaviors of the different methods area similar with the same stabilization parameter, and all these stabilized techniques can obtain a correct and stable solution. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2016

Muon Beam Tracking and Spin-Orbit Correlations for Precision g-2 Measurements

The main goal of the Muon g-2 Experiment (g-2) at Fermilab is to measure the muon anomalous magnetic moment to unprecedented precision. This new measurement will allow to test the completeness of the Standard Model (SM) and to validate other theoretical models beyond the SM. The close interplay of the understanding of particle beam dynamics and the preparation of the beam properties with the experimental measurement is tantamount to the reduction of systematic errors in the determination of the muon anomalous magnetic moment. We describe progress in developing detailed calculations and modeling of the muon beam delivery system in order to obtain a better understanding of spin-orbit correlations, nonlinearities, and more realistic aspects that contribute to the systematic errors of the g-2 measurement. Our simulation is meant to provide statistical studies of error effects and quick analyses of running conditions for when g-2 is taking beam, among others. We are using COSY, a differential algebra solver developed at Michigan State University that will also serve as an alternative to compare results obtained by other simulation teams of the g-2 Collaboration.

43 PARTICLE ACCELERATORS↗
DOE OSTI2016

Sanibel Symposium in the Petascale-Exascale Computational Era

The 56 th Sanibel Symposium was held February 14-19 2016 at the King and Prince Hotel, St. Simons Island, GA. It successfully brought quantum chemists and chemical and condensed matter physicists together in presentations, posters, and informal discussions bridging those two communities. The Symposium has had a significant role in preparing generations of quantum theorists. As computational potency and algorithmic sophistication have grown, the Symposium has evolved to emphasize more heavily computationally oriented method development in chemistry and materials physics, including nanoscience, complex molecular phenomena, and even bio-molecular methods and problems. Given this context, the 56 th Sanibel meeting systematically and deliberately had sessions focused on exascale computation. A selection of outstanding theoretical problems that need serious attention was included. Five invited sessions, two contributed sessions (hot topics), and a poster session were organized with the exascale theme. This was a historic milestone in the evolution of the Symposia. Just as years ago linear algebra, perturbation theory, density matrices, and band-structure methods dominated early Sanibel Symposia, the exascale sessions of the 56 th meeting contributed a transformative influence to add structure and strength to the computational physical science community in an unprecedented way. A copy of the full program of the 56 th Symposium is attached. The exascale sessions were Linear Scaling, Non-Adabatic Dynamics, Interpretive Theory and Models, Computation, Software, and Algorithms, and Quantum Monte Carlo. The Symposium Proceedings will be published in Molecular Physics (2017). Note that the Sanibel proceedings from 2015 and 2014 were published as Molecular Physics vol. 114, issue 3-4 (2016) and vol. 113, issue 3-4 (2015) respectively.

36 MATERIALS SCIENCE↗