Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “velocity steps”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

LS-DYNA User-Defined Internal Ballistic Modeling

- LS-DYNA Explicit Finite Element Analysis software and CADPROG Internal Ballistic Analysis code are integrated as a single modeling tool to analyze pyro-mechanical devices. - Combining these two codes allows the ballistic calculations to be performed without modeling the propellant in an FEA mesh. At each LS-DYNA time step, CADPROG is analytically run using the kinematic data (stroke, velocity, chamber volume) fed back from the LS-DYNA simulation. Calculated pressure is then applied to a piston or any applicable surface mesh in LS-DYNA model interactively. -All external loads such as friction, damping loads, assisting/resisting loads, locking/unlocking loads are modeled in LS-DYNA and coupled with pressure loads from CADPROG to create a fully defined equation of motion -Both codes are written in Fortran 77. CADPROG is inserted in LS-DYNA source code as a user-defined subroutine and re-compiled to generate a custom executable.

pyrotechnics↗
NASA NTRSDate not supplied

LS-DYNA User-Defined Internal Ballistic Modeling

LS-DYNA Explicit Finite Element Analysis software and CADPROG Internal Ballistic Analysis code are integrated as a single modeling tool to analyze pyro-mechanical devices. Coupling these two codes allows the ballistic calculations to be performed without modeling the propellant in an FEA mesh. At each LS-DYNA time step, CADPROG is run analytically using the kinematic data (stroke, velocity, chamber volume) fed back from the LS-DYNA simulation. The calculated pressure is then applied to a piston or any applicable surface mesh in LS-DYNA model interactively. All external loads such as friction, damping loads, assisting/resisting loads, locking/unlocking loads are modeled in LS-DYNA and coupled with pressure loads from CADPROG to create a fully defined equation of motion Both codes are written in Fortran 77. CADPROG is integrated in the LS-DYNA source code as a user-defined subroutine and recompiled to generate a custom executable.

pyrotechnics↗
NASA NTRS1991

Analysis of two-equation turbulence models for recirculating flows

The two-equation kappa-epsilon model is used to analyze turbulent separated flow past a backward-facing step. It is shown that if the model constraints are modified to be consistent with the accepted energy decay rate for isotropic turbulence, the dominant features of the flow field, namely the size of the separation bubble and the streamwise component of the mean velocity, can be accurately predicted. In addition, except in the vicinity of the step, very good predictions for the turbulent shear stress, the wall pressure, and the wall shear stress are obtained. The model is also shown to provide good predictions for the turbulence intensity in the region downstream of the reattachment point. Estimated long time growth rates for the turbulent kinetic energy and dissipation rate of homogeneous shear flow are utilized to develop an optimal set of constants for the two equation kappa-epsilon model. The physical implications of the model performance are also discussed.

Thangam, S.↗
NASA NTRS1993

Second quantization in bit-string physics

Using a new fundamental theory based on bit-strings, a finite and discrete version of the solutions of the free one particle Dirac equation as segmented trajectories with steps of length h/mc along the forward and backward light cones executed at velocity +/- c are derived. Interpreting the statistical fluctuations which cause the bends in these segmented trajectories as emission and absorption of radiation, these solutions are analogous to a fermion propagator in a second quantized theory. This allows us to interpret the mass parameter in the step length as the physical mass of the free particle. The radiation in interaction with it has the usual harmonic oscillator structure of a second quantized theory. How these free particle masses can be generated gravitationally using the combinatorial hierarchy sequence (3,10,137,2(sup 127) + 136), and some of the predictive consequences are sketched.

Noyes, H. Pierre↗
NASA NTRS2005

Growth from Solutions: Kink dynamics, Stoichiometry, Face Kinetics and stability in turbulent flow

