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 901 records · Page 50
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

A sharp interface model for deterministic simulation of dendrite growth

A high order level set model is developed for deterministic simulation of dendritic growth in unstable solidifying systems. The model captures motion of the front implicitly on a structured finite difference grid, enables calculation of its geometric properties and also applies boundary conditions on the immersed interface. Interfacial capillary effect is incorporated in the model through the Gibbs-Thomson condition. Canonical problems for evaluating grid convergence of the numerical method and validation tests for stability of a growing nucleus in the presence of isotropic surface tension are presented. The growth morphology of solidifying nuclei in undercooled metallic melts is quantitatively analyzed. Effects of crystal anisotropy and melt undercooling on the front geometry, propagation speed and formation of branched dendritic structures are examined. In conclusion, the complex morphological changes, such as remelting of secondary perturbations under specific conditions of undercooling, are also captured during the quantitative analysis.

36 MATERIALS SCIENCE↗
DOE OSTI2018

High-order, stable, and conservative boundary schemes for central and compact finite differences

Stable and conservative numerical boundary schemes are constructed such that they do not diminish the overall accuracy of the method for interior schemes of orders 4, 6, and 8 using both explicit (central) and compact finite differences. Previous attempts to develop stable numerical boundary schemes for non-linear problems have resulted in schemes which significantly reduced the global accuracy and/or required some form of artificial dissipation. Thus, the schemes developed in this paper are the first to not require this tradeoff, while also ensuring discrete conservation and allowing for direct boundary condition enforcement. After outlining a general procedure for the construction of conservative boundary schemes of any order, a simple, yet novel, optimization strategy which focuses directly on the compressible Euler equations is presented. Furthermore, the result of this non-linear optimization process is a set of high-order, stable, and conservative numerical boundary schemes which demonstrate excellent stability and convergence properties on an array of linear and non-linear hyperbolic problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Meshless discretization of the discrete-ordinates transport equation with integration based on Voronoi cells

The time-dependent, gray, linear radiation transport equation is discretized using the meshless local Petrov-Galerkin method with reproducing kernels. The integration is performed using a Voronoi tessellation, which creates a partition of unity that only depends on the position and extent of the kernels. The resolution of the integration automatically follows the particles and requires no manual adjustment. The discretization includes streamline-upwind Petrov-Galerkin stabilization to prevent oscillations and improve numerical conditioning. The angular quadrature is selectively refineable to increase angular resolution in chosen directions. The time discretization is done using backward Euler. The transport solve for each direction and the solve for the scattering source are both done using Krylov iterative methods. The results indicate first-order convergence in time and second-order convergence in space for linear reproducing kernels.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Traveling wave solutions in a chain of periodically forced coupled nonlinear oscillators

Artificial perceptions of light called phosphenes were motivated by earlier studies. We analyze traveling wave solutions in a chain of periodically forced coupled nonlinear oscillators modeling this phenomenon. We examine the discrete model problem in its co-traveling frame and systematically obtain the corresponding traveling waves in one spatial dimension. Direct numerical simulations as well as linear stability analysis are employed to reveal the parameter regions where the traveling waves are stable, and these waves are, in turn, connected to the standing waves analyzed in earlier work. We also consider a two-dimensional extension of the model and demonstrate the robust evolution and stability of planar fronts. Moreover, our simulations also suggest the radial fronts tend to either annihilate or expand and flatten out, depending on the phase value inside and the parameter regime. Finally, we observe that solutions that initially feature two symmetric fronts with bulged centers evolve in qualitative agreement with experimental observations of phosphenes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2017

Molecular Retrofitting Adapts a Metal–Organic Framework to Extreme Pressure

Despite numerous studies on chemical and thermal stability of metal-organic frameworks (MOFs), mechanical stability remains largely undeveloped. No strategy exists to control the mechanical deformation of MOFs under ultrahigh pressure, to date. We show that the mechanically unstable MOF-520 can be retrofitted by precise placement of a rigid 4,4'-biphenyldicarboxylate (BPDC) linker as a "girder" to afford a mechanically robust framework: MOF-520-BPDC. This retrofitting alters how the structure deforms under ultrahigh pressure and thus leads to a drastic enhancement of its mechanical robustness. While in the parent MOF-520 the pressure transmitting medium molecules diffuse into the pore and expand the structure from the inside upon compression, the girder in the new retrofitted MOF-520-BPDC prevents the framework from expansion by linking two adjacent secondary building units together. As a result, the modified MOF is stable under hydrostatic compression in a diamond-anvil cell up to 5.5 gigapascal. The increased mechanical stability of MOF-520-BPDC prohibits the typical amorphization observed for MOFs in this pressure range. Direct correlation between the orientation of these girders within the framework and its linear strain was estimated, providing new insights for the design of MOFs with optimized mechanical properties.

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

