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.

SourceResearch connectionFrom yearThrough year

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

At least 649 records · Page 36

All sources NASA NTRS DOE OSTI

NASA NTRS1993

An investigation of small scales of turbulence in a boundary layer at high Reynolds numbers

The assumption that turbulence at large wave-numbers is isotropic and has universal spectral characteristics which are independent of the flow geometry, at least for high Reynolds numbers, has been a cornerstone of closure theories as well as of the most promising recent development in the effort to predict turbulent flows, viz. large eddy simulations. This hypothesis was first advanced by Kolmogorov based on the supposition that turbulent kinetic energy cascades down the scales (up the wave-numbers) of turbulence and that, if the number of these cascade steps is sufficiently large (i.e. the wave-number range is large), then the effects of anisotropies at the large scales are lost in the energy transfer process. Experimental attempts were repeatedly made to verify this fundamental assumption. However, Van Atta has recently suggested that an examination of the scalar and velocity gradient fields is necessary to definitively verify this hypothesis or prove it to be unfounded. Of course, this must be carried out in a flow with a sufficiently high Reynolds number to provide the necessary separation of scales in order unambiguously to provide the possibility of local isotropy at large wave-numbers. An opportunity to use our 12-sensor hot-wire probe to address this issue directly was made available at the 80'x120' wind tunnel at the NASA Ames Research Center, which is normally used for full-scale aircraft tests. An initial report on this high Reynolds number experiment and progress toward its evaluation is presented.

Wallace, James M.↗

NASA NTRS1994

Ignition dynamics of a laminar diffusion flame in the field of a vortex embedded in a shear flow

The role of streamwise-spanwise vorticity interactions that occur in turbulent shear flows on flame/vortex interactions is examined by means of asymptotic analysis and numerical simulation in the limit of small Mach number. An idealized model is employed to describe the interaction process. The model consists of a one-step, irreversible Arrhenius reaction between initially unmixed species occupying adjacent half-planes which are then allowed to mix and react in the presence of a streamwise vortex embedded in a shear flow. It is found that the interaction of the streamwise vortex with shear gives rise to small-scale velocity oscillations which increase in magnitude with shear strength. These oscillations give rise to regions of strong temperature gradients via viscous heating, which can lead to multiple ignition points and substantially decrease ignition times. The evolution in time of the temperature and mass-fraction fields is followed, and emphasis is placed on the ignition time and structure as a function of vortex and shear strength.

Macaraeg, Michele G.↗

NASA NTRS1993

A Fractional-Step Method Of Computing Incompressible Flow

Method of computing time-dependent flow of incompressible, viscous fluid involves numerical solution of Navier-Stokes equations on two- or three-dimensional computational grid based on generalized curvilinear coordinates. Equations of method derived in primitive-variable formulation. Dependent variables are pressure at center of each cell of computational grid and volume fluxes across faces of each cell. Volume fluxes replace Cartesian components of velocity; these fluxes correspond to contravariant components of velocity multiplied by volume of computational cell, in staggered grid. Choice of dependent variables enables simple extension of previously developed staggered-grid approach to generalized curvilinear coordinates and facilitates enforcement of conservation of mass.

Kwak, Dochan↗

DOE OSTI2025

A Fundamentally New Coupled Approach to Contact Mechanics via the Dirichlet‐Neumann Schwarz Alternating Method

Contact phenomena are crucial for understanding the behavior of mechanical systems. However, existing computational approaches for simulating mechanical contact often face numerical challenges, such as inaccurate physical predictions, energy conservation errors, and unwanted oscillations. Here, we introduce an alternative technique for simulating dynamic contact based on the non‐overlapping Schwarz alternating method, originally developed for domain decomposition. In multibody contact scenarios, this method treats each body as a separate, non‐overlapping domain and prevents interpenetration using an alternating Dirichlet–Neumann iterative process. This approach has a strong theoretical foundation, eliminates the need for contact constraints, and offers flexibility, making it ideal for multiscale and multiphysics applications. We conducted a numerical comparison between the Schwarz method and traditional methods, such as the Lagrange multiplier and penalty methods, focusing on a benchmark impact problem. Our results indicate that the Schwarz alternating method outperforms traditional methods in several key areas: it provides more accurate predictions for various measurable quantities and demonstrates exceptional energy conservation capabilities. To address unwanted oscillations in contact velocities and forces, we explored various algorithms and stabilization techniques, ultimately opting for the naïve‐stabilized Newmark scheme for its simplicity and effectiveness. Additionally, we validated the efficiency of the Schwarz method in a three‐dimensional impact problem, highlighting its inherent capacity to accommodate different mesh topologies, time‐integration schemes, and time steps for each interacting body.

