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 667 records · Page 37

All sources NASA NTRS DOE OSTI

DOE OSTI2025

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↗

DOE OSTI2025

TEAMER – Enhanced Flow Measurement for Aquantis Tidal Turbine Test

The AQ10 is a floating, two-bladed, passive yawing tidal turbine developed by Aquantis that has a 10-meterrotor diameter, 160 kW rating, and employs reliable off-the-shelf powertrain and power conversion hardware. Aquantis is planning on-water turbine power performance and loads (blade loading and thrust)testing, where the turbine will be pushed through still water up to 4 knots and placed in a ‘station keeping’ tow in a tidal race up to its rated speed. On-water testing will be conducted using vessels and floating platforms on the sea surface to improve ease of testing and reduce disturbance to the environment. In this TEAMER project, Pacific Northwest National Laboratory (PNNL) will conduct water velocity and turbulence measurements in front of the turbine during on-water testing using acoustic Doppler instrumentation. By measuring both the tidal current flowing past the turbine and the resulting electrical power output, test results will provide a power curve (power vs flow speed) for the turbine up to rated power. Measurements of turbulence and velocity shear in front of the rotor will also provide information to assess the structural response of the rotor blades. With this analysis, Aquantis can use the performance and loads data to validate Tidal Bladed and OpenFAST simulations of the measured operating conditions. Measuring the power performance of a prototype turbine is a valuable step to improving device development and conducting a complete power performance assessment to IEC/TS 62600-200 standards in the future.

16 TIDAL AND WAVE POWER↗

NASA NTRS1982

A new numerical method for the simulation of three dimensional flow in a pipe

A new numerical technique for simulating three dimensional, unsteady, incompressible pipe flows is presented and its utility and accuracy is shown. Each vector function in the expansion of the velocity field is divergence free and satisfies the boundary conditions for viscous flow. Some of the benefits of the expansion technique are that pressure is eliminated from the dynamics, only two unknowns per mesh point are required, implicit treatment of the viscous terms is provided at no extra computational cost, and no fractional time steps are required. The method uses spectral expansions: Fourier series in the azimuthal and streamwise directions, and Jacobi polynominals in the radial direction.

Leonard, A.↗

NASA NTRS1982

A new numerical method for the simulation of three-dimensional flow in a pipe

Leonard, A.↗

NASA NTRS1986

Tether crawler system

A crawler system is designed to move a low-g/variable-g laboratory module along a tether between the Space Station and an attached space platform. An analysis is made of the effects of control law parameter change on the displacement, velocity, and acceleration of the crawler system. The control law is then modified by the addition of a constant-velocity section and the values of distance traveled, velocity, and acceleration are analyzed as a function of time. The power and torque equations are derived for a crawler system moving along a tether in orbit and numerical values of power and torque required for each prescribed movement are calculated versus time for four different cases using the control laws. A two-step control sequence is selected to permit initial location along the tether by distance traveled, followed by a vernier movement to reach the final desired constant net acceleration level. The components for the control system are identified and arranged in a block disgram configuration. The support subsystems are also identified. The sections were integrated to develop a procedure for the determination of crawler system performance requirements and the initial design of tether crawler systems.

Swanson, Frank R.↗

NASA NTRS1996

A viscoelastic higher-order beam finite element

A viscoelastic internal variable constitutive theory is applied to a higher-order elastic beam theory and finite element formulation. The behavior of the viscous material in the beam is approximately modeled as a Maxwell solid. The finite element formulation requires additional sets of nodal variables for each relaxation time constant needed by the Maxwell solid. Recent developments in modeling viscoelastic material behavior with strain variables that are conjugate to the elastic strain measures are combined with advances in modeling through-the-thickness stresses and strains in thick beams. The result is a viscous thick-beam finite element that possesses superior characteristics for transient analysis since its nodal viscous forces are not linearly dependent an the nodal velocities, which is the case when damping matrices are used. Instead, the nodal viscous forces are directly dependent on the material's relaxation spectrum and the history of the nodal variables through a differential form of the constitutive law for a Maxwell solid. The thick beam quasistatic analysis is explored herein as a first step towards developing more complex viscoelastic models for thick plates and shells, and for dynamic analyses. The internal variable constitutive theory is derived directly from the Boltzmann superposition theorem. The mechanical strains and the conjugate internal strains are shown to be related through a system of first-order, ordinary differential equations. The total time-dependent stress is the superposition of its elastic and viscous components. Equations of motion for the solid are derived from the virtual work principle using the total time-dependent stress. Numerical examples for the problems of relaxation, creep, and cyclic creep are carried out for a beam made from an orthotropic Maxwell solid.

