SEARCH · Engineering Papers

Results for “Numerical stability”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 991 records · Page 55

All sources NASA NTRS DOE OSTI

NASA NTRS1999

Detailed Multidimensional Simulations of the Structure and Dynamics of Flames

Numerical simulations in which the various physical and chemical processes can be independently controlled can significantly advance our understanding of the structure, stability, dynamics and extinction of flames. Therefore, our approach has been to use detailed time-dependent, multidimensional, multispecies numerical models to perform carefully designed computational experiments of flames on Earth and in microgravity environments. Some of these computational experiments are complementary to physical experiments performed under the Microgravity Program while others provide a fundamental understanding that cannot be obtained from physical experiments alone. In this report, we provide a brief summary of our recent research highlighting the contributions since the previous microgravity combustion workshop. There are a number of mechanisms that can cause flame instabilities and result in the formation of dynamic multidimensional structures. In the past, we have used numerical simulations to show that it is the thermo-diffusive instability rather than an instability due to preferential diffusion that is the dominant mechanism for the formation of cellular flames in lean hydrogen-air mixtures. Other studies have explored the role of gravity on flame dynamics and extinguishment, multi-step kinetics and radiative losses on flame instabilities in rich hydrogen-air flames, and heat losses on burner-stabilized flames in microgravity. The recent emphasis of our work has been on exploring flame-vortex interactions and further investigating the structure and dynamics of lean hydrogen-air flames in microgravity. These topics are briefly discussed after a brief discussion of our computational approach for solving these problems.

Patnaik, G.↗

NASA NTRS1976

The dynamics of spin stabilized spacecraft with movable appendages, part 2

The dynamics and stability of a spin stabilized spacecraft with a hinged appendage system are treated analytically and numerically. The hinged system consists of a central hub with masses attached to (assumed) massless booms of fixed length whose orientation relative to the main part can change. The general three dimensional deployment dynamics of such a hinged system is considered without any restriction on the location of the hinge points. The equations of motion for the hinged system, with viscous damping at both hinge points, are linearized about the nominal equilibrium position where the booms are orthogonal to the nominal spin axis for the case of two dimensional and three dimensional motion. Analytic stability criteria are obtained from the necessary condition on the sign of all the coefficients in the system characteristic equation.

Bainum, P. M.↗

DOE OSTI2024

Numerical integration in the virtual element method with the scaled boundary cubature scheme

Abstract The virtual element method (VEM) is a stabilized Galerkin method on meshes that consist of arbitrary (convex and nonconvex) polygonal and polyhedral elements. A crucial ingredient in the implementation of low‐ and high‐order VEM is the numerical integration of monomials and nonpolynomial functions over such elements. In this article, we apply the recently proposed scaled boundary cubature (SBC) scheme to compute the weak form integrals in various virtual element formulations over polygonal and polyhedral meshes. In doing so, we demonstrate the flexibility of the approach and the accuracy that it delivers on a broad suite of boundary‐value problems in 2D and 3D over polytopes with affine faces as well as on elements with curved boundaries. In addition, the use of the SBC scheme is exemplified in an enriched Poisson formulation of the VEM in which weakly singular functions are required to be integrated. This study establishes the SBC method as a simple, accurate and efficient integration scheme for use in the VEM.

Chin, Eric B.↗

NASA NTRS1996

The Effects of Dissipation and Coarse Grid Resolution for Multigrid in Flow Problems

The objective of this paper is to investigate the effects of the numerical dissipation and the resolution of the solution on coarser grids for multigrid with the Euler equation approximations. The convergence is accomplished by multi-stage explicit time-stepping to steady state accelerated by FAS multigrid. A theoretical investigation is carried out for linear hyperbolic equations in one and two dimensions. The spectra reveals that for stability and hence robustness of spatial discretizations with a small amount of numerical dissipation the grid transfer operators have to be accurate enough and the smoother of low temporal accuracy. Numerical results give grid independent convergence in one dimension. For two-dimensional problems with a small amount of numerical dissipation, however, only a few grid levels contribute to an increased speed of convergence. This is explained by the small numerical dissipation leading to dispersion. Increasing the mesh density and hence making the problem over resolved increases the number of mesh levels contributing to an increased speed of convergence. If the steady state equations are elliptic, all grid levels contribute to the convergence regardless of the mesh density.

Eliasson, Peter↗

NASA NTRS1981

Hydrostatic calculations of axisymmetric flow and its stability for the AGCE model