Schwarz alternating method↗

DOE OSTI2022

A mass–momentum consistent coupling for mesh-adaptive two-phase flow simulations

Here, we present a novel mass-momentum consistent coupling between a geometric volume-of-fluid scheme and an incompressible flow solver with differing directional-splitting approaches. The advection of the volume fraction is performed using a direction-split algorithm, whereas the momentum advection algorithm uses a traditional unsplit, fractional-step approach. Both algorithms employ finite-volume discretizations based on Cartesian meshes. In solving the mass-momentum consistency problem, momentum fluxes at the cell faces are weighted by the density fluxes based on the already advected volume fraction. The success of our approach lies on introducing a Favre-averaged velocity interpolation at the liquid/gas interface along with a minmod slope limiter. Mesh-convergence studies show that when the minmod slope limiter is used, the two-phase solver retains an accuracy between first and second order, but when a purely upwind scheme is considered, its accuracy drops to first order. Finally, after considering several validation problems, the solver is shown to agree well with reference numerical and experimental data while retaining its robustness and efficiency.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

The effect of non-linear mutual friction on pulsar glitch sizes and rise times

ABSTRACT Observations of pulsar glitches have the potential to provide constraints on the dynamics of the high density interior of neutron stars. However, to do so, realistic glitch models must be constructed and compared to the data. We take a step towards this goal by testing non-linear models for the mutual friction force, which is responsible for the exchange of angular momentum between the neutron superfluid and the observable normal component in a glitch. In particular, we consider a non-linear dependence of the drag force on the relative velocity between superfluid vortices and the normal component, in which the contributions of both kelvin and phonon excitations are included. This non-linear model produces qualitatively new features, and is able to reproduce the observed bimodal distribution of glitch sizes in the pulsar population. The model also suggests that the differences in size distributions in individual pulsars may be due to the glitches being triggered in regions with different pinning strengths, as stronger pinning leads to higher vortex velocities and a qualitatively different mutual friction coupling with respect to the weak pinning case. Glitches in pulsars that appear to glitch quasi-periodically with similar sizes may thus be due to the same mechanisms as smaller events in pulsars that have no preferred glitch size, but simply originate in stronger pinning regions, possibly in the core of the star.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1982

A scheme for tracking and pointing at small celestial bodies during high velocity encounters

A scheme has been developed and verified for closed-loop tracking and pointing space-borne science instruments at small celestial bodies, such as comets and asteroids, during high velocity encounters. To overcome ephemeris uncertainties for these bodies, the scheme involves sequential estimation of flyby model parameters. The design consists of a two-axis gimballed platform mounted on a three-axis stabilized spacecraft. A platform-mounted optical tracker provides closed-loop target measurements and precision micro-step actuators enable high-rate platform slewing. For comet missions which involve dust particle impact disturbances, a dual-mode attitude control scheme is presented for minimizing transient response time.

Burdick, G. M.↗

NASA NTRS1984

Operation of a single-channel, sequential Navstar GPS receiver in a helicopter mission environment

It is pointed out that the future utilization of the Navstar Global Positioning System (GPS) by civil helicopters will provide an enhanced performance not obtainable with current navigations systems. GPS will supply properly equipped users with extremely accurate three-dimensional position and velocity information anywhere in the world. Preliminary studies have been conducted to investigate differential GPS concept mechanizations and cost, and to theoretically predict navigation performance and the impact of degradation of the GPS C/A code for national security considerations. The obtained results are encouraging, but certain improvements are needed. As a second step in the program, a single-channel sequential GPS navigator was installed and operated in the NASA SH-3G helicopter. A series of flight tests were conducted. It is found that performance of the Navstar GPS Z-set is quite acceptable to support area navigation and nonprecision approach operations.

Edwards, F. G.↗

NASA NTRS1946

Effect of compressibility at high subsonic velocities on the lifting force acting on an elliptic cylinder

