Search NASA⌕ Search

SEARCH · Search NASA

Results for “curvature estimation”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS2025

PSE-Based Aerodynamic Flow Transition Prediction Using Automated Unstructured CFD Integration

Accurate, robust, and efficient prediction of transition in viscous flows is a significant challenge in computational fluid dynamics. We present a coupled, high-fidelity, iterative approach that leverages the FUN3D flow solver and the LASTRAC stability code to predict transition in low-disturbance environments, initiated by the linear growth of boundary-layer instability modes. Our method integrates the ability of FUN3D to compute mixed laminar-transitional-turbulent mean flows via transition-sensitized Reynolds-Averaged Navier-Stokes equations with the ability of LASTRAC to perform linear stability analysis, all within an automated framework that requires no intermediate user involvement. Unlike conventional frameworks that rely on classical stability theory or reduced-order metamodels, our approach employs the parabolized stability equations to provide more accurate and reliable estimates of disturbance growth for multiple instability mechanisms, including Tollmien-Schlichting, Kelvin-Helmholtz, and crossflow modes. By accounting for the effects of mean-flow nonparallelism as well as the surface curvature, this approach lays the foundation for improved N-factor correlations for transition onset prediction in a broad class of flows. We apply this method to three distinct flow configurations: 1) flow over a zero-pressure-gradient flat plate, 2) the NLF-0416 airfoil with both natural and separation-induced transitions, and 3) a 6:1 prolate spheroid, where transition is primarily driven by crossflow instability. For the two-dimensional cases, a formulated intermittency distribution is used to model the transition zone in between the laminar and fully turbulent flows. The results include comparisons with experimental measurements, similar numerical approaches, and transport-equations-based models, demonstrating good agreement in surface pressure coefficients, transition onset locations, and skin-friction coefficients for all three configurations. Besides contributing a couple of new insights into boundary-layer transition in these canonical cases, this study provides a powerful tool for transition modeling in both research and design applications in aerodynamics.

Transition Prediction↗
NASA NTRS2005

Estimating Aeroheating of a 3D Body Using a 2D Flow Solver

A method for rapidly estimating the aeroheating, shear stress, and other properties of hypersonic flow about a three-dimensional (3D) blunt body has been devised. First, the geometry of the body is specified in Cartesian coordinates. The surface of the body is then described by its derivatives, coordinates, and principal curvatures. Next, previously relatively simple equations are used to find, for each desired combination of angle of attack and meridional angle, a scaling factor and the shape of an equivalent axisymmetric body. These factors and equivalent shapes are entered as inputs into a previously developed computer program that solves the two-dimensional (2D) equations of flow in a non-equilibrium viscous shock layer (VSL) about an axisymmetric body. The coordinates in the output of the VSL code are transformed back to the Cartesian coordinates of the 3D body, so that computed flow quantities can be registered with locations in the 3D flow field of interest. In tests in which the 3D bodies were elliptic paraboloids, the estimates obtained by use of this method were found to agree well with solutions of 3D, finite-rate-chemistry, thin-VSL equations for a catalytic body.

Scott, Carl D.↗
NASA NTRS1979

Applications of velocity potential function to acoustic duct propagation and radiation from inlets using finite element theory

A finite element velocity potential program has been developed to study acoustic wave propagation in complex geometries. For irrotational flows, relatively low sound frequencies, and plane wave input, the finite element solutions show significant effects of inlet curvature and flow gradients on the attenuation of a given acoustic liner in a realistic variable area turbofan inlet. In addition, as shown in the paper, the velocity potential approach can not be used to estimate the effects of rotational flow on acoustic propagation since the potential acoustic disturbances propagate at the speed of the media in sheared flow. Approaches are discussed that are being considered for extending the finite element solution to include the far field as well as the internal portion of the duct. A new matrix partitioning approach is presented that can be incorporated in previously developed programs to allow the finite element calculation to be marched into the far field. The partitioning approach provides a large reduction in computer storage and running times.

Baumeister, K. J.↗
NASA NTRS1982

Pulsed multiwavelength laser ranging system

