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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Submarine landslides triggered by destabilization of high-saturation hydrate anomalies: LANDSLIDES CAUSED BY HYDRATE ANOMALIES

Submarine landslides occur along continental margins at depths that often intersect the gas hydrate stability zone, prompting suggestions that slope stability may be affected by perturbations that arise from changes in hydrate stability. Here we develop a numerical model to identify the conditions under which the destabilization of hydrates results in slope failure. Specifically, we focus on high-saturation hydrate anomalies at fine-grained to coarse-grained stratigraphic boundaries that can transmit bridging stresses that decrease the effective stress at sediment contacts and disrupt normal sediment consolidation. We evaluate slope stability before and after hydrate destabilization. Hydrate anomalies act to significantly increase the overall slope stability due to large increases in effective cohesion. However, when hydrate anomalies destabilize there is a loss of cohesion and increase in effective stress that causes the sediment grains to rapidly consolidate and generate pore pressures that can either trigger immediate slope failure or weaken the surrounding sediment until the pore pressure diffuses away. In cases where failure does not occur, the sediment can remain weakened for months. In cases where failure does occur, we quantify landslide dynamics using a rate and state frictional model and find that landslides can display either slow or dynamic (i.e., catastrophic) motion depending on the rate-dependent properties, size of the stress perturbation, and the size of the slip patch relative to a critical nucleation length scale. Our results illustrate the fundamental mechanisms through which the destabilization of gas hydrates can pose a significant geohazard.

58 GEOSCIENCES↗
DOE OSTI2018

Changes in tropical cyclones under stabilized 1.5 and 2.0 °C global warming scenarios as simulated by the Community Atmospheric Model under the HAPPI protocols

The United Nations Framework Convention on Climate Change (UNFCCC) invited the scientific community to explore the impacts of a world in which anthropogenic global warming is stabilized at only 1.5°C above preindustrial average temperatures. In this paper, we present a projection of future tropical cyclone statistics for both 1.5 and 2.0°C stabilized warming scenarios with direct numerical simulation using a high-resolution global climate model. As in similar projections at higher warming levels, we find that even at these low warming levels the most intense tropical cyclones become more frequent and more intense, while simultaneously the frequency of weaker tropical storms is decreased. We also conclude that in the 1.5°C stabilization, the effect of aerosol forcing changes complicates the interpretation of greenhouse gas forcing changes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

Analysis of a Stabilized CNLF Method with Fast Slow Wave Splittings for Flow Problems

In this work, we study Crank-Nicolson leap-frog (CNLF) methods with fast-slow wave splittings for Navier-Stokes equations (NSE) with a rotation/Coriolis force term, which is a simplification of geophysical flows. We propose a new stabilized CNLF method where the added stabilization completely removes the method's CFL time step condition. A comprehensive stability and error analysis is given. We also prove that for Oseen equations with the rotation term, the unstable mode (for which u(n+1) + u(n-1) equivalent to 0) of CNLF is asymptotically stable. Numerical results are provided to verify the stability and the convergence of the methods.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1977

Stability analysis for laminar flow control, part 1

The basic equations for the stability analysis of flow over three dimensional swept wings are developed and numerical methods for their solution are surveyed. The equations for nonlinear stability analysis of three dimensional disturbances in compressible, three dimensional, nonparallel flows are given. Efficient and accurate numerical methods for the solution of the equations of stability theory were surveyed and analyzed.

Benney, D. J.↗
DOE OSTI2018

Laser Acceleration of Nuclei for Synthesis of Superheavy Elements in the Stability Island of the Periodic Table

Within a three-dimensionalmodel of interaction of a powerful electromagnetic fieldwith plasma of solid-state density, the effect of acceleration of calcium nuclei ({sub 20}{sup 48}Ca) to high energies that are likely to be sufficient for the synthesis of superheavy elements within the “stability island” of the periodic table of elements is obtained by numerically solving the Vlasov–Maxwell equations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

STABILITY OF SATELLITES IN CLOSELY PACKED PLANETARY SYSTEMS