Baroclinic waves in the atmospherics general circulation experiment (AGCE) apparatus by the use of numerical hydrostatic primitive equation models were determined. The calculation is accomplished by using an axisymmetric primitive equation model to compute, for a given set of experimental parameters, a steady state axisymmetric flow and then testing this axisymmetric flow for stability using a linear primitive equation model. Some axisymmetric flows are presented together with preliminary stability calculations.

Miller, T. L.↗

DOE OSTI2020

First-principles Modeling and Design of Solid-State Interfaces for the Protection and Use of Lithium Metal Anodes

Li-ion batteries are one of the most advanced energy storage technologies in use today. Li-ion batteries are used in a multitude of applications ranging from consumer electronics, medical devices, sensors and grid storage. However, improving the capacity and energy density delivered by current Li-ion technology requires advanced materials research into novel chemical systems. In this project we have focused particularly in the use of solid-state electrolytes with lithium metal electrodes. Research into all solid-state batteries (ASSB) with Li metal electrodes has significantly expanded in recent years, however most studies reported experimental findings, which left substantial room for theoretical and modeling work as a tool to understand and determine design principles allowing reliable and safe use of ASSBs with Li metal. Among the remaining obstacles preventing reliable use of ASSBs with a Li metal electrode, the stability of the interface between the solid electrolyte and Li metal, and the propagation/dendrite formation of Li metal and resulting mechanical degradation of the electrolyte are key phenomenon that are yet to be fully understood. In the current project we have addressed these two coupled phenomena using first principles calculations and mesoscale continuum modeling. We have obtained chemical and electrochemical stability windows for several solid electrolyte materials. Additionally, from mathematical and numerical modeling of Li protrusion and dendrite initiation during plating and stripping we have determined design criteria in terms of chemical, electrochemical, and mechanical properties and operating conditions for which stable deposition can occur. We also considered the effects of mixed electronic-ionic conduction in solid electrolytes, which has more recently been suggested as another important mechanism involved in ASSB failure. Throughout our work we have successfully addressed important questions necessary for the use of ASSB’s. We have determined guiding principles for materials properties and operating conditions necessary to operate ASSB’s. And have proposed novel solid electrolyte materials with predicted chemical stability an ionic conductivity. Although this represents significant progress in our understanding, open questions remain in order to fully develop reliable and safely operate ASSBs with Li metal. Future work, building on this project will require further experimental, theoretical and simulation efforts to address remaining questions.

25 ENERGY STORAGE↗

NASA NTRS1984

A mathematical model of the UH-60 helicopter

This report documents the revisions made to a ten-degree-of-freedom, full-flight envelope, generic helicopter mathematical model to represent the UH-60 helicopter accurately. The major modifications to the model include fuselage aerodynamic force and moment equations specific to the UH-60, a canted tail rotor, a horizontal stabilator with variable incidence, and a pitch bias actuator (PBA). In addition, this report presents a full set of parameters and numerical values which describe the helicopter configuration and physical characteristics. Model validation was accomplished by comparison of trim and stability derivative data generated from the UH-60 math model with data generated from a similar total force and moment math model.

Hilbert, K. B.↗

DOE OSTI2025

An impermeable copper surface monolayer with high-temperature oxidation resistance

Despite numerous efforts involving surface coating, doping, and alloying, maintaining surface stability of metal at high temperatures without compromising intrinsic properties has remained challenging. Here, we present a pragmatic method to address the accelerated oxidation of Cu, Ni, and Fe at temperatures exceeding 200 °C. Inspired by the concept that oxygen (O) itself can effectively obstruct the pathway of O infiltration, this study proposes the immobilization of O on the metal surface. Through extensive calculations considering various elements (C, Al, Si, Ge, Ga, In, and Sn) to anchor O on Cu surfaces, Si emerges as the optimal element. The theoretical findings are validated through systematic sputtering deposition experiments. The introduction of anchoring elements to reinforce Cu–O bonds enables the formation of an atomically thin barrier on the Cu surface, rendering it impermeable to O even at high temperatures (400 °C) while preserving its intrinsic conductivity. This oxidation resistance, facilitated by the impermeable atomic monolayer, opens promising opportunities for researchers and industries to overcome limitations associated with the use of oxidizable metal films.

36 MATERIALS SCIENCE↗

DOE OSTI2025

On resistive interchange, double tearing, and resonant and non-resonant infernal modes in spherical tokamak negative central shear discharges