A pulsed multiwavelength laser ranging system for measuring atmospheric delay was built and tested, and its theoretical performance limits were calculated. The system uses a dye modelocked ND:YAG laser, which transmits 70 psec wide pulses simultaneously at 1064, 532, and 355 nm. The differential delay of the 1064 and 355 nm pulses is measured by a specially calibrated waveform digitizer to estimate the dry atmospheric delay. The delay time of the 532 nm pulse is used to measure the target distance. Static crossed field photomultipliers are used as detectors for all wavelengths. Theoretical analysis shows that path curvature and atmospheric turbulence are fundamental limits to the ranging accuracy of both single and multicolor systems operating over horizontal paths. For two color systems, an additional error is caused by the uncertainty in the path averaged water vapor. The standard deviation of the multicolor instrument's timing measurements is directly proportional to the laser pulse width plus photomultiplier jitter divided by the square root of the received photoelectron number. The prototype system's maximum range is km, which is limited by atmospheric and system transmission losses at 355 nm. System signal detection and false alarm calculations are also presented.

Abshire, J. B.↗
NASA NTRS2012

Coronal Magnetic Field Measurement Using CME-Driven Shock Observations

Collisionless shocks form ahead of coronal mass ejections (CMEs) when the CME speed exceeds the Alfven speed of the ambient plasma in the corona and interplanetary medium. The shock stands at a distance from the CME flux rope that depends on the shock Mach number, the geometry of the driver, and the adiabatic index. While the shock ahead of the CME has been observed for a long time in the in situ data, it has been identified recently near the Sun in the coronagraphic and EUV images. Unlike in situ observations, the imaging observations are two dimensional, so one can better discern the CME-shock relationship near the Sun. Gopalswamy and Yashiro demonstrated that the coronal magnetic field can be derived from the shock standoff distance measured in coronagraphic images. The method involves measuring the standoff distance, the radius of curvature of the flux rope, and assuming the value of the adiabatic index and deriving the Alfvenic Mach number. The next step is to derive the Alfvenic Mach number from the measured shock speed and an estimate of the local solar wind speed. The final step involves deriving the magnetic field from the Alfven speed by measuring the local plasma density either from coronagraphic (polarized brightness) images or from the band-splitting of type II radio bursts. In this paper, we derive the combined magnetic field profile from near the Sun to the edge of the LASCO field of view (1.5 to 30 solar radii) and compare it with the current model profiles.

Gopalswarmy, Nat↗
NASA NTRS1986

A flare-induced cascade model of gamma-ray bursts

An analytical model is developed for the source of gamma ray bursts as a stellar flare in the magnetosphere of a neutron star. It is suggested that the loss of energy through synchrotron radiation experienced by electrons moving through a sufficiently strong magnetic field at a large pitch angle may not be regained. Instead, pulsar theory is applied to show that the acceleration of electrons in an electric field parallel to the magnetic field will rapidly be inhibited by curvature radiation as the loop experiences a reconnection. It is shown that electrons passing through a curvature with a radius of one million with an electric field strength of 10 billion e.s.u. will emit photons with energies of up to 10 to the 12.6 eV by curvature radiation. The photons, gamma rays, would annihilate in the magnetosphere, which they cannot escape. The resulting cascade of electron-positron particles would eventually produce photons of sufficiently low energy to escape. Upper and low bounds are estimated for the resulting emission spectrum, which would vary according to the magnetic field geometry. The model explains the observed 511 keV annihilation line and the optical radiation which at times accompanies gamma-ray bursts.

Sturrock, P. A.↗
NASA NTRS2020

Conditions of Loss Cone Filling by Scattering on the Curved Field Lines for 30 keV Protons During Geomagnetic Storm as Inferred From Numerical Trajectory Tracing

The rate of pitch angle scattering on the curved magnetic field lines is well parameterized by the ratio of the minimum field line curvature radius to the maximum effective particle gyroradius (K = RC/rg). The critical value of this ratio (Kcr) corresponding to the loss cone filling is of special interest since it corresponds to the low altitude isotropic boundaries (IBs). The early theoretical estimates gave Kcr = 8, whereas recent estimations of the K parameter on the field lines corresponding to the observed IBs during the geomagnetic storms revealed KIB values in the range of 3–30. We numerically trace the trajectories of the 30 keV protons in the magnetic field of the global magnetohydrodynamic simulation of the intense storm in order to infer statistical distribution of Kcr. The electric field and effects of nonstationarity are neglected in this study. It is found that although the Kcr values do show some variations during the course of the storm, its range is rather narrow 4 < Kcr < 8. The result suggests that higher KIB values found in the observational studies, if not caused by the magnetosphere-ionosphere mapping error, should be attributed to some other mechanism of pitch angle scattering. The Kcr values tend to be lower (4–6) during the main phase because the region of low K values approaches the Earth and the equatorial loss cone size becomes larger due to a larger equatorial magnetic field in the near earth region. The remaining variation of Kcr is explained by the presence of the guide component of the magnetic field.

