SEARCH · Engineering Papers

Results for “Discrete velocity method”

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 235 records · Page 13

All sources NASA NTRS DOE OSTI

NASA NTRS2010

Adjoint-Based Design of Rotors Using the Navier-Stokes Equations in a Noninertial Reference Frame

Optimization of rotorcraft flowfields using an adjoint method generally requires a time-dependent implementation of the equations. The current study examines an intermediate approach in which a subset of rotor flowfields are cast as steady problems in a noninertial reference frame. This technique permits the use of an existing steady-state adjoint formulation with minor modifications to perform sensitivity analyses. The formulation is valid for isolated rigid rotors in hover or where the freestream velocity is aligned with the axis of rotation. Discrete consistency of the implementation is demonstrated by using comparisons with a complex-variable technique, and a number of single- and multipoint optimizations for the rotorcraft figure of merit function are shown for varying blade collective angles. Design trends are shown to remain consistent as the grid is refined.

Nielsen, Eric J.↗

DOE OSTI2025

Realizability-preserving discontinuous Galerkin method for spectral two-moment radiation transport in special relativity

Here we present a realizability-preserving numerical method for solving a spectral two-moment model to simulate the transport of massless, neutral particles interacting with a steady background material moving with relativistic velocities. The model is obtained as the special relativistic limit of a four-momentum-conservative general relativistic two-moment model. Using a maximum-entropy closure, we solve for the Eulerian-frame energy and momentum. The proposed numerical method is designed to preserve moment realizability, which corresponds to moments defined by a nonnegative phase-space density. The realizability-preserving method is achieved with the following key components: (i) a discontinuous Galerkin phase-space discretization with specially constructed numerical fluxes in the spatial and energy dimensions; (ii) a strong stability-preserving implicit-explicit time-integration method; (iii) a realizability-preserving conserved to primitive moment solver; (iv) a realizability-preserving implicit collision solver; and (v) a realizability-enforcing limiter. Component (iii) is necessitated by the closure procedure, which closes higher order moments nonlinearly in terms of primitive moments. The nonlinear conserved to primitive and the implicit collision solves are formulated as fixed-point problems, which are solved with custom iterative solvers designed to preserve the realizability of each iterate. With a series of numerical tests, we demonstrate the accuracy and robustness of this discontinuous-Galerkin-implicit-explicit method.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1991

Solution of the Navier-Stokes equations for a driven cavity

The flow field in a lid driven cavity is determined by integration of the incompressible Navier-Stokes equations. The numerical integration is accomplished via an operator splitting method known as the theta-scheme. This splitting separates the problem into the solution of a quasi-stokes problem and a nonlinear convection problem. Some details of solution methods used for the two subproblems and results obtained for the driven cavity are described. The schemes developed for the quasi-Stokes problem are more advanced at this stage than those for the nonlinear problem. However, the approaches used for both parts are outlined. As a model problem, a two dimensional square cavity with sides of unit length and a lid moving with unit velocity from left to right is considered. The Navier-Stokes equations are discretized in space on a uniform staggered or MAC mesh. The time discretization is accomplished via the theta-scheme.

Semeraro, B. D.↗

NASA NTRS1991

Computation of rotor aerodynamic loads with a constant vorticity contour free wake model

An analytical method is presented which facilitates the study of isolated rotors with an improved approach to wake simulation. Vortex filaments are simulated along contours of constant sheet strength for the sheet of vorticity resulting from each rotor blade. Curved vortex elements comprise the filaments which can be distorted by the local velocity field. Called the Constant Vorticity Contour wake model, the approach permits the simulation of the blades' wakes corresponding to the full span of the rotor blade. The discretization of the wake of the rotor blade produces spacing and structure that are consistent with the spatial and temporal variations in the loading. A vortex-lattice aerodynamic model of the blade is also included which introduces a finite-element structural model of the blade and consideration of the force and moment trim analysis. Results of the present version of the simulation, called RotorCRAFT, are found to correlate well with H-34 flight-test data.

Quackenbush, Todd R.↗

NASA NTRSDate not supplied

Advancing xEMU Lunar Dust Mitigation Devices