An extended form of the Ackeret iteration method, applicable to arbitrary profiles, is utilized to calculate the compressible flow at high subsonic velocities past an elliptic cylinder. The angle of attack to the direction of the undisturbed stream is small and the circulation is fixed by the Kutta condition at the trailing end of the major axis. The expression for the lifting force on the elliptic cylinder is derived and shows a first-step improvement of the Prandtl-Glauert rule. It is further shown that the expression for the lifting force, although derived specifically for an elliptic cylinder, may be extended to arbitrary symmetrical profiles.

Kaplan, Carl↗

NASA NTRS1952

An empirically derived basis for calculating the area, rate, and distribution of water-drop impingement on airfoils

An empirically derived basis for predicting the area, rate, and distribution of water-drop impingement on airfoils of arbitrary section is presented. The concepts involved represent an initial step toward the development of a calculation technique which is generally applicable to the design of thermal ice-prevention equipment for airplane wing and tail surfaces. It is shown that sufficiently accurate estimates, for the purpose of heated-wing design, can be obtained by a few numerical computations once the velocity distribution over the airfoil has been determined. The calculation technique presented is based on results of extensive water-drop trajectory computations for five airfoil cases which consisted of 15-percent-thick airfoils encompassing a moderate lift-coefficient range. The differential equations pertaining to the paths of the drops were solved by a differential analyzer.

Bergrun, Norman R↗

NASA NTRS1994

The Flux-integral Method for Multidimensional Convection and Diffusion

The flux-integral method is a procedure for constructing an explicit, single-step, forward-in-time, conservative, control volume update of the unsteady, multidimensional convection-diffusion equation. The convective plus diffusive flux at each face of a control-volume cell is estimated by integrating the transported variable and its face-normal derivative over the volume swept out by the convecting velocity field. This yields a unique description of the fluxes, whereas other conservative methods rely on nonunique, arbitrary pseudoflux-difference splitting procedures. The accuracy of the resulting scheme depends on the form of the subcell interpolation assumed, given cell-average data. Cellwise constant behavior results in a (very artificially diffusive) first-order convection scheme. Second-order convection-diffusion schemes correspond to cellwise linear (or bilinear) subcell interpolation. Cellwise quadratic subcell interpolants generate a highly accurate convection-diffusion scheme with excellent phase accuracy. Under constant-coefficient conditions, this is a uniformly third-order polynomial interpolation algorithm (UTOPIA).

Leonard, B. P.↗

NASA NTRS1998

Model Refinement Using Eigensystem Assignment

This paper presented a novel approach for the refinement of finite-element-based analytical models of flexible structures is presented. The proposed approach models the possible refinements in the mass, damping, and stiffness matrices of the finite element model in the form of a constant gain feedback with acceleration, velocity, and displacement measurements, respectively. Once, the free elements of the structural matrices have been defined, the problem of model refinement reduces to obtaining position, velocity, and acceleration gain matrices, which reassign a desired subset of the eigenvalues of the model, along with partial mode shapes, from their baseline values to those obtained from system identification test data. A sequential procedure is used to assign one self-conjugate pair of closed-loop eigenvalues at each step using symmetric output feedback gain matrices, and the closed-loop eigenvectors are partially assigned, while ensuring that the eigenvalues assigned in the previous steps are not disturbed. The procedure can also impose that gain matrices be dissipative in order to guarantee the stability of the refined model. A numerical example, involving finite element model refinement for a structural testbed at NASA Langley (CSI Evolutionary Model) is presented to demonstrate the feasibility of the proposed approach.

Maghami, Pieman G.↗

NASA NTRS2022

Predicting the Penetration of a Shielded TPS Tile