We examine the linear and nonlinear stability of a sequence of reversed shear NSTX equilibria with the same toroidal current density and pressure profile but with different toroidal field strengths. All equilibria have two q = 2 surfaces and have q > 1 everywhere. For this sequence, which all have β P ∼ 0.5, as the minimum value of q decreases below about 1.4, the localized resistive interchange criterion is strongly violated, and the unstable mode with toroidal mode number n=1 changes its dominant poloidal mode number m from a (m,n)=(2,1) double tearing mode to a non-resonant (1,1) infernal mode. This (1,1) mode nonlinearly flattens both the current density and the pressure near the magnetic axis. The higher-n unstable modes are all resistive infernal modes localized around the surface where q has a minimum and the magnetic shear vanishes. They generally saturate nonlinearly at low amplitude but can cause magnetic surface breakup near the axis.

Spitzer resistivity↗

DOE OSTI2020

Ground State and Hidden Symmetry of Magic-Angle Graphene at Even Integer Filling

In magic angle twisted bilayer graphene (TBG), electron-electron interactions play a central role, resulting in correlated insulating states at certain integer fillings. Identifying the nature of these insulators is a central question, and it is potentially linked to the relatively high-temperature superconductivity observed in the same devices. Here, we address this question using a combination of analytical strong-coupling arguments and a comprehensive Hartree-Fock numerical calculation, which includes the effect of remote bands. The ground state we obtain at charge neutrality is an unusual ordered state, which we call the Kramers intervalley-coherent (K-IVC) insulator. In its simplest form, the K-IVC order exhibits a pattern of alternating circulating currents that triples the graphene unit cell, leading to an “orbital magnetization density wave.” Although translation and time-reversal symmetry are broken, a combined “Kramers” time-reversal symmetry is preserved. Our analytic arguments are built on first identifying an approximate U⁡(4) × U⁡(4) symmetry, resulting from the remarkable properties of the TBG band structure, which helps select a low-energy manifold of states that are further split to favor the K-IVC state. This low-energy manifold is also found in the Hartree-Fock numerical calculation. We show that symmetry-lowering perturbations can stabilize other insulators and the semimetallic state, and we discuss the ground state at half-filling and give a comparison with experiments.

Bultinck, Nick↗

DOE OSTI2023

Operating-Envelopes-Aware Decentralized Welfare Maximization for Energy Communities

We propose an operating-envelope-aware, prosumer-centric, and efficient energy community that aggregates individual and shared community distributed energy resources and transacts with a regulated distribution system operator (DSO) under a generalized net energy metering tariff design. To ensure safe network operation, the DSO imposes dynamic export and import limits, known as dynamic operating envelopes, on end-users' revenue meters. Given the operating envelopes, we propose an incentive-aligned community pricing mechanism under which the decentralized optimization of community members' benefit implies the optimization of overall community welfare. The proposed pricing mechanism satisfies the cost-causation principle and ensures the stability of the energy community in a coalition game setting. Numerical examples provide insights into the characteristics of the proposed pricing mechanism and quantitative measures of its performance.

distributed energy resources aggregation↗

DOE OSTI2025

Score-Based Physics-Informed Neural Networks for High-Dimensional Fokker–Planck Equations

The Fokker-Planck (FP) equation is a foundational partial differential equation (PDE) in stochastic processes involving Brownian motions. However, the curse of dimensionality (CoD) poses a formidable challenge when dealing with high-dimensional FP equations. Although Monte Carlo simulation and (vanilla) Physics-Informed Neural Networks (PINNs) have shown the potential to tackle CoD, both methods exhibit significant numerical errors in high dimensions when dealing with the probability density function (PDF) associated with Brownian motion. The point-wise PDF values tend to decrease exponentially as dimensionality increases, surpassing the precision of numerical simulations and resulting in substantial errors. In addition, due to its massive sampling, Monte Carlo fails to offer fast sampling. Modeling the logarithm likelihood (LL) via vanilla PINNs transforms the FP equation into a notoriously difficult Hamilton-Jacobi-Bellman (HJB) equation, which is impractical for PINN learning, whose error grows rapidly with dimension. To this end, we propose a novel approach utilizing a score-based solver to fit the score function in stochastic differential equations (SDEs). The score function, defined as the gradient of the LL, plays a fundamental role in inferring LL and PDF and enables fast SDE sampling, offering an effective means to overcome the CoD. Three fitting methods, Score Matching (SM), Sliced Score Matching (SSM), and Score-PINN, are introduced, each contributing unique advantages in computational complexity, accuracy, and generality. The proposed score-based SDE solver operates in two stages: first, employing score matching or Score-PINN to acquire the score function; and second, solving the LL via an ordinary differential equation (ODE) using the obtained score function. Comparative evaluations across these methods showcase varying trade-offs. The proposed methodology is evaluated across diverse SDEs, including anisotropic Ornstein-Uhlenbeck processes, geometric Brownian motion, and Brownian motion with varying eigenspace. We also test various distributions, including Gaussian, Log-normal, Laplace, and Cauchy distributions. The numerical results demonstrate the score-based SDE solver’s stability, speed, and performance across different experimental settings, solidifying its potential as a solution to CoD for high-dimensional FP equations.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

