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,189 records · Page 66

A real-time guidance algorithm for aerospace plane optimal ascent to low earth orbit

Problems of onboard trajectory optimization and synthesis of suitable guidance laws for ascent to low Earth orbit of an air-breathing, single-stage-to-orbit vehicle are addressed. A multimode propulsion system is assumed which incorporates turbojet, ramjet, Scramjet, and rocket engines. An algorithm for generating fuel-optimal climb profiles is presented. This algorithm results from the application of the minimum principle to a low-order dynamic model that includes angle-of-attack effects and the normal component of thrust. Maximum dynamic pressure and maximum aerodynamic heating rate constraints are considered. Switching conditions are derived which, under appropriate assumptions, govern optimal transition from one propulsion mode to another. A nonlinear transformation technique is employed to derived a feedback controller for tracking the computed trajectory. Numerical results illustrate the nature of the resulting fuel-optimal climb paths.

Calise, A. J.↗

Computerized Generation and Simulation of Meshing and Contact of New Type of Novikov-Wildhaber Helical Gears

In this report, we propose a new geometry for low-noise, increased-strength helical gears of the Novikov-Wildhaber type. Contact stresses are reduced as a result of their convex-concave gear tooth surfaces. The gear tooth surfaces are crowned in the profile direction to localize bearing contact and in the longitudinal direction to obtain a parabolic function of transmission errors. Such a function results in the reduction of noise and vibrations. Methods for the generation of the proposed gear tooth surfaces by grinding and hobbing are considered, and a tooth contact analysis (TCA) computer program to simulate meshing and contact is applied. The report also investigates the influence of misalignment on transmission errors and shift of bearing contact. Numerical examples to illustrate the developed approaches are proposed. The proposed geometry was patented by Ford/UIC (Serial Number 09-340-824, pending) on June 28, 1999.

Litvin, Faydor L.↗

Ultrafast Modulation and Switching of Quantum-Well Lasers using Terahertz Fields

Modulation and switching of semiconductor lasers are important for laser-based information technology. Typically the speed of modulation and switching is limited by interband processes such as stimulated and spontaneous recombinations which occur on a nanosecond time scale. This is why the diode laser modulation has been restricted to tens of GHz. Modulation at higher speed is highly desirable as the information technology enters into the so-called tera-era. In this paper, we study the possibility of utilizing THz-field-induced plasma heating to modulate quantum-well lasers. This is a timely study since, with the advancement of THz solid-state sources and free-electron lasers, THz physics and related technology is currently coming out of its infancy. The investigation of interplaying THz and optical fields is also of intruiging fundamental interest. First, we introduce theoretical plasma heating results for the quantum-well optical amplifier in the presense of an intense half-cycle THz pulse. The heated carrier distributions are then utilized to calculate the THz-pulse-induced change in refractive index and gain profile. Since the electron-hole-plasma is heated using intraband transitions, we circumvent the usual complications due to an overall change in density, and the nonlinear recovery is governed solely by the carrier-LO-phonon interactions, typically 5 ps for a complete recovery. This procedure implies THz and sub-THz switching and recovery rates, respectively; using either gain modulation or index modulation. Plasma heating via steady-state THz fields is also studied. Finally, numerical simulation of a coupled set of equations to investigate the THz modulation based on a simplified model for quantum-well lasers is presented. Our results show that a semiconductor laser can be modulated at up to 1 THz with little distortion with a THz field amplitude at the order of a few kV/cm. Laser responses to a change in THz frequency will be shown. Constraints, practicalities, and applications will be discussed.

Ning, Cun-Zheng↗

Receptivity of Hypersonic Boundary Layers over Straight and Flared Cones

The effects of adverse pressure gradients on the receptivity and stability of hypersonic boundary layers were numerically investigated. Simulations were performed for boundary layer flows over a straight cone and two flared cones. The steady and the unsteady flow fields were obtained by solving the two-dimensional Navier-Stokes equations in axi-symmetric coordinates using the 5th order accurate weighted essentially non-oscillatory (WENO) scheme for space discretization and using third-order total-variation-diminishing (TVD) Runge-Kutta scheme for time integration. The mean boundary layer profiles were analyzed using local stability and non-local parabolized stability equations (PSE) methods. After the most amplified disturbances were identified, two-dimensional plane acoustic waves were introduced at the outer boundary of the computational domain and time accurate simulations were performed. The adverse pressure gradient was found to affect the boundary layer stability in two important ways. Firstly, the frequency of the most amplified second-mode disturbance was increased relative to the zero pressure gradient case. Secondly, the amplification of first- and second-mode disturbances was increased. Although an adverse pressure gradient enhances instability wave growth rates, small nose-tip bluntness was found to delay transition due to the low receptivity coefficient and the resulting weak initial amplitude of the instability waves. The computed and measured amplitude-frequency spectrums in all three cases agree very well in terms of frequency and the shape except for the amplitude.

