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 1,567 records · Page 87

All sources NASA NTRS DOE OSTI

NASA NTRS1957

Handbook of structural stability part II : buckling of composite elements

The local buckling of stiffener sections and the buckling of plates with sturdy stiffeners are reviewed, and the results are summarized in charts and tables. Numerical values of buckling coefficients are presented for longitudinally compressed stiffener sections of various shapes, for stiffened plates loaded in longitudinal compression and in shear, and for stiffened cylinders loaded in torsion. Although the data presented consist primarily of elastic-buckling coefficients, the effects of plasticity are discussed for a few special cases. (author)

Becker, Herbert↗

NASA NTRS2019

Effects of Spin-Orbit Resonance in Stability for Low Altitude Mars Orbits

Orbit stability has been thoughtfully studied in various celestial bodies. The increasing interest in Mars orbiters brings the question of the likelihood of natural decay in low altitude regimes. This paper studies the shape change of low altitude Mars orbits by carrying out large sets of numerical high fidelity simulations. Results showed that various configurations of the orbital elements gave perturbations that resulted in unstable orbits. The paper also studies the potential causes of the observed unstable regions. We computed theoretical spin-orbit resonances to study their implications in the stability at low altitudes. The resonances were tested at different initial Longitudes of the Ascending Node (LAN) and orbit inclinations to check the potential existence of latitude/longitude implications on the stability.

Spin-Orbit Resonance↗

NASA NTRS1987

Flap-lag-torsion aeroelastic stability of a circulation control rotor in forward flight

The aeroelastic stability of a circulation control rotor blade undergoing three degrees of motion (flap, lag, and torsion) is investigated in forward flight. Quasi-steady strip theory is used to evaluate the aerodynamics forces; and the airfoil characteristics are from data tables. The propulsive and the auxiliary power trims are calculated from vehicle and rotor equilibrium equations through the numerical integration of element forces in azimuth as well as in radial directions. The nonlinear time dependent periodic blade response is calculated using an iterative procedure based on Floquet theory. The periodic perturbation equations are solved for stability using Floquet transition matrix theory. The effects of several parameters on blade stability are examined, including advance ratio, collective pitch, thrust level, shaft tilt, structural stiffnesses variation, and propulsive and auxiliary power trims.

Chopra, Inderjit↗

DOE OSTI2020

Direct observation of the formation and stabilization of metallic nanoparticles on carbon supports

Direct formation of ultra-small nanoparticles on carbon supports by rapid high temperature synthesis method offers new opportunities for scalable nanomanufacturing and the synthesis of stable multi-elemental nanoparticles. However, the underlying mechanisms affecting the dispersion and stability of nanoparticles on the supports during high temperature processing remain enigmatic. In this work, we report the observation of metallic nanoparticles formation and stabilization on carbon supports through in situ Joule heating method. We find that the formation of metallic nanoparticles is associated with the simultaneous phase transition of amorphous carbon to a highly defective turbostratic graphite (T-graphite). Molecular dynamic (MD) simulations suggest that the defective T-graphite provide numerous nucleation sites for the nanoparticles to form. Furthermore, the nanoparticles partially intercalate and take root on edge planes, leading to high binding energy on support. This interaction between nanoparticles and T-graphite substrate strengthens the anchoring and provides excellent thermal stability to the nanoparticles. These findings provide mechanistic understanding of rapid high temperature synthesis of metal nanoparticles on carbon supports and the origin of their stability.

36 MATERIALS SCIENCE↗

NASA NTRS1981

HADY-I, a FORTRAN program for the compressible stability analysis of three-dimensional boundary layers

A computer program HADY-I for calculating the linear incompressible or compressible stability characteristics of the laminar boundary layer on swept and tapered wings is described. The eigenvalue problem and its adjoint arising from the linearized disturbance equations with the appropriate boundary conditions are solved numerically using a combination of Newton-Raphson interative scheme and a variable step size integrator based on the Runge-Kutta-Fehlburh fifth-order formulas. The integrator is used in conjunction with a modified Gram-Schmidt orthonormalization procedure. The computer program HADY-I calculates the growth rates of crossflow or streamwise Tollmien-Schlichting instabilities. It also calculates the group velocities of these disturbances. It is restricted to parallel stability calculations, where the boundary layer (meanflow) is assumed to be parallel. The meanflow solution is an input to the program.

El-Hady, N. M.↗

