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 847 records · Page 47

All sources NASA NTRS DOE OSTI

NASA NTRS1989

Effect of wall suction on the stability of compressible subsonic flows over smooth two-dimensional backward-facing steps

The effect of suction on the stability of compressible flows over backward-facing steps is investigated. Mach numbers up to 0.8 are considered. The results show that continuous suction stabilizes the flow outside the separation bubble, but it destabilizes the flow inside it. Nevertheless, the overall N factor decreases as the suction level increases due to the considerable reduction of the separation bubble. For the same suction flow rate, properly distributed suction strips stabilize the flow more than continuous suction. The size of the separation bubble, and hence its effect on the instability can be considerably reduced by placing strips with high suction velocities in the separation region.

Al-Maaitah, Ayman A.↗

NASA NTRS2024

Simplified Structure Models for Premixed n-Alkane Cool Flames

The chemical kinetics of cool flames in mixtures of normal alkanes, oxygen, and nitrogen recently have been reduced to obtain simplified asymptotic flame-structure models that can be used to compare theoretical predictions with experimental measurements. A key step that affects flame structures strongly, namely the abstraction of a hydrogen atom from the alkyl radical by an oxygen molecule, had not been taken into account in those models. This deficiency has now been remedied for diffusion flames. The present contribution reports the corresponding remedy for premixed-flame analyses. In the new model, non-zero chemical reaction rates are restricted to narrow regions in the vicinity of the maximum temperature, whereas previously they extended throughout the preheat zone. Comparisons of the new predictions with experimental burning-velocity measurements seem favorable.

premixed cool flames↗

NASA NTRS1983

The use of dynamical constraints in the analysis of mesoscale rawinsonde data

A process for data interpolation which seeks to maintain dynamical consistency is applied to the analysis of mesoscale rawinsonde data obtained from the National Severe Storms Laboratory. The analysis procedure is in two steps: a successive-corrections interpolation to gridpoints, followed by a variational adjustment which uses a set of forecast equations of horizontal momentum as dynamical constraints. The objective analysis is applied in case studies of conditions on May 22, 1966 and June 8, 1966, and it is shown that the adjusted winds and geopotentials have greatly reduced residuals compared to the initial fields. Adjustments to the wind and geopotential fields have root mean square values of 1-2 m/sec and 30-50 sq m/sec per sec, respectively. In addition, the vertical velocity fields computed from the adjusted winds are found to be in closer agreement with actual weather data.

Bloom, S. C.↗

NASA NTRS1991

Computation of three-dimensional flows using two stream functions

An approach to compute 3-D flows using two stream functions is presented. The method generates a boundary fitted grid as part of its solution. Commonly used two steps for computing the flow fields are combined into a single step in the present approach: (1) boundary fitted grid generation; and (2) solution of Navier-Stokes equations on the generated grid. The presented method can be used to directly compute 3-D viscous flows, or the potential flow approximation of this method can be used to generate grids for other algorithms to compute 3-D viscous flows. The independent variables used are chi, a spatial coordinate, and xi and eta, values of stream functions along two sets of suitably chosen intersecting stream surfaces. The dependent variables used are the streamwise velocity, and two functions that describe the stream surfaces. Since for a 3-D flow there is no unique way to define two sets of intersecting stream surfaces to cover the given flow, different types of two sets of intersecting stream surfaces are considered. First, the metric of the (chi, xi, eta) curvilinear coordinate system associated with each type is presented. Next, equations for the steady state transport of mass, momentum, and energy are presented in terms of the metric of the (chi, xi, eta) coordinate system. Also included are the inviscid and the parabolized approximations to the general transport equations.

Greywall, Mahesh S.↗

NASA NTRS1993

Free-wake computation of helicopter rotor flowfields in forward flight

