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,675 records · Page 93
All sourcesNASA NTRSDOE OSTI
NASA NTRS2015

In-Flight Aeroelastic Stability of the Thermal Protection System on the NASA HIAD, Part II: Nonlinear Theory and Extended Aerodynamics

Conical shell theory and a supersonic potential flow aerodynamic theory are used to study the nonlinear pressure buckling and aeroelastic limit cycle behavior of the thermal protection system for NASA's Hypersonic Inflatable Aerodynamic Decelerator. The structural model of the thermal protection system consists of an orthotropic conical shell of the Donnell type, resting on several circumferential elastic supports. Classical Piston Theory is used initially for the aerodynamic pressure, but was found to be insufficient at low supersonic Mach numbers. Transform methods are applied to the convected wave equation for potential flow, and a time-dependent aerodynamic pressure correction factor is obtained. The Lagrangian of the shell system is formulated in terms of the generalized coordinates for all displacements and the Rayleigh-Ritz method is used to derive the governing differential-algebraic equations of motion. Aeroelastic limit cycle oscillations and buckling deformations are calculated in the time domain using a Runge-Kutta method in MATLAB. Three conical shell geometries were considered in the present analysis: a 3-meter diameter 70 deg. cone, a 3.7-meter 70 deg. cone, and a 6-meter diameter 70 deg. cone. The 6-meter configuration was loaded statically and the results were compared with an experimental load test of a 6-meter HIAD. Though agreement between theoretical and experimental strains was poor, the circumferential wrinkling phenomena observed during the experiments was captured by the theory and axial deformations were qualitatively similar in shape. With Piston Theory aerodynamics, the nonlinear flutter dynamic pressures of the 3-meter configuration were in agreement with the values calculated using linear theory, and the limit cycle amplitudes were generally on the order of the shell thickness. The effect of axial tension was studied for this configuration, and increasing tension was found to decrease the limit cycle amplitudes when the circumferential elastic supports were neglected, but resulted in more complex behavior when the supports were included. The nominal flutter dynamic pressure of the 3.7-meter configuration was significantly lower than that of the 3-meter, and it was found that two sets of natural modes coalesce to flutter modes near the same dynamic pressure. This resulted in a significant drop in the limit cycle frequencies at higher dynamic pressures, where the flutter mode with the lower frequency becomes more critical. Pre-buckling pressure loads and the aerodynamic pressure correction factor were studied for all geometries, and these effects resulted in significantly lower flutter boundaries compared with Piston Theory alone. The maximum dynamic pressure predicted by aerodynamic simulations of a proposed 3.7-meter HIAD vehicle was still lower than any of the calculated flutter dynamic pressures, suggesting that aeroelastic effects for this vehicle are of little concern.

Goldman, Benjamin D.↗
DOE OSTI2025

A Simple, Scalable Large Deformation Solid Mechanics Implementation in the MOOSE Framework

This article describes a large deformation solid mechanics solver implemented as part of the freely available and open source MOOSE finite element simulation framework. The article documents the choices made in developing the solid mechanics framework and describes novel formulations for the gradient operator and constitutive modeling framework made to simplify implementations of different coordinate systems, stabilized gradient operators, and different constitutive model inputs and outputs. In the process, the article describes a new formulation that casts objective integration of the Cauchy stress as a linear transformation of the small stress rate. Finally, the article presents key implementation details and examines the parallel efficiency of the solid mechanics solver implemented in MOOSE. The implementation retains a good weak scaling efficiency beyond 1,000 parallel processes. The article includes a discussion of the factors limiting the parallel efficiency of implicit, large deformation solid mechanics codes on current high-performance computers, with the main current limitation being the scalability of the algebraic multigrid methods used to solve the linearized equilibrium equations.

Applied computing → Computer-aided design↗
DOE OSTI2018

Local gravity theories in conformal superspace