Johnson, Arthur R.↗

NASA NTRS2010

Simulation of Micron-Sized Debris Populations in Low Earth Orbit

The update of ORDEM2000, the NASA Orbital Debris Engineering Model, to its new version ORDEM2010, is nearly complete. As a part of the ORDEM upgrade, this paper addresses the simulation of micro-debris (greater than 10 m and smaller than 1 mm in size) populations in low Earth orbit. The principal data used in the modeling of the micron-sized debris populations are in-situ hypervelocity impact records, accumulated in post-flight damage surveys on the space-exposed surfaces of returned spacecrafts. The development of the micro-debris model populations follows the general approach to deriving other ORDEM2010-required input populations for various components and types of debris. This paper describes the key elements and major steps in the statistical inference of the ORDEM2010 micro-debris populations. A crucial step is the construction of a degradation/ejecta source model to provide prior information on the micron-sized objects (such as orbital and object-size distributions). Another critical step is to link model populations with data, which is rather involved. It demands detailed information on area-time/directionality for all the space-exposed elements of a shuttle orbiter and damage laws, which relate impact damage with the physical properties of a projectile and impact conditions such as impact angle and velocity. Also needed are model-predicted debris fluxes as a function of object size and impact velocity from all possible directions. In spite of the very limited quantity of the available shuttle impact data, the population-derivation process is satisfactorily stable. Final modeling results obtained from shuttle window and radiator impact data are reasonably convergent and consistent, especially for the debris populations with object-size thresholds at 10 and 100 m.

Xu, Y.-L.↗

NASA NTRS2010

Simulation of Micron-Sized Debris Populations in Low Earth Orbit

The update of ORDEM2000, the NASA Orbital Debris Engineering Model, to its new version . ORDEM2010, is nearly complete. As a part of the ORDEM upgrade, this paper addresses the simulation of micro-debris (greater than 10 micron and smaller than 1 mm in size) populations in low Earth orbit. The principal data used in the modeling of the micron-sized debris populations are in-situ hypervelocity impact records, accumulated in post-flight damage surveys on the space-exposed surfaces of returned spacecrafts. The development of the micro-debris model populations follows the general approach to deriving other ORDEM2010-required input populations for various components and types of debris. This paper describes the key elements and major steps in the statistical inference of the ORDEM2010 micro-debris populations. A crucial step is the construction of a degradation/ejecta source model to provide prior information on the micron-sized objects (such as orbital and object-size distributions). Another critical step is to link model populations with data, which is rather involved. It demands detailed information on area-time/directionality for all the space-exposed elements of a shuttle orbiter and damage laws, which relate impact damage with the physical properties of a projectile and impact conditions such as impact angle and velocity. Also needed are model-predicted debris fluxes as a function of object size and impact velocity from all possible directions. In spite of the very limited quantity of the available shuttle impact data, the population-derivation process is satisfactorily stable. Final modeling results obtained from shuttle window and radiator impact data are reasonably convergent and consistent, especially for the debris populations with object-size thresholds at 10 and 100 micron.

Xu, Y.-L.↗

NASA NTRS2012

Frequency Correction for MIRO Chirp Transformation Spectroscopy Spectrum

