Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Domain walls”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

A comparison of time domain boundary conditions for acoustic waves in wave guides

Researchers consider several types of boundary conditions in the context of time domain models for acoustic waves. Experiments with four different duct terminations (hard wall, free radiation, foam, and wedge) were carried out in a wave duct from which reflection coefficients over a wide frequency range were measured. These reflection coefficients were used to estimate parameters in the time domain boundary conditions. A comparison of the relative merits of the models in describing the data is presented. Boundary conditions which yield a good fit of the model to the experimental data were found for all duct terminations except the wedge.

Banks, H. T.↗
NASA NTRSDate not supplied

Computational Optimization of a Rotary Valved Pulse Combustor Concept

A resonant pulse combustor valve concept is introduced that utilizes two slotted, coaxial counterrotating discs for mechanical actuation at the combustor inlet. The intended test article for prototype demonstration is a small, 22 in. long, propane fueled laboratory combustor, flowing approximately 0.006 lbm/s of air. The objective is to develop an externally actuated (i.e., active) valve that yields better performance and longer life than the traditional internally actuated (i.e., passive), reed-type valve found on most pulse combustors. The rotary valve motion is optimized using an axisymmetric, two-dimensional computational fluid dynamic simulation with a domain that includes the valve as a moveable interior wall. Parameters such as slew rate, dwell period in the open position, and total closed period are varied using fuel specific impulse as the figure of merit. Variations in fuel injector location and air fuel ratio are also examined. Additionally, the performance impact of leakage from the rotary valve is quantified since leakage is endemic to the design. The optimized simulation results indicate that the rotary valve concept can deliver the desired performance attributes using disc rotational speeds and stresses that are well within the realm of modern materials. A preliminary mechanical valve design is included in the report.

pressure gain combustion↗
NASA NTRSDate not supplied

Computational Optimization of a Rotary Valved Pulse Combustor Concept

A resonant pulse combustor valve concept is introduced that utilizes two slotted, coaxial counterrotating discs for mechanical actuation at the combustor inlet. The intended test article for prototype demonstration is a small, 22 in. long, propane fueled laboratory combustor, flowing approximately 0.006 lbm/s of air. The objective is to develop an externally actuated (i.e., active) valve that yields better performance and longer life than the traditional internally actuated (i.e., passive), reed-type valve found on most pulse combustors. The rotary valve motion is optimized using an axisymmetric, two-dimensional computational fluid dynamic simulation with a domain that includes the valve as a moveable interior wall. Parameters such as slew rate, dwell period in the open position, and total closed period are varied using fuel specific impulse as the figure of merit. Variations in fuel injector location and air fuel ratio are also examined. Additionally, the performance impact of leakage from the rotary valve is quantified since leakage is endemic to the design. The optimized simulation results indicate that the rotary valve concept can deliver the desired performance attributes using disc rotational speeds and stresses that are well within the realm of modern materials. A preliminary mechanical valve design is included in the report.

pressure gain combustion↗
NASA NTRS1986

Acoustic reflection contamination measurements in the 16-foot NASA Langley Transonic Wind Tunnel

The Propfan Test Assessment (PTA) contract, awarded recently to Lockheed by the NASA-Lewis Research Center, required a comprehesive series of near field acoustic measurements. These were to use transducers mounted on the surfaces on a one-ninth scale model of the Gulfstream G-II aircraft, modified to the PTA single propfan testbed configuration. The 16-foot Transonic Wind Tunnel at the NASA Langley Research Center was chosen as the facility to be used for these model acoustic tests. Since the tunnel was hard-walled, it is not clear to what extent the propfan noise signals, reaching the various transducers, might be contaminated by reflections from the walls. A Multidriver Acoustic Source was built, and using an impulse and time domain averaging technique, the reflection contaminations in the wind tunnel were measured at selected microphone locations. Results of the investigations of near field measurements exists, due to reflections from the hard walls of the wind tunnnel.

Burrin, R. H.↗
NASA NTRS1993

Navier-Stokes simulation of external/internal transonic flow on the forebody/inlet of the AV-8B Harrier II

