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 343 records · Page 19

Fundamental Aerodynamic Investigations for Development of Arrow-Stabilized Projectiles

The numerous patent applications on arrow-stabilized projectiles indicate that the idea of projectiles without spin is not new, but has appeared in various proposals throughout the last decades. As far as projectiles for subsonic speeds are concerned, suitable shapes have been developed for sometime, for example, numerous grenades. Most of the patent applications, though, are not practicable particularly for projectiles with supersonic speed. This is because the inventor usually does not have any knowledge of aerodynamic flow around the projectile nor any particular understanding of the practical solution. The lack of wind tunnels for the development of projectiles made it necessary to use firing tests for development. These are obviously extremely tedious or expensive and lead almost always to failures. The often expressed opinion that arrow-stabilized projectiles cannot fly supersonically can be traced to this condition. That this is not the case has been shown for the first time by Roechling on long projectiles with foldable fins. Since no aerodynamic investigations were made for the development of these projectiles, only tedious series of firing tests with systematic variation of the fins could lead to satisfactory results. These particular projectiles though have a disadvantage which lies in the nature cf foldable fins. They occasionally do not open uniformly in flight, thus causing unsymmetry in flow and greater scatter. The junctions of fins and body are very bad aerodynamically and increase the drag. It must be possible to develop high-performance arrow-stabilized projectiles based on the aerodynamic research conducted during the last few years at Peenemuende and new construction ideas. Thus the final shape, ready for operational use, could be developed in the wind tunnel without loss of expensive time in firing tests. The principle of arrow-stabilized performance has been applied to a large number of caliburs which were stabilized by various means Most promising was the development of a subcaliber wing-stabilized projectile with driving disc (Treibspiegel) where rigid control surfaces extend beyond the caliber of the projectile into the free stream. The stabilized projectiles of full-caliber, wing-stabilized projectiles with fins within the caliber is considerably more difficult. A completely satisfactory solution for the latter has not been found yet.

Kurzweg, Hermann↗

Meso-modeling of Carbon Fiber Composite for Crash Safety Analysis

In the conventional approach, the material properties for crash safety simulations are typically obtained from standard coupon tests, where the test results only provide single layer material properties used in crash simulations. However, the lay-up effects for the failure behaviors of the real structure were not considered in numerical simulations. Hence, there was discrepancy between the crash simulations and experimental tests. Consequently, an intermediate stage is required for accurate predictions. Some component tests are required to correlate the material models in the intermediate stage. In this paper, a Mazda Tube under high-impact velocity is chosen as an example for the crash safety analysis. The tube consists of 24 layers of uni-directional (UD) carbon fiber composite materials, in which 4 layers are perpendicular to, while the other layers are parallel to the impact direction. An LS-DYNA meso-model was constructed with orthotropic material models counting for the single-layer material behaviors. Between layers, a node-based tie-break contact was used for modeling the delamination of the composite material. Since fiber directions are not single-oriented, the lay-up effects could be an important effect. From the first numerical trial, premature material failure occurred due to the use of material parameters obtained directly from the coupon tests. Some parametric studies were conducted to identify the cause of the numerical instability. The finding is that the material failure strength used in the numerical model needs to be enlarged to stabilize the numerical model. Some hypothesis was made to provide the foundation for enlarging the failure strength and the corresponding experiments will be conducted to validate the hypothesis.

Lin, Shih-Po↗

Validation study of RWM stability in DIII-D high- β N plasmas

The n = 1 (n is the toroidal mode number) resistive wall mode (RWM) stability is numerically investigated for two DIII-D high-β N discharges 176440 and 172461, utilizing the MARS-F (Liu et al 2000 Phys. Plasmas 7 3681) and MARS-K (Liu et al 2008 Phys. Plasmas 15 112503) codes. Systematic validation efforts are attempted, for the first time, for discharges with very slow or vanishing toroidal flow for a large fraction of the plasma volume. While gaining physics insights in accessing stable operation regime at β N exceeding the Troyon no-wall limit in these slow-rotation experiments, the predictive capability of fluid and non-perturbative magnetohydrodynamic-kinetic hybrid models for the RWM is further confirmed. The MARS-F fluid model, with a strong but numerically tunable viscosity mimicking ion Landau damping of parallel sound waves, finds complete stabilization of the n = 1 RWM in the considered DIII-D plasmas under the experimental flow conditions. Similarly, either full stabilization (for discharge 176440) or marginal stability (for discharge 172461) of the mode is computed by the MARS-K hybrid model, which is first-principle based without free model parameters. In particular, all drift kinetic resonances, including those of thermal and energetic particles, are found to synergistically act to marginally stabilize the RWM in discharge 172461. These MARS-F/K modeling results explain the experimentally observed stable operational regime in DIII-D, as far as the RWM stability is concerned. Extensive numerical sensitivity studies, with respect to the plasma toroidal flow speed as well as the radial location of the resistive wall, are also carried out to further support the validation study.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A new concept of stability in orbit propagation, useful for quantifying numerical errors