S Dubyagin↗
NASA NTRS1993

Inferring global characteristics of current sheet from local measurements

A new procedure to extract global characteristics of a current sheet based on local measurements of magnetic field, plasma, and current density from a single spacecraft is described. By adopting suitable current density profiles, this method permits an estimate of the following parameters: (1) the asymptotic magnetic field strength outside the current sheet, (2) the north-south thickness of the current sheet, (3) the spacecraft distance from the center of the current sheet, (4) the magnetic field line curvature at the center of the current sheet, (5) the volume current density at the center of the current sheet, (6) the integrated current density across the thickness of the current sheet, and (7) the kappa parameter used by many researchers in the study of particle orbits in a current sheet geometry.

Lui, A. T. Y.↗
NASA NTRS2013

A Hybrid Fiber/Solid-State Regenerative Amplifier with Tunable Pulse Widths for Satellite Laser Ranging

A fiber/solid-state hybrid seeded regenerative amplifier, capable of achieving high output energy with tunable pulse widths, has been developed for satellite laser ranging applications. The regenerative amplifier cavity uses a pair of Nd:YAG zigzag slabs oriented orthogonally to one another in order to make thermal lensing effects symmetrical and simplify optical correction schemes. The seed laser used is a fiber-coupled 1,064-nm narrowband (<0.02 nm) diode laser that is discretely driven in a new short-pulsed mode, enabling continuously tunable seed pulse widths in the 0.2-to-0.4-ns range. The amplifier gain unit consists of a pair of Brewster-cut 6-bounce zigzag Nd:YAG laser slabs, oriented 90deg relative to each other in the amplifier head. This arrangement creates a net-symmetrical thermal lens effect (an opposing singleaxis effect in each slab), and makes thermo-optical corrections simple by optimizing the curvature of the nearest cavity mirror. Each slab is pumped by a single 120-W, pulsed 808-nm laser diode array. In this configuration, the average pump beam distribution in the slabs had a 1-D Gaussian shape, which matches the estimated cavity mode size. A half-wave plate between the slabs reduces losses from Fresnel reflections due to the orthogonal slabs Brewster-cut end faces. Successful "temporal" seeding of the regenerative amplifier cavity results in a cavity Q-switch pulse envelope segmenting into shorter pulses, each having the width of the input seed, and having a uniform temporal separation corresponding to the cavity round-trip time of approx. =10 ns. The pulse energy is allowed to build on successive passes in the regenerative amplifier cavity until a maximum is reached, (when cavity gains and losses are equal), after which the pulse is electro- optically switched out on the next round trip The overall gain of the amplifier is approx. =82 dB (or a factor of 1.26 million). After directing the amplified output through a LBO frequency doubling crystal, approx. = 2.1 W of 532-nm output (>1 mJ) was measured. This corresponds to a nonlinear conversion efficiency of >60%. Furthermore, by pulse pumping this system, a single pulse per laser shot can be created for the SLR (satellite laser ranging) measurement, and this can be ejected into the instrument. This is operated at the precise frequency needed by the measurement, as opposed to commercial short-pulsed, mode-locked systems that need to operate in a continuous fashion, or CW (continuous wave), and create pulses at many MHz. Therefore, this design does not need to throw away or dump 99% of the laser energy to produce what is required; this system can be far smaller, more efficient, cheaper, and readily deployed in the field when packaged efficiently. Finally, by producing custom diode seed pulses electronically, two major advantages over commercial systems are realized: First, this pulse shape is customizable and not affected by the cavity length or gain of the amplifier cavity, and second, it can produce adjustable (selectable) pulse widths by simply adding multiple seed diodes and coupling each into commercial, low-cost fiber-optic combiners.

Coyle, Barry↗
NASA NTRS1973

Reflection and refraction of hydromagnetic waves at the magnetopause

Reflection and transmission coefficients of MHD waves are obtained at a stable, plane interface which separates two compressible, perfectly conducting media in relative motion to each other. The coefficients are evaluated for representative conditions of the quiet-time, near-earth magnetopause. The transmission coefficient averaged over a hemispherical distribution of incident waves is found to be 1-2%. Yet the magnitude of the energy flux deposited into the magnetosphere in a day averaged over a hemispherical distribution of waves having amplitudes of say 2-3 gamma, is estimated to be of the order 10 to the 22-nd power erg. Therefore the energy input of MHD waves must contribute significantly to the energy budget of the magnetosphere. The assumption that the boundary surface is a tangential discontinuity with no curvature limits the present theory to hydromagnetic frequencies higher than about .1 Hz.

