Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Direct numerical simulation”

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 1,549 records · Page 86
All sourcesNASA NTRSDOE OSTI
NASA NTRS2001

A Hybrid Numerical Method for Turbulent Mixing Layers

A hybrid method has been developed for simulations of compressible turbulent mixing layers. Such mixing layers dominate the flows in exhaust systems of modern day aircraft and also those of hypersonic vehicles currently under development. The method configurations in which a dominant structural feature provides an unsteady mechanism to drive the turbulent development in the mixing layer. The hybrid method uses a Reynolds-averaged Navier-Stokes (RANS) procedure to calculate wall bounded regions entering a mixing section, and a Large Eddy Simulation (LES) procedure to calculate the mixing dominated regions. A numerical technique was developed to enable the use of the hybrid RANS-LES method on stretched, non-Cartesian grids. Closure for the RANS equations was obtained using the Cebeci-Smith algebraic turbulence model in conjunction with the wall-function approach of Ota and Goldberg. The wall-function approach enabled a continuous computational grid from the RANS regions to the LES region. The LES equations were closed using the Smagorinsky subgrid scale model. The hybrid RANS-LES method is applied to a benchmark compressible mixing layer experiment. Preliminary two dimensional calculations are used to investigate the effects of axial grid density and boundary conditions. Vortex shedding from the base region of a splitter plate separating the upstream flows was observed to eventually transition to turbulence. The location of the transition, however, was much further downstream than indicated by experiments. Actual LES calculations, performed in three spatial directions, also indicated vortex shedding, but the transition to turbulence was found to occur much closer to the beginning of the mixing section. which is in agreement with experimental observations. These calculations demonstrated that LES simulations must be performed in three dimensions. Comparisons of time-averaged axial velocities and turbulence intensities indicated reasonable agreement with experimental data.

Georgiadis, Nicholas J.↗
NASA NTRS1991

Wake interference behind two flat plates normal to the flow - A finite-element study

A finite-element model of the Navier-Stokes equations is used for numerical simulation of flow past two normal flat plates arranged side by side at Reynolds number 80 and 160. The results from this simulation indicate that when the gap between the plates is twice the width of a single plate, the individual wakes of the plates behave independently, with the antiphase vortex shedding being dominant. At smaller gap sizes, the in-phase vortex shedding, with strong wake interaction, is favored. The gap flow in those cases becomes biased, with one of the wakes engulfing the other. The direction of the biased flow was found to be switching at irregular intervals, with the time histories of the indicative flow parameters and their power spectra resembling those of a chaotic system.

Behr, M.↗
NASA NTRS2011

Analysis of the Radiative Transfer Equation with Highly Asymmetric Phase Function

This paper considers a scalar radiative transfer problem with high scattering anisotropy, Two computational methods are presented based on decomposition of the diffuse light field into a regular and anisotropic part. The first algorithm (DOMAS) singles out the anisotropic radiance in the forward scattering peak using the Small-Angle Modification of RTE. The second algorithm (DOM2+) separates the single scattering radiance as an anisotropic part, which largely defines the fine detail of the total radiance in the backscattering directions. In both cases, the anisotropic part is represented analytically. With anisotropy subtraction, the regular part of the signal. which requires a numerical solution, is essentially smoothed as a function of angles. Further, the transport equation is obtained for the regular part that contains an additional source function from the anisotropic part of the signal. This equation is solved with the discrete ordinates method. A conducted numerical analysis of this work showed that algorithm DOMAS has a strong advantage as compared to the standard discrete ordinates method for simulation of the radiance transmission, and DOM2 + is the best of the three for the reflection computations. Both algorithms offer at least a factor of three acceleration of convergence of the azimuthal series for highly anisotropic phase functions.

Korkin, Sergey V.↗
NASA NTRS1993

A new flux conserving Newton's method scheme for the two-dimensional, steady Navier-Stokes equations

