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 1,027 records · Page 57
All sourcesNASA NTRSDOE OSTI
NASA NTRS1980

The dynamics and control of large flexible space structures. Volume 3, part B: The modelling, dynamics, and stability of large Earth pointing orbiting structures

The dynamics and stability of large orbiting flexible beams, and platforms and dish type structures oriented along the local horizontal are treated both analytically and numerically. It is assumed that such structures could be gravitationally stabilized by attaching a rigid light-weight dumbbell at the center of mass by a spring loaded hinge which also could provide viscous damping. For the beam, the small amplitude inplane pitch motion, dumbbell librational motion, and the anti-symmetric elastic modes are all coupled. The three dimensional equations of motion for a circular flat plate and shallow spherical shell in orbit with a two-degree-of freedom gimballed dumbbell are also developed and show that only those elastic modes described by a single nodal diameter line are influenced by the dumbbell motion. Stability criteria are developed for all the examples and a sensitivity study of the system response characteristics to the key system parameters is carried out.

Bainum, P. M.↗
NASA NTRS1991

Numerical simulations of the nonlinear kink modes in linearly stable supersonic slip surfaces

High-resolution numerical simulations of supersonic slip surfaces were performed using the Cray 2 at the Minnesota Supercomputer Center to confirm earlier analytic nonlinear stability calculations of such structures. The present study presents a numerical solution for the problem solved analytically by Artola and Majda (1987) and follows growing modes well into the large-amplitude regime. The response of an equal-density Mach-4 slip surface to a variety of impinging nonlinear sound wave trains is computed using the piecewise-parabolic method (PPM). The nature of the nonlinear kink modes observed in the simulations was found to be similar to that discussed by Artola and Majda. Most of the speeds moved in either direction, with speeds near the predicted value. It is suggested that the stationary modes will eventually dominate the flow at much later times, and it is argued that the stationary modes are more disruptive than the propagating kink modes.

Pedelty, Jeffrey A.↗
NASA NTRS1989

Instability of jets of arbitrary geometry

This paper describes a calculation technique for determining the stability of jets of arbitrary cross section. In particular, elliptic and rectangular jets are considered. The numerical procedure involves both a conformal transformation between the computational domain and the physical plane and a solution of the transformed stability equation in the computational domain. Modern, efficient, conformal mappings are used for both simply and doubly connected domains. The numerical solution is based on a finite difference/pseudospectral discretization of the stability equation. The technique is verified by comparison with previous calculations for circular and elliptic jets. Calculations are performed for the stability of elliptic and rectangular jets of aspect ratio 2. Growth rates, phase velocities, and pressure eigenfunctions are calculated.

Baty, Roy S.↗
NASA NTRS1995

Wave Journal Bearing. Part 2: Experimental Pressure Measurements and Fractional Frequency Whirl Threshold for Wave and Plain Journal Bearings

A new hydrodynamic bearing concept, the wave journal bearing, is being developed because it has better stability characteristics than plain journal bearings while maintaining similar load capacity. An analysis code to predict the steady state and dynamic performance of the wave journal bearing is also part of the development. To verify numerical predictions and contrast the wave journal bearing's stability characteristics to a plain journal bearing, tests were conducted at NASA Lewis Research Center using an air bearing test rig. Bearing film pressures were measured at 16 ports located around the bearing circumference at the middle of the bearing length. The pressure measurements for both a plain journal bearing and a wave journal bearing compared favorably with numerical predictions. Both bearings were tested with no radial load to determine the speed threshold for self-excited fractional frequency whirl. The plain journal bearing started to whirl immediately upon shaft start-up. The wave journal did not incur self-excited whirl until 800 to 900 rpm as predicted by the analysis. Furthermore, the wave bearing's geometry limited the whirl orbit to less than the bearing's clearance. In contrast, the plain journal bearing did not limit the whirl orbit, causing it to rub.

Walker, James F.↗
DOE OSTI2021

Automatically adaptive stabilized finite elements and continuation analysis for compaction banding in geomaterials

Under compressive creep, viscoplastic solids experiencing internal mass transfer processes can accommodate singular cnoidal wave solutions as material instabilities at the stationary wave limit. These instabilities appear when the loading rate is significantly faster than the material's capacity to diffusive internal perturbations, leading to localized failure features (e.g., cracks and compaction bands). These cnoidal waves, generally found in fluids, have strong nonlinearities that produce periodic patterns. Due to the singular nature of the solutions, the applicability of the theory is currently limited. Additionally, practical simulation tools require proper regularization to overcome the challenges that singularity induces. We focus on the numerical treatment of the governing equation using a nonlinear approach building on a recent adaptive stabilized finite element method. This automatic refinement method provides an error estimate that drives mesh adaptivity, a crucial feature for the problem at hand. We compare the performance of this adaptive strategy against analytical and standard finite element solutions. We then investigate the sensitivity of the diffusivity ratio, the parameter controlling the process, and identify multiple possible solutions with several stress peaks. We also show the evolution of the spacing between peaks for all solutions as a function of that parameter.