Balakumar, Ponnampalam↗

Forebody and base region real gas flow in severe planetary entry by a factored implicit numerical method. II - Equilibrium reactive gas

The factored-implicit, gas-dynamic algorithm has been adapted to the numerical simulation of equilibrium reactive flows. Changes required in the perfect gas version of the algorithm are developed, and the method of coupling gas-dynamic and chemistry variables is discussed. A flow-field solution that approximates a Jovian entry case was obtained by this method and compared with the same solution obtained by HYVIS, a computer program much used for the study of planetary entry. Comparison of surface pressure distribution and stagnation line shock-layer profiles indicates that the two solutions agree well.

Davy, W. C.↗

Inferences about the solar wind dynamics from observed distributions of electrons and ions

Recent measurements, mostly from Helios and Voyager, were reviewed. The observational material was compiled into a series of radial profiles for various macroscopic plasma parameters. Several examples for experimentally determined three-dimensional distribution functions of protons are shown. The physical properties of the solar wind electrons were surveyed. The still unresolved question of the driving mechanism of the solar wind is considered. The problem of the energy balance in pure, high velocity streams emanating from solar coronal holes was outlined. It is suggested that the suprathermal electrons may play an important role in driving the wind. Numerical estimates, referring to the lower corona, are offered in support of this hypothesis.

Olbert, S.↗

The growth of the tearing mode - Boundary and scaling effects

A numerical model of resistive magnetic tearing is developed in order to verify and relate the results of the principal approximations used in analytic analyses and to investigate the solutions and their growth-rate scalings over a large range of primary parameters which include parametric values applicable to the solar atmosphere. The computations cover the linear behavior for a variety of boundary conditions, emphasizing effects which differentiate magnetic tearing in astrophysical situations from that in laboratory devices. Eigenfunction profiles for long and short wavelengths are computed and the applicability of the 'constant psi' approximation is investigated. The growth rate is computed for values of the magnetic Reynolds number up to a trillion and of the dimensionless wavelength parameter down to 0.001. The analysis predicts significant effects due to differing values of the magnetic Reynolds number.

Steinolfson, R. S.↗

Error analysis applied to several inversion techniques used for the retrieval of middle atmospheric constituents from limb-scanning MM-wave spectroscopic measurements

The formal retrieval error analysis of Rodgers (1990) allows the quantitative determination of such retrieval properties as measurement error sensitivity, resolution, and inversion bias. This technique was applied to five numerical inversion techniques and two nonlinear iterative techniques used for the retrieval of middle atmospheric constituent concentrations from limb-scanning millimeter-wave spectroscopic measurements. It is found that the iterative methods have better vertical resolution, but are slightly more sensitive to measurement error than constrained matrix methods. The iterative methods converge to the exact solution, whereas two of the matrix methods under consideration have an explicit constraint, the sensitivity of the solution to the a priori profile. Tradeoffs of these retrieval characteristics are presented.

Puliafito, E.↗

Thermophysical Property Measurements of Molten Semiconductors in 1-g and Reduced-g Condition