The conformal method for the initial value formulation of GR is one of the most powerful and widely usedtools in both analytical studies and numerical simulations of solutions. As is well-known, it exploits a hidden spatial conformal symmetry that can be extracted from the equations. But why conformal symmetry? Motivated by these questions and by well-known obstacles to quantum gravity, I look for the most general geometrodynamical symmetries compatible with a reduced physical configuration space for metric gravity. I argue that they are indeed spatial conformal diffeomorphisms. The next question is: what sort of theories most naturally incorporate this symmetry? Demanding locality for an action for metric gravity compatible with these principles determines the allowed operators, both for the purely gravitational part as well as matter couplings. The symmetries guarantee that there are two gravitational propagating physical degrees of freedom, but no explicit refoliation-invariance. The simplest such system has a geometric interpretation as a geodesic model in infinite-dimensional conformal superspace. One example of solution to the equations of motion corresponds to a static Bianchi IX spatial ansatz. The unique coupling to electromagnetism forces the electromagnetic equations to be hyperbolic, enabling us to “build” a standard space–time causal structure. There are, however, deviations from the standard Maxwell equations when space–time anisotropies become too large. Regarding quantization, with the geometric interpretation and the lack of refoliation invariance, the path integral treatment of the symmetries becomes much less involved than the similar approaches to GR. The symmetries form an (infinite-dimensional) Lie algebra, and no BFV treatment is necessary. Moreover, one can use gauge-invariant variational principles for selecting the boundary conditions of the path integral. We find that the propagator around the homogeneous solution has up to 6th order spatial derivatives, giving it plausible regularization properties as in Horava–Lifschitz.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Dynamic mode decomposition with core sketch

With the increase in collected data volumes, either from experimental measurements or high fidelity simulations, there is an ever-growing need to develop computationally efficient tools to process, analyze, and interpret these datasets. Modal analysis techniques have gained great interest due to their ability to identify patterns in the data and extract valuable information about the system being considered. Dynamic mode decomposition (DMD) relies on elements of the Koopman approximation theory to compute a set of modes, each associated with a fixed oscillation frequency and a decay/growth rate. Extracting these details from large datasets can be computationally expensive due to the need to implement singular value decomposition of the input data matrix. Sketching algorithms have become popular in numerical linear algebra where statistical theoretic approaches are utilized to reduce the cost of major operations. A sketch of a matrix is another matrix, which is significantly smaller, but still sufficiently approximates the original system. We put forth an efficient DMD framework, SketchyDMD, based on a core sketching algorithm that captures information about the range and corange (their mutual relationship) of input data. The proposed sketching-based framework can accelerate various portions of the DMD routines, compared to classical methods that operate directly on the raw input data. We conduct numerical experiments using the spherical shallow water equations as a prototypical model in the context of geophysical flows. In conclusion, we show that the proposed SketchyDMD is superior to existing randomized DMD methods that are based on capturing only the range of the input data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1999

Quasi-Static Viscoelastic Finite Element Model of an Aircraft Tire

An elastic large displacement thick-shell mixed finite element is modified to allow for the calculation of viscoelastic stresses. Internal strain variables are introduced at the element's stress nodes and are employed to construct a viscous material model. First order ordinary differential equations relate the internal strain variables to the corresponding elastic strains at the stress nodes. The viscous stresses are computed from the internal strain variables using viscous moduli which are a fraction of the elastic moduli. The energy dissipated by the action of the viscous stresses is included in the mixed variational functional. The nonlinear quasi-static viscous equilibrium equations are then obtained. Previously developed Taylor expansions of the nonlinear elastic equilibrium equations are modified to include the viscous terms. A predictor-corrector time marching solution algorithm is employed to solve the algebraic-differential equations. The viscous shell element is employed to computationally simulate a stair-step loading and unloading of an aircraft tire in contact with a frictionless surface.

Johnson, Arthur R.↗
DOE OSTI2015

A transmission line model for propagation in elliptical core optical fibers

The calculation of mode propagation constants of elliptical core fibers has been the purpose of extended research leading to many notable methods, with the classic step index solution based on Mathieu functions. This paper seeks to derive a new innovative method for the determination of mode propagation constants in single mode fibers with elliptic core by modeling the elliptical fiber as a series of connected coupled transmission line elements. We develop a matrix formulation of the transmission line and the resonance of the circuits is used to calculate the mode propagation constants. The technique, used with success in the case of cylindrical fibers, is now being extended for the case of fibers with elliptical cross section. The advantage of this approach is that it is very well suited to be able to calculate the mode dispersion of arbitrary refractive index profile elliptical waveguides. The analysis begins with the deployment Maxwell’s equations adjusted for elliptical coordinates. Further algebraic analysis leads to a set of equations where we are faced with the appearance of harmonics. Taking into consideration predefined fixed number of harmonics simplifies the problem and enables the use of the resonant circuits approach. According to each case, programs have been created in Matlab, providing with a series of results (mode propagation constants) that are further compared with corresponding results from the ready known Mathieu functions method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2012

