Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numeric profiling”

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,315 records · Page 73

Numerical Simulations For the F-16XL Aircraft Configuration

Numerical simulations of flow around the F-16XL are presented as a contribution to the Cranked Arrow Wing Aerodynamic Project International II (CAWAPI-II). The NASA Tetrahedral Unstructured Software System (TetrUSS) is used to perform numerical simulations. This CFD suite, developed and maintained by NASA Langley Research Center, includes an unstructured grid generation program called VGRID, a postprocessor named POSTGRID, and the flow solver USM3D. The CRISP CFD package is utilized to provide error estimates and grid adaption for verification of USM3D results. A subsonic high angle-of-attack case flight condition (FC) 25 is computed and analyzed. Three turbulence models are used in the calculations: the one-equation Spalart-Allmaras (SA), the two-equation shear stress transport (SST) and the kε turbulence models. Computational results, and surface static pressure profiles are presented and compared with flight data. Solution verification is performed using formal grid refinement studies, the solution of Error Transport Equations, and adaptive mesh refinement. The current study shows that the USM3D solver coupled with CRISP CFD can be used in an engineering environment in predicting vortex-flow physics on a complex configuration at flight Reynolds numbers.

Elmiligui, Alaa A.↗

Status of the NASA GMAO Observing System Simulation Experiment

An Observing System Simulation Experiment (OSSE) is a pure modeling study used when actual observations are too expensive or difficult to obtain. OSSEs are valuable tools for determining the potential impact of new observing systems on numerical weather forecasts and for evaluation of data assimilation systems (DAS). An OSSE has been developed at the NASA Global Modeling and Assimilation Office (GMAO, Errico et al 2013). The GMAO OSSE uses a 13-month integration of the European Centre for Medium- Range Weather Forecasts 2005 operational model at T511/L91 resolution for the Nature Run (NR). Synthetic observations have been updated so that they are based on real observations during the summer of 2013. The emulated observation types include AMSU-A, MHS, IASI, AIRS, and HIRS4 radiance data, GPS-RO, and conventional types including aircraft, rawinsonde, profiler, surface, and satellite winds. The synthetic satellite wind observations are colocated with the NR cloud fields, and the rawinsondes are advected during ascent using the NR wind fields. Data counts for the synthetic observations are matched as closely as possible to real data counts, as shown in Figure 2. Errors are added to the synthetic observations to emulate representativeness and instrument errors. The synthetic errors are calibrated so that the statistics of observation innovation and analysis increments in the OSSE are similar to the same statistics for assimilation of real observations, in an iterative method described by Errico et al (2013). The standard deviations of observation minus forecast (xo-H(xb)) are compared for the OSSE and real data in Figure 3. The synthetic errors include both random, uncorrelated errors, and an additional correlated error component for some observational types. Vertically correlated errors are included for conventional sounding data and GPS-RO, and channel correlated errors are introduced to AIRS and IASI (Figure 4). HIRS, AMSU-A, and MHS have a component of horizontally correlated error. The forecast model used by the GMAO OSSE is the Goddard Earth Observing System Model, Version 5 (GEOS-5) with Gridpoint Statistical Interpolation (GSI) DAS. The model version has been updated to v. 5.13.3, corresponding to the current operational model. Forecasts are run on a cube-sphere grid with 180 points along each edge of the cube (approximately 0.5 degree horizontal resolution) with 72 vertical levels. The DAS is cycled at 6-hour intervals, with 240 hour forecasts launched daily at 0000 UTC. Evaluation of the forecasting skill for July and August is currently underway. Prior versions of the GMAO OSSE have been found to have greater forecasting skill than real world forecasts. It is anticipated that similar forecast skill will be found in the updated OSSE.

OSSE↗

Numerical study of nozzle wall cooling for nuclear thermal rockets

The flowfields and performance of nuclear thermal rockets, which utilize radiation and film-cooling to cool the nozzle extension, are studied by solving the Navier-Stokes equations and species equations. The thrust level of the rocket for the present study is about 75,000 lb(f) for a chamber pressure of 68 atm(l,000 psi) and a chamber temperature of 2700 K. The throat radius of the nozzle is 0.0936 m and the area ratios of the nozzles are 300 and 500. It is assumed that the flow is chemically frozen and the turbulence is simulated by the modified Baldwin-Lomax turbulence model. The calculated results for various area ratios and film mass-flow rates are presented as Mach number contours, variations of nozzle wall temperature, exit profiles, and vacuum specific impulses. The present study shows that by selecting the flow rate of the film-cooling hydrogen and area ratio of the nozzle correctly, high area ratio nozzle extensions can be cooled effectively with radiation and film-cooling without significant penalty in performance.

Kim, Suk C.↗

