Large deflection and stability analysis by the direct stiffness method
Matrix numerical procedure by digital computer equipment, and application to elastic stability and deflection problems
ELASTIC STABILITY↗
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.
Dynamics and Statics of Nonaxisymmetric Liquid Bridges
Theoretical and experimental investigation of the stability of nonaxisymmetric and nonaxisymmetric bridges contained between equal and unequal radii disks as a function of Bond and Weber number with emphasis on the transition from unstable axisymmetric to stable nonaxisymmetric shapes, are conducted. Numerical analysis of the stability of nonaxisymmetric bridges between unequal disks for various orientations of the gravity vector is performed. Experimental and theoretical investigation of large (nonaxisymmetric) oscillations and breaking of liquid bridges are also conducted.
Alexander, J. Iwan D.↗
Edge flame speed analysis of a swirl-stabilized spray flame using sustainable aviation fuels: A direct numerical simulation study
Not Available
Soriano, Bruno S. (ORCID:0000000331566499)↗
Finite volume computation of unsteady inviscid rotational transonic flows past airfoils in rigid body motion
Unsteady inviscid transonic flow over airfoils in arbitrary rigid body motion is analyzed numerically by solving the two-dimensional unsteady Euler equations in integral form using a finite volume scheme. The solution procedure is based on an explicit Runge-Kutta time-stepping scheme wherein the spatial terms are central-differenced and a combination of second- and fourth-differences in the flow variables are used to form the numerical dissipation terms to stabilize the scheme. Unsteady calculations are started from converged steady-state solutions as initial conditions. Nonreflective boundary conditions are imposed on the far-field boundaries. Results are presented and, where possible, validated against available numerical and experimental data for airfoils subjected to a step change in angle of attack, airfoils oscillating and plunging in transonic flow, and airfoils immersed in a time-varying free stream.
Damodaran, Murali↗
A numerical and analytical study of nonlinear bifurcations associated with the morphological stability of two-dimensional single crystals
The nonlinear stability of a two-dimensional single crystal of pure material in an undercooled melt is studied both analytically and numerically. The quasi-steady state approximation is used for the thermal fields, and the effects of different solid and liquid thermal conductivities, isotropic interfacial growth kinetics, and isotropic surface tension are included. The bifurcation analysis is performed by calculating the instantaneous value of the fundamental component of the local normal growth speed for an interface perturbed by a single Fourier shape component. Numerically, the fundamental component of the interfacial growth speed is found by Fourier analysis of the solution to an integrodifferential equation obeyed at the interface. Analytically, an expansion technique is used to derive a solvability condition defining each of these bifurcation points. The analytical and numerical results are in very close agreement. Almost all of the bifurcations are subcritical, and the results are presented by giving values of the Landau coefficient as a function of the different dimensionless parameters used in the model.
Brush, L. N.↗
Libration point staging options for SEI lunar missions - Station keeping implications
Libration point staging is one alternative to lunar orbit rendezvous which is being examined for the Space Exploration Initiative lunar missions. The libration point staging strategy being considered would eliminate the lunar mid-latitude accessibility constraint for long duration missions. Staging options for such a strategy are discussed with particular emphasis on the impact of the libration points' stability. The equations of motion for a satellite at the libration points are numerically integrated to investigate the stability and to determine the station-keeping requirements. These requirements are presented, and overall requirements for libration point staging briefly summarized. The results indicate that the slight eccentricity of the Moon's orbit would result in excessive station-keeping requirements for long duration missions.
Haynes, Davy A.↗
Numerical Simulation of an Enclosed Laminar Jet Diffusion Flame in Microgravity Environment: Comparison with ELF Data
Enclosed diffusion flames are commonly found in practical combustion systems, such as the power-plant combustor, gas turbine combustor, and jet engine after-burner. In these systems, fuel is injected into a duct with a co-flowing or cross-flowing air stream. The diffusion flame is found at the surface where the fuel jet and oxygen meet, react, and consume each other. In combustors, this flame is anchored at the burner (i.e., fuel jet inlet) unless adverse conditions cause the flame to lift off or blow out. Investigations of burner stability study the lift off, reattachment, and blow out of the flame. Flame stability is strongly dependent on the fuel jet velocity. When the fuel jet velocity is sufficiently low, the diffusion flame anchors at the burner rim. When the fuel jet velocity is increased, the flame base gradually moves downstream. However, when the fuel jet velocity increases beyond a critical value, the flame base abruptly jumps downstream. When this "jump" occurs, the flame is said to have reached its lift-off condition and the critical fuel jet velocity is called the lift-off velocity. While lifted, the flame is not attached to the burner and it appears to float in mid-air. Flow conditions are such that the flame cannot be maintained at the burner rim despite the presence of both fuel and oxygen. When the fuel jet velocity is further increased, the flame will eventually extinguish at its blowout condition. In contrast, if the fuel jet velocity of a lifted flame is reduced, the flame base moves upstream and abruptly returns to anchor at the burner rim. The fuel jet velocity at reattachment can be much lower than that at lift off, illustrating the hysteresis effect present in flame stability. Although there have been numerous studies of flame stability, the controlling mechanisms are not well understood. This uncertainty is described by Pitts in his review of various competing theories of lift off and blow out in turbulent jet diffusion flames. There has been some research on the stability of laminar flames, but most studies have focused on turbulent flames. It is also well known that the airflow around the fuel jet can significantly alter the lift off, reattachment and blow out of the jet diffusion flame. Buoyant convection is sufficiently strong in 1-g flames that it can dominate the flow-field, even at the burner rim. In normal-gravity testing, it is very difficult to delineate the effects of the forced airflow from those of the buoyancy-induced flow. Comparison of normal-gravity and microgravity flames provides clear indication of the influence of forced and buoyant flows on the flame stability. The overall goal of the Enclosed Laminar Flames (ELF) investigation (STS-87/USMP-4 Space Shuttle mission, November to December 1997) is to improve our understanding of the effects of buoyant convection on the structure and stability of co-flow diffusion flame, e.g., see http://zeta.lerc.nasa.gov/expr/elf.htm. The ELF hardware meets the experiment hardware limit of the 35-liter interior volume of the glovebox working area, and the 180x220-mm dimensions of the main door. The ELF experiment module is a miniature, fan-driven wind tunnel, equipped with a gas supply system. A 1.5-mm diameter nozzle is located on the duct's flow axis. The cross section of the duct is nominally a 76-mm square with rounded corners. The forced air velocity can be varied from about 0.2 to 0.9 m/s. The fuel flow can be set as high as 3 std. cubic centimeter (cc) per second, which corresponds to a nozzle exit velocity of up to 1.70 m/s. The ELF hardware and experimental procedure are discussed in detail in Brooker et al. The 1-g test results are repeated in several experiments following the STS-87 Mission. The ELF study is also relevant to practical systems because the momentum-dominated behavior of turbulent flames can be achieved in laminar flames in microgravity. The specific objectives of this paper are to evaluate the use reduced model for simulation of flame lift-off and blowout.
Jia, Kezhong↗
The stability of thick, self-gravitating disks in protostellar systems
Using numerical hydrodynamics techniques, we perform a nonlinear stability analysis of accretion disk systems that contain thick, self-gravitating disks. The systems are initially represented by a point mass M(sub c) at the center and a geometrically thick, axisymmetric disk of mass M(sub d) that supports uniform specific angular momentum and obeys an n = 3/2, polytropic equation of state. The equilibrium disk structure is uniquely defined upon the specification of two key dimensionless system parameters: M(sub d)/M(sub c) and T/absolute value of W (the ratio of rotational kinetic energy if the disk to the gravitational potential energy of the system). The focus of this work is on the identification of systems within this two-dimensional parameter space that are marginally unstable toward the development of nonaxisymmetric distortions. The geometric form and relative pattern speed of the disk's distortion as well as the likelihood of disk fragmentation as a result of such instabilities is examined, particularly in the context of protostellar systems. The value of T/absolute value of W at which thick disks first become dynamically unstable to nonaxisymmetric distortions is found to vary significantly with the mass ratio of the accertion disk system. Nonaxisymmetric eigenmodes with four distinctly different characters are identified in systems with mass ratios in the range 0.2 less than or equal to (M(sub d)/M(sub c)) less than or equal to 5.
Woodward, John W.↗
A variational multiscale immersed meshfree method for heterogeneous materials
Abstract We introduce an immersed meshfree formulation for modeling heterogeneous materials with flexible non-body-fitted discretizations, approximations, and quadrature rules. The interfacial compatibility condition is imposed by a volumetric constraint, which avoids a tedious contour integral for complex material geometry. The proposed immersed approach is formulated under a variational multiscale based formulation, termed the variational multiscale immersed method (VMIM). Under this framework, the solution approximation on either the foreground or the background can be decoupled into coarse-scale and fine-scale in the variational equations, where the fine-scale approximation represents a correction to the residual of the coarse-scale equations. The resulting fine-scale solution leads to a residual-based stabilization in the VMIM discrete equations. The employment of reproducing kernel (RK) approximation for the coarse- and fine-scale variables allows arbitrary order of continuity in the approximation, which is particularly advantageous for modeling heterogeneous materials. The effectiveness of VMIM is demonstrated with several numerical examples, showing accuracy, stability, and discretization efficiency of the proposed method.
36 MATERIALS SCIENCE↗
Physics basis for the Wisconsin HTS Axisymmetric Mirror (WHAM)
The Wisconsin high-temperature superconductor axisymmetric mirror experiment (WHAM) will be a high-field platform for prototyping technologies, validating interchange stabilization techniques and benchmarking numerical code performance, enabling the next step up to reactor parameters. A detailed overview of the experimental apparatus and its various subsystems is presented. WHAM will use electron cyclotron heating to ionize and build a dense target plasma for neutral beam injection of fast ions, stabilized by edge-biased sheared flow. At 25 keV injection energies, charge exchange dominates over impact ionization and limits the effectiveness of neutral beam injection fuelling. This paper outlines an iterative technique for self-consistently predicting the neutral beam driven anisotropic ion distribution and its role in the finite beta equilibrium. Beginning with recent work by Egedal et al. ( Nucl. Fusion , vol. 62, no. 12, 2022, p. 126053) on the WHAM geometry, we detail how the FIDASIM code is used to model the charge exchange sources and sinks in the distribution function, and both are combined with an anisotropic magnetohydrodynamic equilibrium solver method to self-consistently reach an equilibrium. We compare this with recent results using the CQL3D code adapted for the mirror geometry, which includes the high-harmonic fast wave heating of fast ions.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Energy stability of thermocapillary convection in a model of the float-zone crystal-growth process
Energy stability theory has been applied to a basic state of thermocapillary convection occurring in a cylindrical half-zone of finite length to determine conditions under which the flow will be stable. Because of the finite length of the zone, the basic state must be determined numerically. Instead of obtaining stability criteria by solving the related Euler-Lagrange equations, the variational problem is attacked directly by discretization of the integrals in the energy identity using finite differences. Results of the analysis are values of the Marangoni number below which axisymmetric disturbances to the basic state will decay, for various values of the other parameters governing the problem.
Shen, Y.↗
Constrained model calibration of grain structure dependent spall dynamics in shock-loaded tantalum
Here, we perform a gas gun experiment by shock loading tantalum samples of varying grain structures to assess the suitability of a numerical model for simulating spall behavior. The observed differences in spall strength, as well spallation and re-compression history, are not captured in uncalibrated hydrodynamic simulations. An optimization is performed on the Johnson spall model to determine the best parameters that fit the observed trends. Linear stability analysis is employed to motivate bounds on those parameters. Herein, optimized simulations agree well with the experimental results, reproducing pullback depth and recompression timescales across the different samples tested. Further, the observed pullback time of the single crystal sample was found to imply, via the stability analysis, a percolation threshold in good agreement with the theoretical value for a body centered cubic lattice. Therefore, the combined linear stability and percolation analysis shows promise and may be applied to other materials with diverse microstructures. Collectively, the findings demonstrate that the model is suitable for reproducing spall-induced free surface behavior across various microstructures, but also points to caution in using model coefficients for uncalibrated microstructures.
58 GEOSCIENCES↗
Transition and laminar instability
The linear stability theory was applied to the problem of boundary layer transition in incompressible flow. The theory was put into a form suitable for three-dimensional boundary layers; both the temporal and spatial theories were examined; and a generalized Gaster relation for three-dimensional boundary layers was derived. Numerical examples include the stability characteristics of Falkner-Skan boundary layers, the accuracy of the two-dimensional Gaster relation for these boundary layers, and the magnitude and direction of the group velocity for oblique waves in the Blasius boundary layer. Available experiments which bear on the validity of stability theory and its relation to transition are reviewed and the stability theory is applied to transition prediction. The amplitude method is described in which the wide band disturbance amplitude in the boundary layer is estimated from stability theory and an interaction relation for the initial amplitude density of the most unstable frequency.
Mack, L. M.↗
Structure preserving transport stabilized compatible finite element methods for magnetohydrodynamics
Here, we present compatible finite element space discretizations for the ideal compressible magnetohydrodynamic equations. The magnetic field is considered both in div- and curl-conforming spaces, leading to a strongly or weakly preserved zero-divergence condition, respectively. The equations are discretized in space such that transfers between the kinetic, internal, and magnetic energies are consistent, leading to a preserved total energy. We also discuss further adjustments to the discretization required to additionally achieve magnetic helicity preservation. Finally, we describe new transport stabilization methods for the magnetic field equation which maintain the zero-divergence and energy conservation properties, including one method which also preserves magnetic helicity. The methods' preservation and improved stability properties are confirmed numerically using a steady state and a magnetic dynamo test case.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Analytic stability boundaries for compressional and global Alfvén eigenmodes driven by fast ions. II. Interaction via Landau resonance
We report that conditions for net fast ion drive are derived for beam-driven, co-propagating, sub-cyclotron compressional (CAE) and global (GAE) Alfvén eigenmodes driven by the Landau resonance with super-Alfvénic fast ions. Approximations applicable to realistic neutral beam distributions and mode characteristics observed in spherical tokamaks enable the derivation of marginal stability conditions for these modes. Such conditions successfully reproduce the stability boundaries found from numerical integration of the exact expression for local fast ion drive/damping. Coupling between the CAE and GAE branches of the dispersion due to finite ω/ω ci and |k ∥ /k ⊥ | is retained and found to be responsible for the existence of the GAE instability via this resonance. Finally, encouraging agreement is demonstrated between the approximate stability criterion, simulation results, and a database of NSTX observations of co-CAEs.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Regularizing the linearly extrapolated BDF2 scheme for incompressible flows with time relaxation
This paper presents a highly-efficient finite element scheme for the time relaxation model (TRM). The efficiency is achieved through the second-order BDF2 time-stepping scheme with linear extrapolation (BDF2LE). The accuracy of the scheme is also greatly enhanced through the use of the divergence-free Scott-Vogeulis finite elements, and van Cittert approximate deconvolution. A complete finite element analysis is provided, which includes rigorous proofs for the stability, well-possessedness, and convergence of both velocity and pressure solutions. Furthermore, we also demonstrate that the inclusion of the linear time relaxation term preserves the long-time stability of the unregularized BDF2LE scheme. Finally, numerical experiments are presented that demonstrate the added stability and accuracy that time relaxation can provide.
97 MATHEMATICS AND COMPUTING↗
High Order Finite Difference Methods for Multiscale Complex Compressible Flows
The classical way of analyzing finite difference schemes for hyperbolic problems is to investigate as many as possible of the following points: (1) Linear stability for constant coefficients; (2) Linear stability for variable coefficients; (3) Non-linear stability; and (4) Stability at discontinuities. We will build a new numerical method, which satisfies all types of stability, by dealing with each of the points above step by step.
Sjoegreen, Bjoern↗
Linear Wave Motion from Concentrated Harmonic Sources in Blasius Flow
The motion of individual linear instability waves in shear flows is well described by existing theoretical and numerical methods. However, naturally occuring sources produce coherent wave motions with broadband spanwise wavenumber and frequency spectra, and the different spectral components interact both linearly and nonlinearily. This paper describes a series of calculations for the parameters of three different experiments using locally parallel linear stability theory (LST), the Parabolized Stability Equations (PSE), and Direct Numerical Simulation (DNS). The calculations illustrate the strengths and weaknesses of the different methods, the extent to which the methods agree or disagree, and, finally the extent to which agreement with the measurements can be attained, given that the experiments also have their own difficulties.
linear↗