We perform numerical integrations of four-body (star, planet, planet, satellite) systems to investigate the stability of satellites in planetary systems with tightly packed inner planets (STIPs). We find that the majority of closely spaced stable two-planet systems can stably support satellites across a range of parameter-space which is only slightly decreased compared to that seen for the single-planet case. In particular, circular prograde satellites remain stable out to ∼0.4 R{sub H} (where R{sub H} is the Hill radius) as opposed to 0.5 R{sub H} in the single-planet case. A similarly small restriction in the stable parameter-space for retrograde satellites is observed, where planetary close approaches in the range 2.5-4.5 mutual Hill radii destabilize most satellites orbits only if a ∼ 0.65 R{sub H} . In very close planetary pairs (e.g., the 12:11 resonance) the addition of a satellite frequently destabilizes the entire system, causing extreme close approaches and the loss of satellites over a range of circumplanetary semi-major axes. The majority of systems investigated stably harbored satellites over a wide parameter-space, suggesting that STIPs can generally offer a dynamically stable home for satellites, albeit with a slightly smaller stable parameter-space than the single-planet case. As we demonstrate that multi-planet systems are not a priori poor candidates for hosting satellites, future measurements of satellite occurrence rates in multi-planet systems versus single-planet systems could be used to constrain either satellite formation or past periods of strong dynamical interaction between planets.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1985

Three-dimensional baroclinic instability of a Hadley cell for small Richardson number

A three-dimensional, linear stability analysis of a baroclinic flow for Richardson number, Ri, of order unity is presented. The model considered is a thin horizontal, rotating fluid layer which is subjected to horizontal and vertical temperature gradients. The basic state is a Hadley cell which is a solution of the complete set of governing, nonlinear equations and contains both Ekman and thermal boundary layers adjacent to the rigid boundaries; it is given in a closed form. The stability analysis is also based on the complete set of equations; and perturbation possessing zonal, meridional, and vertical structures were considered. Numerical methods were developed for the stability problem which results in a stiff, eighth-order, ordinary differential eigenvalue problem. The previous work on three-dimensional baroclinic instability for small Ri was extended to a more realistic model involving the Prandtl number, sigma, and the Ekman number, E, and to finite growth rates and a wider range of the zonal wavenumber.

Antar, B. N.↗
NASA NTRS2015

Calculated Low-Speed Steady and Time-Dependent Aerodynamic Derivatives for Some Airfoils Using a Discrete Vortex Method

This paper contains a collection of some results of four individual studies presenting calculated numerical values for airfoil aerodynamic stability derivatives in unseparated inviscid incompressible flow due separately to angle-of-attack, pitch rate, flap deflection, and airfoil camber using a discrete vortex method. Both steady conditions and oscillatory motion were considered. Variables include the number of vortices representing the airfoil, the pitch axis / moment center chordwise location, flap chord to airfoil chord ratio, and circular or parabolic arc camber. Comparisons with some experimental and other theoretical information are included. The calculated aerodynamic numerical results obtained using a limited number of vortices provided in each study compared favorably with thin airfoil theory predictions. Of particular interest are those aerodynamic results calculated herein (such as induced drag) that are not readily available elsewhere.

Riley, Donald R.↗
DOE OSTI2025

Kelvin–Helmholtz instability under stabilizing parallel magnetic field in nonhomogeneous compressible MHD flows

We study the Kelvin–Helmholtz instability (KHI) for the general case of a compressible, nonhomogeneous, magnetized plasma flow. The study is limited to a vortex sheet interface with an imposed parallel magnetic field. We introduce a new formalism based on a convective Mach number M c , a convective Alfvénic Mach number M Ac , and a total convective Mach number that combines the two. We derive an analytic expression of the KHI growth rate for a homogeneous flow (i.e., zero Atwood number, A=0) that converges toward both the expression for unmagnetized compressible flow and Chandrasekhar's expression for magnetized incompressible flow. Otherwise, the dispersion relation is solved numerically and allows deriving general stability diagrams of magnetized KHI for the triplet (A, M c , β −plasma) parameters. We show these parameters uniquely define all configurations for a parallel magnetic field. We also construct diagrams with respect to the convective Alfvénic Mach number, the β − plasma parameter, or the magnetic field showing which magnetic field strength is required for stabilizing a given shear flow. The theoretical growth rates are compared with 18 simulations made with the GAMERA code, currently used for 3D magnetospheric simulations. Finally, we apply our results to the analysis of a past KHI experiment performed at the OMEGA laser facility, showing linear theory succeeds to provide accurate estimates of the growth rate at early times. We further discuss how our results can inform future experiments in the high-Mach magnetized regime at the National Ignition Facility. Possible limitations of the study due to resistive, mixing, or turbulence effects are discussed.