42 ENGINEERING↗
DOE OSTI2016

A second-order ensemble method based on a blended backward differentiation formula timestepping scheme for time-dependent Navier-Stokes equations

Here, we present a second-order ensemble method based on a blended three-step backward differentiation formula (BDF) timestepping scheme to compute an ensemble of Navier–Stokes equations. Compared with the only existing second-order ensemble method that combines the two-step BDF timestepping scheme and a special explicit second-order Adams–Bashforth treatment of the advection term, this method is more accurate with nominal increase in computational cost. We give comprehensive stability and error analysis for the method. Numerical examples are also provided to verify theoretical results and demonstrate the improved accuracy of the method.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

First‐Order Empirical Interpolation Method for Real‐Time Solution of Parametric Time‐Dependent Nonlinear PDEs

ABSTRACT We present a model reduction approach for the real‐time solution of time‐dependent nonlinear partial differential equations (PDEs) with parametric dependencies. A major challenge in constructing efficient and accurate reduced‐order models for nonlinear PDEs is the efficient treatment of nonlinear terms. We address this by unifying the implementation of hyperreduction methods to deal with nonlinear terms. Furthermore, we introduce a first‐order empirical interpolation method (EIM) to provide an efficient approximation of the nonlinear terms in time‐dependent PDEs. We demonstrate the effectiveness of our approach on the Allen–Cahn equation, which models phase separation, and the Buckley–Leverett equation, which describes two‐phase fluid flow in porous media. Numerical results highlight the accuracy, efficiency, and stability of the proposed method compared with both the Galerkin–Newton approach and hyper‐reduced models using the standard EIM.