A new numerical method is developed for the solution of the two-dimensional, steady Navier-Stokes equations. The method that is presented differs in significant ways from the established numerical methods for solving the Navier-Stokes equations. The major differences are described. First, the focus of the present method is on satisfying flux conservation in an integral formulation, rather than on simulating conservation laws in their differential form. Second, the present approach provides a unified treatment of the dependent variables and their unknown derivatives. All are treated as unknowns together to be solved for through simulating local and global flux conservation. Third, fluxes are balanced at cell interfaces without the use of interpolation or flux limiters. Fourth, flux conservation is achieved through the use of discrete regions known as conservation elements and solution elements. These elements are not the same as the standard control volumes used in the finite volume method. Fifth, the discrete approximation obtained on each solution element is a functional solution of both the integral and differential form of the Navier-Stokes equations. Finally, the method that is presented is a highly localized approach in which the coupling to nearby cells is only in one direction for each spatial coordinate, and involves only the immediately adjacent cells. A general third-order formulation for the steady, compressible Navier-Stokes equations is presented, and then a Newton's method scheme is developed for the solution of incompressible, low Reynolds number channel flow. It is shown that the Jacobian matrix is nearly block diagonal if the nonlinear system of discrete equations is arranged approximately and a proper pivoting strategy is used. Numerical results are presented for Reynolds numbers of 100, 1000, and 2000. Finally, it is shown that the present scheme can resolve the developing channel flow boundary layer using as few as six to ten cells per channel width, depending on the Reynolds number.

Scott, James R.↗
NASA NTRS1992

Vorticity generation by contoured wall injectors

A class of contoured wall fuel injectors was designed to enable shock-enhancement of hypervelocity mixing for supersonic combustion ramjet applications. Previous studies of these geometries left unresolved questions concerning the relative importance of various axial vorticity sources in mixing the injectant with the freestream. The present study is a numerical simulation of two generic fuel injectors which is aimed at elucidating the relative roles of axial vorticity sources including: baroclinic torque through shock-impingement, cross-stream shear, turning of boundary layer vorticity, shock curvature, and diffusive flux. Both the magnitude of the circulation, and the location of vorticity with respect to the mixing interface were considered. Baroclinic torque and cross-stream shear were found to be most important in convectively mixing the injectant with the freestream, with the former providing for deposition of vorticity directly on the fuel/air interface.

Waitz, Ian A.↗
NASA NTRS1993

An evaluation of excitation techniques for time domain based flutter data processing

Numerically simulated, wind tunnel and flight test experimental data were used to evaluate various types of system excitation sources for subcritical flutter testing. Emphasis was placed upon the determination of modal frequencies and damping from relatively short time series records. Parameter identification based upon digital time series models was used to predict modal characteristics from subcritical test data. Various types of excitation were considered. These were: 1) impulsive loading and the resulting transient free vibration response, 2) random forced response, in which the 'mechanical' forcing function could be directly measured, 3) sine sweep forced response, in which the 'mechanical' forcing function could be directly measured, and 4) forced response from unsteady aerodynamic excitation or turbulence, in which the excitation could only be indirectly inferred from other measurements.

Pinkelman, J. K.↗
NASA NTRS1993

A local dynamic model for large eddy simulation

The dynamic model is a method for computing the coefficient C in Smagorinsky's model for the subgrid-scale stress tensor as a function of position from the information already contained in the resolved velocity field rather than treating it as an adjustable parameter. A variational formulation of the dynamic model is described that removes the inconsistency associated with taking C out of the filtering operation. This model, however, is still unstable due to the negative eddy-viscosity. Next, three models are presented that are mathematically consistent as well as numerically stable. The first two are applicable to homogeneous flows and flows with at least one homogeneous direction, respectively, and are, in fact, a rigorous derivation of the ad hoc expressions used by previous authors. The third model in this set can be applied to arbitrary flows, and it is stable because the C it predicts is always positive. Finally, a model involving the subgrid-scale kinetic energy is presented which attempts to model backscatter. This last model has some desirable theoretical features. However, even though it gives results in LES that are qualitatively correct, it is outperformed by the simpler constrained variational models. It is suggested that one of the constrained variational models should be used for actual LES while theoretical investigation of the kinetic energy approach should be continued in an effort to improve its predictive power and to understand more about backscatter.