We present the concept of topological stability in the numerical propagation of orbits, and show how it results in a useful new method for measuring the global numerical error of an orbit propagation. The concept applies to any problem in orbital dynamics. Moreover, it can be extended to any three-dimensional system of di erential equations of second order. In order to assess the topological stability of a given integration a special metric is introduced, which can be used to estimate the numerical errors robustly. The method is particularly well suited for dealing with strongly perturbed and chaotic systems. The construction is based on the constraint imposed by the Hopf map that supports the Kustaanheimo-Stiefel transformation. Generic concepts of stability are translated to KS space.

Pelaez, Jesus↗

Numerical experiments in climate stability.

Two semiempirical climate models are formulated as time-dependent problems in order to study the stability of their asymptotic steady state equilibrium solutions to perturbations in internal (in this case, initial) conditions. For fixed external conditions the time-dependent versions of these models exhibit fully transitive behavior to positive perturbations in initial conditions, slight intransivity to negative perturbations up to -18 K, and an ice-covered earth regime is obtained for extremely large negative perturbations in initial temperatures (below -18 K). The parametrization found to be most critical in these models is the albedo-temperature coupling, especially in tropical regions. The temperature albedo parametrization from an intransitive climate model used in a time-dependent form of one of the preceding models is found to be highly sensitive to negative perturbations in both internal and external conditions. Numerical experimentation with these semiempirical models shows the important role of the tropics in maintaining the equilibrium climate and suggests that the radiation balance in equatorial latitudes might have a controlling influence on the equatorward extent of the polar ice cap.

Schneider, S. H.↗

Solution of the minimum time-to-climb problem by matched asymptotic expansions

Application of singular perturbation techniques to trajectory optimization problems of flight mechanics is discussed. The method of matched asymptotic expansions is used to obtain an approximate solution to the aircraft minimum time-to-climb problem. Outer, boundary-layer, and composite solutions are obtained to zeroth and first orders. A stability criterion is derived for the zeroth-order boundary-layer solutions (the theory requires a form of boundary-layer stability). A numerical example is considered for which it is shown that the stability criterion is satisfied and a useful numerical solution is obtained. The zeroth-order solution proves to be a poor approximation, but the first-order solution gives a good approximation for both the trajectory and the minimum time-to-climb. The computational cost of the singular perturbation solution is considerably less than that of a steepest descent solution. Thus singular perturbation methods appear to be promising for the solution of optimal control problems.

Ardema, M. D.↗

Experimental and numerical investigation of flame stabilization and pollutant formation in matrix stabilized ammonia-hydrogen combustion

Ammonia (NH 3 ) is a carbon-free fuel that offers an attractive alternative for reducing greenhouse gas emissions. However, the slow flame speed, low heating value, and emissions of nitrogen-containing pollutants present significant issues for practical combustion applications. To address these issues, we investigate the use of matrix stabilized combustion. In this type of burner, combustion is performed within an inert porous ceramic foam, heat is recirculated by solid conduction and radiation, which enhances flame speed and combustion stabilization, thereby permitting combustion over a wide range of equivalence ratio conditions. We present a new porous media burner (PMB) capable of stabilizing NH 3 /air flames at ambient conditions. An extensive experimental characterization of the stability of this burner is conducted with up to 30% by volume of hydrogen (H 2 ) in the fuel stream. A 15:1 turndown ratio is demonstrated, with a high thermal power density of 62 MW m -3 . Concentrations of NO, unburnt NH 3 , and H 2 in the exhaust stream are measured. Two regimes are identified for low NO operation: rich and very lean. For rich conditions, NO emissions decrease with increasing equivalence ratio and decreasing H 2 blending. Unburnt NH 3 emissions follow opposite trends. These measurements are complemented by simulations in which the burner is represented by a coupled solid-gas reactor network. In conclusion, this model captures the burner’s pollutant emissions to good accuracy and is used to analyze the mechanisms of pollutant formation.