This software processes the flyby spectra of the Chirp Transform Spectrometer (CTS) of the Microwave Instrument for Rosetta Orbiter (MIRO). The tool corrects the effect of Doppler shift and local-oscillator (LO) frequency shift during the flyby mode of MIRO operations. The frequency correction for CTS flyby spectra is performed and is integrated with multiple spectra into a high signal-to-noise averaged spectrum at the rest-frame RF frequency. This innovation also generates the 8 molecular line spectra by dividing continuous 4,096-channel CTS spectra. The 8 line spectra can then be readily used for scientific investigations. A spectral line that is at its rest frequency in the frame of the Earth or an asteroid will be observed with a time-varying Doppler shift as seen by MIRO. The frequency shift is toward the higher RF frequencies on approach, and toward lower RF frequencies on departure. The magnitude of the shift depends on the flyby velocity. The result of time-varying Doppler shift is that of an observed spectral line will be seen to move from channel to channel in the CTS spectrometer. The direction (higher or lower frequency) in the spectrometer depends on the spectral line frequency under consideration. In order to analyze the flyby spectra, two steps are required. First, individual spectra must be corrected for the Doppler shift so that individual spectra can be superimposed at the same rest frequency for integration purposes. Second, a correction needs to be applied to the CTS spectra to account for the LO frequency shifts that are applied to asteroid mode.

Lee, Seungwon↗

DOE OSTI2025

Integrated phononic waveguide on thin-film lithium niobate on diamond

Here, we demonstrate wavelength-scale phononic waveguides formed by transfer-printed thin-film lithium niobate (LN) on bulk diamond, a material stack that combines the strong piezoelectricity of LN with the high acoustic velocity and color-center compatibility of diamond. We characterize a delay line based on a 100 μm long phononic waveguide at room and cryogenic temperatures. The total insertion loss through the device at 4 K is −5.8 dB, corresponding to a $>$ 50 % transducer efficiency, at a frequency of 2.8 GHz. Our work represents a step toward phonon-mediated hybrid quantum systems consisting of strain-sensitive color centers in diamond.

Malik, Sultan [Stanford Univ., CA (United States)]↗

NASA NTRS1975

Unified nonlinear approach to both weak and strong-interaction problems

The numerical results are obtained for heat transfer, skin-friction, and viscous interaction induced pressure for a step-wise accelerated flat plate in hypersonic flow. In the unified approach here the results are presented for both weak and strong-interaction problems without employing any linearization scheme. With the help of the numerical method used in this work an accurate prediction of wall shear can be made for the problems with plate velocity changes of 1% or larger. The obtained results indicate that the transient contribution to the induced pressure for helium is greater than that for air.

Gupta, R. N.↗

NASA NTRS1981

Analytical description of charged particle transport along arbitrary guiding field configurations

Aspects of charged particle transport theory have been studied by Earl (1973, 1974, 1976) in a number of investigations. The investigations were conducted under the assumption that particles propagate along a DC guiding field with negligible perpendicular transport and that pitch-angle scattering by random fields can be described by the Fokker-Planck formalism. The present study provides a generalized and simplified treatment of focused transport. The first step is taken with respect to a perturbation analysis by which the Boltzmann equation is solved in terms of successively improved approximations to its exact solution. The first approximation is the sum of isotropic pseudodiffusive and anisotropic supercoherent components whose coefficients must satisfy two coupled partial differential equations. The normalization constant K and the pseudodiffusive velocity, which appear in the first-order transport equations, are given by simple integrals over the pitch angle.

Earl, J. A.↗

NASA NTRS1987

An approach to multivariable control of manipulators

The paper presents simple schemes for multivariable control of multiple-joint robot manipulators in joint and Cartesian coordinates. The joint control scheme consists of two independent multivariable feedforward and feedback controllers. The feedforward controller is the minimal inverse of the linearized model of robot dynamics and contains only proportional-double-derivative (PD2) terms - implying feedforward from the desired position, velocity and acceleration. This controller ensures that the manipulator joint angles track any reference trajectories. The feedback controller is of proportional-integral-derivative (PID) type and is designed to achieve pole placement. This controller reduces any initial tracking error to zero as desired and also ensures that robust steady-state tracking of step-plus-exponential trajectories is achieved by the joint angles. Simple and explicit expressions of computation of the feedforward and feedback gains are obtained based on the linearized model of robot dynamics. This leads to computationally efficient schemes for either on-line gain computation or off-line gain scheduling to account for variations in the linearized robot model due to changes in the operating point. The joint control scheme is extended to direct control of the end-effector motion in Cartesian space. Simulation results are given for illustration.