In this work, the computation of combined external/internal transonic flow on the complex forebody/inlet configuration of the AV-8B Harrier II is performed. The actual aircraft has been measured and its surface and surrounding domain, in which the fuselage and inlet have a common wall, have been described using structured grids. The 'thin-layer' Navier-Stokes equations were used to model the flow along with the Chimera embedded multi-block technique. A fully conservative, alternating direction implicit (ADI), approximately factored, partially fluxsplit algorithm was employed to perform the computation. Comparisons to some experimental wind tunnel data yielded good agreement for flow at zero incidence and angle of attack. The aim of this paper is to provide a methodology or computational tool for the numerical solution of complex external/internal flows.

Mysko, Stephen J.↗
NASA NTRS2003

Field Effects of Buoyancy on Lean Premixed Turbulent Flames

The study of field effects of buoyancy on premixed turbulent flames is directed towards the advancement of turbulent combustion theory and the development of cleaner combustion technologies. Turbulent combustion is considered the most important unsolved problem in combustion science and laboratory studies of turbulence flame processes are vital to theoretical development. Although buoyancy is dominant in laboratory flames, most combustion models are not yet capable to consider buoyancy effects. This inconsistency has impeded the validation of theories and numerical simulations with experiments. Conversely, the understanding of buoyancy effects is far too limited to help develop buoyant flame models. Our research is also relevant to combustion technology because lean premixed combustion is a proven method to reduce the formation of oxides of nitrogen (NOx). In industrial lean premixed combustion systems, their operating conditions make them susceptible to buoyancy thus affecting heat distribution, emissions, stability, flashback and blowoff. But little knowledge is available to guide combustion engineers as to how to avoid or overcome these problems. Our hypothesis is that through its influence on the mean pressure field, buoyancy has direct and indirect effects on local flame/turbulence interactions. Although buoyancy acts on the hot products in the farfield the effect is also felt in the nearfield region upstream of the flame. These changes also influence the generation and dissipation of turbulent kinetic energy inside the flame brush and throughout the flowfield. Moreover, the plume of an open flame is unstable and the periodic fluctuations make additional contributions to flame front dynamics in the farfield. Therefore, processes such as flame wrinkling, flow acceleration due to heat release and flame- generated vorticity are all affected. Other global flame properties (e.g. flame stabilization limits and flame speed) may all be coupled to buoyancy. This problem poses major challenges to combustion modeling due to its need for a computation domain extending into the farfield and full specifications of upstream, wall and downstream boundary conditions.

Cheng, R. K.↗
NASA NTRS2017

A New Formulation of Time Domain Boundary Integral Equation for Acoustic Wave Scattering in the Presence of a Uniform Mean Flow

It has been well-known that under the assumption of a constant uniform mean flow, the acoustic wave propagation equation can be formulated as a boundary integral equation, in both the time domain and the frequency domain. Compared with solving partial differential equations, numerical methods based on the boundary integral equation have the advantage of a reduced spatial dimension and, hence, requiring only a surface mesh. However, the constant uniform mean flow assumption, while convenient for formulating the integral equation, does not satisfy the solid wall boundary condition wherever the body surface is not aligned with the uniform mean flow. In this paper, we argue that the proper boundary condition for the acoustic wave should not have its normal velocity be zero everywhere on the solid surfaces, as has been applied in the literature. A careful study of the acoustic energy conservation equation is presented that shows such a boundary condition in fact leads to erroneous source or sink points on solid surfaces not aligned with the mean flow. A new solid wall boundary condition is proposed that conserves the acoustic energy and a new time domain boundary integral equation is derived. In addition to conserving the acoustic energy, another significant advantage of the new equation is that it is considerably simpler than previous formulations. In particular, tangential derivatives of the solution on the solid surfaces are no longer needed in the new formulation, which greatly simplifies numerical implementation. Furthermore, stabilization of the new integral equation by Burton-Miller type reformulation is presented. The stability of the new formulation is studied theoretically as well as numerically by an eigenvalue analysis. Numerical solutions are also presented that demonstrate the stability of the new formulation.

Hu, Fang↗
NASA NTRS2002