Verzariu, P.↗
NASA NTRS1992

Experimental measurement and theoretical modeling of microwave scattering and the structure of the sea surface influencing radar observations from space

The electromagnetic (EM) bias 'epsilon' is an error present in radar altimetry of the ocean due to the nonuniform reflection from wave troughs and crests. The EM bias is defined as the difference between the mean reflecting surface and the mean sea surface. A knowledge of the EM bias is necessary to permit error reduction in mean sea level measurements by satellite radar altimeters. Direct measurements of the EM bias were made from a Shell Offshore oil production platform in the Gulf of Mexico for a six month period during 1989 and 1990. Measurements of the EM bias were made at 5 and 14 Ghz. During the EM bias experiments by Melville et al., a wire wave gauge was used to obtain the modulation of the high frequency waves by the low frequency waves. It became apparent that the EM bias was primarily caused by the modulation of the short waves. This was reported by Arnold et al. The EM bias is explained using physical optics scattering and an empirical model for the short wave modulation. Measurements of the short wave modulation using a wire wave gauge demonstrated a linear dependence of the normalized bias on the short wave modulation strength, M. The theory accurately predicts this dependence by the relation epsilon = -alphaMH sub 1/3. The wind speed dependence of the normalized bias is explained by the dependence of the short wave modulation strength on the wind speed. While other effects such as long wave tilt and curvature will have an effect on the bias, the primary cause of the bias is shown to be due to the short wave modulation. This report will present a theory using physical optics scattering and an empirical model of the short wave modulation to estimate the EM bias. The estimated EM bias will be compared to measurements at C and Ku bands.

Arnold, David↗
NASA NTRS1977

Influence of transverse surface waves on turbulent boundary layers

Wavy wall experiments using solid waves and progressive waves have been reported. For this paper, the major effects of waviness of the wall on the flow are identified as due to oscillatory curvature (convex-concavity) and oscillatory acceleration/deceleration of the flow, which imposes a highly nonequilibrium influence upon the turbulence structure. The theoretical analysis in this presentation takes into account proper turbulence modeling (including the nonequilibrium effects) for the wavy wall problem. The analysis proceeds in three stages: (1) inviscid solution for induced pressure due to the physical wall, (2) solution of a turbulent boundary layer with pressure gradients and curvature effects in the modeling from which the profile correction is computed, and (3) induced pressure computations for the corrected profile. The phase shift of pressure perturbations with respect to the physical wall can be predicted, and pressure drag and skin friction drag can be estimated, with nonlinear viscous effects included. Comparison of the theoretical estimates with experimental data are also presented.

Balasubramanian, R.↗
NASA NTRS1986

Viscous pumping and the spin-down of thermospheric gyres and jets

Strong gyres and jets can be generated at auroral latitudes in the thermosphere by enhanced electric fields during geomagnetic substorms. Typical height profiles of ion density suggest that the ion drag force should generate large curvature in the vertical profile of the winds in the highly viscous region of the thermosphere above about 200 km. It is proposed that the poststorm spin-down of these gyres and jets proceeds via Ekman circulations driven by the curvatures in the height profiles of the winds. Analytic and numerical calculations of the ageostrophic winds forced by curvature in model geostrophic wind profiles show that the ageostrophic wind speeds and directions depend mainly on the kinematic viscosity in the region of curvature and the total change in shear in the geostrophic wind. Ageostrophic wind speeds for typical thermospheric jets can exceed 200 m/s (about 50 percent of the jet winds). Spin-down times of thermospheric jets and cyclonic gyres by the Ekman pumping mechanism are estimated at less than about 6 hours.

Walterscheid, R. L.↗
NASA NTRS1987

Latitudinal and radial gradients of anomalous and galactic cosmic rays in the outer heliosphere