Leading-edge effects on boundary-layer receptivity

Numerical calculations are presented for the incompressible flow over a parabolic cylinder. The computational domain extends from a region upstream of the body downstream to the region where the Blasius boundary-layer solution holds. A steady mean flow solution is computed and the results for the scaled surface vorticity, surface pressure and displacement thickness are compared to previous studies. The unsteady problem is then formulated as a perturbation solution starting with and evolving from the mean flow. The response to irrotational time harmonic pulsation of the free-stream is examined. Results for the initial development of the velocity profile and displacement thickness are presented. These calculations will be extended to later times to investigate the initiation of instability waves within the boundary-layer.

Gatski, Thomas B.↗

Theoretical Analysis of Stationary Potential Flows and Boundary Layers at High Speed

The present report consists of two parts. The first part deals with the two-dimensional stationary flow in the presence of local supersonic zones. A numerical method of integration of the equation of gas dynamics is developed. Proceeding from solutions at great distance from the body the flow pattern is calculated step by step. Accordingly the related body form is obtained at the end of the calculation. The second part treats the relationship between the displacement thickness of laminar and turbulent boundary layers and the pressure distribution at high speeds. The stability of the boundary layer is investigated, resulting in basic differences in the behavior of subsonic and supersonic flows. Lastly, the decisive importance of the boundary layer for the pressure distribution, particularly for thin profiles, is demonstrated.

Oswaititsch, K.↗

Numerical Study of Lander Engine Plume Impingement on the Surface of Europa

Plume exhaust from lander engines would be of concern when landing on the surface of Europa since it could have detrimental effects on both the landing surface and powered descent vehicle. The plume could also entrain particles and redirect them up toward the landing vehicle, as well as erode and contaminate the surface where science would be conducted. In this work, a numerical methodology is developed and validated to conduct a first-order assessment of individual engine plumes of a potential Europa Lander vehicle. Computational Fluid Dynamics (CFD) is used to solve the flow field inside and immediately downstream of the nozzle, while the Direct Simulation Monte Carlo (DSMC) method is applied further downstream where the flow becomes rarefied. An interface between the two domains is defined where macroscopic flow data is passed from the CFD domain to the DSMC domain, establishing a one-way coupling. Plume pressure and velocity fields, as well as ground heating, pressure and density flux profiles, are obtained at altitudes from 25 m down to 10 m, spanning the final stages of landing. The codes and methodologies used in this study are successfully validated with simulations conducted by Morris et al. of the Apollo Lunar Module Descent Engine (LMDE) exhaust plume impinging onto the lunar surface.

Lam, Rebekah↗

Numerical study of lander engine plume impingement on the surface off Europa

Plume exhaust from lander engines would be of concern when landing on the surface of Europa since it could have detrimental effects on both the landing surface and powered descent vehicle. The plume could also entrain particles and redirect them up toward the landing vehicle, as well as erode and contaminate the surface where science would be conducted. In this work, a numerical methodology is developed and validated to conduct a first-order assessment of individual engine plumes of a potential Europa Lander vehicle. Computational Fluid Dynamics (CFD) is used to solve the flow field inside and immediately downstream of the nozzle, while the Direct Simulation Monte Carlo (DSMC) method is applied further downstream where the flow becomes rarefied. An interface between the two domains is defined where macroscopic flow data is passed from the CFD domain to the DSMC domain, establishing a one-way coupling. Plume pressure and velocity fields, as well as ground heating, pressure and density flux profiles, are obtained at altitudes from 25 m down to 10 m, spanning the final stages of landing. The codes and methodologies used in this study are successfully validated with simulations conducted by Morris et al. [1-3] of the Apollo Lunar Module Descent Engine (LMDE) exhaust plume impinging onto the lunar surface.

Hoey, William↗

Global Reference Atmospheric Model (GRAM) Suite Overview and Current Status

• Engineering-oriented atmospheric models that estimate mean values and statistical variations of atmospheric properties for numerous planetary destinations • Currently available for Earth, Mars, Venus, Titan, Neptune, Uranus, and Jupiter • Outputs include atmospheric density, temperature, pressure, chemical composition, radiative fluxes (for Mars-GRAM), and wind components along a user-defined path – Includes seasonal, diurnal, geographic, and altitude variations • Widely used by the engineering community because of their ability to create realistic atmospheric dispersions • Can be integrated into high fidelity flight dynamic simulations of launch, entry, descent and landing (EDL), aerobraking and aerocapture • Optional trajectory input file consisting of time, height, latitude, and longitude can be used to provide the GRAM trajectory path • Optional auxiliary profile consisting of height, latitude, longitude, temperature, pressure, density, eastward wind, and northward wind may be used to replace model data in the GRAMs • Not a forecast model