NASA NTRS1993

On the nonlinear interfacial instability of rotating core-annular flow

The interfacial stability of rotating core-annular flows is investigated. The linear and nonlinear effects are considered for the case when the annular region is very thin. Both asymptotic and numerical methods are used to solve the flow in the core and film regions which are coupled by a difference in viscosity and density. The long-term behavior of the fluid-fluid interface is determined by deriving its nonlinear evolution in the form of a modified Kuramoto-Sivashinsky equation. We obtain a generalization of this equation to three dimensions. The flows considered are applicable to a wide array of physical problems where liquid films are used to lubricate higher or lower viscosity core fluids, for which a concentric arrangement is desired. Linearized solutions show that the effects of density and viscosity stratification are crucial to the stability of the interface. Rotation generally destabilizes non-axisymmetric disturbances to the interface, whereas the centripetal forces tend to stabilize flows in which the film contains the heavier fluid. Nonlinear affects allow finite amplitude helically travelling waves to exist when the fluids have different viscosities.

Coward, Aidrian V.↗

DOE OSTI2022

Orbital stability analysis and photometric characterization of the second Earth Trojan asteroid 2020 XL5

Trojan asteroids are small bodies orbiting around the L 4 or L 5 Lagrangian points of a Sun-planet system. Due to their peculiar orbits, they provide key constraints to the Solar System evolution models. Despite numerous dedicated observational efforts in the last decade, asteroid 2010 TK 7 has been the only known Earth Trojan thus far. Here we confirm that the recently discovered 2020 XL 5 is the second transient Earth Trojan known. To study its orbit, we used archival data from 2012 to 2019 and observed the object in 2021 from three ground-based observatories. Our study of its orbital stability shows that 2020 XL 5 will remain in L 4 for at least 4 000 years. With a photometric analysis we estimate its absolute magnitude to be ${H}_{r}=18.5{8}_{-0.15}^{+0.16}$, and color indices suggestive of a C-complex taxonomy. Assuming an albedo of 0.06 ± 0.03, we obtain a diameter of 1.18 ± 0.08 km, larger than the first known Earth Trojan asteroid.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

A novel defoamer for processing nuclear waste: Testing and performance

Abstract Legacy radioactive waste from nuclear weapons material processing is the biggest environmental challenge in the state of South Carolina, and one of the biggest in the United States. Although substantial progress has been made in processing and vitrifying high‐level radioactive waste at the Savannah River site, approximately 35 million gallons remains to be treated and dispositioned. In this article, we show the development of a new defoamer for use in the processing of high‐level radioactive waste. The application of the new defoamer is not only more effective at controlling foam but will shorten the batch processing time leading to safer and more efficient processing. This has the potential to help shorten the overall site mission lifetime, saving the federal government hundreds of millions of dollars. The novel aspect of the new defoamer is its effectiveness and stability in the harsh conditions needed for processing high‐level radioactive waste: high temperatures, pH 4–13, strong oxidizing and reducing agents, and numerous metal catalysts. The defoamer will replace an existing antifoam agent that decomposes to form multiple flammable gases. Because the existing antifoam agent decomposes quickly, processing changes increased processing time. Statement of novelty We demonstrate that a new defoamer has been developed for use in processing high‐level radioactive waste. This novel defoamer is effective and stable in the harsh processing conditions with no detectable flammable by‐products. Because the existing antifoam agent decomposes quickly, processing changes, including lowering acid addition rates and evaporation rates, increased processing time. The new defoamer is not only more effective at controlling foam but will shorten the batch processing time and lead to safer and more efficient processing by eliminating the production of flammable gases. Statement of industrial relevance Foaming is of great concern in many industrial processes involving three‐phase gas/liquid/fine‐solids systems, like in water evaporation, mineral floatation, air sparging in situ remediation techniques, and in those found in the paper industry.

Lambert, Dan P.↗

DOE OSTI2023

Gaussian process hydrodynamics