Ghosal, Sandip↗
NASA NTRSDate not supplied

Extending mARC II Arc-Jet Test Duration via Design and Implementation of Model System Cooling Sleeve

The mARC II is a 30 kW arc-jet facility at NASA Ames Research Center used to generate high-enthalpy flows for low-cost thermal protection system (TPS) technology development. Sustained operation of downstream instrumentation and material samples is constrained by thermal loading transmitted through the arc-jet test environment, limiting achievable run times and experimental throughput. This work presents the design, integration, and validation of a cooling sleeve implemented on the sweep arm drive motor feedthrough to mitigate thermal accumulation during testing. The addition of the cooling sleeve is a simple, robust upgrade that translates directly into enhanced facility capability by supporting longer run durations, reduced turnaround time, and higher throughput.

sting arm↗
NASA NTRS1988

NAS: The first year

Discussed are the capabilities of NASA's Numerical Aerodynamic Simulation (NAS) Program and its application as an advanced supercomputing system for computational fluid dynamics (CFD) research. First, the paper describes the NAS computational system, called the NAS Processing System Network, and the advanced computational capabilities it offers as a consequence of carrying out the NAS pathfinder objective. Second, it presents examples of pioneering CFD research accomplished during NAS's first operational year. Examples are included which illustrate CFD applications for predicting fluid phenomena, complementing and supplementing experimentation, and aiding in design. Finally, pacing elements and future directions for CFD and NAS are discussed.

Bailey, F. R.↗
NASA NTRS1988

NAS - The first year

Discussed are the capabilities of NASA's Numerical Aerodynamic Simulation (NAS) Program and its application as an advanced supercomputing system for computational fluid dynamics (CFD) research. First, the paper describes the NAS computational system, called the NAS Processing System Network, and the advanced computational capabilities it offers as a consequence of carrying out the NAS Pathfinder objective. Second, it presents examples of pioneering CFD research accomplished during NAS's first operational year. Examples are included which illustrate CFD applications for predicting fluid phenomena, complementing and supplementing experimentation, and aiding in design. Finally, pacing elements and future directions for CFD and NAS are discussed.

Bailey, F. R.↗
NASA NTRS1990

Conditions at the downstream boundary for simulations of viscous incompressible flow

The proper specification of boundary conditions at artificial boundaries for the simulation of time-dependent fluid flows has long been a matter of controversy. A general theory of asymptotic boundary conditions for dissipative waves is applied to the design of simple, accurate conditions at downstream boundary for incompressible flows. For Reynolds numbers far enough below the critical value for linear stability, a scaling is introduced which greatly simplifies the construction of the asymptotic conditions. Numerical experiments with the nonlinear dynamics of vortical disturbances to plane Poiseuille flow are presented which illustrate the accuracy of our approach. The consequences of directly applying the scalings to the equations are also considered.

Hagstrom, Thomas↗
NASA NTRS1989

Aeroassist flight experiment aerodynamics and aerothermodynamics

The problem is to determine the transitional flow aerodynamics and aerothermodynamics, including the base flow characteristics, of the Aeroassist Flight Experiment (AFE). The justification for the computational fluid dynamic (CFD) Application stems from MSFC's system integration responsibility for the AFE. To insure that the AFE objectives are met, MSFC must understand the limitations and uncertainties of the design data. Perhaps the only method capable of handling the complex physics of the rarefied high energy AFE trajectory is Bird's Direct Simulation Monte Carlo (DSMC) technique. The 3-D code used in this analysis is applicable only to the AFE geometry. It uses the Variable Hard Sphere (VHS) collision model and five specie chemistry model available from Langley Research Center. The code is benchmarked against the AFE flight data and used as an Aeroassisted Space Transfer Vehicle (ASTV) design tool. The code is being used to understand the AFE flow field and verify or modify existing design data. Continued application to lower altitudes is testing the capability of the Numerical Aerodynamic Simulation Facility (NASF) to handle 3-D DSMC and its practicality as an ASTV/AFE design tool.