atmospheric models↗

Solution of the comoving-frame equation of transfer in spherically symmetric flows. IV - Frequency-dependent source functions for scattering by atoms and electrons

A numerical method is presented of solving the radiative transfer equation in the comoving frame of a spherically symmetric expanding atmosphere in which both the line and the electron-scattering source function can depend on frequency (i.e., when there is partial frequency redistribution in the scattering process). This method is used to assess the adequacy of various assumptions regarding frequency redistribution in the comoving frame and to discuss the effects of electron scattering more accurately than previously possible. The methods developed here can be used in realistic model atmospheres to account for the (major) effects of electron scattering upon emergent flux profiles.

Mihalas, D.↗

Experimental studies on two dimensional shock boundary layer interactions

Experiments have been performed on the interaction of oblique shock waves with flat plate boundary layers in the 30.48 cm x 30.48 cm (1 ft. x 1 ft.) supersonic wind tunnel at NASA Lewis Research Center. High accuracy measurements of the plate surface static pressure and shear stress distributions as well as boundary layer velocity profiles were obtained through the interaction region. Documentation was also performed of the tunnel test section flow field and of the two-dimensionality of the interaction regions. The findings provide detailed description of two-dimensional interaction with initially laminar boundary layers over the Mach number range 2.0 to 4.0. Additional information with regard to interactions involving initially transitional boundary layers is presented over the Mach number range 2.0 to 3.0 and those for initially turbulent boundary layers at Mach 2.0. These experiments were directed toward providing well documented information of high accuracy useful as test cases for analytic and numerical calculations. Flow conditions encompassed a Reynolds number range of 4.72E6 to 2.95E7 per meter. The shock boundary layer interaction results were found to be generally in good agreement with the experimental work of previous authors both in terms of direct numerical comparison and in support of correlations establishing laminar separation characteristics.

Skebe, S. A.↗

Hypersonic chemically reacting viscous shock layers over sphere-cones and cylinder-wedges

Hypersonic, nonequilibrium viscous flow over nonanalytic blunt bodies is considered. The equations which govern the viscous shock-layer flow are presented and the method by which the equations are solved is discussed. The predictions of the present finite-difference method are compared with other numerical predictions as well as with experimental data. Three flow conditions are considered; the experimental, wind tunnel conditions of Pappas and Lee for a 7.5 deg sphere-cone at Mach 13 and two cases considered by Kang and Dunn, a 9 deg sphere-cone at 233,000 ft and a 20 deg sphere-cone at 280 and 310 Kft. The predictions of the present method agreed well with the experimental heat-transfer data, but substantial differences were found between the present predictions and the more approximate predictions of Kang and Dunn for heat-transfer distributions and temperature profiles.

Miner, E. W.↗

Methods for computing comet core temperatures

The temperature profile within the comet nucleus provides the key to an understanding of the history of the volatiles within a comet. Certain difficulties arise in connection with current cometary temperature models. It is shown that the constraint of zero net heat flow can be used to derive general analytical expressions which will allow for the determination of comet core temperature for a spherically symmetric comet, taking into account information about the surface temperature and the thermal conductivity. The obtained results are compared with the expression for comet core temperatures considered by Klinger (1981). Attention is given to analytical results, an example case, and numerical models. The formalization developed makes it possible to determine the core temperature on the basis of the numerical models of the surface temperature.

Mckay, C. P.↗

Tethered satellite system

The Tethered Satellite System (TSS) will operate from the Space Shuttle as an earth-orbiting facility that will permit tethered deployment of numerous types of payloads to altitudes both above and below that of the Shuttle. The TSS has evolved from initial studies, beginning in the sixties, to a cooperative development program that is currently being carried out by NASA and the National Space Plan of Italy's CNR. This cooperative agreement between NASA and CNR includes system development and science instrumentation for the first mission, and planning for two additional missions. This paper will (1) discuss predevelopment activities, (2) review the development approach and management relationships between NASA and CNR, (3) describe the current hardware/software configuration and functional interfaces, (4) discuss the mission profile and flight operations planning, and (5) overview science experiment status/plans for the first TSS mission.

Laue, Jay H.↗

Theory of Influence of a Low-Volatility, Soluble Impurity on Spherically-Symmetric Combustion of Fuel Droplets

Analyses are given for the evolution of liquid-phase mass fraction profiles and temperature during combustion of fuel droplets composed of binary miscible mixtures of low-volatility and high-volatility constituents, with small initial mass fractions of the low-volatility material, or slow introduction of low-volatility material into the liquid phase by absorption from the gas phase. The gas phase is assumed to remain quasi-steady and the liquid temperature spatially uniform. The ratio of the liquid-phase diffusion coefficient to the initial burning-rate constant is treated as a small parameter. Asymptotic analyses in this parameter are developed, and the conservation equations are integrated numerically in obtaining descriptions of the combustion history. It is shown that at the surface of the droplet a boundary layer arises in which the mass fraction of the low-volatility component increases with time. The results are used to explain qualitatively some observed conditions of flame contraction and liquid disruption in droplet combustion.