Abstract We present a Gaussian process (GP) approach, called Gaussian process hydrodynamics (GPH) for approximating the solution to the Euler and Navier-Stokes (NS) equations. Similar to smoothed particle hydrodynamics (SPH), GPH is a Lagrangian particle-based approach that involves the tracking of a finite number of particles transported by a flow. However, these particles do not represent mollified particles of matter but carry discrete/partial information about the continuous flow. Closure is achieved by placing a divergence-free GP prior ξ on the velocity field and conditioning it on the vorticity at the particle locations. Known physics (e.g., the Richardson cascade and velocity increment power laws) is incorporated into the GP prior by using physics-informed additive kernels. This is equivalent to expressing ξ as a sum of independent GPs ξ l , which we call modes, acting at different scales (each mode ξ l self-activates to represent the formation of eddies at the corresponding scales). This approach enables a quantitative analysis of the Richardson cascade through the analysis of the activation of these modes, and enables us to analyze coarse-grain turbulence statistically rather than deterministically. Because GPH is formulated by using the vorticity equations, it does not require solving a pressure equation. By enforcing incompressibility and fluid-structure boundary conditions through the selection of a kernel, GPH requires significantly fewer particles than SPH. Because GPH has a natural probabilistic interpretation, the numerical results come with uncertainty estimates, enabling their incorporation into an uncertainty quantification (UQ) pipeline and adding/removing particles (quanta of information) in an adapted manner. The proposed approach is suitable for analysis because it inherits the complexity of state-of-the-art solvers for dense kernel matrices and results in a natural definition of turbulence as information loss. Numerical experiments support the importance of selecting physics-informed kernels and illustrate the major impact of such kernels on the accuracy and stability. Because the proposed approach uses a Bayesian interpretation, it naturally enables data assimilation and predictions and estimations by mixing simulation data and experimental data.

Mathematics↗

DOE OSTI2023

Machine learning assisted derivation of minimal low-energy models for metallic magnets

Abstract We consider the problem of extracting a low-energy spin Hamiltonian from a triangular Kondo Lattice Model (KLM). The non-analytic dependence of the effective spin-spin interactions on the Kondo exchange excludes the use of perturbation theory beyond the second order. We then introduce a Machine Learning (ML) assisted protocol to extract effective two- and four-spin interactions. The resulting spin model reproduces the phase diagram of the original KLM as a function of magnetic field and single-ion anisotropy and reveals the effective four-spin interactions that stabilize the field-induced skyrmion crystal phase. Moreover, this model enables the computation of static and dynamical properties with a much lower numerical cost relative to the original KLM. A comparison of the dynamical spin structure factor in the fully polarized phase computed with both models reveals a good agreement for the magnon dispersion even though this information was not included in the training data set.

Sharma, Vikram (ORCID:0000000189105519)↗

DOE OSTI2022

Atomically dispersed iron sites with a nitrogen–carbon coating as highly active and durable oxygen reduction catalysts for fuel cells

Nitrogen-coordinated single atom iron sites (FeN 4 ) embedded in carbon (Fe–N–C) are the most active platinum group metal-free oxygen reduction catalysts for proton-exchange membrane fuel cells. Still, current Fe–N–C catalysts lack sufficient long-term durability and are not yet viable for practical applications. Here we report a highly durable and active Fe–N–C catalyst synthesized using heat treatment with ammonia chloride followed by high-temperature deposition of a thin layer of nitrogen-doped carbon on the catalyst surface. We propose that catalyst stability is improved by converting defect-rich pyrrolic N-coordinated FeN 4 sites into highly stable pyridinic N-coordinated FeN 4 sites. The stability enhancement is demonstrated in membrane electrode assemblies using accelerated stress testing and a long-term steady-state test (>300 h at 0.67 V), approaching a typical Pt/C cathode (0.1 mg Pt cm -2 ). The encouraging stability improvement represents a critical step in developing viable Fe–N–C catalysts to overcome the cost barriers of hydrogen fuel cells for numerous applications.

08 HYDROGEN↗

DOE OSTI2024

Fast ion relaxation in ITER mediated by Alfvén instabilities

Abstract We address the critical issue for future burning plasmas of whether high-energy fusion products or auxiliary heating-beam ions will be confined for a sufficiently long time to compensate for thermal plasma energy losses. This issue can be mitigated by one of the most deleterious collective phenomena—the instability of low, sub-cyclotron frequency Alfvén eigenmodes (AEs), such as toroidicity-induced AEs and reversed-shear AEs in the ITER steady-state scenario. Using a revised quasi-linear (QL) theory applied to energetic particle (EP) relaxation in the presence of AEs, we find that the AE instabilities can affect both neutral beam ions and alpha particles, although the resulting fast ion transport is expected to be modest if classical particle slowing down is assumed. On the other hand, the QL theory predicts that the AE amplitudes will be enhanced by the background microturbulence, although this topic remains outside our scope due to the significant numerical effort required to evaluate these effects. We report our results for EP relaxation dynamics obtained utilizing several tools: (i) a comprehensive linear stability study of the sub-cyclotron Alfvénic spectrum as computed by ideal magnetohydrodynamic NOVA simulations for the AE eigenproblem, (ii) drift kinetic NOVA-C calculations for wave–particle interaction and AE growth/damping rates, and (iii) predictive QL modeling coupled with the global transport code TRANSP to assess the EP relaxation on the equilibrium timescale.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1976