42 ENGINEERING↗

Numerical modeling of pedestal stability and broadband turbulence of wide-pedestal QH-mode plasmas on DIII-D

The wide-pedestal quiescent high confinement mode discovered on DIII-D in recent years is a stationary and quiescent H-mode (QH-mode) with the pedestal width exceeding EPED prediction by at least 25%. Its characteristics, such as low rotation, high energy confinement and edge localized mode-free operation, make it an attractive operation mode for future reactors. Linear and nonlinear simulations using BOUT++ reduced two fluid MHD models and awere carried out to investigate the bursty broadband turbulence often observed in the edge of wide-pedestal QH-mode plasmas. Two kinds of MHD-scale instabilities in different spatial locations within the pedestal were found in the simulations: one mild peeling–ballooning (PB) mode γ PB < 0.04ω A ) located near the minimum in E r</:sub> well propagating in ion diamagnetic drift direction; and one drift-Alfvén wave locates at smaller radius compared to E r</:sub> well propagating in the electron diamagnetic drift direction and unstable only when the parallel electron dynamics is included in the simulation. The coupling between drift wave and shear Alfvén wave provides a possible cause of the experimentally observed local profile flattening in the upper-pedestal. The rotation direction, mode location, as well as the wavenumber of these two modes from BOUT++ simulations agree reasonably well with the experimental measurements, while the lack of quantitative agreement is likely due to the lack of trapped electron physics in current fluid model. This work presents improved physics understanding of the pedestal stability and turbulence dynamics for wide-pedestal QH-mode.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Evaluation of Thin Plate Hydrodynamic Stability through a Combined Numerical Modeling and Experimental Effort

An experimental and computational effort was undertaken in order to evaluate the capability of the fluid-structure interaction (FSI) simulation tools to describe the deflection of a Missouri University Research Reactor (MURR) fuel element plate redesigned for conversion to lowenriched uranium (LEU) fuel due to hydrodynamic forces. Experiments involving both flat plates and curved plates were conducted in a water flow test loop located at the University of Missouri (MU), at conditions and geometries that can be related to the MURR LEU fuel element. A wider channel gap on one side of the test plate, and a narrower on the other represent the differences that could be encountered in a MURR element due to allowed fabrication variability. The difference in the channel gaps leads to a pressure differential across the plate, leading to plate deflection. The induced plate deflection the pressure difference induces in the plate was measured at specified locations using a laser measurement technique. High fidelity 3-D simulations of the experiments were performed at MU using the computational fluid dynamics code STAR-CCM+ coupled with the structural mechanics code ABAQUS. Independent simulations of the experiments were performed at Argonne National Laboratory (ANL) using the STAR-CCM+ code and its built-in structural mechanics solver. The simulation results obtained at MU and ANL were compared with the corresponding measured plate deflections.

97 MATHEMATICS AND COMPUTING↗

A numerical study of capillary stability in a circular cylindrical container with a concave spheroidal bottom

The paper describes a theoretical investigation of the stability under gravitational and surface forces of a liquid in a circular cylindrical container with a concave spheroidal bottom for the case in which the volume of the liquid is sufficiently small so that the bottom is not covered completely. The gravitational field is assumed to be directed along the symmetry axis of the container, and for a specific container shape the critical Bond number is calculated as a function of liquid volume for contact angles of gamma = 0, 1, 2, and 4 deg. For gamma = 0 deg, some critical equilibrium configurations and corresponding perturbation modes are presented.

Concus, P.↗

Numerical experiments on the stability of controlled shear flows

Three-dimensional, nonlinear simulations are presented for the parallel water boundary layer subjected to pressure gradient, suction, or heating control. The focus is on the secondary instability responsible for the K-type transition. Heating control is found to be the least effective in delaying the secondary instability. Flow-field details (including temperature profiles) are presented for both the uncontrolled boundary layer and the heated boundary layer.

Zang, T. A.↗

A numerical study of flare stabilized projectiles at Mach 7.0