Understanding and controlling the formation kinetics of varieties of crystal imperfections such as point defects, non uniform distribution of doping atoms, and impurity atoms in growing crystals are very important. Theoretical (numerical) modeling of the crystal growth process is an essential step to achieving these objectives. In order to obtain reliable modeling results, input parameters, i.e. various thermophysical parameters, must be accurate. The importance of accurate thermophysical properties of semiconductors in crystal growth cannot be overly emphasized. The total hemispherical emissivity, for instance, has a dramatic impact on the thermal environment. It determines the radiative emission from the surface of the melt which determines to a large extent the profile of the solidified crystal. In order to understand the convection and the turbulence in a melt, viscosity becomes an important parameter. The liquid surface tension determines the shape of the liquid-atmosphere interface near the solid-liquid-atmosphere triple point. Currently used values for these parameters are rather inaccurate, and this program intends to provide more reliable measurements of these thermophysical properties. Thus, the objective of this program is in the accurate measurements of various thermophysical properties which can be reliably used in the modeling of various crystal growth processes. In this program, thermophysical properties of molten semiconductors, such as Si, Ge, Si-Ge, and InSb will be measured as a function of temperature using the High Temperature Electrostatic Levitator at JPL. Each material will be doped by different kinds of impurities at various doping levels. Thermophysical properties which will be measured include: density, thermal expansion coefficient, surface tension, viscosity, specific heat, hemispherical total emissivity, and perhaps electrical and thermal conductivities. Many molten semiconductors are chemically reactive with crucibles. As a result, these dispersed impurities in the melts tend to substantially modify the properties of pure semiconductors. Sample levitation done in a vacuum clearly helps maintain the sample purity. However, in the 1-g environment, all gravity caused effects such as convection, sedimentation and buoyancy are still present in the sample. In addition, large forces needed to levitate a sample in the presence of the gravity can cause additional flows in the melt. The use of the High Temperature Electrostatic Levitator (HTESL) for the present research is a recent development and little is known about the flows induced by the electrostatic forces. In this ground base program, we will define the limits of HTESL technology as various thermophysical properties of molten semiconductors are measured.

Rhim, Won-Kyu↗

Simulation and Modeling of Hypersonic Turbulent Boundary Layers Subject to Favorable Pressure Gradients due to Streamline Curvature

Direct numerical simulations (DNS) of favorable-pressure-gradient turbulent boundary layers are presented for a nominal freestream Mach number of 5, with the objective of assessing the limitations of the currently available Reynolds-averaged Navier-Stokes (RANS) models. The favorable pressure gradient is induced by the streamwise curvature of the two-dimensional, planar, convex measurement surface used during experiments at the Texas A&M University. The DNS data shows good comparison with the measured velocity profiles, strain rates, and some, but not all, of the Reynolds-stress components. The discrepancies between the predicted and the measured wall-normal as well as shear stress components are primarily attributed to the lower than actual values inferred from typical PIV measurements of turbulent boundary layers. The DNS data shows a zero or slightly negative Reynolds shear stress in the outer part of the boundary-layer, which is indicative of the decaying turbulent motion under a strong favorable pressure gradient. The DNS data is also compared with the results of RANS computations based on commonly used zero, one, and two equation eddy-viscosity models. The RANS models yield reasonable comparisons with the DNS-based skin friction under zero and weak pressure gradients, but significant discrepancies under a strong pressure gradient. The k-w SST model provided the best overall predictions of skin friction, except in the region where the flow transitions from a favorable to an adverse pressure gradient. While the RANS models examined herein also give good predictions of the Reynolds shear stress under a sufficiently weak pressure gradient, none of those models are able to appropriately capture the reduction in the Reynolds stresses when the flow was subjected to a strong pressure gradient. An a priori assessment of the turbulent heat-flux prediction based on the assumption of a constant turbulent Prandtl number with the DNS data shows that while the constant turbulent Prandtl number model is effective in predicting the wall-normal component of turbulent heat flux, it does not capture the turbulent heat transfer in the streamwise direction for all the pressure gradient cases. The failure of the constant turbulent Prandtl number model highlights a requirement for more advanced models of the turbulent heat flux.

turbulent boundary layers↗

Coarse, Intermediate and High Resolution Numerical Simulations of the Transition of a Tropical Wave Critical Layer to a Tropical Storm