compressible flows↗
NASA NTRS2018

Hypersonic Chemically Reacting Boundary-Layer Stability using LASTRAC

Hypersonic boundary layer transition is critical to the design of all hypersonic vehicles due to its effect on the heat transfer into the vehicle surface and potential drag enhancement or reduction during reentry. Boundary layer transition and boundary layer stability analysis under hypersonic conditions has been studied for decades, yet there is ample room for improved accuracy and further investigations into the relevant phenomena. In this work, we present a recent implementation of chemical equilibrium, finite-rate chemistry, and thermochemical nonequilibrium capabilities into LASTRAC, an existing well-established boundary-layer stability analysis code. Verification against existing numerical results in the literature are presented. LASTRAC was previously able to address calorically perfect flows. By using solutions of the Parabolized Stability Equations (PSE) with chemical and thermal nonequilibrium, we are able to investigate the effects of chemical and thermal nonequilibrium on a variety of phenomena including stationary crossflow instability on a swept wing and 2nd mode instabilities over a wedge.

Kline, H. L.↗
NASA NTRS1997

Development of Curved-Plate Elements for the Exact Buckling Analysis of Composite Plate Assemblies Including Transverse-Shear Effects

The analytical formulation of curved-plate non-linear equilibrium equations including transverse-shear-deformation effects is presented. The formulation uses the principle of virtual work. A unified set of non-linear strains that contains terms from both physical and tensorial strain measures is used. Linearized, perturbed equilibrium equations (stability equations) that describe the response of the plate just after buckling occurs are then derived after the application of several simplifying assumptions. These equations are then modified to allow the reference surface of the plate to be located at a distance z(sub c) from the centroidal surface. The implementation of the new theory into the VICONOPT exact buckling and vibration analysis and optimum design computer program is described as well. The terms of the plate stiffness matrix using both Classical Plate Theory (CPT) and first-order Shear-Deformation Plate Theory (SDPT) are presented. The necessary steps to include the effects of in-plane transverse and in-plane shear loads in the in-plane stability equations are also outlined. Numerical results are presented using the newly implemented capability. Comparisons of results for several example problems with different loading states are made. Comparisons of analyses using both physical and tensorial strain measures as well as CPT and SDPF are also made. Results comparing the computational effort required by the new analysis to that of the analysis currently in the VICONOPT program are presented. The effects of including terms related to in-plane transverse and in-plane shear loadings in the in-plane stability equations are also examined. Finally, results of a design-optimization study of two different cylindrical shells subject to uniform axial compression are presented.

McGowan, David Michael↗
NASA NTRS2016

Direct Numerical Simulation of Transition Due to Traveling Crossflow Vortices

Previous simulations of laminar breakdown mechanisms associated with stationary crossflow instability over a realistic swept-wing configuration are extended to investigate the alternate scenario of transition due to secondary instability of traveling crossflow modes. Earlier analyses based on secondary instability theory and parabolized stability equations have shown that this alternate scenario is viable when the initial amplitude of the most amplified mode of the traveling crossflow instability is greater than approximately 0.03 times the initial amplitude of the most amplified stationary mode. The linear growth predictions based on the secondary instability theory and parabolized stability equations agree well with the direct numerical simulation. Nonlinear effects are initially stabilizing but subsequently lead to a rapid growth followed by the onset of transition when the amplitude of the secondary disturbance exceeds a threshold value. Similar to the breakdown of stationary vortices, the transition zone is rather short and the boundary layer becomes completely turbulent across a distance of less than 15 times the boundary layer thickness at the completion of transition.

Li, Fei↗
NASA NTRS1983

Transformation of two and three-dimensional regions by elliptic systems

Finite difference methods for composite grids were analyzed. It was observed that linear interpolation between grids would suffice only where low order accuracy was required. In the context of fluid flow, this would be in regions where the flow was essentially free stream. Higher order interpolation schemes were also investigated. The well known quadratic and cubic interpolating polynomials would increase the formal accuracy of the overall numerical algorithm. However, it can also be shown that the stability of the algorithm may be adversely affected. Further numerical results are needed in order to assess the nature of this instability induced by the interpolation procedure. Finally, error analysis and the order of difference expressions on general curvilinear coordinates are discussed.