Brewer, Edwin B.↗
NASA NTRS2007

Characteristics of a Linearly Tapered Slot Antenna (LTSA) Conformed Longitudinally Around a Cylinder

The family of tapered slot antennas (TSA s) is suitable for numerous applications. Their ease of fabrication, wide bandwidth, and high gain make them desirable for military and commercial systems. Fabrication on thin, flexible substrates allows the TSA to be conformed over a given body, such as an aircraft wing or a piece of clothing for wearable networks. Previously, a Double Exponentially Tapered Slot Antenna (DETSA) was conformed around an exponential curvature, which showed that the main beam skewed towards the direction of curvature. This paper presents a Linearly Tapered Slot Antenna (LTSA) conformed longitudinally around a cylinder. Measured and simulated radiation patterns and the direction of maximum H co-polarization (Hco) as a function of the cylinder radius are presented.

Jordan, Jennifer L.↗
NASA NTRS2007

Energetic Neutral Atom Emissions From Venus: VEX Observations and Theoretical Modeling

Venus has almost no intrinsic magnetic field to shield itself from its surrounding environment. The solar wind thus directly interacts with the planetary ionosphere and atmosphere. One of the by-products of this close encounter is the production of energetic neutral atom (ENA) emissions. Theoretical studies have shown that significant amount of ENAs are emanated from the planet. The launch of the Venus Express (VEX) in 2005 provided the first light ever of the Venus ENA emissions. The observed ENA flux level and structure are in pretty good agreement with the theoretical studies. In this paper, we present VEX ENA data and the comparison with numerical simulations. We seek to understand the solar wind interaction with the planet and the impacts on its atmospheres.

Fok, M.-C.↗
NASA NTRS2008

Marshall Space Flight Center's Impact Testing Facility Capabilities

Marshall Space Flight Center's (MSFC) Impact Testing Facility (ITF) serves as an important installation for space and missile related materials science research. The ITF was established and began its research in spacecraft debris shielding in the early 1960% then played a major role in the International Space Station debris shield development. As NASA became more interested in launch debris and in-flight impact concerns, the ITF grew to include research in a variety of impact genres. Collaborative partnerships with the DoD led to a wider range of impact capabilities being relocated to MSFC as a result of the closure of Particle Impact Facilities in Santa Barbara, California. The Particle Impact Facility had a 30 year history in providing evaluations of aerospace materials and components during flights through rain, ice, and solid particle environments at subsonic through hypersonic velocities. The facility's unique capabilities were deemed a "National Asset" by the DoD. The ITF now has capabilities including environmental, ballistic, and hypervelocity impact testing utilizing an array of air, powder, and two-stage light gas guns to accommodate a variety of projectile and target types and sizes. Relocated test equipment was dated and in need of upgrade. Numerous upgrades including new instrumentation, triggering circuitry, high speed photography, and optimized sabot designs have been implemented. Other recent research has included rain drop demise characterization tests to obtain data for inclusion in on-going model development. Future ITF improvements will be focused on continued instrumentation and performance enhancements. These enhancements will allow further, more in-depth, characterization of rain drop demise characterization and evaluation of ice crystal impact. Performance enhancements also include increasing the upper velocity limit of the current environmental guns to allow direct environmental simulation for missile components. The current and proposed ITF capabilities range from rain to micrometeoroids allowing the widest test parameter range possible for materials investigations in support of space, atmospheric, and ground environments. These test capabilities including hydrometeor, single/multi-particle, ballistic gas grins, exploding wire gun, and light gas guns combined with Smooth Particle Hydrodynamics Code (SPHC) simulations represent the widest range of impact test capabilities in the country.

Evans, Steve↗
NASA NTRS2008