Computing Aerodynamic Performance of a 2D Iced Airfoil: Blocking Topology and Grid Generation

The ice accrued on airfoils can have enormously complicated shapes with multiple protruded horns and feathers. In this paper, several blocking topologies are proposed and evaluated on their ability to produce high-quality structured multi-block grid systems. A transition layer grid is introduced to ensure that jaggedness on the ice-surface geometry do not to propagate into the domain. This is important for grid-generation methods based on hyperbolic PDEs (Partial Differential Equations) and algebraic transfinite interpolation. A 'thick' wrap-around grid is introduced to ensure that grid lines clustered next to solid walls do not propagate as streaks of tightly packed grid lines into the interior of the domain along block boundaries. For ice shapes that are not too complicated, a method is presented for generating high-quality single-block grids. To demonstrate the usefulness of the methods developed, grids and CFD solutions were generated for two iced airfoils: the NLF0414 airfoil with and without the 623-ice shape and the B575/767 airfoil with and without the 145m-ice shape. To validate the computations, the computed lift coefficients as a function of angle of attack were compared with available experimental data. The ice shapes and the blocking topologies were prepared by NASA Glenn's SmaggIce software. The grid systems were generated by using a four-boundary method based on Hermite interpolation with controls on clustering, orthogonality next to walls, and C continuity across block boundaries. The flow was modeled by the ensemble-averaged compressible Navier-Stokes equations, closed by the shear-stress transport turbulence model in which the integration is to the wall. All solutions were generated by using the NPARC WIND code.

Chi, X.↗
NASA NTRS1988

Velocity-coupled flow oscillations in a simulated solid-propellant rocket environment

A comprehensive numerical analysis has been carried out to study the unsteady flowfields in a simulated rocket-motor environment. The model is based on the time-dependent compressible Navier-Stokes equations with a two-equation turbulence closure scheme. Various important aspects of the coupling between acoustic oscillations and mean flowfields, including flow reversal, modification of transport properties, etc., are addressed. Results indicate that multi-dimensional effects play important roles in determining local flow structures and wave characteristics. In much of the domain, acoustic velocity nodal points are observed in the near-wall region. The classical one-dimensional theory fails to describe several important mechanisms associated with velocity-induced flow instabilities.

Yang, Vigor↗
NASA NTRS1996

Time-Dependent Parabolic Finite Difference Formulation for Harmonic Sound Propagation in a Two-Dimensional Duct with Flow

An explicit finite difference real time iteration scheme is developed to study harmonic sound propagation in aircraft engine nacelles. To reduce storage requirements for future large 3D problems, the time dependent potential form of the acoustic wave equation is used. To insure that the finite difference scheme is both explicit and stable for a harmonic monochromatic sound field, a parabolic (in time) approximation is introduced to reduce the order of the governing equation. The analysis begins with a harmonic sound source radiating into a quiescent duct. This fully explicit iteration method then calculates stepwise in time to obtain the 'steady state' harmonic solutions of the acoustic field. For stability, applications of conventional impedance boundary conditions requires coupling to explicit hyperbolic difference equations at the boundary. The introduction of the time parameter eliminates the large matrix storage requirements normally associated with frequency domain solutions, and time marching attains the steady-state quickly enough to make the method favorable when compared to frequency domain methods. For validation, this transient-frequency domain method is applied to sound propagation in a 2D hard wall duct with plug flow.

Kreider, Kevin L.↗
NASA NTRS2018

Wall-Resolved Large-Eddy Simulations of Transonic Shock-Induced Flow Separation

Wall-resolved large-eddy simulations of shock-induced flow separation over an axisymmet- ric bump for a Mach number of 0.875 and a chord-based Reynolds number of 2.763 million are performed. The incoming boundary layer has a momentum-thickness Reynolds number of 6600 at 1.5 chords upstream of the bump. The calculations, which employ up to 24 billion grid points, simulate the experimental model of Bachalo and Johnson (AIAA Journal, Vol. 24, No. 3, 1986), except that the tunnel walls are ignored and free air is assumed. The effects of domain span and grid resolution are examined along with the main flowfield features. The predicted shock position as well as separation and reattachment locations agree well with the experiment. Grid convergence is observed in the attached region well upstream of separation. Two-point azimuthal correlations suggest that a span of at least 20 degrees is needed in the attached region, while a 120-degree span might be barely large enough in the separated region. Computed root-mean-square surface pressure fluctuations at the shock foot reach about 5 percent of ambient pressure and the level at reattachment location drops to about half of the primary peak. Simulations reveal evidence of the low-frequency shock unsteadiness; however, a much longer statistical sample is needed to investigate this phenomenon.

