Engineering Papers⌕ Search

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.

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

At least 901 records · Page 50

Development and Breakdown of Goertler Vortices in High Speed Boundary Layers

The nonlinear development of G rtler instability over a concave surface gives rise to a highly distorted stationary flow in the boundary layer that has strong velocity gradients in both spanwise and wall-normal directions. This distorted flow is susceptible to strong, high frequency secondary instability that leads to the onset of transition. For high Mach number flows, the boundary layer is also subject to the second mode instability. The nonlinear development of G rtler vortices and the ensuing growth and breakdown of secondary instability, the G rtler vortex interactions with second mode instabilities as well as oblique second mode interactions are examined in the context of both internal and external hypersonic configurations using nonlinear parabolized stability equations, 2-D eigenvalue analysis and direct numerical simulation. For G rtler vortex development inside the Purdue Mach 6 Ludwieg tube wind tunnel, multiple families of unstable secondary eigenmodes are identified and their linear and nonlinear evolution is examined. The computation of secondary instability is continued past the onset of transition to elucidate the physical mechanisms underlying the laminar breakdown process. Nonlinear breakdown scenarios associated with transition over a Mach 6 compression cone configuration are also explored.

Li, Fei↗

Instability Wave–Streak Interactions in a High Mach Number Boundary Layer at Flight Conditions

The interaction of stationary streaks undergoing nonmodal growth with modally unstable instability waves in a hypersonic boundary-layer flow is studied using numerical computations. The geometry and flow conditions are selected to match a relevant trajectory location from the ascent phase of the HIFiRE-1 flight experiment; namely, a 7 degree half-angle, circular cone with 2.5 mm nose radius, freestream Mach number equal to 5.30, unit Reynolds number equal to 13.42 m -1 , and wall-to-adiabatic temperature ratio of approximately 0.35 over most of the vehicle. This paper investigates the nonlinear evolution of initially linear optimal disturbances that evolve into finite-amplitude streaks, followed by an analysis of the modal instability characteristics of the perturbed, streaky boundary-layer flow. The investigation is performed with stationary direct numerical simulations (DNS) and plane-marching parabolized stability equations (PSE), in conjunction with partial-differential-equation-based planar eigenvalue analysis. The overall effect of streaks is to reduce the peak amplification factors of instability waves, indicating a possible downstream shift in the onset of laminar-turbulent transition. The present study confirms previous findings that the mean flow distortion of the nonlinear streak perturbation reduces the amplification rates of the Mack-mode instability. More importantly, however, the present results demonstrate that the spanwise varying component of the streak can produce a larger effect on the Mack-mode amplification. The study with selected azimuthal wavenumbers for the stationary streaks reveals that a wavenumber of approximately 1.4 times larger than the optimal wavenumber is more effective in stabilizing the planar Mack-mode instabilities. In the absence of unstable first-mode waves for the present cold-wall condition, transition onset is expected to be delayed until the peak streak amplitude increases to nearly 35 percent of the freestream velocity, when intrinsic instabilities of the boundary-layer streaks begin to dominate the transition process. For streak amplitudes below that limit a significant net stabilization is achieved, yielding a potential transition delay that can exceed 100 percent of the length of the laminar region in the uncontrolled case.

Boundary layer stability↗

Fast Solution of Fully Implicit Runge--Kutta and Discontinuous Galerkin in Time for Numerical PDEs, Part II: Nonlinearities and DAEs