A new method has been developed for computing advancing rotor flows. This method uses the Vorticity Embedding technique, which has been developed and validated over the last several years for hovering rotor problems. In this work, the unsteady full potential equation is solved on an Eulerian grid with an embedded vortical velocity field. This vortical velocity accounts for the influence of the wake. Dynamic grid changes that are required to accommodate prescribed blade motion and deformation are included using a novel grid blending method. Free wake computations have been performed on a two-bladed AH-1G rotor at low advance ratios including blade motion. Computed results are compared with experimental data. The sudden variations in airloads due to blade-vortex interactions on the advancing and retreating sides are well captured. The sensitivity of the computed solution to various factors like core size, time step and grids has been investigated. Computed wake geometries and their influence on the aerodynamic loads at these advance ratios are also discussed.

Ramachandran, K.↗

DOE OSTI2026

Fragmentation dynamics of CS 2 dications and trications following S 2p ionization

Here, we present the results from a detailed study of the fragmentation dynamics of CS$_2^{2+}$ and CS$_2^{3+}$, formed in intense femtosecond soft x-ray pulses above the sulfur 2p edge, primarily through single core photoionization from the S 2p site, and subsequent Auger–Meitner decay(s). By combining three-dimensional velocity map imaging with covariance analysis, we determine the relative momenta of the ions produced in each two- and three-body fragmentation channel, at significantly higher ion count rates than conventional coincidence measurements. We shed new light on the wide range of fragmentation channels observed from the CS 2 dication and trication, including channels that involve ionization-induced bond formation and fragmentations producing undetected neutral cofragments. In the latter case, a “native frames” approach is used to isolate contributions from concerted and sequential fragmentations and extract dynamical information about each step of a concerted fragmentation process. While dications often fragment sequentially, the trication is dominated by concerted fragmentation. The main trication fragmentation channel into S + + C + + S + can be well-approximated by classical Coulombic simulations of the ground-state geometry distribution, reflecting both the nature of the trication potential energy surface and the rapid multiple ionization prior to substantial structural dynamics. This study demonstrates ways in which fundamental insights into the fragmentation dynamics of polycations following x-ray ionization may be extracted, which will be beneficial to future studies that employ time-resolved x-ray Coulomb explosion imaging to study ultrafast photochemistry.