Nguyen, Ngoc Cuong [Center for Computational Engin↗
DOE OSTI2019

The Additive Approximation for Heat Transfer and for Stability Calculations in a Multi-filamentary Superconductor - part A

This paper addresses the applicability of the “additive approximation” of total thermal conductivity in heat transfer and in superconductor stability calculations. If cases (a), (b) and (c) denote total (conductive plus radiative), or only conductive or only radiative heat flux, respectively, each flux q-dot calculated with Fourier’s conduction law using the corresponding thermal conductivity (λ{sub a}, λ{sub b}, λ{sub c}), the additive approximation would be confirmed if the heat flux difference Δ q-dot = q-dot {sub a}− q-dot {sub b}− q-dot {sub c}, at any position of an investigated object, and at any time, converges to zero. This is not trivial because of the strong, non-linear temperature dependence of the radiation component. Heat transfer calculations including radiative transfer are presented in this paper, first for simple, homogeneous, thin film test samples and later for a multi-filamentary BSCCO 2223 superconductor. The simulated heat sources either result from a sudden increase of conductor boundary temperature or from flux flow and Ohmic resistances in the superconductor under a disturbance (like transport current exceeding critcal current density). The conductors, though very thin, are non-transparent to mid-IR radiation. Validity of the additive approximation is critical for superconductor stability against quench. Based on the applied numerical scheme, a hypothesis is suggested concerning correlation of the results of the simulation (the “numerical space”) with the experimental situation (the “physical reality”): Non-convergence of the numerical scheme might tightly be correlated with occurrence of a quench in the simulated superconductor.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Processing and Thermoelectric Properties of TiNiSn Materials: A Review

The half-Heusler alloy TiNiSn is a promising material for high-temperature thermoelectric applications due to its good thermal stability and semiconductor-like electrical properties. Numerous synthesis techniques have been utilized to make TiNiSn, but a comparative study on its thermoelectric properties with respect to processing parameters has not been reported. Based on published experimental data, primary melting methods, such as arc melting or induction levitation melting, do not have a noticeable effect on the final thermoelectric properties, although choice of densification technique and annealing parameters correlate with each thermoelectric property. Thermal conductivity (2.47-6.08 W/m K), governed by lattice scattering effects, was maximized with high density from the spark plasma sintering (SPS) technique as well with low inclusions of metallic impurities and interstitial Ni defects from low-temperature, non-densified processes. Electrical resistivity (4.75-30 µΩ m), inversely related to thermal conductivity, is minimized when density is maximized and chemical defects are favorable from using SPS. All reported TiNiSn alloys contain intermetallic impurities, imparting a global trend of lowered electrical resistivity. Seebeck coefficient (− 273.20 to − 50.71 µV/K), proportional to electrical resistivity, is maximized when charge carrier concentration is minimized by eliminating phase impurities and Ni defects. Porosity has been observed to scatter low-energy electrons, which can increase Seebeck coefficient and electrical resistivity, although no statistical significance is found. An ideal TiNiSn alloy should have intermediate values of all interconnected thermoelectric properties to maximize ZT. Future research work should strive to include quantitative phase analysis to better characterize the contributions of Ni defects and impurity phases toward thermoelectric properties.

36 MATERIALS SCIENCE↗
DOE OSTI2018

A front-fixing finite element method for pricing American options under regime-switching jump-diffusion models

A front-fixing finite element method is applied to solve the partial integro-differential equations (PIDEs) arising in pricing American options under Markov-modulated jump-diffusion models with the free boundaries feature. For this purpose, we used a front-fixing method to transform the pricing problem into a nonlinear parabolic integro-differential equation on a fixed domain. Then the variational form of the resulting problem is solved by a finite element method. Under some appropriate assumptions, we establish the stability of the method and illustrate some numerical results to examine the rate of convergence of the proposed method for the pricing problem and compare its accuracy to some recent works on pricing American options under regime-switching jump-diffusion models.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Plastic additives in the ocean: Use of a comprehensive dataset for meta-analysis and method development

In excess of 13,000 chemicals are added to plastics (‘additives’) to improve performance, durability, and production of plastic products. They are categorized into numerous chemical classes including flame retardants, light stabilizers, antioxidants, and plasticizers. While research on plastic additives in the marine environment has increased over the past decade, there is a lack of methodological standardization. To direct future measurement of plastic additives, we compiled a first-of-its-kind dataset of literature assessing plastic additives in marine environments, delineated by sample type (plastic debris, seawater, sediment, biota). Using this dataset, we performed a meta-analysis to summarize the state of the science. Currently, our dataset includes 217 publications published between 1978 and May 2023. The majority of publications analyzed plastic additives in biota collected from Europe and Asia. Analyses concentrated on plasticizers, brominated flame retardants, and bisphenols. Common sample preparation techniques included Solvent - Agitation extraction for plastic, sediment, and biota samples, and Solid Phase Extraction for seawater samples with dichloromethane and solvent mixtures including dichloromethane as the organic extraction solvent. Finally, most analyses were performed utilizing gas chromatography/mass spectrometry. There are a variety of data gaps illuminated by this meta-analysis, most notably the small number of compounds that have been targeted for detection compared to the large number of additives used in plastic production. The provided dataset facilitates future investigation of trends in plastic additive concentration data in the marine environment (allowing for comparison to toxicity thresholds) and acts as a starting point for optimizing and harmonizing plastic additive analytical methods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

High-order dimensionally-split Cartesian embedded boundary method for non-dissipative schemes

Centered finite-difference schemes are commonly used for high-fidelity turbulent flow simulations in canonical configurations because of their non-dissipative property and computational efficiency. However, their use in flow simulations over complex geometries is limited by the requirements of a structured grid and a stable boundary treatment in the absence of artificial (numerical) dissipation. Cartesian embedded boundary (EB) approaches provide an efficient structured-grid framework to apply difference schemes over complex domains. However, they are often restricted to low orders of accuracy because of numerical instabilities at the embedded boundaries and the issues of small-cell problem that are difficult to address with high-order accuracy. The present work discusses a systematic approach to obtain high-order EB methods with non-dissipative centered schemes in the interior. This approach, based on satisfying the primary and secondary conservation conditions, is employed to derive EB schemes that are up to sixth-order accurate in the interior and fourth-order accurate globally for hyperbolic, parabolic as well as incompletely parabolic problems. The proposed finite-difference discretization is, by construction, dimensionally split and addresses the small-cell problem without any cell/geometry transformations, thus, highly simplifying implementation in a flow solver. Various linear and non-linear numerical tests are performed to evaluate the stability and the accuracy of the proposed EB schemes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

A face-upwinded spectral element method

Here we present a new high-order accurate discretisation on unstructured meshes of quadrilateral elements. Our Face Upwinded Spectral Element (FUSE) method uses the same node distribution as a high-order continuous Galerkin (CG) method, but with a particular choice of node locations within each element and an upwinded stencil on the face nodes. This results in a number of benefits, including fewer degrees of freedom and straight-forward integration with CG. We present the derivation of the scheme and the analysis of its properties, in particular showing stability using von Neumann analysis. We show numerical evidence for its accuracy and efficiency on multiple classes of problems including convection-dominated flows, Poisson's equation, and the incompressible Navier-Stokes equations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

A high-order finite difference method for moving immersed domain boundaries and material interfaces

Here, we present a high-order sharp treatment of immersed moving domain boundaries and material interfaces, and apply it to the advection-diffusion equation in two and three dimensions. The spatial discretization combines dimension-split finite difference schemes with an immersed boundary treatment based on a weighted least-squares reconstruction of the solution, providing stable discretizations with up to sixth order accuracy for diffusion terms and third order accuracy for advection terms. The temporal discretization relies on a novel strategy for maintaining high-order temporal accuracy in problems with moving boundaries that minimizes implementation complexity and allows arbitrary explicit or diagonally-implicit Runge-Kutta schemes. The approach is broadly compatible with popular PDE-specialized Runge-Kutta time integrators, including low-storage, strong stability preserving, and diagonally implicit schemes. Through numerical experiments we demonstrate that the full discretization maintains high-order spatial and temporal accuracy in the presence of complex 3D geometries and for a range of boundary conditions, including Dirichlet, Neumann, and flux conditions with large jumps in coefficients.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Swelling properties of montmorillonite and beidellite clay minerals from molecular simulation: Comparison of temperature interlayer cation, and charge location effects

In this study, the swelling properties of smectite clay minerals are relevant to many engineering applications including environmental remediation, repository design for nuclear waste disposal, borehole stability in drilling operations, and additives for numerous industrial processes and commercial products. We used molecular dynamics and grand canonical Monte Carlo simulations to study the effects of layer charge location, interlayer cation, and temperature on intracrystalline swelling of montmorillonite and beidellite clay minerals. For a beidellite model with layer charge exclusively in the tetrahedral sheet, strong ion–surface interactions shift the onset of the two-layer hydrate to higher water contents. In contrast, for a montmorillonite model with layer charge exclusively in the octahedral sheet, weaker ion–surface interactions result in the formation of fully hydrated ions (two-layer hydrate) at much lower water contents. Clay hydration enthalpies and interlayer atomic density profiles are consistent with the swelling results. Water adsorption isotherms from grand canonical Monte Carlo simulations are used to relate interlayer hydration states to relative humidity, in good agreement with experimental findings.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Nonlinear periodic waves solutions of the nonlinear self-dual network equations

The new classes of periodic solutions of nonlinear self-dual network equations describing the breather and soliton lattices, expressed in terms of the Jacobi elliptic functions have been obtained. The dependences of the frequencies on energy have been found. Numerical simulations of soliton lattice demonstrate their stability in the ideal lattice and the breather lattice instability in the dissipative lattice. However, the lifetime of such structures in the dissipative lattice can be extended through the application of ac driving terms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Stability of a pinned magnetic domain wall as a function of its internal configuration

It is shown that there are many stable configurations for a domain wall pinned by a notch along a magnetic stripe. The stability of several of these configurations is investigated numerically as a function of the thickness of the magnetic film. The depinning mechanism depends on the structure of the domain wall and on the thickness of the magnetic film. In the case of a spin-valve structure, it appears that the stray fields emerging from the hard layer at the notch location influence the stability of the micromagnetic configuration. Different depinning mechanisms are thus observed for the same film thickness depending on the magnetization orientation of the propagating domain. This conclusion qualitatively explains experimental magnetoresistance measurements.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Effect of flexibility on the growth of concentration fluctuations in a suspension of sedimenting fibers: Particle simulations

Three-dimensional numerical simulations are performed to study the stability of a sedimenting suspension of weakly flexible fibers. It is well known that a suspension of rigid rods sedimenting under gravity at low Reynolds number is unstable to concentration fluctuations owing to hydrodynamic interactions. Flexible fibers, however, reorient while settling and even weak flexibility can alter their collective dynamics. In our recent work [Manikantan et al., “The instability of a sedimenting suspension of weakly flexible fibres,” J. Fluid Mech. 756, 935–964 (2014)], we developed a mean-field theory to predict the linear stability of such a system. Here, we verify these predictions using accurate and efficient particle simulations based on a slender-body model. We also demonstrate the mechanisms by which flexibility-induced reorientation alters suspension microstructure, and through it, its stability. Specifically, we first show that the anisotropy of the base state in the case of a suspension of flexible fibers has a destabilizing effect compared to a suspension of rigid rods. Second, a conflicting effect of flexibility is also shown to suppress particle clustering and slow down the growth of the instability. The relative magnitude of filament flexibility and rotational Brownian motion dictates which effect dominates, and our simulations qualitatively follow theoretically predicted trends. The mechanism for either effects is tied to the flexibility-induced reorientation of particles, which we illustrate using velocity and orientation statistics from our simulations. Finally, we also show that, in the case of an initially homogeneous and isotropic suspension, flexibility always acts to suppress the growth of the instability.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