Search NASA⌕ Search

SEARCH · Search NASA

Results for “numerical stabilization”

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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

A non-Boltzmannian behavior of the energy distribution for quasi-stationary regimes of the Fermi–Pasta–Ulam β system

In a recent paper [M. Leo, R.A. Leo, P. Tempesta, C. Tsallis, Phys. Rev. E 85 (2012) 031149], the existence of quasi-stationary states for the Fermi–Pasta–Ulam β system has been shown numerically, by analyzing the stability properties of the N/4-mode exact nonlinear solution. Here we study the energy distribution of the modes N/4, N/3 and N/2, when they are unstable, as a function of N and of the initial excitation energy. We observe that the classical Boltzmann weight is replaced by a different weight, expressed by a q-exponential function. -- Highlights: ► New statistical properties of the Fermi–Pasta–Ulam beta system are found. ► The energy distribution of specific observables are studied: a deviation from the standard Boltzmann behavior is found. ► A q-exponential weight should be used instead. ► The classical exponential weight is restored in the large particle limit (mesoscopic nature of the phenomenon)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Robust Approaches to Handling Complex Geometries with Galerkin Difference Methods

The Galerkin difference (GD) basis is a set of continuous, piecewise polynomials defined using a finite-difference-like grid of degrees of freedom. The one dimensional GD basis functions are naturally extended to multiple dimensions using a tensor product construction on quadrilateral elements. The GD basis can be used to define the solution space for a discontinuous Galerkin finite element discretization of partial differential equations. In this work we propose two approaches to handling complex geometries within this setting: (1) using nonconforming, curvilinear GD elements and (2) coupling affine GD elements with curvilinear simplicial elements. In both cases the (semidiscrete) discontinuous Galerkin method is provably energy stable even when variational crimes are committed. Additionally, for both element types a weight-adjusted mass matrix is used, which ensures that only the reference mass matrix must be inverted. We also present sufficient conditions on the treatment of metric terms for the curvilinear, nonconforming GD elements to ensure that the scheme is both constant preserving and conservative. Numerical experiments confirm the stability results and demonstrate the accuracy of the coupled schemes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Tunable Atomic Layer Deposition into Ultra-High-Aspect-Ratio (>60000:1) Aerogel Monoliths Enabled by Transport Modeling

Atomic layer deposition (ALD) modification of ultra-high-aspect-ratio structures (>10,000:1) is a powerful platform with applications in catalysis, filtration, and energy conversion. However, the deposition of conformal and tunable ALD coatings at these aspect ratios remains challenging, resulting in empirical tradeoffs between precursor utilization and reaction time. Here, we demonstrate tunable control of the ALD infiltration depth into an aerogel monolith (AM) and develop a reaction-diffusion model to accurately describe the coating process. Specifically, we investigate the ALD exposure time and precursor dose needed to conformally coat a silica AM with pore sizes of ~20 nm, a monolith thickness of ~2.5 mm, and aspect ratios exceeding 60,000:1. We demonstrate complete infiltration into the AM, which is quantified by elemental mapping. A reaction-diffusion model is developed, which accounts for multiple doses and precursor depletion in the ALD chamber during an exposure step. The experimentally validated model enables the prediction and tuning of infiltration depth into a tortuous, high-aspect-ratio structure such as an AM, allowing for the synthesis of rationally designed material architectures. Additionally, the model allows for co-optimization of the total deposition time and percentage of unreacted precursor, which are important for the manufacturability and sustainability of ALD processing. Lastly, we demonstrate that ultra-thin ALD Al 2 O 3 coatings can be used to stabilize silica AMs against structural degradation under high-temperature annealing conditions (700-800 °C), by limiting changes in surface area and monolith volume. As a result, this improved high-temperature stability has implications for numerous aerogel applications, including catalysis and thermal insulation.

14 SOLAR ENERGY↗
DOE OSTI2018

Metal-Free Aqueous Flow Battery with Novel Ultrafiltered Lignin as Electrolyte