1. Kink dynamics. The first segment of a polygomized dislocation spiral step measured by AFM demonstrates up to 60% scattering in the critical length l*- the length when the segment starts to propagate. On orthorhombic lysozyme, this length is shorter than that the observed interkink distance. Step energy from the critical segment length based on the Gibbs-Thomson law (GTL), l* = 20(omega)alpha/(Delta)mu is several times larger than the energy from 2D nucleation rate. Here o is tine building block specific voiume, a is the step riser specific free energy, Delta(mu) is the crystallization driving force. These new data support our earlier assumption that the classical Frenkel, Burton -Cabrera-Frank concept of the abundant kink supply by fluctuations is not applicable for strongly polygonized steps. Step rate measurements on brushite confirms that statement. This is the1D nucleation of kinks that control step propagation. The GTL is valid only if l* <Dk/vk, the diffusion path of a kink that has diffusivity Dk and average growth velocity vk. This is equivalent to supersaturations sigma less than approx. alpha/2l*, where alpha is the building block size. For lysozyme, sigma much less than (1%). Conventionally used interstep distance generated by screw dislocation, 19(omega)alpha/Delta(mu) should be replaced by the very different real one, approx.4l*. 2. Stoichiometry. Kink, and thus step and face rates of a non-Kossel complex molecular monocomponent or any binary, AB, lattice was found theoretically to be proportional to 1/(zeta(sup 1/2) + zeta(sup - 1/2)), where zeta = [B]/[A] is the stoichiometry ratio in solution. The velocities reach maxima at zeta = 1. AFM studies of step rates on CaOxalate monohydrate (kidney stones) from aqueous solution was found to obey the law mentioned above. Generalization for more complex lattice will be discussed. 3. Turbulence. In agreement with theory, high precision in-situ laser interferometry of the (101) KDP crystal face shows step bunching if solution flows parallel to the step flow. The bunch height increases with the distance the bunch travels, i.e. with the face size. However, when the flow rate, u, increases, at u greater than approx. 1 m / s , the average step bunch height decreases as 1/u. The pheonomenon is attributed to the turbulent rather than laminar viscous boundary layer where diffusivity Dt = 0.5u(sub tau),y, i.e. increases linearly with the distance y from the solid face. Friction velocity, u(sub tau) approx. u(sup 7/8). Dramatically larger rate of the mass/heat transport within the turbulent, as compared to the laminar, viscous layer will be discussed.

Chernov, A. A.↗
NASA NTRS2001

Horizontal vestibuloocular reflex evoked by high-acceleration rotations in the squirrel monkey. IV. Responses after spectacle-induced adaptation

The horizontal angular vestibuloocular reflex (VOR) evoked by sinusoidal rotations from 0.5 to 15 Hz and acceleration steps up to 3,000 degrees /s(2) to 150 degrees /s was studied in six squirrel monkeys following adaptation with x2.2 magnifying and x0.45 minimizing spectacles. For sinusoidal rotations with peak velocities of 20 degrees /s, there were significant changes in gain at all frequencies; however, the greatest gain changes occurred at the lower frequencies. The frequency- and velocity-dependent gain enhancement seen in normal monkeys was accentuated following adaptation to magnifying spectacles and diminished with adaptation to minimizing spectacles. A differential increase in gain for the steps of acceleration was noted after adaptation to the magnifying spectacles. The gain during the acceleration portion, G(A), of a step of acceleration (3,000 degrees /s(2) to 150 degrees /s) increased from preadaptation values of 1.05 +/- 0.08 to 1.96 +/- 0.16, while the gain during the velocity plateau, G(V), only increased from 0.93 +/- 0.04 to 1.36 +/- 0.08. Polynomial fits to the trajectory of the response during the acceleration step revealed a greater increase in the cubic than the linear term following adaptation with the magnifying lenses. Following adaptation to the minimizing lenses, the value of G(A) decreased to 0.61 +/- 0.08, and the value of G(V) decreased to 0.59 +/- 0.09 for the 3,000 degrees /s(2) steps of acceleration. Polynomial fits to the trajectory of the response during the acceleration step revealed that there was a significantly greater reduction in the cubic term than in the linear term following adaptation with the minimizing lenses. These findings indicate that there is greater modification of the nonlinear as compared with the linear component of the VOR with spectacle-induced adaptation. In addition, the latency to the onset of the adapted response varied with the dynamics of the stimulus. The findings were modeled with a bilateral model of the VOR containing linear and nonlinear pathways that describe the normal behavior and adaptive processes. Adaptation for the linear pathway is described by a transfer function that shows the dependence of adaptation on the frequency of the head movement. The adaptive process for the nonlinear pathway is a gain enhancement element that provides for the accentuated gain with rising head velocity and the increased cubic component of the responses to steps of acceleration. While this model is substantially different from earlier models of VOR adaptation, it accounts for the data in the present experiments and also predicts the findings observed in the earlier studies.