All spacecraft are subject to the possibility of high-speed particle impacts during their mission life. In low earth orbit, those impacts could be the result of collisions with pieces of orbital debris or with meteorites. Beyond LEO, and especially beyond GEO, those impacts will likely be caused by meteorites. Such high-speed impacts on spacecraft surfaces create debris clouds that travel towards and eventually impact other downstream spacecraft components. In addition to the impulsive load that such debris clouds would impart to the spacecraft elements with which they subsequently collide, the largest fragment in these debris clouds poses a significant threat on its own to those spacecraft elements. In order to be able to assess the severity of the threat posed by such a fragment, it is important to be able to predict the extent of damage sustained by the impacted spacecraft element. In this paper, we present a new process for determining the penetration depth in a shielded TPS tile system. The requirements for the process were that it be applicable across a full spectrum of impact low velocities (i.e. from ~ 7 km/s to above 50 km/s), that it be applicable over a wide range of projectile materials densities (i.e. ranging from water to steel, for example), and that it be applicable over a wide range of trajectory obliquities (i.e. not just normal impacts). Furthermore, the process developed should be sufficiently flexible so that it can be used for shield design parameter and / or system configuration trade studies (i.e. be written in terms of as many shield and TPS tile system parameters as possible). Following the initial impact of the projectile on the outer wall of the dual-wall shield, a debris cloud (i.e. the primary debris cloud) is created that travels towards and impacts the inner wall of the shield. Within this debris cloud is a combination of solid, liquid, and vaporized material, depending on the impact velocity and the impedance mismatch of the projectile and outer wall materials. The impact of this debris cloud on the inner wall creates another debris cloud (i.e. the secondary debris cloud) that then travels towards and impacts the TPS tile. Within the particulate distribution of the material in each of these debris clouds is a so-called largest fragment. The process used to calculate the penetration depth in a shielded TPS tile system developed herein consists of the following three steps: 1. Does the initial projectile perforate the outer wall of the shield? If NO, the TPS tile remains undamaged. If YES, calculate the size and speed of the largest fragment exiting the outer wall of the shield and traveling towards the shield’s inner wall. 2. Does the largest fragment in the debris cloud exiting the outer wall of the shield perforate the inner wall of the shield? If NO, the TPS tile again remains undamaged. If YES, calculate the size and speed of the largest fragment exiting the rear of the shield inner wall. 3. Calculate the penetration depth of the largest fragment in the debris cloud exiting the inner wall into the TPS. Perforation of the inner and outer shield walls is predicted using the Cour-Palais penetration depth equation and an associated failure criterion [1]. The largest debris cloud fragment diameters and associated velocities are calculated using regression equations developed from a mix of hydrocode and experimental diameter and velocity information [2]. Finally, TPS penetration depths are calculated using an empirical penetration depth predictor equation based on high-speed impact test data and are calibrated using a limited number of hydrocode penetration depth predictions for water, nylon, aluminum, Dunite, and steel projectiles The predictions of the penetration depth calculation process are compared against the predictions of more than 60 SPHC hydrocode runs. These comparisons show that the process used to calculate TPS tile penetration depths usually yields values that are within 1 mm of the predictions of hydrocode simulations. If a maximum allowable TPS penetration depth is known, the process developed herein can now be used to develop a failure limit equation that would predict whether or not an impacting particle would penetrate deeply enough into the TPS to result in a “failed” or “not failed” end state. References [1] E.L. Christiansen, “Design and performance equations for advanced meteoroid and debris shields”, International Journal of Impact Engineering, Vol. 14, 1993, pp. 145-156. [2] W.P. Schonberg, “Predicting the size of the largest particle fragment in a debris cloud created by an orbital debris impact and its associated velocity”, 2021 Applied Space Environments Conference, NASA Jet Propulsion Laboratory, Pasadena, California, November, 2021.

Meteroid↗

NASA NTRS1982

Scout trajectory error propagation computer program

Since 1969, flight experience has been used as the basis for predicting Scout orbital accuracy. The data used for calculating the accuracy consists of errors in the trajectory parameters (altitude, velocity, etc.) at stage burnout as observed on Scout flights. Approximately 50 sets of errors are used in Monte Carlo analysis to generate error statistics in the trajectory parameters. A covariance matrix is formed which may be propagated in time. The mechanization of this process resulted in computer program Scout Trajectory Error Propagation (STEP) and is described herein. Computer program STEP may be used in conjunction with the Statistical Orbital Analysis Routine to generate accuracy in the orbit parameters (apogee, perigee, inclination, etc.) based upon flight experience.

Myler, T. R.↗

NASA NTRS1989

Estimation of the two-dimensional ocean current shear field with a synthetic aperture radar

A double inverse method to convert SAR data to estimates of current shear was developed theoretically. Problems of applying the technique to real SAR data were examined. The first step of the inversion process involved elimination of noise in the (x-vector, k-vector) wave spectra deduced from SAR images. The second step, estimation of current shear, required the existence of well defined swell and a favorable single of interaction between swell and ocean current. Application of the method in the Antarctic circumpolar current encountered sea conditions that both satisfied and violated these requirements. Where the wave-current geometry was favorable, realistic surface current velocity fields were obtained, suggesting the presence of a substantial mesoscale eddy. Another attempt to apply the method in the Gulf Stream produced wave spectral estimates that showed no clear peak and hence produced no current estimates.