The effect of nose shape on two flare stabilized projectiles was studied using a Parabolized Navier-Stokes code. Pressure coefficients, forces and moments, skin friction coefficients, and Stanton number calculations are presented for the hemisphere-cylinder-flare and the cone-cylinder-flare configurations. Pitching moment and static margin plots versus flare angles and flare lengths are presented in a parametric study to show aerodynamic stability effects. Distinctly different flow field values of pressure, local Mach number, and dynamic pressure were generated by the two different nose configurations. These flow field values just upstream of the flare will be examined. The code demonstrated its value as a design tool by making a clear distinction between aerodynamically stable characteristics for this variety of nose shapes, flare angles, and flare lengths.

Hughson, Montgomery C.↗

CONDIF - A modified central-difference scheme with unconditional stability and very low numerical diffusion

This paper describes a method, called CONDIF, which retains the essential nature of the central difference scheme (CDS) but eliminates the over- and under-shoots. It modifies the CDS by introducing a controlled amount of numerical diffusion based on local gradients. The results obtained for a number of test problems show that the CONDIF is significantly more accurate than the hybrid scheme when the Peclet number is very high and the flow is at large angles to the grid. In the worst case, it yields results similar to those of the hybrid scheme. Though the CONDIF is unconditionally stable, under some conditions the rate of convergence deteriorates as the grid Peclet number is increased. Efforts are currently underway to develop a faster converging variation of CONDIF.

Runchal, Akshai K.↗

Numerical experiments on the stability of controlled boundary layers

Nonlinear simulations are presented for instability and transition in parallel water boundary layers subjected to pressure gradient, suction, or heating control. In the nonlinear regime, finite amplitude, 2-D Tollmein-Schlichting waves grow faster than is predicted by linear theory. Moreover, this discrepancy is greatest in the case of heating control. Likewise, heating control is found to be the least effective in delaying secondary instabilities of both the fundamental and subharmonic type. Flow field details (including temperature profiles) are presented for both the uncontrolled boundary layer and the heated boundary layer.

Zang, Thomas A.↗

A numerical study of the nonlinear thermal stability of solar loops

A time-dependent numerical model is used to investigate the nonlinear thermal stability of static loops of various heights. Simulations show that the instability of a hot state with loop heights of less than about 1000 km is physically significant, with an initially hot atmosphere in low-lying compact loops evolving to an extended atmosphere with temperatures far below 100,000 K. Results also show that high-lying loops are stable to all reasonable perturbations, including those of large initial amplitude and long wavelength. The simulation results suggest that low-lying compact loops should not be common to the sun, and that cool loops with temperatures near 100,000 K must be formed in the cool state initially and cannot evolve from preexisiting loops.

Klimchuk, J. A.↗

Thermal stability during rotation in space - A scaling and numerical analysis

An analysis of thermal stability in a fluid during rotation is presented. The conditions resulting in thermal stability are examined and numerically verified. It is shown that at least one thermal field exists for any rotating fluid, with or without a gravitational background field, where convection does not occur. The numerical model used is fully nonlinear and 3D. Results are presented on the basis of both steady-state and fully transient models. The thermally stable condition is further investigated in the context of the average net acceleration and the inhomogeneous part of the acceleration field (gradient of acceleration). It is shown that in the thermally stable condition the flow that would result from the average net acceleration is equal and opposite to the flow that would result from the inhomogeneous part of the acceleration field. It is concluded that the thermally stable condition is actually a flow mode transition at Ad = 1.

Arnold, William A.↗

Numerical solution of wave equations for the stability of the inner cometo-sheath

Numerical solution of the MHD wave equations for stability of the cometary sheath determined by the balance between the inward Lorentz body force and the outward ion-neutral drag force is obtained by using a two-point boundary value method. The eigenvalues and the eigenfunctions are obtained numerically by treating the cometary inner sheath as a layer of finite thickness, bounded by the contact surface, i.e., the diamagnetic cavity boundary. The magnetic field structure discovered in the ionosphere of Comets Halley and Giacobini-Zinner is found to be unstable. The effects of finite plasma pressure, dissociative recombination, and mass loading due to photoionization are found to be stabilizing but are unable to quench the instability completely. It is also found that the higher the neutral production rate the lesser is the growth rate for the instability.

Srivastava, Krishna M.↗