Recent work has hypothesized that tropical cyclones in the deep Atlantic and eastern Pacific basins develop from within the cyclonic Kelvin cat's eye of a tropical easterly wave critical layer located equatorward of the easterly jet axis. The cyclonic critical layer is thought to be important to tropical cyclogenesis because its cat's eye provides (i) a region of cyclonic vorticity and weak deformation by the resolved flow, (ii) containment of moisture entrained by the developing flow and/or lofted by deep convection therein, (iii) confinement of mesoscale vortex aggregation, (iv) a predominantly convective type of heating profile, and (v) maintenance or enhancement of the parent wave until the developing proto-vortex becomes a self-sustaining entity and emerges from the wave as a tropical depression. This genesis sequence and the overarching framework for describing how such hybrid wave-vortex structures become tropical depressions/storms is likened to the development of a marsupial infant in its mother's pouch, and for this reason has been dubbed the "marsupial paradigm". Here we conduct the first multi-scale test of the marsupial paradigm in an idealized setting by revisiting the Kurihara and Tuleya problem examining the transformation of an easterly wave-like disturbance into a tropical storm vortex using the WRF model. An analysis of the evolving winds, equivalent potential temperature, and relative vertical vorticity is presented from coarse (28 km), intermediate (9 km) and high resolution (3.1 km) simulations. The results are found to support key elements of the marsupial paradigm by demonstrating the existence of rotationally dominant region with minimal strain/shear deformation near the center of the critical layer pouch that contains strong cyclonic vorticity and high saturation fraction. This localized region within the pouch serves as the "attractor" for an upscale "bottom up" development process while the wave pouch and proto-vortex move together.

Montgomery, M. T.↗

Intermittency and Cessation of Vortex Shedding in the Wake of a Thin Flat Plate with a Circular Trailing Edge

The near and very near wake of thin flat plates with both sharp and circular trailing edges (TEs) are investigated with direct numerical simulations (DNSs). The TE is circular in two of the cases (IN & NS) and sharp in one of them (ST). The separating boundary layers are turbulent in all cases. The objectives of this study are twofold. The first is to explore the effect of significantly reducing Re(sub D) (Reynolds number based on circular TE diameter, D) on the flow in the TE region, and the shedding process (Cases IN and NS). The second is to better understand the reasons underlying the findings of an earlier experimental wake investigation (sharp TE) where (1) the center-line values in normal intensity, and the peak in shear stress profiles in the cross-stream direction, were found to first increase in the streamwise direction (x), from that obtained at the TE, before diminishing further downstream, and (2) a broadband peak was observed in centerline cross-stream velocity (v) spectra (indicating quasi-periodicity, possibly due to vortices or wave-like motions). Case ST from the present study showed a near wake instability resulting in spanwise vortices (with a streamwise component). The instability is intermittent and contributes to both the broadband peak in the v spectrum and the initial increases in normal intensity and shear stress (as in the experiment). Case NS, with the lower value of ReD is an "essentially" non-shedding case where the flow in the TE region continually changes direction (upward/downward) because of turbulence. Case IN, with twice the value of ReD as Case NS, also exhibits a swaying motion in the TE region. In addition, vortex shedding is initiated during periods when the flow direction changes rapidly. Shedding in this case is intermittent. It results in a peak in the v spectrum obtained at the centerline (x/D = 1.0).

direct numerical simulations↗

The Global Ocean Observing System

A Global Ocean Observing System (GOOS) should be established now with international coordination (1) to address issues of global change, (2) to implement operational ENSO forecasts, (3) to provide the data required to apply global ocean circulation models, and (4) to extract the greatest value from the one billion dollar investment over the next ten years in ocean remote sensing by the world's space agencies. The objectives of GOOS will focus on climatic and oceanic predictions, on assessing coastal pollution, and in determining the sustainability of living marine resources and ecosystems. GOOS will be a complete system including satellite observations, in situ observations, numerical modeling of ocean processes, and data exchange and management. A series of practical and economic benefits will be derived from the information generated by GOOS. In addition to the marine science community, these benefits will be realized by the energy industries of the world, and by the world's fisheries. The basic oceanic variables that are required to meet the oceanic and predictability objectives of GOOS include wind velocity over the ocean, sea surface temperature and salinity, oceanic profiles of temperature and salinity, surface current, sea level, the extent and thickness of sea ice, the partial pressure of CO2 in surface waters, and the chlorophyll concentration of surface waters. Ocean circulation models and coupled ocean-atmosphere models can be used to evaluate observing system design, to assimilate diverse data sets from in situ and remotely sensed observations, and ultimately to predict future states of the system. The volume of ocean data will increase enormously over the next decade as new satellite systems are launched and as complementary in situ measuring systems are deployed. These data must be transmitted, quality controlled, exchanged, analyzed, and archived with the best state-of-the-art computational methods.

Kester, Dana↗

Comparison of linear inversion methods by examination of the duality between iterative and inverse matrix methods