Shaw, B. D.↗

Near-critical spherical accretion by neutron stars - General relativistic treatment

Time-independent spherical accretion by a neutron star is studied using general relativistic radiation hydrodynamics. Numerical integrations of the flow equations are presented. These show that when the luminosity is sufficiently close to (but below) the Eddington limit, the flow velocity increases with decreasing radius far from the neutron star, reaches a maximum at an intermediate radius, and decreases at small radii. A large fraction of the binding energy of the flow is transferred to the radiation through scattering before the flow strikes the surface of the neutron star. Following Miller's treatment of accretion at luminosites near the Eddington limit (which neglected general relativistic effects), analytic approximations for the decelerating phase of the flow's velocity profile are derived. The dependence of the solutions on the variable Eddington factor prescription chosen to close the radiation moment equations is also examined.

Park, Myeong-Gu↗

Passive microwave remote sensing of precipitation from space - Impact of cloud microphysics

Numerical experiments to investigate the impact of time-dependent cloud microphysical structure on the transfer to space of passive microwave radiation are carried out, using a vertically and angularly detailed, plane-parallel radiative transfer model for unpolarized microwave radiation. The object is to explore the feasibility of multichannel passive-microwave retrieval techniques to estimate precipitation from space. Results on the time-dependent influence of large graupel-size ice particles on top-of-atmosphere brightness temperature - rainfall rate relationships at 6.6, 10.7, 18, 19.35, 22.235, 37, 89.5, 130, and 231 GHz are given. Frequency-dependent vertical profiles of a generalized brightness-temperature weighting function, introduced to vertically resolve the contributions by individual cloud and precipitation layers to the brightness temperatures, are also shown.

Mugnai, Alberto↗

Genesis of Pre-Hurricane Felix (2007): Warm Core Formation, Precipitation Evolution, and Predictability - Part 2

This is the second of a two-part study examining the simulated formation of Atlantic Hurricane Felix (2007) in a cloud-representing framework. Here several open issues are addressed concerning the formation of the storm's warm core, the evolution and respective contribution of stratiform versus convective precipitation within the parent wave's pouch, and the sensitivity of the development pathway reported in Part I to different model physics options and initial conditions. All but one of the experiments include ice microphysics as represented by one of several parameterizations, and the partition of convective versus stratiform precipitation is accomplished using a standard numerical technique based on the high-resolution control experiment. The transition to a warm-core tropical cyclone from an initially cold-core, lower tropospheric wave disturbance is analyzed first. As part of this transformation process, it is shown that deep moist convection is sustained near the pouch center. Both convective and stratiform precipitation rates increase with time. While stratiform precipitation occupies a larger area even at the tropical storm stage, deep moist convection makes a comparable contribution to the total rain rate at the pregenesis stage, and a larger contribution than stratiform processes at the storm stage. The convergence profile averaged near the pouch center is found to become dominantly convective with increasing deep moist convective activity there. Low-level convergence forced by interior diabatic heating plays a key role in forming and intensifying the near-surface closed circulation, while the midlevel convergence associated with stratiform precipitation helps to increase the midlevel circulation and thereby contributes to the formation and upward extension of a tropospheric-deep cyclonic vortex. Sensitivity tests with different model physics options and initial conditions demonstrate a similar pregenesis evolution. These tests suggest that the genesis location of a tropical storm is largely controlled by the parent wave's critical layer, whereas the genesis time and intensity of the protovortex depend on the details of the mesoscale organization, which is less predictable. Some implications of the findings are discussed.

Wang, zhuo↗

Roughness Induced Transition in a Supersonic Boundary Layer

Direct numerical simulation is used to investigate the transition induced by threedimensional isolated roughness elements in a supersonic boundary layer at a free stream Mach number of 3.5. Simulations are performed for two different configurations: one is a square planform roughness and the other is a diamond planform roughness. The mean-flow calculations show that the roughness induces counter rotating streamwise vortices downstream of the roughness. These vortices persist for a long distance downstream and lift the low momentum fluid from the near wall region and place it near the outer part of the boundary layer. This forms highly inflectional boundary layer profiles. These observations agree with recent experimental observations. The receptivity calculations showed that the amplitudes of the mass-flux fluctuations near the neutral point for the diamond shape roughness are the same as the amplitude of the acoustic disturbances. They are three times smaller for the square shape roughness.

Balakumar, Ponnampalam↗