High Fidelity Simulations of a Methane-Oxygen Rotating Detonation Rocket Engine

The rotating detonation combustor (RDC) is an important realization of pressure gain combustion for practical applications. The RDC system is characterized by a highly unsteady flow field, with multiple reflected pressure waves following the strong detonations and the existence of partially-burnt gases entrained within the detonation annulus during the injector recovery process. While experimental efforts have provided macroscopic properties of an RDC with a given fuel-oxidizer mixture and geometry, limited optical access and reduced lifespan of flow-field diagnostic equipment constrain understanding of the mechanisms that lend to wave stability, controllability, and sustainability. To this end, high-fidelity numerical simulations of the Air Force Research Laboratory (AFRL) methane-oxygen rotating detonation rocket engine (RDRE) are performed to provide detailed insight into the detonation physics and anomalous behavior within the combustor. The results highlight the merit of high-fidelity numerical studies in analyzing the combustion processes and injector dynamics of a complex combustion system. The simulations capture the general wave behavior observed within the experimental studies. Three different operating conditions, with varying global equivalence ratio and system mass flow rate, are considered: (1) baseline configuration, (2) high equivalence ratio, and (3) high mass flow rate. The high mass flow rate case results in increased detonation wave strength whereas the high equivalence ratio results in increased deflagrative heat release. The high equivalence ratio case results in chaotic behavior with two pairs of counter-rotating detonation waves, whereas a single pair of co-rotating waves is observed in the baseline and high mass flow rate cases. Further, parasitic combustion of the reactant mixture ahead of the detonation wave diminishes its strength and reduces propagation speed to approximately 50-60\% of the Chapman-Jouguet (CJ) condition. In this work, a detailed analysis of the detonation wave behavior, injector dynamics, and the statistics of combustion parameters are presented.

Prakash, Supraj↗

NASA NTRS1973

Stability analysis of cylinders with circular cutouts

The stability of axially compressed cylinders with circular cutouts is analyzed numerically. An extension of the finite-difference method is used which removes the requirement that displacement components be defined in the directions of the grid lines. The results of this nonlinear analysis are found to be in good agreement with earlier experimental results.

Almroth, B. O.↗

NASA NTRS1973

Partial implicitization

The steady-state solution to the full Navier-Stokes equations for complicated flows is generally difficult to obtain. The Burgers (1948) equation is used as a model of the Navier-Stokes equations. The steady-state solution is obtained by a one-step explicit technique resulting from a partial implicitization of the difference equation. Stability analysis shows that the technique is unconditionally stable, and numerical tests show the technique to be accurate.

Graves, R. A., Jr.↗

NASA NTRS1976

Dynamics of spin stabilized spacecraft during deployment of telescoping appendages

The dynamics of a spin stabilized spacecraft with telescoping appendages are treated both analytically and numerically. Two types of telescoping appendages are considered: (a) where an end mass is mounted at the end of an (assumed) massless boom; and (b) where the appendage is assumed to consist of a uniformly distributed homogeneous mass throughout its length. For the telescoping system Eulerian equations of motion are developed. During all deployment sequences it is assumed that the transverse component of angular momentum is much smaller than the component along the major spin axis. Closed form analytical solutions for the time response of the transverse components of angular velocities are obtained when the spacecraft hub has a nearly spherical mass distribution. For the more general case, a series solution is obtained and this solution is limited by its radius of convergence. The comparison of the different approximate analytical methods with numerical integration results are studied and it is observed that the oscillatory nature of the responses of the transverse angular velocity components reduces rapidly with faster extension rates.

Sellappan, R.↗

NASA NTRS1977

A numerical analysis for a comparative study of the pointing accuracy of flexible satellites

The analytical models for flexible satellites have been reported earlier. Numerical results comparing the pointing accuracy of three-axes stabilized, spinning, and dual-spin satellites are presented here. The responses are compared for deterministic and random models with and without external torques. The performances of a uniform control system have been compared to a time-optimal control system for the same base-line configurations.

Huang, T. C.↗

NASA NTRS1981

Spherical models: Basic states and stability

Steady axisymmetric basic states in a spherical cap and their stability to azimuthally varying perturbations of different wave numbers were numerically investigated. The main features of the planned computer code and numerical methods are outlined.

Roberts, G. O.↗