Fully implicit Runge--Kutta (IRK) methods have many desirable accuracy and stability properties as time integration schemes, but high-order IRK methods are not commonly used in practice with large-scale numerical PDEs because of the difficulty of solving the stage equations. This paper introduces a theoretical and algorithmic framework for solving the nonlinear equations that arise from IRK methods (and discontinuous Galerkin discretizations in time) applied to nonlinear numerical PDEs, including PDEs with algebraic constraints. Several new linearizations of the nonlinear IRK equations are developed, offering faster and more robust convergence than the often-considered simplified Newton, as well as an effective preconditioner for the true Jacobian if exact Newton iterations are desired. Inverting these linearizations requires solving a set of block 2 x 2 systems. Under quite general assumptions, it is proven that the preconditioned 2 x 2 operator's condition number is bounded by a small constant close to one, independent of the spatial discretization, spatial mesh, and time step, and with only weak dependence on the number of stages or integration accuracy. Moreover, the new method is built using the same preconditioners needed for backward Euler-type time stepping schemes, so can be readily added to existing codes. The new methods are applied to several challenging fluid flow problems, including the compressible Euler and Navier--Stokes equations, and the vorticity-streamfunction formulation of the incompressible Euler and Navier--Stokes equations. Up to 10th-order accuracy is demonstrated using Gauss IRK, while in all cases fourth-order Gauss IRK requires roughly half the number of preconditioner applications as required by standard Singly diagonally implicit Runge--Kutta methods.

97 MATHEMATICS AND COMPUTING↗

Diffused morphotropic phase boundary in relaxor-PbTiO 3 crystals: High piezoelectricity with improved thermal stability

Solid solution ferroelectrics are the most widely used piezoelectric material for numerous electromechanical applications, including sensors, actuators, and transducers. A milestone in ferroelectric research is the discovery of the morphotropic phase boundary that was first reported in Pb(Zr x Ti 1–x )O 3 , which has been extensively solicited to improve the performance of various solid solution ferroelectrics, including those having the highest piezoelectricity known today. However, due to the inherent correlation between phase transition and thermodynamic imbalance, the efforts of building the phase boundary encounter the challenge that high performance materials are generally accompanied by property instability. In this paper, we report a comprehensive study on the crystalline symmetry and polar nano-regions in relaxor-PbTiO 3 crystals by synchrotron measurements and property characterizations. In contrast to the common belief that the morphotropic phase boundary is a narrow composition region, the morphotropic phase boundary in ternary Pb(In 1/2 Nb 1/2 )O 3 -Pb(Mg 1/3 Nb 2/3 )O 3 -PbTiO 3 is demonstrated as an extensive region covering at least 9% of the PbTiO 3 content. Such a diffused morphotropic phase boundary is associated with the intensive interaction of polar nano-regions, leading to high piezoelectricity (>1500 pC/N) with greatly improved thermal stability, where the piezoelectric variation is ~90% over the temperature range of 273–373 K, which is about a factor of 3 lower compared to its binary counterpart Pb(Mg 1/3 Nb 2/3 )O 3 -PbTiO 3 . This work sheds light on the fundamental understanding of nanoscale inhomogeneity and macroscopic symmetry in relaxor-PbTiO 3 , and successfully interlinks the structure and properties in complex solid solutions. The existence of a diffused morphotropic phase boundary is also expected to benefit other ferroic systems beyond ferroelectrics, such as ferromagnetic and multiferroic materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

On the stability of Kuzmin/Toomre disks

A simple two-component system consiting of a Kuzmin/Toomre disk embedded in a spherically symmetrical 'halo' is studied with Plummer's density distribution. The bounds on the mass fraction of the disk component which result from applying the semiempirical criterion against bar instabilities due to Ostriker and Peebles are studied analytically and the results are compared with recent numerical simulations. It is found that the locus t = 0.14, calculated using analytic approximations for a variety of rotation laws, agrees reasonably well with the neutral stability line derived from numerical results.

Frank, Juhan↗

Oscillatory Thermocapillary Convection