Catalytic production of tetrahydropyran (THP): a biomass-derived, economically competitive solvent with demonstrated use in plastic dissolution

Tetrahydropyran (THP) is a five-carbon heterocyclic ether that is non-carcinogenic, non-peroxide forming, biodegradable, and economically competitive with tetrahydrofuran (THF) as a solvent. In this work, THP has been synthesized from renewable biomass at >99.8% selectivity and 98% yield via hydrogenation of furfural-derived 3,4-dihydropyran (DHP) over Ni/SiO 2 in a continuous flow reactor at 150–200 °C. The apparent activation energy of THP formation is 31 kJ mol -1 , and the reaction orders with respect to the partial pressures of H 2 , DHP, and THP are: 2, 1, and -0.3. The kinetic data has been fitted to a Hougen–Watson model where the rate limiting step is the hydrogenation of adsorbed DHP. Ni/SiO 2 is shown to have a low deactivation rate constant of 0.012 h -1 over 100 h time on stream and can be regenerated in situ. As a performance advantage, THP is shown to be resistant to ring opening polymerization under strongly acidic conditions that THF is not, revealing it to be a superior solvent. Further, the minimum selling price of THP is competitive with the market price of THF ($$900 – 1400 per ton ) at a DHP feedstock cost of $1000 per ton. Conductor-like Screening Model for Real Solvents (COSMO-RS) and molecular dynamics (MD) simulations with 1008 solvents and 8 common plastics have demonstrated that THP can serve as an alternative solvent to THF, 2-methyltetrahydrofuran (MeTHF), and cyclopentyl methyl ether (CPME) for plastic dissolution, especially low-density polyethylene (LDPE), polypropylene (PP), polystyrene (PS) and polyvinyl chloride (PVC). This work establishes THP as a green solvent with excellent thermal, chemical and peroxidative stability that can be used for numerous applications, including waste plastic recycling.

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

The direct criterion of Newcomb for the ideal MHD stability of an axisymmetric toroidal plasma

A method is presented for determining the ideal magnetohydrodynamic stability of an axisymmetric toroidal plasma, based on a toroidal generalization of the method developed by Newcomb for fixed-boundary modes in a cylindrical plasma. For toroidal mode number n≠0, the stability problem is reduced to the numerical integration of a high-order complex system of ordinary differential equations, the Euler-Lagrange equation for extremizing the potential energy, for the coupled amplitudes of poloidal harmonics m as a function of the radial coordinate ψ in a straight-fieldline flux coordinate system. Unlike the cylindrical case, different poloidal harmonics couple to each other, which introduces coupling between adjacent singular intervals. A boundary condition is used at each singular surface, where m = nq and q(ψ) is the safety factor, to cross the singular surface and continue the solutions beyond it. Fixed-boundary instability is indicated by the vanishing of a real determinant of a Hermitian complex matrix constructed from the fundamental matrix of solutions, the generalization of Newcomb's crossing criterion. In the absence of fixed-boundary instabilities, an M × M plasma response matrix W P with M the number of poloidal harmonics used, is constructed from the Euler-Lagrange solutions at the plasma-vacuum boundary. This is added to a vacuum response matrix W V to form a total response matrix W T . Finally, the existence of negative eigenvalues of W T indicates the presence of free-boundary instabilities. The method is implemented in the fast and accurate DCON code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

The roll-up and merging of coherent structures in shallow mixing layers

The current study seeks a fundamental explanation to the development of two-dimensional coherent structures (2DCSs) in shallow mixing layers. A nonlinear numerical model based on the depth-averaged shallow water equations is used to investigate the temporal evolution of shallow mixing layers, where the mapping from temporal to spatial results is made using the velocity at the center of the mixing layers. The flow is periodic in the streamwise direction. Transmissive boundary conditions are used in the cross-stream boundaries to prevent reflections. Numerical results are compared to linear stability analysis, mean-field theory, and secondary stability analysis. Results suggest that the onset and development of 2DCS in shallow mixing layers are the result of a sequence of instabilities governed by linear theory, mean-field theory, and secondary stability theory. The linear instability of the shearing velocity gradient gives the onset of 2DCS. When the perturbations reach a certain amplitude, the flow field of the perturbations changes from a wavy shape to a vortical (2DCS) structure because of nonlinearity. The development of the vertical 2DCS does not appear to follow weakly nonlinear theory; instead, it follows mean-field theory. After the formation of 2DCS, separate 2DCSs merge to form larger 2DCS. In this way, 2DCSs grow and shallow mixing layers develop and grow in scale. The merging of 2DCS in shallow mixing layers is shown to be caused by the secondary instability of the 2DCS. Eventually 2DCSs are dissipated by bed friction. The sequence of instabilities can cause the upscaling of the turbulent kinetic energy in shallow mixing layers.