As the number of generation sources from intermittent renewable technologies on the electric grid increases, the need for large-scale energy storage devices is becoming essential to ensure grid stability. Flow batteries offer numerous advantages over conventional sealed batteries for grid storage. In this work, for the first time, we investigated lignin, the second most abundant wood derived biopolymer, as an anolyte for the aqueous flow battery. Lignosulfonate, a water-soluble derivative of lignin, is environmentally benign, low cost and abundant as it is obtained from the byproduct of paper and biofuel manufacturing. The lignosulfonate utilizes the redox chemistry of quinone to store energy and undergoes a reversible redox reaction. Here, we paired lignosulfonate with Br2/Br-, and the full cell runs efficiently with high power density. Also, the large and complex molecular structure of lignin considerably reduces the electrolytic crossover, which ensures very high capacity retention. The flowcell was able to achieve current densities of up to 20 mA/cm2 and charge polarization resistance of 15 ohm cm2. This technology presents a unique opportunity for a low-cost, metal-free flow battery capable of large-scale sustainable energy storage.

25 ENERGY STORAGE↗
DOE OSTI2013

Open-loop control of noise amplification in a separated boundary layer flow

Linear optimal gains are computed for the subcritical two-dimensional separated boundary-layer flow past a bump. Very large optimal gain values are found, making it possible for small-amplitude noise to be strongly amplified and to destabilize the flow. The optimal forcing is located close to the summit of the bump, while the optimal response is the largest in the shear layer. The largest amplification occurs at frequencies corresponding to eigenvalues which first become unstable at higher Reynolds number. Nonlinear direct numerical simulations show that a low level of noise is indeed sufficient to trigger random flow unsteadiness, characterized here by large-scale vortex shedding. Next, a variational technique is used to compute efficiently the sensitivity of optimal gains to steady control (through source of momentum in the flow, or blowing/suction at the wall). A systematic analysis at several frequencies identifies the bump summit as the most sensitive region for control with wall actuation. Based on these results, a simple open-loop control strategy is designed, with steady wall suction at the bump summit. Linear calculations on controlled base flows confirm that optimal gains can be drastically reduced at all frequencies. Nonlinear direct numerical simulations also show that this control allows the flow to withstand a higher level of stochastic noise without becoming nonlinearly unstable, thereby postponing bypass transition. In the supercritical regime, sensitivity analysis of eigenvalues supports the choice of this control design. Full restabilization of the flow is obtained, as evidenced by direct numerical simulations and linear stability analysis.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Energetic ions in ITER plasmas

This paper discusses the behaviour and consequences of the expected populations of energetic ions in ITER plasmas. It begins with a careful analytic and numerical consideration of the stability of Alfvén Eigenmodes in the ITER 15 MA baseline scenario. The stability threshold is determined by balancing the energetic ion drive against the dominant damping mechanisms and it is found that only in the outer half of the plasma (r/a>0.5) can the fast ions overcome the thermal ion Landau damping. This is in spite of the reduced numbers of alpha-particles and beam ions in this region but means that any Alfvén Eigenmode-induced redistribution is not expected to influence the fusion burn process. The influence of energetic ions upon the main global MHD phenomena expected in ITER's primary operating scenarios, including sawteeth, neoclassical tearing modes and Resistive Wall Modes, is also reviewed. Fast ion losses due to the non-axisymmetric fields arising from the finite number of toroidal field coils, the inclusion of ferromagnetic inserts, the presence of test blanket modules containing ferromagnetic material, and the fields created by the Edge Localised Mode (ELM) control coils in ITER are discussed. The greatest losses and associated heat loads onto the plasma facing components arise due to the use of the ELM control coils and come from neutral beam ions that are ionised in the plasma edge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Quenching behavior of molten pool with different strategies – A review

After the major severe accident in nuclear reactor, there has been lot of concerns regarding long term core melt stabilization following a severe accident in nuclear reactors. Numerous strategies have been though for quenching and stabilization of core melt like top flooding, bottom flooding, indirect cooling, etc. However, the effectiveness of these schemes is yet to be determined properly, for which, lot of experiments are needed. Several experiments have been performed for coolability of melt pool under bottom flooding as well as for indirect cooling. Besides these tests are very scattered because they involve different simulants material initial temperatures and masses of melt, which makes it very complex to judge the effectiveness of a particular technique and advantage over the other. In this review paper, a study has been carried on different cooling techniques of simulant materials with same mass. Three techniques have been compared here and the results are discussed. Under top flooding technique it took several hours to cool the melt under without decay heat condition. In bottom flooding technique was found to be the best technique among in indirect cooling technique, top flooded technique, and bottom flooded technique.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2017