Stability analysis of thermocapillary convection in rectangular cavities is performed using direct numerical simulation. The influence of the Reynolds number (Re), the fluid Prandtl number (Pr) and cavity aspect ratio (Ar) on the motion is investigated. Neutral stability curves for transition to time-dependent convection are delineated in the Re-Ar plane for fluids with Pr=1.0, 4.4, 6.78, and 10. Several interesting features of these diagrams are discussed. One important conclusion is that Ar(epsilon tau) increases as Pr decreases. Thus large values of both Ar and Re are necessary to induce thermocapillary oscillations for small Pr fluids such as liquid metals and semiconductor melts. Energy analysis is also performed for the oscillatory flow in the neighborhood of critical points in order to gain insight into the mechanisms leading to instability.

Xu, Jieyong↗

Long-term risk analysis associated with nuclear waste disposal in space

An assessment and verification of previous analytic results on the long term risk of earth reentry for hazardous payloads is presented. The two areas were studied: (1) stability of nominal, near-circular storage orbits in the regions between Venus and earth and between earth and Mars, and (2) probability of earth reentry for off-nominal planet-crossing orbits resulting from deployment system failures. In the first area, numerical integrations of the equations of motion are compared with stability predications based on secular perturbation theory. The agreement is good in terms of the heliocentric distances covered and the general behavior of the orbital history, although certain near-resonance situations can lead to difficulty. In the second area, a Monte Carlo simulation of orbital evolution is used and the results compared with Opik's analytic theory of planetary encounters and collision statistics, with data verified to within a close order-of-magnitude.

Friedlander, A. L.↗

Stability analysis of a one-dimensional multiphysics model of a molten salt fast reactor

Reactor designs with a fast neutron spectrum have recently moved into the research focus for possible Generation IV reactor concepts - a representative of which is the molten salt fast reactor (MSFR). In order to assess the safety and reliability of this reactor concept, it is necessary to obtain knowledge of the system behaviour. An important part of this research concerns the dynamical stability of the MSFR. Research regarding the MSFR stability so far has either been using linearized models for stability criteria or only been examining certain transients at fixed parameter values. This work delivers a comprehensive stability analysis for a non-linearized MSFR model, and a wide range of values for all relevant parameters. For this purpose, a one-dimensional MSFR model was set up, taking neutron kinetics and thermal hydraulics in the form of a system of differential equations into account. The stability of this model was investigated by means of the numerical tool MATCONT, which was used to monitor the evolution of a so-called fixed-point solution, here referring to the steady state at operating conditions of the system, while system parameters got varied. The MATCONT analyses showed no loss of stability for any of the considered parameter variations and no solutions that might co-exist in parallel with the steady-state fixed-point solution. Therefore these results indicate a stable fixed point to which all solution transients converge, making the considered MSFR model insusceptible to deviations from the equilibrium. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

The properties of close multiple stars

A reexamination of the properties of close multiple systems is conducted, taking into account recent numerical, analytical, and observational results. The orbital elements of 43 spectroscopic multiple systems are presented in a table. Tables with photometry data of the systems and data describing their properties are also provided. Attention is given to period ratios, mass ratios, orbital coplanarity, stability, lithium abundances and age estimates, and prospects of future observations. Numerical studies indicate that fragmentation produces binary components with mass ratios usually close to 1.0, while fission usually produces mass ratios in the range 0.1-0.5. Using analytical stability criteria, all 27 close multiple systems considered are stable if they are corotating. If they are counter-rotating, four systems are possibly unstable. However, other factors suggest that even these four are probably stable. The solar-type close multiple systems ranged in age from 100 million years to 2,000 million years.

Fekel, F. C., Jr.↗

Combined Experimental and Computational Efforts to Establish Ion Mobility, Solubility and Stability of Functional Liquids for Electrochemical Energy Storage