Barnett, Tim P.↗

NASA NTRS1990

The diffuse interstellar medium

The last 20 years of the efforts to understand the diffuse ISM are reviewed, with recent changes of fundamental aspects being highlighted. Attention is given to the interstellar pressure and its components, the weight of the ISM, the midplane pressure contributions, and pressure contributions at 1 kpc. What velocity dispersions, cosmic ray pressure, and magnetic field pressure that can be expected for a gas in a high magnetic field environment is addressed. The intercloud medium is described, with reference to the work of Cox and Slavin (1989). Various caveats are discussed and a number of areas for future investigation are identified. Steps that could be taken toward a successful phase segregation model are discussed.

Cox, Donald P.↗

NASA NTRS2016

Small Scale Motions Observed by Aircraft in the Tropical Tropopause Layer - Convective and Non-Convective Environments

Indirect evidence indicates a role for vertical mixing in the Tropical Tropopause Layer (TTL). In particular, detailed model studies suggest that such vertical mixing may be required to explain the value of the water vapor minimum in the TTL. There have been previous observations during the STEP Tropical aircraft campaign (1987) of bursts of high frequency activity associated with convectively generated gravity waves in the tropical western Pacific. Higher frequency, higher quality measurements from NASA high altitude aircraft (ER-2, WB-57, and Global Hawk) have been made available in the last 20 years. These include measurements of vertical velocity and other meteorological parameters. Most recently, during the ATTREX Global Hawk aircraft mission (Airborne Tropical TRopopause EXperiment), there have been extensive measurements at all altitudes of the TTL in both convective (winter western Pacific) and less convective (winter eastern Pacific) regions. This presentation represents an initial analysis of high frequency small scale (a few km max) meteorological measurements from the ATTREX dataset. We obtain some basic information about the distribution and character of high frequency activity in vertical velocity in the TTL. In particular, we focus on relating the high frequency activity to nearby tropical convection and to vertical shears associated with gravity and inertia-gravity waves.

convection↗

DOE OSTI2026

Sap Velocity Data for Urban Trees in Chicago, Illinois (2024-2025)

This dataset contains uncorrected sap velocity measurements using the heat ratio method (HRM) collected using ICT International SFM1x sensors at five urban sites in Chicago, Illinois, as part of the DOE CROCUS project. The data includes continuous monitoring of sap velocity from various tree species, including Maples (Acer spp.): Sugar Maple (Acer saccharum), Silver Maple (Acer saccharinum), and Red Maple (Acer rubrum); Oaks (Quercus spp.): Swamp White Oak (Quercus bicolor); American Elm (Ulmus americana); Honey Locust (Gleditsia triacanthos); Cottonwood (Populus deltoides); and Tree of Heaven (Ailanthus altissima) across Chicago State University (CSU), Northeastern Illinois University (NEIU), Northwestern University (NU), University of Illinois Chicago (UIC), and West Woodlawn "Blacks in Green" (BIG). These include both street trees and those in urban park locations. Measurements were collected at 15-20 minute intervals, depending on the sensor, and transmitted via Long Range Wide Area Network (LoRaWAN) protocols. The wireless data was collected by Sage Network (https://sagecontinuum.org/) nodes. The dataset includes sensor ID, Global Positioning System (GPS) coordinates, tree species (common and scientific names), tree identification number, diameter at breast height (DBH in cm), uncorrected sap velocity measurements (cm/hr) from both inner and outer probes, and Sage Node identifiers so the data can be mapped to related variables such as air quality and wind speed that were collected on the Sage nodes. All timestamps are in local Chicago time (CDT/CST). Quality control flags are provided using a 3-bit binary system indicating physical range violations (< -10 or > 60 cm/hr), step spikes (absolute difference > 36 cm/hr), and stuck sensor conditions (> 10 consecutive identical values). These are raw data, not corrected for wood anatomy or species-specific characteristics. Data is provided in comma separated (CSV) format. This dataset is part of a larger collection of CROCUS environmental monitoring data, including linked datasets from Air Quality Transmitter (AQT) sensors, Weather Transmitter (WXT) sensors, and Multi-Function Research LoRaWAN (MFR) Nodes. DOIs for the supporting data are provided as part of this data package.

Chicago↗