NASA Discipline Neuroscience↗
DOE OSTI2020

Numerical integrators for Lagrangian oceanography

Abstract. A common task in Lagrangian oceanography is to calculate a large number of drifter trajectories from a velocity field precalculated with an ocean model. Mathematically, this is simply numerical integration of an ordinary differential equation (ODE), for which a wide range of different methods exist. However, the discrete nature of the modelled ocean currents requires interpolation of the velocity field in both space and time, and the choice of interpolation scheme has implications for the accuracy and efficiency of the different numerical ODE methods. We investigate trajectory calculation in modelled ocean currents with 800 m, 4 km, and 20 km horizontal resolution, in combination with linear, cubic and quintic spline interpolation. We use fixed-step Runge–Kutta integrators of orders 1–4, as well as three variable-step Runge–Kutta methods (Bogacki–Shampine 3(2), Dormand–Prince 5(4) and 8(7)). Additionally, we design and test modified special-purpose variants of the three variable-step integrators, which are better able to handle discontinuous derivatives in an interpolated velocity field. Our results show that the optimal choice of ODE integrator depends on the resolution of the ocean model, the degree of interpolation, and the desired accuracy. For cubic interpolation, the commonly used Dormand–Prince 5(4) is rarely the most efficient choice. We find that in many cases, our special-purpose integrators can improve accuracy by many orders of magnitude over their standard counterparts, with no increase in computational effort. Equivalently, the special-purpose integrators can provide the same accuracy as standard methods at a reduced computational cost. The best results are seen for coarser resolutions (4 and 20 km), thus the special-purpose integrators are particularly advantageous for research using regional to global ocean models to compute large numbers of trajectories. Our results are also applicable to trajectory computations on data from atmospheric models.

58 GEOSCIENCES↗
NASA NTRS2023

Study of Pairwise Deconfliction Metrics to Analyze Air Traffic Complexity in Upper Class E Airspace

Upper Class E Traffic Management (ETM) is envisioned to cooperatively facilitate operations of a diverse set of aerial vehicles, such as high-altitude long-endurance fixed-wing unmanned aircraft (low-speed and high-speed), high-altitude platforms, airships, stratospheric balloons, supersonic unmanned and commercial aircraft, etc., with a wide variety of mission types, performance characteristics, communication, navigation and surveillance capabilities, maneuverability, and on-board avionics in the National Airspace System (NAS) ’above’ 60,000 feet above mean sea level, without an active and direct control from human air traffic controllers. A diverse mixture of aerial vehicle types creates significant challenges in understanding air traffic complexity, which may not correlate strongly with air traffic density. One key step for determining air traffic complexity in upper class E airspace is to first understand pairwise deconfliction metrics such as reachability, reserve area, and reserve flight time for each pair of unique aerial vehicle types under potential conflict. Therefore, pairwise deconfliction metrics are first defined, and analytical equations are derived for conflict resolution using the heading change maneuver. Next, case studies are performed to analyze deconfliction metrics to avoid secondary conflicts in upper class E airspace. The study shows that pairwise deconfliction metrics are functions of maneuverability, performance characteristics, uncertainty in position and velocity, heading angle change, and conflict angle of aerial vehicles. The next step for this research is to build a mathematical model for air traffic complexity using pairwise deconfliction metrics and validate it in an upper Class E simulation environment.

Airspace Complexity↗
NASA NTRSDate not supplied

Study of Pairwise Deconfliction Metrics to Analyze Air Traffic Complexity in Upper Class E Airspace

Upper Class E Traffic Management (ETM) is envisioned to cooperatively facilitate operations of a diverse set of aerial vehicles, such as high-altitude long-endurance fixed-wing unmanned aircraft (low-speed and high-speed), high-altitude platforms, airships, stratospheric balloons, supersonic unmanned and commercial aircraft, etc., with a wide variety of mission types, performance characteristics, communication, navigation and surveillance capabilities, maneuverability, and on-board avionics in the National Airspace System (NAS) ’above’ 60,000 feet above mean sea level, without an active and direct control from human air traffic controllers. A diverse mixture of aerial vehicle types creates significant challenges in understanding air traffic complexity, which may not correlate strongly with air traffic density. One key step for determining air traffic complexity in upper class E airspace is to first understand pairwise deconfliction metrics such as reachability, reserve area, and reserve flight time for each pair of unique aerial vehicle types under potential conflict. Therefore, pairwise deconfliction metrics are first defined, and analytical equations are derived for conflict resolution using the heading change maneuver. Next, case studies are performed to analyze deconfliction metrics to avoid secondary conflicts in upper class E airspace. The study shows that pairwise deconfliction metrics are functions of maneuverability, performance characteristics, uncertainty in position and velocity, heading angle change, and conflict angle of aerial vehicles. The next step for this research is to build a mathematical model for air traffic complexity using pairwise deconfliction metrics and validate it in an upper Class E simulation environment.