Measurements from instruments on Voyagers 1 and 2 and Pioneer 10 are used to derive simultaneous radial and latitudinal gradients of anomalous cosmic-ray oxygen during the latter part of 1985. It is found that in the energy interval from 7.1-30.6 MeV/nucleon the latitudinal gradient is -3 to -4 percent/deg. The sign of the latitudinal gradient is opposite to that reported for the last solar cycle when the solar magnetic field polarity was reversed, as predicted by propagation models in which curvature and gradient drifts are significant. The ratios of the radial and latitudinal gradients of anomalous oxygen, anomalous helium, and galactic cosmic rays are similar, also in agreement with drift theory predictions. Based on drift theory, these observed ratios, which are near unity, lead to estimates of the perpendicular diffusion coefficient/particle velocity ratio in the range 3-8 x 10 to the 22nd sq cm/sec at 24.9 AU for particles with rigidities of about 2-3 GV.

Cummings, A. C.↗
NASA NTRS1990

Magnetospheric interaction with Triton's ionosphere

The large electron densities measured by the Voyager radio occultation experiment are attributed to the precipitation of magnetospheric electrons with energy above 10 keV. Because the ionospheric electric Pedersen conductivity of Triton is about 10,000-20,000 mho and the Alfven conductance is about 3.5 mho, direct convective flow of plasma into the essentially infinitely conducting ionosphere is negligible. Magnetospheric electrons are transported to Triton's ionopause by curvature drift as a result of weak magnetic field line draping in a sub-Alfvenic plasma interaction with Triton. At the ionopause energetic electrons have a high probability of elastic and inelastic scattering and precipitate into the upper atmosphere. The average power dissipation is estimated to be about (2 - 3) x 10 to the 8th W.

Strobel, Darrell F.↗
NASA NTRS1993

On the Accuracy of CFD-Based Pressure Drop Predictions for Right-Angle Ducts

The predictive capability of computational fluid dynamics (CFD) codes for turbulent flow through curved ducts is of significant importance to the design and performance analysis of modern rocket engine flowpaths. Code calibration and validation studies for this class of flow are desireable to estimate the performance margin and operating range of components designed using Navier-Stokes methods. Parametric experimental studies such as that of Weske (NACA ARR W-39) provided a wealth of performance data for the design of single- and compound elbow configurations with various cross-sections, curvature and aspect ratios at varying Reynolds numbers. In that work, the majority of data is presented in the form of loss coefficients, characterizing pressure losses due to duct curvature, and including losses due to wall friction. Using measured friction coefficients, losses of equivalent straight lengths of duct are subtracted, resulting in performance curves useful for design computations. These data are currently used in a CFD-based parametric study covering a broad range of operating conditions. Of particular interest for the accuracy of CFD predictions are the effects on pressure loss due to inlet boundary layer thickness (dependent on upstream development length), and the wall treatment for the turbulence equations (conventional wall functions vs. wall integration using a two-layer model). The experimental data are reassessed in the form of an error analysis, and are compared with CFD predictions for 18 computational cases. Grid-independence, grid spacing, and convergence requirements of the cases are discussed. Conclusions regarding the relative importance of the parametric variables will be presented.

Brankovic, Andreja↗
NASA NTRS1979

Design and calibration of slotted walls for transonic airfoil wind tunnels

The traditional procedure for estimating the performance of slotted walls for airfoil wind tunnels is reviewed, and a modification which improves the accuracy of this procedure is described. Unlike the traditional procedure, the modified procedure indicates that the design of airfoil wind-tunnel walls which induce minimal blockage and streamline-curvature effects is feasible. The design and testing of such a slotted wall is described. It is shown experimentally that the presence of a model can affect the plenum pressure and thus make the use of the plenum pressure as a calibration reference questionable. Finally, an ONERA experiment which shows the effect of the sidewall boundary layer on the measured model normal force is discussed.

Barnwell, R. W.↗
NASA NTRS1979

Fundamental mechanisms that influence the estimate of heat transfer to gas turbine blades

Heat tranfer problems in aircraft gas turbines required for improved prediction of turbine blade or vane gas-side heat transfer are examined. Estimates of the heat transfer from the gas to vanes or rotating blades are uncertain due to the complexity of the heat transfer processes, since the gas flow is three dimensional with complex secondary viscous flow patterns that interact with the endwalls and blade surfaces. In addition, upstream disturbances, stagnation flow, curvature effects, and flow acceleration complicate the thermal transport mechanisms in the boundary layers. The thermal state and flow characteristics of the hot gases that enter the turbine blade row, analytical methods for calculating the gas-side heat transfer to turbine blades, and flow phenomena such as stagnation, curvature effects, acceleration, secondary flows, and transition that influences local heat transfer rates are discussed.

Graham, R. W.↗