Thompson, J. F.↗
NASA NTRS2006

Geometrically Nonlinear Static Analysis of 3D Trusses Using the Arc-Length Method

Rigorous analysis of geometrically nonlinear structures demands creating mathematical models that accurately include loading and support conditions and, more importantly, model the stiffness and response of the structure. Nonlinear geometric structures often contain critical points with snap-through behavior during the response to large loads. Studying the post buckling behavior during a portion of a structure's unstable load history may be necessary. Primary structures made from ductile materials will stretch enough prior to failure for loads to redistribute producing sudden and often catastrophic collapses that are difficult to predict. The responses and redistribution of the internal loads during collapses and possible sharp snap-back of structures have frequently caused numerical difficulties in analysis procedures. The presence of critical stability points and unstable equilibrium paths are major difficulties that numerical solutions must pass to fully capture the nonlinear response. Some hurdles still exist in finding nonlinear responses of structures under large geometric changes. Predicting snap-through and snap-back of certain structures has been difficult and time consuming. Also difficult is finding how much load a structure may still carry safely. Highly geometrically nonlinear responses of structures exhibiting complex snap-back behavior are presented and analyzed with a finite element approach. The arc-length method will be reviewed and shown to predict the proper response and follow the nonlinear equilibrium path through limit points.

Hrinda, Glenn A.↗
DOE OSTI2022

Current state and future trends in boundary layer control on lifting surfaces

Successful flow control may bring numerous benefits, such as flow stabilization, flow reattachment, separation delay, drag reduction, lift increase, aerodynamic performance improvement, energy efficiency increase, shock delay or weakening, noise reduction, etc. For these purposes, many different flow control devices, which can be classified as passive, semi-active and active, have been designed and tested. This review paper aims to highlight the most promising and commonly employed boundary layer control methods as well as outline their potential in specific applications in aerospace and energy engineering. Referenced studies, performed on various geometries (flat plates, channels, airfoils, wings, blades, cylinders), are primarily numerical or experimental. Although enhanced aerodynamic performance is achieved in many cases, further research is required to draw general conclusions. This paper aims to demonstrate that, in the future, we may expect further developments of flow control actuators, as well as their increased application.

42 ENGINEERING↗
NASA NTRS1984

Development of numerical procedures for analysis of complex structures

The paper is concerned with the development of novel numerical procedures for the solution of static, stability, free vibration and dynamic response analysis of large, complex practical structures. Thus, details of numerical algorithms evolved for dynamic analysis of usual non-rotating and also rotating structures as well as finite dynamic elements are presented in the paper. Furthermore, the article provides some description of a general-purpose computer program STARS specifically developed for efficient analysis of complex practical structures.

Gupta, K. K.↗
NASA NTRS1987

Stability of axisymmetric boundary layers on sharp cones at hypersonic Mach numbers

The stability of sharp-cone boundary layers at zero angle of attack is investigated. Standard linear stability theory is used to perform a numerical study at an edge Mach number M(e) of 6.8 of normal-mode stability characteristics on an adiabatic-wall cone with special reference to the conditions of the stability experiment of Stetson et al. (1983). Comparisons of the calculations with experimental measurements bring out major areas of disagreement which remain to be resolved even in this simple case. Finally, a series of calculations of both cone and flat-plate N factors at M(e) of 4.5, 5.8, and 6.8 is used to study the relation between transition on a cone and flat plate based on stability theory.

Mack, Leslie M.↗
NASA NTRS1982

The dependence of the shape and stability of captive rotating drops on multiple parameters

Asymptotic and numerical techniques in bifurcation theory are applied to the Young-Laplace equation governing meniscus shape in order to analyze the dependence of the shape and stability of rigidly rotating drops held captive between corotating solid faces on multiple parameters. Asymptotic analysis of the evolution of drop shape from the cylindrical as a function of distance between the solid faces, drop volume, rotational Bond number and gravitational Bond number shows that some shape bifurcations from cylinders to wavy, axisymmetric menisci are ruptured by small changes in drop volume or gravity. Computer calculations of axisymmetric drop shapes based on a finite element representation of the interface and numerical algorithms for tracking shape families and singular points are then used to map drop stability for the four-dimensional parameter space. The results of the asymptotic and numerical analyses are shown to agree well within the limited range of parameters where the asymptotic analysis is valid.

Ungar, L. H.↗