Allum, Felix [Deutsches Elektronen-Synchrotron (DE↗

NASA NTRS1983

Mass and momentum turbulent transport experiments with confined swirling coaxial jets

Swirling coaxial jets mixing downstream, discharging into an expanded duct was conducted to obtain data for the evaluation and improvement of turbulent transport models currently used in a variety of computational procedures throughout the combustion community. A combination of laser velocimeter (LV) and laser induced fluorescence (LIF) techniques was employed to obtain mean and fluctuating velocity and concentration distributions which were used to derive mass and momentum turbulent transport parameters currently incorporated into various combustor flow models. Flow visualization techniques were also employed to determine qualitatively the time dependent characteristics of the flow and the scale of turbulence. The results of these measurements indicated that the largest momentum turbulent transport was in the r-z plane. Peak momentum turbulent transport rates were approximately the same as those for the nonswirling flow condition. The mass turbulent transport process for swirling flow was complicated. Mixing occurred in several steps of axial and radial mass transport and was coupled with a large radial mean convective flux. Mixing for swirling flow was completed in one-third the length required for nonswirling flow.

Roback, R.↗

NASA NTRS1990

Computation of three-dimensional turbulent boundary layers with heat transfer in a plane of symmetry using embedded wall-layer functions

In the calculation of turbulent boundary layers, a large number of mesh points are required to adequately resolve the intense variation in the velocity and enthalpy in the near-wall region. A substantial reduction in computational effort may be realized by representing the velocity and enthalpy profiles in the wall layer by analytical embedded functions. The effectively inviscid flow in the outer part of the boundary layer may then be resolved by employing a relatively coarser mesh. To obtain complete profiles, the outer numerical solution is matched asymptotically to the inner wall-layer analytical solution. To date, this approach has been restricted to two-dimensional flows; in the present study, a method which may be utilized for turbulent boundary layers with heat transfer in a plane of symmetry is developed as a first step in the application of the embedded-function method to full three-dimensional flows. The present method uses only about half as many mesh points as that required in a conventional procedure, which calculates the flow all the way to the wall, but there is no degradation in accuracy of the computed results.

Degani, A. T.↗

NASA NTRS1993

Moisture and heat budgets of a cirrus cloud from aircraft measurements during FIRE

This study is based on the NCAR King Air aircraft and radiosonde observations on 31 October 1986 during the FIRE in Wisconsin over Oshkosh. Aircraft step-up and spiral descent flights are used to obtain kinematic and thermodynamic data. In the step-up maneuver, six different penetrations were made between 1528 and 1616 UTC. Each penetration was about 30 km long separated in the vertical by about 300 m. The time difference between the two spiral soundings was about 43 min. The aircraft descended at a rate of 1.5 m/s during these spiral soundings. Kinematic, cloud physical, and radiometric observations from various instruments are used to estimate the different terms in the moisture- and heat-budget equations. The results show that the advection terms, estimated using the mean longitudinal wind and vertical velocities, and radiative fluxes are important in forming budgets for the cirrus layers. Ice-crystal growth is significant in the upper layers. The maintenance of cirrus can be attributed to relatively warm and moist air advection, radiative cooling at upper levels, and moisture advection in the vertical. Turbulent heat and moisture fluxes are found to be significant in the low levels of cirrus.

Gultepe, Ismail↗

DOE OSTI2024

Liquid lithium divertor analysis using coupled plasma material interaction model

A liquid lithium divertor can improve performance of future fusion devices by creating efficient power exhaust and improving the energy confinement via pumping of the hydrogen isotopes. In addition, significantly higher heat fluxes can be handled if controlled vapor shielding is used to redistribute the divertor heat flux over a wider area. Design and optimization of such a system calls for an analysis model which includes a strong two-way coupling between the plasma and divertor material. The incoming plasma heat and particle flux will affect the divertor surface temperature, which is a defining factor of the lithium evaporative and sputtered flux going into the plasma. Results of the coupled model based on the plasma edge code SOLPS-ITER and the computational fluid dynamics (CFD) code ANSYS-CFX will be presented for different configurations. An analytical slab flow model is used as a heat transfer boundary condition for SOLPS, defining particle flux from the wall via calculation of the surface temperature. At the final step, results of the SOLPS analysis are verified using a 3D CFD magnetohydrodynamics (MHD) analysis which uses heat and particle flux from SOLPS as a boundary condition. In addition to plasma heat flux, both analytical and CFD temperature models include several plasma material interaction effects, such as lithium evaporation, condensation and sputtering based on deuterium target flux. New adatom sputtering model based on the available experimental data is presented. Analytical model is expanded to include free surface axisymmetric configurations. Results of parametric studies of the divertor configurations with different lithium inlet temperature and velocity will be presented leading to the optimal design resulting in the lowest possible lithium contamination in the core.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1986

The 60 GHz IMPATT diode development

The objective is to develop 60 GHz IMPATT diodes suitable for communications applications. The performance goals of the 60 GHz IMPATT is 1W CW output power with a conversion efficiency of 15 percent and 10-year lifetime. The final design of the 60 GHz IMPATT structure evolved from computer simulations performed at the University of Michigan. The initial doping profile, involving a hybrid double-drift (HDD) design, was derived from a drift-diffusion model that used the static velocity-field characteristics for GaAs. Unfortunately, the model did not consider the effects of velocity undershoot and delay of the avalanche process due to energy relaxation. Consequently, the initial devices were oscillating at a much lower frequency than anticipated. With a revised simulation program that included the two effects given above, a second HDD profile was generated and was used as a basis for fabrication efforts. In the area of device fabrication, significant progress was made in epitaxial growth and characterization, wafer processing, and die assembly. The organo-metallic chemical vapor deposition (OMCVD) was used. Starting with a baseline X-Band IMPATT technology, appropriate processing steps were modified to satisfy the device requirements at V-Band. In terms of efficiency and reliability, the device requirements dictate a reduction in its series resistance and thermal resistance values. Qualitatively, researchers were able to reduce the diodes' series resistance by reducing the thickness of the N+ GaAs substrate used in its fabrication.

Dat, Rovindra↗

NASA NTRS1995

Numerical Simulations of High Enthalpy Pulse Facilities

Axisymmetric flows within shock tubes and expansion tubes are simulated including the effects of finite rate chemistry and both laminar and turbulent boundary layers. The simulations demonstrate the usefulness of computational fluid dynamics for characterizing the flows in high enthalpy pulse facilities. The modeling and numerical requirements necessary to simulate these flows accurately are also discussed. Although there is a large body of analysis which explains and quantifies the boundary layer growth between the shock and the interface in a shock tube, there is a need for more detailed solutions. Phenomena such as thermochemical nonequilibrium. or turbulent transition behind the shock are excluded in the assumptions of Mirels' analysis. Additionally there is inadequate capability to predict the influence of the boundary layer on the expanded gas behind the interface. Quantifying the gas in this region is particularly important in expansion tubes because it is the location of the test gas. Unsteady simulations of the viscous flow in shock tubes are computationally expensive because they must follow features such as a shock wave over the length of the facility and simultaneously resolve the small length scales within the boundary layer. As a result, efficient numerical algorithms are required. The numerical approach of the present work is to solve the axisymmetric gas dynamic equations using an finite-volume formulation where the inviscid fluxes are computed with a upwind TVD scheme. Multiple species equations are included in the formulation so that finite-rate chemistry can be modeled. The simulations cluster grid points at the shock and interface and translate this clustered grid with these features to minimize numerical errors. The solutions are advanced at a CFL number of less than one based on the inviscid gas dynamics. To avoid limitations on the time step due to the viscous terms, these terms are treated implicitly. This requires a block tri-diagonal matrix inversion along each line of cells normal to the wall. The cost of this inversion is more than offset by the larger allowable time step. The source terms representing the finite-rate chemical kinetics are also treated implicitly. An algebraic turbulence model for compressible flow is used. The flow in a low pressure shock tube is computed and the results are compared with Mirels'analysis. The driven gas is nitrogen at 70 Pa, and the incident shock speed is approximately 2.9 km/sec so that there is little dissociation. The simulations include a laminar boundary layer and are run until the limiting flow regime is achieved. At this limit, the shock and interface travel at the same velocity because the amount of driven gas between these two features remains the same: the mass flow across the shock is equal to the mass of gas being entrained at the interface by the boundary layer. Simulations with several grids are presented to establish the grid independence of the solution, Good agreement is achieved between Mirels' correlations and the computations. This is expected since the flow conditions are chosen to be consistent with the assumptions used in Mirels' analysis. This comparison adds credibility to the numerical approach and highlights some of the differences between the theory and the detailed simulations. In addition, simulations of the HYPULSE expansion tube are presented for two operating conditions and the computations are compared to experimental data. The operating gas for both cases is nitrogen. One test condition is at a total enthalpy of 15.2 MJ/Kg and a relatively low pressure of 2 kPa. This case is characterized by a laminar boundary layer and significant chemical nonequilibrium. in the acceleration gas. The second test condition is at a total enthalpy of 10.2 MJ/Kg and a pressure of 38 kPa and is characterized by a turbulent boundary layer. The simulations compare well with experiment and reveal that the nonuniformity in pressure observed during the test time is related to variations in the boundary layer displacement thickness.

Wilson, Gregory J.↗

DOE OSTI2023

Shadow energy functionals and potentials in Born–Oppenheimer molecular dynamics

In Born–Oppenheimer molecular dynamics (BOMD) simulations based on the density functional theory (DFT), the potential energy and the interatomic forces are calculated from an electronic ground state density that is determined by an iterative self-consistent field optimization procedure, which, in practice, never is fully converged. The calculated energies and forces are, therefore, only approximate, which may lead to an unphysical energy drift and instabilities. Here, we discuss an alternative shadow BOMD approach that is based on backward error analysis. Instead of calculating approximate solutions for an underlying exact regular Born–Oppenheimer potential, we do the opposite. Instead, we calculate the exact electron density, energies, and forces, but for an underlying approximate shadow Born–Oppenheimer potential energy surface. In this way, the calculated forces are conservative with respect to the approximate shadow potential and generate accurate molecular trajectories with long-term energy stabilities. We show how such shadow Born–Oppenheimer potentials can be constructed at different levels of accuracy as a function of the integration time step, δt, from the constrained minimization of a sequence of systematically improvable, but approximate, shadow energy density functionals. For each energy functional, there is a corresponding ground state Born–Oppenheimer potential. These pairs of shadow energy functionals and potentials are higher-level generalizations of the original “zeroth-level” shadow energy functionals and potentials used in extended Lagrangian BOMD. The proposed shadow energy functionals and potentials are useful only within this extended dynamical framework, where also the electronic degrees of freedom are propagated as dynamical field variables together with the atomic positions and velocities. The theory is quite general and can be applied to MD simulations using approximate DFT, Hartree–Fock, or semi-empirical methods, as well as to coarse-grained flexible charge models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1995

Sensitivity analysis for aeroacoustic and aeroelastic design of turbomachinery blades

A new method for computing the effect that small changes in the airfoil shape and cascade geometry have on the aeroacoustic and aeroelastic behavior of turbomachinery cascades is presented. The nonlinear unsteady flow is assumed to be composed of a nonlinear steady flow plus a small perturbation unsteady flow that is harmonic in time. First, the full potential equation is used to describe the behavior of the nonlinear mean (steady) flow through a two-dimensional cascade. The small disturbance unsteady flow through the cascade is described by the linearized Euler equations. Using rapid distortion theory, the unsteady velocity is split into a rotational part that contains the vorticity and an irrotational part described by a scalar potential. The unsteady vorticity transport is described analytically in terms of the drift and stream functions computed from the steady flow. Hence, the solution of the linearized Euler equations may be reduced to a single inhomogeneous equation for the unsteady potential. The steady flow and small disturbance unsteady flow equations are discretized using bilinear quadrilateral isoparametric finite elements. The nonlinear mean flow solution and streamline computational grid are computed simultaneously using Newton iteration. At each step of the Newton iteration, LU decomposition is used to solve the resulting set of linear equations. The unsteady flow problem is linear, and is also solved using LU decomposition. Next, a sensitivity analysis is performed to determine the effect small changes in cascade and airfoil geometry have on the mean and unsteady flow fields. The sensitivity analysis makes use of the nominal steady and unsteady flow LU decompositions so that no additional matrices need to be factored. Hence, the present method is computationally very efficient. To demonstrate how the sensitivity analysis may be used to redesign cascades, a compressor is redesigned for improved aeroelastic stability and two different fan exit guide vanes are redesigned for reduced downstream radiated noise. In addition, a framework detailing how the two-dimensional version of the method may be used to redesign three-dimensional geometries is presented.

Lorence, Christopher B.↗

DOE OSTI2023

Large-Eddy Simulation of Flow Over Boeing Gaussian Bump Using Multiagent Reinforcement Learning Wall Model: Preprint

We develop a wall model for large-eddy simulation (LES) that takes into account various pressure-gradient effects using multi-agent reinforcement learning. The model is trained using low-Reynolds-number flow over periodic hills with agents distributed on the wall at various computational grid points. It utilizes a wall eddy-viscosity formulation as the boundary condition to apply the modeled wall shear stress. Each agent receives states based on local instantaneous flow quantities at an off-wall location, computes a reward based on the estimated wall-shear stress, and provides an action to update the wall eddy viscosity at each time step. The trained wall model is validated in wall-modeled LES of flow over periodic hills at higher Reynolds numbers, and the results show the effectiveness of the model on flow with pressure gradients. The analysis of the trained model indicates that the model is capable of distinguishing between the various pressure gradient regimes present in the flow. To further assess the robustness of the developed wall model, simulations of flow over the Boeing Gaussian bump are conducted at a Reynolds number of 2 x 10^6, based on the free-stream velocity and the bump width. The results of mean skin friction and pressure on the bump surface, as well as the velocity statistics of the flow field, are compared to those obtained from equilibrium wall model (EQWM) simulations and published experimental data sets. The developed wall model is found to successfully capture the acceleration and deceleration of the turbulent boundary layer on the bump surface, providing better predictions of skin friction near the bump peak and exhibiting comparable performance to the EQWM with respect to the wall pressure and velocity field. We also conclude that the subgrid-scale model is crucial to the accurate prediction of the flow field, in particular the prediction of separation.

boundary layer↗

DOE OSTI2021

A physics-based model for frost buildup under turbulent flow using direct numerical simulations

We present a new model for frost buildup under turbulent (and laminar) flow using direct numerical simulations. The physical model consists of two layers, the air and the frost. The air layer is fully resolved and consists of solving for the velocity, temperature, and vapor mass fraction fields. The frost layer thickness is resolved using conservation of mass and energy. Both phases are dynamically coupled using the immersed boundary method. Three-dimensional simulations are conducted in an open-channel configuration. A number of challenges need to be overcome to make these simulations feasible. First, to enforce far-field conditions of zero gradient and prescribed mean temperature and humidity, a source term is added to the energy and transport equations in the flow solver. Second, the mean frost thickness is subtracted after each time step to ensure a constant mean flow thickness and level of turbulence in the numerical domain. Third, a slow-time acceleration approach, which accelerates the frost buildup by a predetermined factor, is employed to bridge the gap between the fast turbulent and slow frost buildup time scales. Finally, a frost densification scheme is used to overcome the difficulties of vertically varying frost properties. The model is validated by comparing the frost thickness and frost thickness buildup rate over a period of one hour from a cooled flat plate experiment. As a result, both quantities compare favorably with experiments.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2023

On the design of stable, consistent, and conservative high-order methods for multi-material hydrodynamics

Obtaining stable and high-order numerical solutions for multi-material hydrodynamics is an open challenge. Although slope limiters are widely used to maintain monotonicity near discontinuities, typical limiting procedures violate closure laws at the discrete level when applied to multi-material hydrodynamics equations. Due to this, the high-order expansions of quantities related by the closure laws are no longer consistent. The commonly observed symptom of this consistency-violation is that the numerical method fails to maintain constant pressure and velocity across material interfaces. This leads to sub-optimal convergence rates for smooth multi-material problems as well. Specialized limiting procedures that satisfy consistency while maintaining conservation need to be developed for such equations. A novel procedure that re-instates consistency into slope-limited high-order discretizations applied to the multi-material hydrodynamics equations is presented here. Using simple examples, it is demonstrated that the presented method satisfies closure laws at the discrete level, while maintaining conservative properties of the high-order method. Furthermore, this procedure involves a projection step which relies on the compact basis of the underlying spatial discretization, i.e. for discontinuous schemes (viz. DG and FV) the projection is local, and does not involve global matrix solves. Comparisons with conventional approaches emphasizes the necessity of the consistent closure-law preserving limiting approach, in order to maintain design order of accuracy for smooth multi-material problems.

36 MATERIALS SCIENCE↗

NASA NTRS2005

Conceptual Design of a Communications Relay Satellite for a Lunar Sample Return Mission

In 2003, NASA solicited proposals for a robotic exploration of the lunar surface. Submissions were requested for a lunar sample return mission from the South Pole-Aitken Basin. The basin is of interest because it is thought to contain some of the oldest accessible rocks on the lunar surface. A mission is under study that will land a spacecraft in the basin, collect a sample of rock fragments, and return the sample to Earth. Because the Aitken Basin is on the far side of the Moon, the lander will require a communications relay satellite (CRS) to maintain contact with the Earth during its surface operation. Design of the CRS's orbit is therefore critical. This paper describes a mission design which includes potential transfer and mission orbits, required changes in velocity, orbital parameters, and mission dates. Several different low lunar polar orbits are examined to compare their availability to the lander versus the distance over which they must communicate. In addition, polar orbits are compared to a halo orbit about the Earth-Moon L2 point, which would permit continuous communication at a cost of increased fuel requirements and longer transmission distances. This thesis also examines some general parameters of the spacecraft systems for the mission under study. Mission requirements for the lander dictate the eventual choice of mission orbit. This mission could be the first step in a period of renewed lunar exploration and eventual human landings.

Brunner, Christopher W.↗