Electric Dipole Transitions Within The Ab initio No-Core Shell Model With Continuum

The no-core shell model with continuum (NCSMC) offers a unified approach to nuclear bound and continuum states based on the coupling of the no-core shell model (NCSM), a bound-state tech nique, with the no-core shell model/resonating group method (NCSM/RGM), a nuclear scattering technique. This new ab initio method leads to convergence properties superior to either NCSM or NCSM/RGM while providing a balanced approach to different classes of states. The NCSMC many nucleon wave function contains: i) a pre-diagonalized A-nucleon NCSM component expanded in the harmonic oscillator basis; ii) a binary-cluster NCSM/RGM component with asymptotic boundary conditions that properly describe weakly bound states, resonances and scattering. Here, we review the main aspects of the NCSM formalism, and provide algebraic expression for the reduced matrix elements of the electric dipole operator, needed to calculate E1 transition probabilities, radiative captures and photodisintegration cross sections.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1984

Advanced numerics for multi-dimensional fluid flow calculations

In recent years, there has been a growing interest in the development and use of mathematical models for the simulation of fluid flow, heat transfer and combustion processes in engineering equipment. The equations representing the multi-dimensional transport of mass, momenta and species are numerically solved by finite-difference or finite-element techniques. However despite the multiude of differencing schemes and solution algorithms, and the advancement of computing power, the calculation of multi-dimensional flows, especially three-dimensional flows, remains a mammoth task. The following discussion is concerned with the author's recent work on the construction of accurate discretization schemes for the partial derivatives, and the efficient solution of the set of nonlinear algebraic equations resulting after discretization. The present work has been jointly supported by the Ramjet Engine Division of the Wright Patterson Air Force Base, Ohio, and the NASA Lewis Research Center.

Vanka, S. P.↗
NASA NTRS1994

A high-speed linear algebra library with automatic parallelism

Parallel or distributed processing is key to getting highest performance workstations. However, designing and implementing efficient parallel algorithms is difficult and error-prone. It is even more difficult to write code that is both portable to and efficient on many different computers. Finally, it is harder still to satisfy the above requirements and include the reliability and ease of use required of commercial software intended for use in a production environment. As a result, the application of parallel processing technology to commercial software has been extremely small even though there are numerous computationally demanding programs that would significantly benefit from application of parallel processing. This paper describes DSSLIB, which is a library of subroutines that perform many of the time-consuming computations in engineering and scientific software. DSSLIB combines the high efficiency and speed of parallel computation with a serial programming model that eliminates many undesirable side-effects of typical parallel code. The result is a simple way to incorporate the power of parallel processing into commercial software without compromising maintainability, reliability, or ease of use. This gives significant advantages over less powerful non-parallel entries in the market.

Boucher, Michael L.↗
DOE OSTI2021

A Segregated Approach for Modeling the Electrochemistry in the 3-D Microstructure of Li-Ion Batteries and Its Acceleration Using Block Preconditioners

Abstract Battery performance is strongly correlated with electrode microstructure. Electrode materials for lithium-ion batteries have complex microstructure geometries that require millions of degrees of freedom to solve the electrochemical system at the microstructure scale. A fast-iterative solver with an appropriate preconditioner is then required to simulate large representative volume in a reasonable time. In this work, a finite element electrochemical model is developed to resolve the concentration and potential within the electrode active materials and the electrolyte domains at the microstructure scale, with an emphasis on numerical stability and scaling performances. The block Gauss-Seidel (BGS) numerical method is implemented because the system of equations within the electrodes is coupled only through the nonlinear Butler–Volmer equation, which governs the electrochemical reaction at the interface between the domains. The best solution strategy found in this work consists of splitting the system into two blocks—one for the concentration and one for the potential field—and then performing block generalized minimal residual preconditioned with algebraic multigrid, using the FEniCS and the Portable, Extensible Toolkit for Scientific Computation libraries. Significant improvements in terms of time to solution (six times faster) and memory usage (halving) are achieved compared with the MUltifrontal Massively Parallel sparse direct Solver. Additionally, BGS experiences decent strong parallel scaling within the electrode domains. Last, the system of equations is modified to specifically address numerical instability induced by electrolyte depletion, which is particularly valuable for simulating fast-charge scenarios relevant for automotive application.

25 ENERGY STORAGE↗
DOE OSTI2021