42 ENGINEERING↗
DOE OSTI2019

Design of semiconductor surface pits for fabrication of regular arrays of quantum dots and nanorings

We report a systematic computational study on the formation of complex nanostructures consisting of quantum dots and nanorings on surfaces of coherently strained thin films grown epitaxially on pit-patterned substrates. The analysis is based on self-consistent dynamical simulations according to a film surface evolution model that has been validated experimentally by comparison of its predictions with experimental observations on Ge films grown on Si pit-patterned substrates and is supported by linear stability theory that explains the film surface nanopattern formation as the outcome of a Stranski-Krastanow instability. Emphasis is placed on the design of conical pit patterns and the effects on the resulting film surface nanopattern of varying geometrical design parameters including film thickness, pit-pattern period, pit depth, pit opening diameter, and pit wall inclination. We demonstrate that varying the pit opening diameter and the pit wall slope leads to formation of complex nanostructures inside the pits of a regular pit pattern on the film surface, which include quantum dots, as well as single nanorings and multiple concentric nanorings that may or may not surround a central quantum dot inside each pit. Our simulation predictions demonstrate that the ordered nanostructure patterns forming on the film surface can be controlled precisely by tuning the geometrical parameters of the pits on the pit-patterned substrate. Our findings have important implications for designing optimal semiconductor surface patterns toward enabling future nanofabrication technologies.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

High order explicit local time stepping methods for hyperbolic conservation laws

In this paper we present and analyze a general framework for constructing high order explicit local time stepping (LTS) methods for hyperbolic conservation laws. In particular, we consider the model problem discretized by Runge-Kutta discontinuous Galerkin (RK-DG) methods and design LTS algorithms based on the strong stability preserving Runge-Kutta (SSP-RK) schemes, that allow spatially variable time step sizes to be used for time integration in different regions of the computational domain. The proposed algorithms are of predictor-corrector type, in which the interface information along the time direction is first predicted based on the SSP-RK approximations and Taylor expansions, and then the fluxes over the region of the interface are corrected to conserve mass exactly at each time step. Following the proposed framework, we detail the corresponding LTS schemes with accuracy up to the fourth order, and prove their conservation property and nonlinear stability for the scalar conservation laws. Numerical experiments are also presented to demonstrate excellent performance of the proposed LTS algorithms.

58 GEOSCIENCES↗
DOE OSTI2014

Quantifying the uncertainty of kinetic-theory predictions of clustering. Final Report covering 21 September 2011 - 20 September 2014

Previous work has indicated that inelastic grains undergoing homogeneous cooling may be unstable, giving rise to the formation of velocity vortices and particle clusters for sufficiently large systems. Such instabilities are observed in industrial coal and biomass gasifiers and are known to influence gas-solid contact area, mixing dynamics, and heat/mass transfer rates. However, the driving mechanisms that lead to vortices and clusters are not well understood. Discrete-particle simulations provide a well-established method for understanding such mechanisms but are not a feasible technique for predicting the behavior of large-scale systems. Kinetic-theory-based continuum models offer an effective means of describing such flows, and instabilities present a stringent test of such models due to the transient, three-dimensional nature of instabilities and the large range of time and length scales over which these mechanisms occur.This work begins with the study, via a combination of continuum models and discrete- particle simulations, of a relatively simple flow and includes additional complexities in a stepwise manner to assess various driving mechanisms. Comparisons with discrete-particle simulations, which offer detailed, well-established (but computationally limited) descriptions of particle flows, indicate the ability of continuum models to accurately incorporate each mechanism. Specifically, the critical length scale for velocity vortices and/or particle clusters are studied via direct numerical simulation, molecular dynamics simulations, linear stability analyses of the continuum model, and transient simulations of the continuum model in a range of flow complexities, including moderate dissipation and particle concentration, frictional particles collisions, high gradients, and gas-solid flows. Strong agreement between kinetic-theory-based continuum models and discrete-particle simulations is found for a range for conditions. Furthermore, discrete-particle simulations offer insight into the relative importance of various mechanisms for instability in both granular and gas-solid flows. Normal inelasticity and friction are both found to be important in granular flows. Particle collisions, thermal drag, and mean drag are found to be important in gas-solid flows.