Uzun, Ali↗
NASA NTRS1986

Wall reflection of a viscous vortex ring

The behavior of a viscous axisymmetric vortex ring being reflected from a wall is investigated. The incompressible Navier-Stokes equations formulated in terms of the vorticity function and vector potential are numerically integrated by implicit finite difference methods. To specify the vector potential at a far boundary from the wall, the existing integral method used so far only for an unbounded domain is modified by a kind of image method. The trajectory of the vortex ring calcualted as a result closely resembles that observable from the experiment.

Sa, J. Y.↗
NASA NTRS2000

Mineralogy, Petrology, Chemistry, and Ar-39 - Ar-40 and Ages of the Caddo County IAB Iron: Evidence for Early Partial Melt Segregation of a Gabbro Area Rich in Plagioclase-Diopside

We found coarse-grained gabbroic material rich in plagioclase and diopside in the Caddo County IAB iron meteorite. The polished thin sections studied were made from areas rich in Al and Ca detected by a micro-focus X-ray fluorescence (XRF) mapping technique. The gabbro is not a clast within a breccia, but rather this area is located mainly at silicate-metal boundaries only a few cm away from an area with fine-grained, ultramafic silicate similar to winonaites. Medium-grained orthopyroxene and olivine are found in transitional areas showing no disturbance of their crystalline textures. A vein-like region, starting at the area rich in fine-grained mafic silicate, extends towards the gabbroic area with a gradual increase in abundance of plagioclase and diopside. This texture and our accumulated knowledge of the formation mechanism of IAB/winonaltes meteorites, suggest that the gabbroic materials were formed by inhomogeneous segregation of partial melts of chondritic source materials. Compositional data on two mineralogically distinct samples of the gabbro-rich portion of the inclusion were obtained by INAA. Compared to an average of LAB silicate inclusions or winonaites, the Caddo County gabbro is enriched in the incompatible lithophile elements Na, Ca, Sc, REE and Hf, which is consistent with a melt origin for the gabbro. The cosmogenic space exposure age of Caddo County (511 Ma) is significantly younger than exposure ages of some other IAB meteorites, An 39Ar-40Ar age determination of the gabbroic material indicates a series of upward steps in age from 4.516 Ga to 4.523 Ga, with a few high temperature ages up to 4.54 Ga. The older age could approximate the primary recrystallization age of silicates. The stepped Ar age spectrum may indicate differences in Ar closure temperatures during slow cooling of -2-20'C/Myr in the parent body. Alternatively, the younger Ar-Ar ages may date a shock event which occurred while Caddo County was hot and which also created textures such as a linear metal vein with rounded zigzag walls in a diopside crystal and slightly miss-oriented rounded plagioclase domain.

Takeda, Hiroshi↗
NASA NTRS1996

Finite Difference Time Marching in the Frequency Domain: A Parabolic Formulation for the Convective Wave Equation

An explicit finite difference iteration scheme is developed to study harmonic sound propagation in ducts. To reduce storage requirements for large 3D problems, the time dependent potential form of the acoustic wave equation is used. To insure that the finite difference scheme is both explicit and stable, time is introduced into the Fourier transformed (steady-state) acoustic potential field as a parameter. Under a suitable transformation, the time dependent governing equation in frequency space is simplified to yield a parabolic partial differential equation, which is then marched through time to attain the steady-state solution. The input to the system is the amplitude of an incident harmonic sound source entering a quiescent duct at the input boundary, with standard impedance boundary conditions on the duct walls and duct exit. The introduction of the time parameter eliminates the large matrix storage requirements normally associated with frequency domain solutions, and time marching attains the steady-state quickly enough to make the method favorable when compared to frequency domain methods. For validation, this transient-frequency domain method is applied to sound propagation in a 2D hard wall duct with plug flow.