Tusas: A fully implicit parallel approach for coupled phase-field equations

In this study, we develop a fully-coupled, fully-implicit approach for phase-field modeling of solidification in metals and alloys. Predictive simulation of solidification in pure metals and metal alloys remains a significant challenge in the field of materials science, as microstructure formation during the solidification process plays a critical role in the properties and performance of the solid material. Our simulation approach consists of a finite element spatial discretization of the fully-coupled nonlinear system of partial differential equations at the microscale, which is treated implicitly in time with a preconditioned Jacobian-free Newton-Krylov method. The approach is algorithmically scalable as well as efficient due to an effective preconditioning strategy based on algebraic multigrid and block factorization. We implement this approach in the open-source Tusas framework, which is a general, flexible tool developed in C++ for solving coupled systems of nonlinear partial differential equations. The performance of our approach is analyzed in terms of algorithmic scalability and efficiency, while the computational performance of Tusas is presented in terms of parallel scalability and efficiency on emerging heterogeneous architectures. We demonstrate that modern algorithms, discretizations, and computational science, and heterogeneous hardware provide a robust route for predictive phase-field simulation of microstructure evolution during additive manufacturing.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1992

A non-linearly stable implicit finite element algorithm for hypersonic aerodynamics

A generalized curvilinear coordinate Taylor weak statement implicit finite element algorithm is developed for the two-dimensional and axisymmetric compressible Navier-Stokes equations for ideal and reacting gases. For accurate hypersonic simulation, air is modeled as a mixture of five perfect gases, i.e., molecular and atomic oxygen and nitrogen as well as nitric oxide. The associated pressure is then determined via Newton solution of the classical chemical equilibrium equation system. The directional semidiscretization is achieved using an optimal metric data Galerkin finite element weak statement, on a developed 'companion conservation law system', permitting classical test and trial space definitions. Utilizing an implicit Runge-Kutta scheme, the terminal algorithm is then nonlinearly stable, and second-order accurate in space and time on arbitrary curvilinear coordinates. Subsequently, a matrix tensor product factorization procedure permits an efficient numerical linear algebra handling for large Courant numbers. For ideal- and real-gas hypersonic flows, the algorithm generates essentially nonoscillatory numerical solutions in the presence of strong detached shocks and boundary layer-inviscid flow interactions.

Iannelli, G. S.↗
DOE OSTI2015

Superscaling in electron-nucleus scattering and its link to CC and NC QE neutrino-nucleus scattering

The superscaling approach (SuSA) to neutrino-nucleus scattering, based on the assumed universality of the scaling function for electromagnetic and weak interactions, is reviewed. The predictions of the SuSA model for bot CC and NC differential and total cross sections are presented and compared with the MiniBooNE data. The role of scaling violations, in particular the contribution of meson exchange currents in the two-particle two-hole sector, is explored.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1985

A remark about pointed bubbles

The polymer expansion is a formal algebraic identity between a partition function and logarithm in statistical physics problems. The expansion gives a systematic method to control the free energy or to establish exponential tree-graph decay of connected correlations. Here, the convergence properties of the polymer expansion are analyzed in connection with three practical examples, including: intersecting bonds in chemical polymer chains; a connected closed hypersurface built from the (d-1)-faces of the d-dimensional unit cubes; and the set of Feynamn diagrams in the perturbation series of the Euclidean field theory partition function Z. The example of connected polymer chains is generalized to apply to other lattice models, including n-state Ising models at high temperature; short range lattice gases at high temperature; and weak coupling lattice field and gauge theories.

Garabedian, P. R.↗
NASA NTRS1989

Stable Kalman filters for processing clock measurement data

Kalman filters have been used for some time to process clock measurement data. Due to instabilities in the standard Kalman filter algorithms, the results have been unreliable and difficult to obtain. During the past several years, stable forms of the Kalman filter have been developed, implemented, and used in many diverse applications. These algorithms, while algebraically equivalent to the standard Kalman filter, exhibit excellent numerical properties. Two of these stable algorithms, the Upper triangular-Diagonal (UD) filter and the Square Root Information Filter (SRIF), have been implemented to replace the standard Kalman filter used to process data from the Deep Space Network (DSN) hydrogen maser clocks. The data are time offsets between the clocks in the DSN, the timescale at the National Institute of Standards and Technology (NIST), and two geographically intermediate clocks. The measurements are made by using the GPS navigation satellites in mutual view between clocks. The filter programs allow the user to easily modify the clock models, the GPS satellite dependent biases, and the random noise levels in order to compare different modeling assumptions. The results of this study show the usefulness of such software for processing clock data. The UD filter is indeed a stable, efficient, and flexible method for obtaining optimal estimates of clock offsets, offset rates, and drift rates. A brief overview of the UD filter is also given.