Airspace Complexity↗
NASA NTRS2003

Study of Unsteady Flows with Concave Wall Effect

This paper presents computational fluid dynamic studies of the inlet turbulence and wall curvature effects on the flow steadiness at near wall surface locations in boundary layer flows. The time-stepping RANS numerical solver of the NASA Glenn-HT RANS code and a one-equation turbulence model, with a uniform inlet turbulence modeling level of the order of 10 percent of molecular viscosity, were used to perform the numerical computations. The approach was first calibrated for its predictabilities of friction factor, velocity, and temperature at near surface locations within a transitional boundary layer over concave wall. The approach was then used to predict the velocity and friction factor variations in a boundary layer recovering from concave curvature. As time iteration proceeded in the computations, the computed friction factors converged to their values from existing experiments. The computed friction factors, velocity, and static temperatures at near wall surface locations oscillated periodically in terms of time iteration steps and physical locations along the span-wise direction. At the upstream stations, the relationship among the normal and tangential velocities showed vortices effects on the velocity variations. Coherent vortices effect on the velocity components broke down at downstream stations. The computations also predicted the vortices effects on the velocity variations within a boundary layer flow developed along a concave wall surface with a downstream recovery flat wall surface. It was concluded that the computational approach might have the potential to analyze the flow steadiness in a turbine blade flow.

Wang, Chi R.↗
NASA NTRS1972

Strongly coupled stress waves in heterogeneous plates.

Consideration of coupled stress waves generated by an impulsive load applied at one end of a semiinfinite plate. For the field equations governing the one-dimensional coupled waves a hyperbolic system of equations is obtained in which a strong coupling in the second derivatives exists. The method of characteristics described by Chou and Mortimer (1967) is extended to cover the case of strong coupling, and a study is made of the transient stress waves in a semiinfinite plate subjected to an initial step input. Coupled discontinuity fronts are found to propagate at different velocities. The normal plate stress and the bending moment at different time regimes are illustrated by graphs.

Wang, A. S. D.↗
NASA NTRS1973

Precise calculation of the magnetosphere surface for a tilted dipole.

The shape of the magnetosphere has been calculated self-consistently for inclinations of the earth's magnetic dipole from perpendicular to the solar wind. Inclination angles of 0-35 deg have been chosen in steps of 5 deg and various smooth trends in the surface characteristics with increasing inclination angle noted. The surface points and the complete field at the surface points have been calculated for the entire surfaces. The neutral point region has been given precise study in one degree steps and is found to be tangent to the solar wind velocity and to have a smooth continuous curvature.

Choe, J. Y.↗
NASA NTRS1977

Evaluation of subgrid-scale turbulence models using a fully simulated turbulent flow

An exact turbulent flow field was calculated on a three-dimensional grid with 64 points on a side. The flow simulates grid-generated turbulence from wind tunnel experiments. In this simulation, the grid spacing is small enough to include essentially all of the viscous energy dissipation, and the box is large enough to contain the largest eddy in the flow. The method is limited to low-turbulence Reynolds numbers, in our case R sub lambda = 36.6. To complete the calculation using a reasonable amount of computer time with reasonable accuracy, a third-order time-integration scheme was developed which runs at about the same speed as a simple first-order scheme. It obtains this accuracy by saving the velocity field and its first-time derivative at each time step. Fourth-order accurate space-differencing is used.

Clark, R. A.↗
NASA NTRS1993

The turbulent thermal boundary layer with an abrupt change from a rough to a smooth wall