42 ENGINEERING↗
DOE OSTI2016

Mechanical behavior of nanotwinned materials – experimental and computational approaches

Nanotwinned materials exhibit high strength combined with excellent thermal stability, making them potentially attractive for numerous applications. When deposited on cold substrates at high rates, for example, silver films can be prepared with a high-density of growth twins with an average twin boundary spacing of less than 10 nm. These films show a very strong {111} texture, with the twin boundaries being perpendicular to the growth direction. The origins of superior mechanical and thermal properties of nanotwinned materials, however, are not yet fully understood and need further improvements.

36 MATERIALS SCIENCE↗
DOE OSTI2019

High-Order, Stable, and Conservative Boundary Schemes for Central and Compact Finite Differences

Stable and conservative numerical boundary schemes are constructed such that they do not diminish the overall accuracy of the method for interior schemes of orders 4, 6, and 8 using both explicit (central) and compact finite differences. Previous attempts to develop stable numerical boundary schemes for non-linear problems have resulted in schemes which significantly reduced the global accuracy and/or required some form of artificial dissipation. Thus, the schemes developed are the first to not require this tradeoff, while also ensuring discrete conservation and allowing for direct boundary condition enforcement. After outlining a general procedure for the construction of conservative boundary schemes of any order, a simple, yet novel, optimization strategy which focuses directly on the compressible Euler equations is applied. The result of this non-linear optimization process is a set of high-order, stable, and conservative numerical boundary schemes which demonstrate excellent stability and convergence properties on an array of linear and non-linear hyperbolic problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Impact of Vibrational Nonequilibrium on the Supersonic Mode in Hypersonic Boundary Layers

The supersonic mode in hypersonic boundary layers has recently been shown to achieve a larger peak growth rate than the traditional second mode in particular circumstances relevant to free flight. In this environment, high levels of aerodynamic heating are encountered such that incorporation of chemical and vibrational nonequilibrium effects can become critical. The impact of thermal nonequilibrium (TNE) on the supersonic mode in particular has not been thoroughly and systematically investigated. This work uses thermochemical nonequilibrium (TCNE) direct numerical simulation (DNS) and TCNE linear stability theory (LST) to obtain a more complete investigation of the supersonic mode using both nonequilibrium and frozen gas models. The simulation is Mach 5 flow over a 1-mm-nose-radius axisymmetric cone 1 m in length, representative of experimental conditions in ground test facilities. Both LST and DNS results indicate that TNE effects are destabilizing compared with frozen flow and can lead to an N-factor difference of approximately +1 at the end of the cone. LST results also indicate that TNE effects are stabilizing to the second and supersonic modes in comparison with thermal equilibrium. The supersonic mode appeared to be slightly more affected by TNE than the second mode. In conclusion, collectively, the findings in this work suggest that even in ground test facilities where TCNE effects are relatively small compared with free flight, incorporating nonequilibrium models in both the mean flow and stability solvers can have a practical effect on hypersonic boundary-layer transition prediction analyses.

42 ENGINEERING↗
NASA NTRS1978

Modification of the Douglas Neumann program to improve the efficiency of predicting component interference and high lift characteristics

The Douglas Neumann method for low-speed potential flow on arbitrary three-dimensional lifting bodies was modified by substituting the combined source and doublet surface paneling based on Green's identity for the original source panels. Numerical studies show improved accuracy and stability for thin lifting surfaces, permitting reduced panel number for high-lift devices and supercritical airfoil sections. The accuracy of flow in concave corners is improved. A method of airfoil section design for a given pressure distribution, based on Green's identity, was demonstrated. The program uses panels on the body surface with constant source strength and parabolic distribution of doublet strength, and a doublet sheet on the wake. The program is written for the CDC CYBER 175 computer. Results of calculations are presented for isolated bodies, wings, wing-body combinations, and internal flow.

Bristow, D. R.↗
NASA NTRS1981

Weighted-mean scheme for solving incompressible viscous flow

The problem of how a boundary layer responds to the motion of a convexed vortex on a porous wall was investigated. The wall velocity is approximately given by Darcy's law. The vorticity-stream function approach was adopted for solving Navier-Stokes equations of two dimensional incompressible viscous flows. The weighted-mean scheme was used for constructing finite difference approximations of spatial derivatives. Several test problems were solved and numerical results demonstrate clearly the accuracy, stability, and efficiency of the scheme. The weighted mean scheme then can be applied to the vortical flow problem.

Huynh, Q. Q.↗