Baumeister, K. J.↗
NASA NTRS1996

Finite Difference Time Marching in the Frequency Domain: A Parabolic Formulation for Aircraft Acoustic Nacelle Design

An explicit finite difference iteration scheme is developed to study harmonic sound propagation in aircraft engine nacelles. To reduce storage requirements for large 3D problems, the time dependent potential form of the acoustic wave equation is used. To insure that the finite difference scheme is both explicit and stable, time is introduced into the Fourier transformed (steady-state) acoustic potential field as a parameter. Under a suitable transformation, the time dependent governing equation in frequency space is simplified to yield a parabolic partial differential equation, which is then marched through time to attain the steady-state solution. The input to the system is the amplitude of an incident harmonic sound source entering a quiescent duct at the input boundary, with standard impedance boundary conditions on the duct walls and duct exit. The introduction of the time parameter eliminates the large matrix storage requirements normally associated with frequency domain solutions, and time marching attains the steady-state quickly enough to make the method favorable when compared to frequency domain methods. For validation, this transient-frequency domain method is applied to sound propagation in a 2D hard wall duct with plug flow.

Baumeister, Kenneth J.↗
NASA NTRS1991

One-way nesting for a primitive equation ocean model

Prognostic numerical models for atmospheric and oceanic circulations require initial fields, boundary conditions, and forcing functions in addition to a consistent set of partial differential equations, including a state relation and equations expressing conservation of mass, momentum, and energy. Depending on the horizontal domain to be modeled, the horizontal boundary conditions are either physically obvious or extremely difficult to specify consistently. If the entire atmosphere is modeled, periodic horizontal boundary conditions are appropriate. On the other hand, the physical horizontal boundaries on the entire ocean are solid walls. Obviously, the normal velocity at a solid wall is zero while the specification of the tangential velocity depends on the mathematical treatment of the horizontal viscous terms. Limitations imposed by computer capacity and cost, as well as research interests, have led to the use of limited area models to study flows in the atmosphere and ocean. The limited area models do not have physical horizontal boundaries, merely numerical ones. Correctly determining these open boundary conditions for limited-area numerical models has both intrigued and frustrated numerical modelers for decades. One common approach is to use the closed or solid wall boundary conditions for a limited-area model. The argument given for this approach is that the boundary conditions affect flow near the walls but that none of these effects are propagated into the interior. Therefore, one chooses a big enough domain that the central region of interest is not corrupted by the boundary flow. Research in progress to model the North Atlantic circulation vividly illustrates the pitfalls of this approach. Two model runs are compared: (1) the southern boundary at 20S between latitudes 0 and 40W is artificially closed; and (2) the same boundary is specified as open with an inward transport of 15 Sv (determined from a global model with the same physics) uniformly spread across the boundary. A comparison of both runs is presented.

Blake, D. W.↗
NASA NTRS1993

Flow resolution and domain of influence in rarefied hypersonic blunt-body flows

The study assesses the effects of the upstream domain size and grid resolution upon flow properties and body aerodynamics computed for rarefied flows over cold blunt bodies with a direct simulation Monte Carlo (DSMC) particle method. Empirical correlations are suggested for aerodynamic coefficients for two-dimensional flows past a perpendicular flat plate. Free-stream parameters which were varied in the study include the Mach number, Knudsen number, surface temperature, and intermolecular potential. Insufficient grid resolution leads to overprediction of aerodynamic heating and forces in the DSMC method. Solution accuracy correlates well with the Reynolds number defined at the wall temperature and the stagnation mean free path relative to the cell dimension. Insufficient upstream domain size in the DSMC method leads to overprediction of heating and drag. Errors in aerodynamic coefficients correlate well with the distance ahead of the body where flow temperature reaches half of its peak value. Simulation of a hard-sphere gas is more sensitive to grid resolution, while simulation of a Maxwell gas is more sensitive to upstream domain size.

Haas, Brian L.↗