Seraji, H.↗

NASA NTRS1988

A velocity-pressure integrated, mixed interpolation, Galerkin finite element method for high Reynolds number laminar flows

A velocity-pressure integrated, mixed interpolation, Galerkin finite element method for the Navier-Stokes equations is presented. In the method, the velocity variables were interpolated using complete quadratic shape functions and the pressure was interpolated using linear shape functions. For the two dimensional case, the pressure is defined on a triangular element which is contained inside the complete biquadratic element for velocity variables; and for the three dimensional case, the pressure is defined on a tetrahedral element which is again contained inside the complete tri-quadratic element. Thus the pressure is discontinuous across the element boundaries. Example problems considered include: a cavity flow for Reynolds number of 400 through 10,000; a laminar backward facing step flow; and a laminar flow in a square duct of strong curvature. The computational results compared favorable with those of the finite difference methods as well as experimental data available. A finite elememt computer program for incompressible, laminar flows is presented.

Kim, Sang-Wook↗

NASA NTRS1992

Geomechanical models of impact cratering: Puchezh-Katunki structure

Impact cratering is a complex natural phenomenon that involves various physical and mechanical processes. Simulating these processes may be improved using the data obtained during the deep drilling at the central mound of the Puchezh-Katunki impact structure. A research deep drillhole (named Vorotilovskaya) has been drilled in the Puchezh-Katunki impact structure (European Russia, 57 deg 06 min N, 43 deg 35 min E). The age of the structure is estimated at about 180 to 200 m.y. The initial rim crater diameter is estimated at about 40 km. The central uplift is composed of large blocks of crystalline basement rocks. Preliminary study of the core shows that crystalline rocks are shock metamorphosed by shock pressure from 45 GPa near the surface to 15-20 GPa at a depth of about 5 km. The drill core allows the possibility of investigating many previously poorly studied cratering processes in the central part of the impact structure. As a first step one can use the estimates of energy for the homogeneous rock target. The diameter of the crater rim may be estimated as 40 km. The models elaborated earlier show that such a crater may be formed after collapse of a transient cavity with a radius of 10 km. The most probable range of impact velocities from 11.2 to 30 km/s may be inferred for the asteroidal impactor. For the density of a projectile of 2 g/cu cm the energy of the impact is estimated as 1E28 to 3E28 erg. In the case of vertical impact, the diameter of an asteroidal projectile is from 1.5 to 3 km for the velocity range from 11 to 30 km/s. For the most probable impact angle of 45 deg, the estimated diameter of an asteroid is slightly larger: from 2 to 4 km. Numerical simulation of the transient crater collapse has been done using several models of rock rheology during collapse. Results show that the column at the final position beneath the central mound is about 5 km in length. This value is close to the shock-pressure decay observed along the drill core. Further improvement of the model needs to take into account the blocky structure of target rocks revealed by drilling.

Ivanov, B. A.↗

NASA NTRS1993

A normal stress subgrid-scale eddy viscosity model in large eddy simulation

The Smagorinsky subgrid-scale eddy viscosity model (SGS-EVM) is commonly used in large eddy simulations (LES) to represent the effects of the unresolved scales on the resolved scales. This model is known to be limited because its constant must be optimized in different flows, and it must be modified with a damping function to account for near-wall effects. The recent dynamic model is designed to overcome these limitations but is compositionally intensive as compared to the traditional SGS-EVM. In a recent study using direct numerical simulation data, Horiuti has shown that these drawbacks are due mainly to the use of an improper velocity scale in the SGS-EVM. He also proposed the use of the subgrid-scale normal stress as a new velocity scale that was inspired by a high-order anisotropic representation model. The testing of Horiuti, however, was conducted using DNS data from a low Reynolds number channel flow simulation. It was felt that further testing at higher Reynolds numbers and also using different flows (other than wall-bounded shear flows) were necessary steps needed to establish the validity of the new model. This is the primary motivation of the present study. The objective is to test the new model using DNS databases of high Reynolds number channel and fully developed turbulent mixing layer flows. The use of both channel (wall-bounded) and mixing layer flows is important for the development of accurate LES models because these two flows encompass many characteristic features of complex turbulent flows.