Plasma-induced flow instabilities in atmospheric pressure plasma jets

Pulsed plasma excitation of rare gases flowing into air has been shown to impact the stability of the flow in non-equilibrium atmospheric pressure plasma jets (APPJs). In this paper, the results from a numerical modeling investigation of the stability of a round He APPJ with a powered electrode exposed to the gas flow are discussed. Localized gas heating at the powered electrode occurs on the time scale of the voltage pulse, tens to 100 ns, which is short compared to the fluid timescales. An acoustic wave propagates from this heated, expanding gas and exits the jet. The wave disturbs the shear layer between the He and surrounding humid air, exciting a shear instability which grows downstream with the flow and increases the mixing of the humid air into the He. The effects of the eddy-dominated flow on ionization wave (IW) propagation in an APPJ were investigated. The IW followed the regions of the highest helium concentration, resulting in an increased production of NO, HO 2 , and NO 2 .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Modeling MHD accretion-ejection—from the launching area to propagation scales

We present the results of axisymmetric magnetohydrodynamic simulations investigating the launching of jets and outflows from a magnetically diffusive accretion disk. The time evolution of the disk structure is self-consistently taken into account. In contrast to previous works, we have applied spherical coordinates for the numerical grid, implying substantial benefits concerning the numerical resolution and the stability of the simulation. Due to the new setup, we were able to run simulations for more than 150,000 dynamical times on a domain extending 1500 inner disk radii with a resolution of up to 24 cells per disk height in the inner disk. Depending on the disk magnetization, jet launching occurs in two different but complementary regimes—jets driven predominantly by centrifugal or magnetic forces. These regimes differ in the ejection efficiency concerning mass, energy, and angular momentum. We show that it is the actual disk magnetization and not so much the initial magnetization that best describes the disk-jet evolution. Considering the actual disk magnetization, we also find that simulations starting with different initial magnetization evolve in a similar—typical—way, due to advection and diffusion, as the magnetic flux in the disk evolves in time. Exploring a new, modified diffusivity model, we confirm the self-similar structure of the global jet-launching disk, obtaining power laws for the radial profiles of the disk physical variables such as density, magnetic field strength, or accretion velocity.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Properties of Nb x Ti (1–x) N thin films deposited on 300 mm silicon wafers for upscaling superconducting digital circuits

Scaling superconducting digital circuits requires fundamental changes in the current material set and fabrication process. The transition to 300 mm wafers and the implementation of advanced lithography are instrumental in facilitating mature CMOS processes, ensuring uniformity, and optimizing the yield. Here, this study explores the properties of Nb x Ti (1–x) N films fabricated by magnetron DC sputtering on 300 mm Si wafers. As a promising alternative to traditional Nb in device manufacturing, Nb x Ti (1–x) N offers numerous advantages, including enhanced stability and scalability to smaller dimensions, in both processing and design. As a ternary material, Nb x Ti (1–x) N allows engineering material parameters by changing deposition conditions. The engineered properties can be used to modulate device parameters through the stack and mitigate failure modes. We report characterization of Nb x Ti (1–x) N films at less than 2% thickness variability, 2.4% T c variability and 3% composition variability. Film resistivity (140–375 Ωcm) shows a strong correlation with the film oxygen content, while the critical temperature T c (4.6 K–14.1 K) is strongly affected by film stoichiometry and its microstructure has only a moderate effect on modifying T c . Our results offer insights about the interplay between film stoichiometry, film microstructure and critical temperature.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Discrete solitons and vortices in anisotropic hexagonal and honeycomb lattices

We consider the self-focusing discrete nonlinear Schrödinger equation on hexagonal and honeycomb lattice geometries. Our emphasis is on the study of the effects of anisotropy, motivated by the tunability afforded in recent optical and atomic physics experiments. We find that multi-soliton and discrete vortex states undergo destabilizing bifurcations as the relevant anisotropy control parameter is varied. Furthermore, we quantify these bifurcations by means of explicit analytical calculations of the solutions, as well as of their spectral linearization eigenvalues. Finally, we corroborate the relevant stability picture through direct numerical computations. In the latter, we observe the prototypical manifestation of these instabilities to be the spontaneous rearrangement of the solution, for larger values of the coupling, into localized waveforms typically centered over fewer sites than the original unstable structure. In weak coupling, the instability appears to result in a robust breathing of the relevant waveforms.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Reaction-Infiltration Instabilities in Fractured and Porous Rocks