Fine, electrically charged, glass like dust particles caused significant damage to the Apollo EMU9 during lunar EVAs, identified as one of the greatest challenges to future exploration. Passive Lunar Dust Mitigation Devices (LDMD) were developed, within a SBIR Phase II, to prohibit this dust from interrupting venting space suit component operation. A Computational Fluid-Dynamics and Discrete Element Method (CFD-DEM) Simulation Tool was developed at the University of Colorado, Boulder to predict venting gas flow ability to self-clean adhered dust particles from LDMD surfaces. Lunar dust properties (i.e., adhesion and cohesion strengths) required to complete Simulation Tool analysis are relatively unknown due to considerable differences between the Earth and the Moon (i.e., gravity, humidity) and due to an absence of dust particles in their native state. Analytically determining gas flow velocity, density and direction within the fluid Boundary Layer, microns from LDMD surfaces presented a second challenge. Dusty Plasma Chamber testing is being performed at Auburn University, Auburn to observe electrostatically charged dust behavior as it adheres to LDMD prototypes and specific geometry and is then blown away by metered gas flow. Test articles were developed to offer insight into the impact of different flow geometries, surface roughness and dust removal within the gas Boundary-Layer. Observed dust behavior is currently being developed to support the CFD-DEM analysis. Many Simulation Tool analytical cases have been processed to support the intention of completing sensitivity studies to assess how different dust adherence values and Boundary Layer fluid properties impact LDMD self-cleaning effectivity.

Thomas J Stapleton↗

NASA NTRS1986

A computer code for three-dimensional incompressible flows using nonorthogonal body-fitted coordinate systems

In this report, a numerical method for solving the equations of motion of three-dimensional incompressible flows in nonorthogonal body-fitted coordinate (BFC) systems has been developed. The equations of motion are transformed to a generalized curvilinear coordinate system from which the transformed equations are discretized using finite difference approximations in the transformed domain. The hybrid scheme is used to approximate the convection terms in the governing equations. Solutions of the finite difference equations are obtained iteratively by using a pressure-velocity correction algorithm (SIMPLE-C). Numerical examples of two- and three-dimensional, laminar and turbulent flow problems are employed to evaluate the accuracy and efficiency of the present computer code. The user's guide and computer program listing of the present code are also included.

Chen, Y. S.↗

DOE OSTI2026

Onset of Fluidization in MP-PIC Simulations using MFIX-Exa

Fluidized bed reactors are used across a variety of industries, including for energy processes like pyrolysis that result in low-cost energy products. Design and scale-up of fluidized beds is de-risked by modeling and simulation, utilizing tools like NETL’s MFIX-Exa High-Performance Computing (HPC) code for reacting multiphase flow. This report summarizes an investigation into the breadth of problems to which MFIX-Exa may be applied, specifically with regard to low fluid velocities and the onset of fluidization. A simple fluidization study is conducted both experimentally and numerically for particles of interest, then reactor simulations are compared to cold flow experiments for uniform distributor plates. Approaches for modeling bubble caps are also presented.

discrete particle method↗

DOE OSTI2024

Weak boundary conditions for Lagrangian shock hydrodynamics: A high-order finite element implementation on curved boundaries

Here, we propose a new Nitsche-type approach for weak enforcement of normal velocity boundary conditions for a Lagrangian discretization of the compressible shock-hydrodynamics equations using high-order finite elements on curved boundaries. Specifically, the variational formulation is appropriately modified to enforce free-slip wall boundary conditions, without perturbing the structure of the function spaces used to represent the solution, with a considerable simplification with respect to traditional approaches. Total energy is conserved and the resulting mass matrices are constant in time. The robustness and accuracy of the proposed method are validated with an extensive set of tests involving nontrivial curved boundaries.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS2001

Response of the Alliance 1 Proof-of-Concept Airplane Under Gust Loads

This report presents the work performed by Lockheed Martin's Langley Program Office in support of NASA's Environmental Research Aircraft and Sensor Technology (ERAST) program. The primary purpose of this work was to develop and demonstrate a gust analysis method which accounts for the span-wise variation of gust velocity. This is important because these unmanned aircraft having high aspect ratios and low wing loading are very flexible, and fly at low speeds. The main focus of the work was therefore to perform a two-dimensional Power Spectrum Density (PSD) analysis of the Alliance 1 Proof-of-Concept Unmanned Aircraft, As of this writing, none of the aircraft described in this report have been constructed. They are concepts represented by analytical models. The process first involved the development of suitable structural and aeroelastic Finite Element Models (FEM). This was followed by development of a one-dimensional PSD gust analysis, and then the two-dimensional (PSD) analysis of the Alliance 1. For further validation and comparison, two additional analyses were performed. A two-dimensional PSD gust analysis was performed on a simplet MSC/NASTRAN example problem. Finally a one-dimensional discrete gust analysis was performed on Alliance 1. This report describes this process, shows the relevant comparisons between analytical methods, and discusses the physical meanings of the results.