Shielding, the bulk chemical potential, and cohesion in alloys

It is shown that the bulk chemical potential in alloys is intimately related to the spatial dependence of the shielding cloud that results when the electronic charge rearranges itself as one atom type is replaced by another at a given site. Such a relationship fixes the relative energy scale between the alloy and its pure constituents, important in determining the stability of alloys. A correct treatment of shielding is thus essential to quantitative calculations of alloy stability. A model calculation of the bulk chemical potential and cohesion of alloys in the tight-binding approximation is presented as a numerical example. In the course of this investigation a general invariant of an integral over the shielding cloud is derived.

Stern, E. A.↗

NASA NTRS1980

Note on the eigensolution of a homogeneous equation with semi-infinite domain

The 'variation-iteration' method using Green's functions to find the eigenvalues and the corresponding eigenfunctions of a homogeneous Fredholm integral equation is employed for the stability analysis of fluid hydromechanics problems with a semiinfinite (infinite) domain of application. The objective of the study is to develop a suitable numerical approach to the solution of such equations in order to better understand the full set of equations for 'real-world' flow models. The study involves a search for a suitable value of the length of the domain which is a fair finite approximation to infinity, which makes the eigensolution an approximation dependent on the length of the interval chosen. In the examples investigated y = 1 = a seems to be the best approximation of infinity; for y greater than unity this method fails due to the polynomial nature of Green's functions.

Wadia, A. R.↗

NASA NTRS1980

A perturbation solution for transverse pressure-sensitive combustion instability in an annulus

A method for analysis of combustion instability in rocket motors based on a combination of the Galerkin method and the two-variable (or multiple-scales) perturbation method is developed. The method is illustrated by applying it to the problem of pressure-sensitive combustion instability in a liquid-fuel annular combustion chamber. To the order of approximation inherent in the method it is found that the complete linear stability analysis and most of the nonlinear stability analysis can be carried out in closed form. The results thus obtained are used to clarify several aspects of the interpretation of numerical solutions obtained by previous investigators.

Googerdy, A.↗

NASA NTRS1987

Vibrations in a moving flexible robot arm

The vibration in a flexible robot arm modeled by a moving slender prismatic beam is considered. It is found that the extending and contracting motions have destabilizing and stabilizing effects on the vibratory motions, respectively. The vibration analysis is based on a Galerkin approximation with time-dependent basis functions. Typical numerical results are presented to illustrate the qualitative features of vibrations.

Wang, P. K. C.↗

NASA NTRS1988

The dynamics of coorbital satellite systems

The dynamical behavior of N coorbital satellites moving with the same average mean motion around a primary has been studied both analytically and by numerical integrations for N from 2 to 9 satellites. Simplified dynamical equations are used to determine the different stationary configurations and their local stability against infinitesimal perturbations. The ring of equally spaced identical satellites is found to be locally unstable for N of six or less, while for N from two to eight there exists another stable compact stationary configuration with separations of 60 deg or less between adjacent satellites. For N of seven or more the equally spaced configuration becomes locally stable, and for N of nine or more it is the only stationary configuration. The motion becomes chaotic for large-amplitude perturbations. The chaotic motion fills a restricted region of phase space whose outer boundary is determined by the maximum velocity curve.

Salo, H.↗

NASA NTRS1990

Prominence formation in a coronal loop

A model is presented which depends on the preferential deposition of heating in the legs of a coronal loop and which produces a stable prominence-scale condensation at the loop top. Dynamic stability is attained by the subsequent adjustment of local parallel gravity by a magnetic inversion at the loop (or arcade) apex. A nonlinear numerical simulation of this process, which includes a deep chromosphere, a heating rate with a fixed dissipation length, and full solar gravity is described.

Mok, Y.↗