In this project we are developing a multiscale analysis of the evolution of fracture permeability, using numerical simulations and linear stability analysis. Our simulations include fully three-dimensional simulations of the fracture topography, fluid flow, and reactant transport, two-dimensional simulations based on aperture models, and linear stability analysis.

58 GEOSCIENCES↗
DOE OSTI2020

The virtual element method for linear elastodynamics models: Design, analysis, and implementation

We design the conforming virtual element method for the numerical simulation of two dimensional time-dependent elastodynamics problems. We investigate the performance of the method both theoretically and numerically. We prove the stability and the convergence of the semi-discrete approximation in the energy norm and derive optimal error estimates. We also show the convergence in the L 2 norm. The performance of the virtual element method is assessed on a set of different computational meshes, including non-convex cells up to order four in the h-refinement setting. Exponential convergence is also experimentally seen in the p-refinement setting.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1962

GENERAL INSTABILITY OF ORTHOGONALLY STIFFENED CYLINDRICAL SHELLS

Earlier research at the National Aeronautics Research Institute (N.L.R.), Amsterdam, which forms the basis of recent work is reviewed. This early work refers to 2 schemes: the orthotropic shell and, in view of buckling modes where the half wave length is of the order of the ring distance, the shell with continuously distributed stringers and discrete rings. Linear theory is considered to be adequate for these structures, where the imperfections are small in comparison to the height of the ring sections. Recent developments account for pressure difference in addition to axial compression, for the correct stiffness matrix of skin panels in the post-buckling stage and for stringer bending due to hoop stresses in the skin, which are of importance as has been shown by the investigation of the post-buckling behaviour. Numerical data for the stiffness matrix of skin panels have been established. Numerical evaluation of the stability equation has not been performed as yet.

SHELL STABILITY↗
NASA NTRS1973

Numerical studies of unsteady two dimensional subsonic flows using the ICE method

A numerical program was developed to compute transient compressible and incompressible laminar flows in two dimensions with multicomponent mixing and chemical reaction. The algorithm used the Los Alamos Scientific Laboratory ICE (Implicit Continuous-Fluid Eulerian) method as its base. The program can compute both high and low speed compressible flows. The numerical program incorporating the stabilization techniques was quite successful in treating both old and new problems. Detailed calculations of coaxial flow very close to the entry plane were possible. The program treated complex flows such as the formation and downstream growth of a recirculation cell. An implicit solution of the species equation predicted mixing and reaction rates which compared favorably with the literature.

Wieber, P. R.↗
NASA NTRS1975

Upwind second-order difference schemes and applications in unsteady aerodynamic flows

Explicit second-order upwind difference schemes in combination with spatially symmetric schemes can produce larger stability bounds and better numerical resolution than symmetric schemes alone. However, if conservation form is essential, a special operator is required for transition between schemes. An operational approach has been devised for deriving transition operators so that strict conservation and local consistency are maintained. Various aspects of hybrid schemes are studied numerically for model linear and nonlinear equations. To demonstrate the utility of combining two different algorithms, MacCormack's explicit, noncentered, second-order method is combined with a completely upwind version, and numerical solutions of the Euler equations are obtained for two-dimensional, transonic flows with embedded supersonic regions and shock waves. The general utility of the operational approach for combining schemes is emphasized by deriving a second-order conservative scheme for the steady transonic small-disturbance potential equation.

Warming, R. F.↗
NASA NTRS1980

Heat- and chemical-resistant oxdiazole elastomers

Heat and chemical resistant polymers with triazine crosslinks are prepared by thermal condensation reactions to form 1,2,4-oxdiazole linkages. They are compounded with variety of fillers, extenders, and modifiers for numerous applications in which stability, impermiability to liquids and gases, good plasticity, and elasticity or rigidity are important.

Rosser, R. W.↗