The work reported here was motivated by concern over the use of smooth heat flux gages for heat transfer measurements on the otherwise rough turbine blades. Stanton number distributions and boundary layer profiles of mean temperature, mean velocity, and turbulence intensity are reported for a surface with a step change from a rough to a smooth surface. In most cases, the Stanton number immediately downstream of the change in roughness drops below the all-smooth-wall data at the same x-Reynolds number. The alignment of the smooth surface between the bases and crests of the roughness elements is shown to have only a weak effect on the Stanton number distribution. It is concluded that the use of smooth heat flux gages on otherwise rough surfaces can cause large errors. It is recommended that heat transfer data collected in this manner be used with caution.

Taylor, Robert P.↗
NASA NTRS2012

Simulation of a Wall-Bounded Flow using a Hybrid LES/RAS Approach with Turbulence Recycling

Simulations of a supersonic recessed-cavity flow are performed using a hybrid large-eddy/ Reynolds-averaged simulation approach utilizing an inflow turbulence recycling procedure and hybridized inviscid flux scheme. Calorically perfect air enters the three-dimensional domain at a free stream Mach number of 2.92. Simulations are performed to assess grid sensitivity of the solution, efficacy of the turbulence recycling, and effect of the shock sensor used with the hybridized inviscid flux scheme. Analysis of the turbulent boundary layer upstream of the rearward-facing step for each case indicates excellent agreement with theoretical predictions. Mean velocity and pressure results are compared to Reynolds-averaged simulations and experimental data for each case, and these comparisons indicate good agreement on the finest grid. Simulations are repeated on a coarsened grid, and results indicate strong grid density sensitivity. The effect of turbulence recycling on the solution is illustrated by performing coarse grid simulations with and without inflow turbulence recycling. Two shock sensors, one of Ducros and one of Larsson, are assessed for use with the hybridized inviscid flux reconstruction scheme.

Quinlan, Jesse R.↗
NASA NTRS2013

Hybrid Large-Eddy/Reynolds-Averaged Simulation of a Supersonic Cavity Using VULCAN

Simulations of a supersonic recessed-cavity flow are performed using a hybrid large-eddy/Reynolds-averaged simulation approach utilizing an inflow turbulence recycling procedure and hybridized inviscid flux scheme. Calorically perfect air enters a three-dimensional domain at a free stream Mach number of 2.92. Simulations are performed to assess grid sensitivity of the solution, efficacy of the turbulence recycling, and the effect of the shock sensor used with the hybridized inviscid flux scheme. Analysis of the turbulent boundary layer upstream of the rearward-facing step for each case indicates excellent agreement with theoretical predictions. Mean velocity and pressure results are compared to Reynolds-averaged simulations and experimental data for each case and indicate good agreement on the finest grid. Simulations are repeated on a coarsened grid, and results indicate strong grid density sensitivity. Simulations are performed with and without inflow turbulence recycling on the coarse grid to isolate the effect of the recycling procedure, which is demonstrably critical to capturing the relevant shear layer dynamics. Shock sensor formulations of Ducros and Larsson are found to predict mean flow statistics equally well.

Quinlan, Jesse↗
NASA NTRS2017

Time Resolved Stereo Particle Image Velocimetry Measurements of the Instabilities Downstream of a Backward-Facing Step in a Swept-Wing Boundary Layer

Time-resolved particle image velocimetry (TRPIV) measurements are performed down-stream of a swept backward-facing step, with a height of 49% of the boundary-layer thickness. The results agree well qualitatively with previously reported hotwire measurements, though the amplitudes of the fluctuating components measured using TRPIV are higher. Nonetheless, the low-amplitude instabilities in the flow are fairly well resolved using TR- PIV. Proper orthogonal decomposition is used to study the development of the traveling cross flow and Tollmien-Schlichting (TS) instabilities downstream of the step and to study how they interact to form the large velocity spikes that ultimately lead to transition. A secondary mode within the traveling cross flow frequency band develops with a wavelength close to that of the stationary cross flow instability, so that at a certain point in the phase, it causes an increase in the spanwise modulation initially caused by the stationary cross flow mode. This increased modulation leads to an increase in the amplitude of the TS mode, which, itself, is highly modulated through interactions with the stationary cross flow. When the traveling cross flow and TS modes align in time and space, the large velocity spikes occur. Thus, these three instabilities, which are individually of low amplitude when the spikes start to occur (U'rms/Ue <0.03), interact and combine to cause a large flow disturbance that eventually leads to transition.

Eppink, Jenna L.↗