Naser, A. S.↗

DOE OSTI2025

A moment-conserving discontinuous Galerkin representation of the relativistic Maxwellian distribution

Kinetic simulations of relativistic gases and plasmas are critical for understanding diverse astrophysical and terrestrial systems, but the accurate construction of the relativistic Maxwellian, the Maxwell–Jüttner distribution, on a discrete simulation grid is challenging. Difficulties arise from the finite velocity bounds of the domain, which may not capture the entire distribution function, as well as errors introduced by projecting the function onto a discrete grid. Here, we present a novel scheme for iteratively correcting the moments of the projected distribution applicable to all grid-based discretizations of the relativistic kinetic equation. In addition, we describe how to compute the needed nonlinear quantities, such as Lorentz boost factors, in a discontinuous Galerkin scheme through a combination of numerical quadrature and weak operations. The resulting method accurately captures the distribution function and ensures that the moments match the desired values to machine precision.

astrophysical plasmas↗

NASA NTRS1990

Interfacial wave theory of dendrite growth - Global mode solution and quantum condition

The signal feedback process at the edge of the dendrite tip is investigated, and the global instability mechanism of the system is derived. A mechanism is developed to describe a discrete set of unstable global modes for the system. Called WEASR, the method considers the wave emission at the turning point and signal reflections between the turning point and the front edge of the tip. It is shown that the method can obtain the asymptotic solutions for the unstable global modes and the quantum condition for the corresponding eigenvalues. A turning point called the pattern formation condition is shown to be crucial in the formation of dendritic structure and the choice of the final tip velocity. The wave emission process is outlined, and the importance of a signal feedback process at the edge of the dendrite tip is demonstrated. Parameters such as stability and surface tension can be analyzed in terms of their effects on WEASR modes.

Xu, Jian-Jun↗

NASA NTRS2009

Liquid-Propellant Rocket Engine Throttling: A Comprehensive Review

Liquid-Propellant Rocket Engines (LREs) are capable of on-command variable thrust or thrust modulation, an operability advantage that has been studied intermittently since the late 1930s. Throttleable LREs can be used for planetary entry and descent, space rendezvous, orbital maneuvering including orientation and stabilization in space, and hovering and hazard avoidance during planetary landing. Other applications have included control of aircraft rocket engines, limiting of vehicle acceleration or velocity using retrograde rockets, and ballistic missile defense trajectory control. Throttleable LREs can also continuously follow the most economical thrust curve in a given situation, compared to discrete throttling changes over a few select operating points. The effects of variable thrust on the mechanics and dynamics of an LRE as well as difficulties and issues surrounding the throttling process are important aspects of throttling behavior. This review provides a detailed survey of LRE throttling centered around engines from the United States. Several LRE throttling methods are discussed, including high-pressure-drop systems, dual-injector manifolds, gas injection, multiple chambers, pulse modulation, throat throttling, movable injector components, and hydrodynamically dissipative injectors. Concerns and issues surrounding each method are examined, and the advantages and shortcomings compared.

Casiano, Matthew↗

NASA NTRS1971

Development and application of a method for predicting rotor free wake positions and resulting rotor blade air loads. Volume 1: Model and results

Rotor wake geometries are predicted by a process similar to the startup of a rotor in a free stream. An array of discrete trailing and shed vortices is generated with vortex strengths corresponding to stepwise radial and azimuthal blade circulations. The array of shed and trailing vortices is limited to an arbitrary number of azimuthal steps behind each blade. The remainder of the wake model of each blade is an arbitrary number of trailing vortices. Vortex element end points were allowed to be transported by the resultant velocity of the free stream and vortex-induced velocities. Wake geometry, wake flow, and wake-induced velocity influence coefficients are generated by this program for use in the blade loads portion of the calculations. Blade loads computations include the effects of nonuniform inflow due to a free wake, nonlinear airfoil characteristics, and response of flexible blades to the applied loads. Computed wake flows and blade loads are compared with experimentally measured data. Predicted blade loads, response and shears and moments are obtained for a model rotor system having two independent rotors. The effects of advance ratio, vertical separation of rotors, different blade radius ratios, and different azimuthal spacing of the blades of one rotor with respect to the other are investigated.

Sadler, S. G.↗

NASA NTRS2018

Using Adaptive Mesh Refinement to Study Grid Resolution Effects for Shock-Boundary Layer Interactions