This work provides a computation-driven investigation of the stability of organic electrolytes for lithium-air batteries. Electrolyte instability is currently a key challenge that limits practical use of aprotic Li-air batteries, and the chemical processes that cause this instability are often kinetically-driven. Computational screening for kinetic stability involves the determination of reaction barriers for the numerous potential reaction mechanisms, barriers that are challenging to calculate due to the difficulty of locating transition state structures. Here we screen a broad set of substituted electrolytes for susceptibility to nucleophilic attack by superoxide. We find that carbonates are not typically expected to be stable and that sulfones are generally stable, validating literature trends. We study the effects of chemical functionalization with electron-donating and withdrawing groups and their interplay with steric factors, identifying functional groups and other chemical modifications that increase stability in these groups. User-input driven transition state identification is used for these initial calculations, and an automated computational pipeline is subsequently presented and validated as a means to perform further high-throughput searches across mechanisms and chemistries. The pipeline integrates cheminformatics-based reaction encoding, relaxed potential energy scans, and nudged elastic band calculations for an end-to-end approach to barrier calculations. We review this automated search approach and its current limitations, and discuss challenges and further work.

25 ENERGY STORAGE↗

Ameliorating the Courant-Friedrichs-Lewy condition in spherical coordinates: A double FFT filter method for general relativistic MHD in dynamical spacetimes

Numerical simulations of merging compact objects and their remnants form the theoretical foundation for gravitational wave and multimessenger astronomy. While Cartesian-coordinate-based adaptive mesh refinement is commonly used for simulations, spherical-like coordinates are more suitable for nearly spherical remnants and azimuthal flows due to lower numerical dissipation in the evolution of fluid angular momentum, as well as requiring fewer numbers of computational cells. However, the use of spherical coordinates to numerically solve hyperbolic partial differential equations can result in severe Courant-Friedrichs-Lewy (CFL) stability condition time step limitations, which can make simulations prohibitively expensive. This paper addresses this issue for the numerical solution of coupled spacetime and general relativistic magnetohydrodynamics evolutions by introducing a double fast Fourier transform (FFT) filter and implementing it within the fully message passing interface (mpi)-parallelized sphericalnr framework in the einstein toolkit. In conclusion, we demonstrate the effectiveness and robustness of the filtering algorithm by applying it to a number of challenging code tests, and show that it passes these tests effectively, demonstrating convergence while also increasing the time step significantly compared to unfiltered simulations.

79 ASTRONOMY AND ASTROPHYSICS↗

The influence of dynamic inflow and torsional flexibility on rotor damping in forward flight from symbolically generated equations

The combined effects of blade torsion and dynamic inflow on the aeroelastic stability of an elastic rotor blade in forward flight are studied. The governing sets of equations of motion (fully nonlinear, linearized, and multiblade equations) used in this study are derived symbolically using a program written in FORTRAN. Stability results are presented for different structural models with and without dynamic inflow. A combination of symbolic and numerical programs at the proper stage in the derivation process makes the obtainment of final stability results an efficient and straightforward procedure.

Reddy, T. S. R.↗

Comparison between stability limits for satellite motion

Three methods of obtaining stability information on satellite motion are compared by means of numerical and analytical computations. The model of the restricted problem of three bodies is used to describe the motion. Kuiper's (1961) approximate results, Szebehely's (1978) approximate results, and computer solutions obtained by successive iterations show close agreement regarding the maximum values of the orbital radii for stability. The lowest value, i.e., the most conservative estimate, is provided by the simplified form of Szebehely's formula.

Szebehely, V.↗

Numerical simulation of non-Newtonian free shear flows

Free shear flows, like those of mixing layers, are encountered in aerodynamics, in the atmosphere, and in the ocean as well as in many industrial applications such as flow reactors or combustion chambers. It is, therefore, crucial to understand the mechanisms governing the process of transition to turbulence in order to predict and control the evolution of the flow. Delaying transition to turbulence as far downstream as possible allows a gain in energy expenditure while accelerating the transition can be of interest in processes where high mixing is desired. Various methods, including the use of polymer additives, can be effective in controlling fluid flows. The drag reduction obtained by the addition of small amounts of high polymers has been an active area of research for the last three decades. It is now widely believed that polymer additives can affect the stability of a large variety of flows and that dilute solutions of these polymers have been shown to produce drag reductions of over 80 percent in internal flows and over 60 percent in external flows under a wide range of conditions. The major thrust of this work is to study the effects of polymer additives on the stability of the incompressible mixing layer through large scale numerical simulations. In particular, we focus on the two dimensional flow and examine how the presence of viscoelasticity may affect the typical structures of the flow, namely roll-up and pairing of vortices.