Marshall Space Flight Center's Impact Testing Facility Capabilities

Marshall Space Flight Center's (MSFC) Impact Testing Facility (ITF) serves as an important installation for space and missile related materials science research. The ITF was established and began its research in spacecraft debris shielding in the early 1960s, then played a major role in the International Space Station debris shield development. As NASA became more interested in launch debris and in-flight impact concerns, the ITF grew to include research in a variety of impact genres. Collaborative partnerships with the DoD led to a wider range of impact capabilities being relocated to MSFC as a result of the closure of Particle Impact Facilities in Santa Barbara, California, The Particle Impact Facility had a 30 year history in providing evaluations of aerospace materials and components during flights through rain, ice, and solid particle environments at subsonic through hypersonic velocities. The facility's unique capabilities were deemed a 'National Asset' by the DoD, The ITF now has capabilities including environmental, ballistic, and hypervelocity impact testing utilizing an array of air, powder, and two-stage light gas guns to accommodate a variety of projectile and target types and sizes. Relocated test equipment was dated and in need of upgrade. Numerous upgrades including new instrumentation, triggering circuitry, high speed photography, and optimized sabot designs have been implemented. Other recent research has included rain drop demise characterization tests to obtain data for inclusion in on-going model development. Future ITF improvements will be focused on continued instrumentation and performance enhancements. These enhancements will allow further, more in-depth, characterization of rain drop demise characterization and evaluation of ice crystal impact. Performance enhancements also include increasing the upper velocity limit of the current environmental guns to allow direct environmental simulation for missile components. The current and proposed ITF capabilities range from rain to micrometeoroids allowing the widest test parameter range possible for materials investigations in support of space, atmospheric, and ground environments. These test capabilities including hydrometeor, single/multi-particle, ballistic gas guns, exploding wire gun, and light gas guns combined with Smooth Particle Hydrodynamics Code (SPHC) simulations represent the widest range of impact test capabilities in the country.

Evans, Steve↗
NASA NTRS2007

Post-Newtonian Black-Hole Inspiral Initial Data with Waves for Puncture Simulations

We present post-Newtonian-inspired initial data for non-spinning black-hole binaries, suitable for numerical evolution with punctures. We revisit the work of Tichey et al. and explicitly calculate remaining integral terms.Thereby we improve the accuracy in the far zone, by including realistic gravitational waves in the initial data. We investigate the behavior of this data both at the center of mass and in the far zone, demonstrating agreement of the transverse-traceless parts of the new metric with quadrupole-approximation waveforms. An advantage of these data is that they can be used for numerical evolutions to make a direct connection between the merger waveforms and the post-Newtonian inspiral waveforms.

Kelly, B. J.↗
NASA NTRS1999

Numerical Modeling and In-Situ Observations of the Dynamics of the Solid/Liquid Interface Morphology During Directional Solidification of Alloys

The departure from interface planarity and the subsequent evolution to a periodic array of cells or dendrites is a fundamental process that characterizes most microstructures in solidified alloys. The growing demand for high quality alloys and semiconductor crystals requires a precise methodology to predict and subsequently control both the interface morphology and the distribution of impurities, additives, and phases in the grown crystal. Apart from its practical significance, the study of morphological evolution has also been viewed as a means to unearth a general paradigm for pattern formation in nature. A previously developed 2D numerical model for the solid/liquid interface tracking has been further refined and used to simulate the time-evolution of the perturbations on the interface. The dynamics of the local growth velocity, interface undercooling and solute concentration at the interface has been theoretically predicted by means of the numerical model for Al-Cu and Pb-Sn alloys. The model shows that perturbations with a wavelengths, lambda greater than a critical wavelength lambda(sub c) continue to grow in time whereas perturbations with lambda < lambda(sub c) cease to propagate. The model further predicts that under certain conditions perturbation can also propagate along the interface. Comparison of these predictions with existing theories of pattern formation and experimental results will be discussed.

Catalina, Adrian V.↗