Adaptive Mesh Refinement (AMR) promises a much more computationally efficient means to obtain a discrete approximation to a continuous boundary value problem of a specified accuracy than classic isotropic grid refinement. The AMR capability of OVERFLOW (a computational fluid dynamics (CFD) code) is utilized to provide estimates of the exact analytical solutions to problems of interest to turbulence modeling. Predictions of surface pressure and skin friction, essentially the state of stress at the surface, shows little difference with grids believed to be "grid resolved." Velocity profiles, on the other hand, show marked differences in flows with shocks. The AMR method, as implemented in OVERFLOW 2.2k, appears to provide the ability to produce arbitrarily accurate solutions at a predictable cost much smaller than classic uniform mesh refinement.

Turbulence Modeling↗

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.↗

NASA NTRS1985

Coherent large-scale structures in high Reynolds number supersonic jets

The flow structure of a 50.8 mm (2 in) diameter jet operated at a full expanded Mach number of 1.37, with Reynolds numbers in the range 1.7 to 2.35 million, was examined for the first 20 jet diameters. To facilitate the study of the large scale structure, and determine any coherence, a discrete tone acoustic excitation method was used. Phase locked flow visualization as well as laser velocimeter quantitative measurements were made. The main conclusions derived from this study are: (1) large scale coherent like turbulence structures do exist in large Reynolds number supersonic jets, and they prevail even beyond the potential core; (2) the most preferential Strouhal number for these structures is in the vicinity of 0.4; and (3) quantitatively, the peak amplitudes of these structures are rather low, and are about 1% of the jet exit velocity. Finally, since a number of unique problems related to LV measurements in supersonic jets were encountered, a summary of these problems and lessons learned therefrom are also reported.

Lepicovsky, J.↗

NASA NTRS1992

Discrete-vortex simulation of pulsating flow on a turbulent leading-edge separation bubble

Studies are made of the turbulent separation bubble in a two-dimensional semi-infinite blunt plate aligned to a uniform free stream with a pulsating component. The discrete-vortex method is applied to simulate this flow situation because this approach is effective for representing the unsteady motions of the turbulent shear layer and the effect of viscosity near the solid surface. The numerical simulation provides reasonable predictions when compared with the experimental results. A particular frequency with a minimum reattachment is related to the drag reduction. The most effective frequency is dependent on the amplified shedding frequency. The turbulent flow structure is scrutinized. This includes the time-mean and fluctuations of the velocity and the surface pressure, together with correlations between the fluctuating components. A comparison between the pulsating flow and the non-pulsating flow at the particular frequency of the minimum reattachment length of the separation bubble suggests that the large-scale vortical structure is associated with the shedding frequency and the flow instabilities.

Sung, Hyung Jin↗

NASA NTRS1996

Wind Variability of B Supergiants: The Two-component Stellar Wind of gamma Arae - No. 2

The stellar wind of the rapidly rotating early-B supergiant, gamma Ara, is studied using time series, high-resolution IUE spectroscopy secured over approx. 6 days in 1993 March. Results are presented based on an analysis of several line species, including N(N), C(IV), Si(IV), Si(III), C(II), and Al(III). The wind of this star is grossly structured, with evidence for latitude-dependent mass loss which reflects the role of rapid rotation. Independent, co-existing time variable features are identified at low-velocity (redward of approx. -750 km/s) and at higher-speeds extending to approx. -1500 km/s. The interface between these structures is 'defined' by the appearance of a discrete absorption component which is extremely sharp (in velocity space). The central velocity of this 'Super DAC' changes only gradually, over several days, between approx. -400 and -750 km/s in most of the ions. However, its location is shifted redward by almost 400 km/s in Al(III) and C(II), indicating that the physical structure giving rise to this feature has a substantial velocity and ionization jump. Constraints on the relative ionization properties of the wind structures are discussed, together with results based on SEI line-profile-fitting methods. The overall wind activity in gamma Ara exhibits a clear ion dependence, such that low-speed features are promoted in low-ionization species, including Al(III), C(II), and Si(III). We also highlight that - in contrast to most OB stars - there are substantial differences in the epoch-to-epoch time-averaged wind profiles of gamma Ara. We interpret the results in terms of a two-component wind model for gamma Ara, with an equatorially compressed low ionization region, and a high speed, higher-ionization polar outflow. This picture is discussed in the context of the predicted bi-stability mechanism for line-driven winds in rapidly rotating early-B type stars, and the formation of compressed wind regions in rapidly rotating hot stars. The apparent absence of a substantial shift in the wind ionization mixture of gamma Ara, and the normal nature of its photospheric spectrum, suggests wind-compression as the likely dominant cause for the observed equatorial density enhancements.

Prinja, R. K.↗