Homsy, G. M.↗

Advanced Characterisation of Aerosol Size Properties from Measurements of Spectral Optical Depth Using the GRASP Algorithm

This study evaluates the potential of using aerosol optical depth (τ (sub a)) measurements to characterise the microphysical and optical properties of atmospheric aerosols. With this aim, we used the recently developed GRASP (Generalized Retrieval of Aerosol and Surface Properties) code for numerical testing of six different aerosol models with different aerosol loads. The direct numerical simulations (self-consistency tests) indicate that the GRASP-AOD retrieval provides modal aerosol optical depths (fine and coarse) to within 0.01 of the input values. The retrieval of the fine-mode radius, width and volume concentration are stable and precise if the real part of the refractive index is known. The coarse-mode properties are less accurate, but they are significantly improved when additional a priori information is available. The tests with random simulated errors show that the uncertainty in the bimodal log-normal size distribution parameters increases as the aerosol load decreases. Similarly, the reduction in the spectral range diminishes the stability of the retrieved parameters. In addition to these numerical studies, we used optical depth observations at eight AERONET locations to validate our results with the standard AERONET inversion products. We found that bimodal log-normal size distributions serve as useful input assumptions, especially when the measurements have inadequate spectral coverage and/or limited accuracy, such as moon photometry. Comparisons of the mode median radii between GRASP-AOD and AERONET indicate average differences of 0.013µm for the fine mode and typical values of 0.2–0.3µm for the coarse mode. The dominant mode (i.e. fine or coarse) indicates a 10 % difference in mode radii between the GRASP-AOD and AERONET inversions, and the average of the difference in volume concentration is around 17 % for both modes. The retrieved values of the fine-mode τ(sub a)(500) using GRASP-AOD are generally between those values obtained by the standard AERONET inversion and the values obtained by the AERONET spectral deconvolution algorithm (SDA), with differences typically lower than 0.02 between GRASP-AOD and both algorithms. Finally, we present some examples of application of GRASP-AOD inversion using moon photometry and the airborne PLASMA sun photometer during the ChArMEx summer 2013 campaign in the western Mediterranean.

Torres, Benjamin↗

ON THE ROLE OF INVOLUTIONS IN THE DISCONTINUOUS GALERKIN DISCRETIZATION OF MAXWELL AND MAGNETOHYDRODYNAMIC SYSTEMS

The role of involutions in energy stability of the discontinuous Galerkin (DG) discretization of Maxwell and magnetohydrodynamic (MHD) systems is examined. Important differences are identified in the symmetrization of the Maxwell and MHD systems that impact the construction of energy stable discretizations using the DG method. Specifically, general sufficient conditions to be imposed on the DG numerical flux and approximation space are given so that energy stability is retained These sufficient conditions reveal the favorable energy consequence of imposing continuity in the normal component of the magnetic induction field at interelement boundaries for MHD discretizations. Counterintuitively, this condition is not required for stability of Maxwell discretizations using the discontinuous Galerkin method.

Barth, Timothy↗

Comparison of approximate and numerical analyses of nonlinear combustion instability

At the present time, there are three general analytical techniques available to study problems of unsteady motions in rocket motors: linear stability analysis; approximate nonlinear analysis, founded on examining the behavior of coupled normal modes; and numerical calculations based on the conservation equations for one-dimensional flows. The last two yield the linear results as a limit. It is the main purpose of this paper to check the accuracy of the approximate analysis against the numerical analysis for some special cases. The results provide some justification for using the approximate analysis to study three-dimensional problems.

Culick, F. E. C.↗