Horiuti, K.↗

NASA NTRS1997

Effect of 17 days of bed rest on peak isometric force and unloaded shortening velocity of human soleus fibers

The purpose of this study was to examine the effect of prolonged bed rest (BR) on the peak isometric force (P0) and unloaded shortening velocity (V0) of single Ca(2+)-activated muscle fibers. Soleus muscle biopsies were obtained from eight adult males before and after 17 days of 6 degrees head-down BR. Chemically permeabilized single fiber segments were mounted between a force transducer and position motor, activated with saturating levels of Ca2+, and subjected to slack length steps. V0 was determined by plotting the time for force redevelopment vs. the slack step distance. Gel electrophoresis revealed that 96% of the pre- and 87% of the post-BR fibers studied expressed only the slow type I myosin heavy chain isoform. Fibers with diameter > 100 microns made up only 14% of this post-BR type I population compared with 33% of the pre-BR type I population. Consequently, the post-BR type I fibers (n = 147) were, on average, 5% smaller in diameter than the pre-BR type I fibers (n = 218) and produced 13% less absolute P0. BR had no overall effect on P0 per fiber cross-sectional area (P0/CSA), even though half of the subjects displayed a decline of 9-12% in P0/CSA after BR. Type I fiber V0 increased by an average of 34% with BR. Although the ratio of myosin light chain 3 to myosin light chain 2 also rose with BR, there was no correlation between this ratio and V0 for either the pre- or post-BR fibers. In separate fibers obtained from the original biopsies, quantitative electron microscopy revealed a 20-24% decrease in thin filament density, with no change in thick filament density. These results raise the possibility that alterations in the geometric relationships between thin and thick filaments may be at least partially responsible for the elevated V0 of the post-BR type I fibers.

STS-78 Shuttle Project↗

NASA NTRSDate not supplied

Expanding the Capability of A Legacy Combustion Flametube to Test High Temperature Engine Materials in Relevant Environments

New materials and component designs are needed to advance gas turbine engine technology and provide the performance and efficiency needs for future applications. In order to advance these materials, testing in combustion environments is a critical step prior to engine testing. In this work, we detail the design and the fabrication of a materials test sector in a flametube combustor facility. The facility simulates a combustion environment similar to that experienced by components in gas turbine engines. The flow regime is characterized by a combination of high-temperature, high-velocity, high-heat flux, and high-velocity that components experience in gas turbine engines. Exposure of components in this facility allows for the study of combined environmental effects and the impact on both coating and substrate durability. The test facility may operate across a wide range of pressures from 275-400 psig (1,896-2,758 kPa) and an air flow rate of 5 lb/s (2.27 kg/s). While combustion gas temperature is expected in excess of 3,000°F (1,649°C), 900°F (482°C) cooling air may be supplied to the backside of components or test articles. The flametube combustor was previously used to evaluate fuel injectors and combustion products, and the new test configuration will also allow for materials exposure to complex, engine-like conditions. The interior of the Test Section was additively manufactured from GRCop-84 and cryogenically fit and brazed to a stainless steel 304 housing. The use of a copper liner minimizes welds and with active cooling, is expected to provide better durability over traditional hardware using stainless steel or Inconel with a ceramic liner. The Test Section has two opposing removable windows approximately 230 mm x 80 mm that can accommodate articles up to 85 mm tall. This modular design allows for custom platforms to hold coupons, panels, or airfoil shapes to be tested with minimal re-engineering or fabrication. The bolted joint and sealing remains consistent, so any new testing only needs to work within the existing design footprint. This paper will provide an overview of the facility capabilities, design considerations, as well as thermal and structural analysis of the hardware. Future testing of ceramic matrix composite (CMC) airfoils and advanced environmental barrier coatings (EBCs) will also be discussed.

Combustion↗