Clements, P. A.↗
DOE OSTI2022

Nucleon-pair truncation of the shell model for medium-heavy nuclei

Background: Establishing computationally tractable models of atomic nuclei is a long-time goal of nuclear structure physics. A flexible framework which easily includes excited states and many-body correlations is the configuration-interaction shell model (SM), but the exponential growth of the basis means one needs an efficient truncation scheme, ideally one that includes both deformation and pairing correlations. Purpose: Here, we propose an efficient truncation scheme of the SM: starting from a pair condensate variationally defined by Hartree-Fock single-particle states and the particle-number conserved Bardeen-Cooper-Schrieffer (NBCS) approximation, we carry out projection of states with good angular momentum. Methods: After generating Hartree-Fock single-particle states with Kramers degeneracy in a SM space, we optimize the pair amplitudes in the NBCS by minimizing the energy, and then use linear algebra projection (LAP) of states with good angular momentum. Both NBCS and LAP are computationally fast. Results: Our calculations yield good agreement with full configuration-interaction SM calculations for low-lying states of transitional and rotational nuclei with axially symmetric and triaxial deformation in medium- and heavy-mass regions: 44,46,48 Ti, 48,50 Cr, 52 Fe, 60,62,64 Zn, 66,68 Ge, 68 Se, and 108,110 Xe. We predict low-lying states of 112–114 Ba and 116–120 Ce, nuclei difficult to reach by large-scale SM calculations. Conclusions: Both pair correlation and the configuration mixing between different intrinsic states play a key role in reproducing collectivity and shape coexistence, demonstrating the utility of this truncation scheme of the SM to study transitional and deformed nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

Numerical approach based on fractional-order Lagrange polynomials for solving a class of fractional differential equations

In this study, we propose a new set of fractional functions based on the Lagrange polynomials to solve a class of fractional differential equations. Fractional differential equations are the best tools for modelling natural phenomenon that are elaborated by fractional calculus. Therefore, we need an accurate and efficient technique for solving them. The main purpose of this article is to generalize new functions based on Lagrange polynomials to the fractional calculus. At first, we present a new representation of Lagrange polynomials and in continue, we propose a new set of fractional-order functions which are called fractional-order Lagrange polynomials (FLPs). Besides, a general formulation for operational matrices of fractional integration and derivative of FLPs on arbitrary nodal points are extracted. These matrices are obtained using Laplace transform. The initial value problems is reduced to the system of algebraic equations using the operational matrix of fractional integration and collocation method. Also, we find the upper bound of error vector for the fractional integration operational matrix and we indicate convergence of approximations of FLPs. Illustrative examples are included to demonstrate the validity and applicability of the proposed technique.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Lagrangian conditional statistics and flow topology in edge plasma turbulence

Lagrangian statistics and particle transport in edge plasma turbulence are investigated using the Hasegawa–Wakatani model and its modified version. The latter shows the emergence of pronounced zonal flows. Different values of the adiabaticity parameter are considered. The main goal is to characterize the role of coherent structures, i.e., vortices and zonal flows, and their impact on the Lagrangian statistics of particles. Computationally intensive long time simulations following ensembles of test particles over hundreds of eddy turnover times are considered in statistically stationary turbulent flows. The flow topology is characterized using the Lagrangian Okubo–Weiss criterion in order to split the flow into topologically different domains. In elliptic and hyperbolic regions, the probability density functions (PDFs) of the residence time have self-similar algebraic decaying tails. However, in the intermediate regions, the PDFs exhibit exponentially decaying tails. Topologically conditioned PDFs of the Lagrangian velocity, and acceleration and density fluctuations are likewise computed. The differences between the classical Hasegawa–Wakatani system and its modified version are assessed, and the role of zonal flows is highlighted. Here, the density flux spectrum, which characterizes the contributions of different length scales, is studied, and its inertial scaling is found to be in agreement with predictions based on dimensional arguments. Analyzing the angular change of particle tracers at different time scales, corresponding to coarse grained curvature, completes the study, and these multiscale geometric statistics quantify the directional properties of the particle motion in different flow regimes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