Linear numerical inversion methods applied to atmospheric remote sounding generally can be categorized in two ways: (1) iterative, and (2) inverse matrix methods. However, these two categories are not unrelated; a duality exists between them. In other words, given an iterative scheme, a corresponding inverse matrix method exists, and conversely. This duality concept is developed for the more familiar linear methods. The iterative duals are compared with the classical linear iterative approaches and their differences analyzed. The importance of the initial profile in all methods is stressed. Calculations using simulated data are made to compare accuracies and to examine the dependence of the solution on the initial profile.

Fleming, H. E.↗

Computation of viscous flows in turbomachinery cascades

Analytic and numerical experiments involving the use of the space-marching method for the prediction of the viscous flow field in turbomachinery cascades indicate that the grid lines in the physical domain must be nearly orthogonal in order to arrive at accurate predictions. For staggered cascades, a nonperiodic grid system was used and appropriate approximations were applied in place of the periodicity boundary conditions upstream and downstream of the cascade. Predicted drag coefficient, turning angle, boundary layer momentum thicknesses, and velocity profiles are compared with experimental data. Good agreement is obtained in most cases.

Pouagare, M.↗

Numerical simulation of a turbulent flame stabilized behind a rearward-facing step

Flow of combustible mixtures in a plane channel past a smooth contraction followed by an abrupt expansion, in a typical dump combustor configuration, is modeled by a two-dimensional numerical technique based on the random vortex method. Both the inert and the reacting case are considered. In the latter, the flame is treated as an interface, self-advancing at a prescribed normal burning speed, while the dynamic effects of expansion due to the exothermicity of combustion are expressed by volumetric source lines delineated by its front. Solutions are shown to be in satisfactory agreement with experimental results, especially with respect to global properties such as the average velocity profiles and the reattachment length. The stochastic turbulent velocity components manifest interesting differences, especially near the walls where three-dimensional effects of turbulence are expected to be of importance.

Hsiao, C. C.↗

Effect of Blade Cutout on Power Required by Helicopters Operating at High Tip-Speed Ratios

A numerical study was made of the effects of blade cutout on the power required by a sample helicopter rotor traveling at tip-speed ratios of 0.3, 0.4, and 0.5. The amount of cutout varied from 0 to 0.5 of the rotor radius and the calculations were carried out for a thrust coefficient-solidity ratio of 0.04. In these calculations the blade within the cutout radius was assumed to have zero chord. The effect of such cutout on profile-drag power ranged from almost no effect at a tip-speed ratio of 0.3 to as much as a 60 percent reduction at a tip-speed ratio of 0.5. Optimum cutout was about 0.3 of the rotor radius. Part of the large power reduction at a tip-speed ratio of 0.5 resulted from a reduction in tip-region stall, brought about by cutout. For tip-speed ratios greater than 0.3, cutout also effected a significant increase in the ability of the rotor to overcome helicopter parasite drag. It is thus seen that the adverse trends (at high tip-speed ratios) indicated by the uniform-chord theoretical charts are caused in large measure by the center portion of the rotor. The extent to which a modified-design rotor can actually be made more efficient at high speeds than a uniform-chord rotor will depend in practice on the degree of success in minimizing the blade plan form near the center and on special modifications in center-section profiles. A few suggestions and estimates in regard to such modifications are included herein.

Gessow, Alfred↗

External-Compression Supersonic Inlet Design Code

A computer code named SUPIN has been developed to perform aerodynamic design and analysis of external-compression, supersonic inlets. The baseline set of inlets include axisymmetric pitot, two-dimensional single-duct, axisymmetric outward-turning, and two-dimensional bifurcated-duct inlets. The aerodynamic methods are based on low-fidelity analytical and numerical procedures. The geometric methods are based on planar geometry elements. SUPIN has three modes of operation: 1) generate the inlet geometry from a explicit set of geometry information, 2) size and design the inlet geometry and analyze the aerodynamic performance, and 3) compute the aerodynamic performance of a specified inlet geometry. The aerodynamic performance quantities includes inlet flow rates, total pressure recovery, and drag. The geometry output from SUPIN includes inlet dimensions, cross-sectional areas, coordinates of planar profiles, and surface grids suitable for input to grid generators for analysis by computational fluid dynamics (CFD) methods. The input data file for SUPIN and the output file from SUPIN are text (ASCII) files. The surface grid files are output as formatted Plot3D or stereolithography (STL) files. SUPIN executes in batch mode and is available as a Microsoft Windows executable and Fortran95 source code with a makefile for Linux.

Slater, John W.↗