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,711 records · Page 95

Volumetric pattern analysis of fuselage-mounted airborne antennas

A volumetric pattern analysis of fuselage-mounted airborne antennas at high frequencies was investigated. The primary goal of the investigation was to develop a numerical solution for predicting radiation patterns of airborne antennas in an accurate and efficient manner. An analytical study of airborne antenna pattern problems is presented in which the antenna is mounted on the fuselage near the top or bottom. Since this is a study of general-type commercial aircraft, the aircraft was modeled in its most basic form. The fuselage was assumed to be an infinitely long perfectly conducting elliptic cylinder in its cross-section and a composite elliptic cylinder in its elevation profile. The wing, cockpit, stabilizers (horizontal and vertical) and landing gear are modeled by "N" sided bent or flat plates which can be arbitrarily attached to the fuselage. The volumetric solution developed utilizes two elliptic cylinders, namely, the roll plane and elevation plane models to approximate the principal surface profile (longitudinal and transverse) at the antenna location. With the belt concept and the aid of appropriate coordinate system transformations the solution can be used to predict the volumetric patterns of airborne antennas in an accurate and efficient manner. Applications of this solution to various airborne antenna problems show good agreement with scale model measurements. Extensive data are presented for a microwave landing antenna system.

Yu, C. L.↗

Radiative acceleration in outflows from broad absorption line quasi-stellar objects. 2: Wind models

We investigate the dynamics of radiatively driven broad absorption-line (BAL) outflows in quasi-stellar objects (QSOs) by developing radial and time-independent numerical models. Two limits are explored. The first assumes that the absorbing matter is not forced to comove with the substrate, which provides pressure confinement. This assumption allows us to explore in detail a case in which the acceleration is entirely due to radiation pressure. Using the parameters inferred from observations, we find that under these conditions radiative acceleration (mainly due to resonance line scattering) can readily accelerate the flow to the observed velocities. An important feature of the noncoupled flow is that the line profiles tend to stay relatively flat throughout the velocity interval covered by the line. We discuss how relaxing the assumptions of radial symmetry and time independent may help to explain the structures observed in BALs. In the second class of models, the absorbing flow is assumed to be completely coupled to the substrate in which it is embedded. Aside from being more plausible physically, these models produce line profiles that trail off at higher velocities, a behavior observed in some BALs. We show that, even if the substrate is massless, we have to assume a starting radius very close to the inferred radius of the broad emission-line region (approximately 0.1 pc) in order to obtain a significant contribution from radiative acceleration, given a typical active galactic nucleus (AGN) spectrum. The reason is that the energy input needed to pressurize the substrate, allowing the flow to become supersonic and to retain a reasonable ionization equilibrium, at the same time contributes appreciably to the acceleration. A way to relax the small starting radius constraint is to use a softer ionizing spectrum.

Arav, Nahum↗

1DVAR Analysis of Temperature and Humidity using GPS Radio Occulation Data

The Global Positioning System (GPS) enables positioning with a very small receiver. The signals transmitted by the GPS satellites are sensitive to the atmosphere and can be used to perform soundings with the radio occultation technique (e.g., Kursinski et al., 1997). The GPS signal can be converted to refractivity N via the Abel transform. The refractivity can then be related to atmospheric pressure P, temperature T, and water vapor partial pressure P (sub w) the GPS measurement, (between 0.5 and 1.5 km), its self-calibration, and it's nearly all-weather capabilities make it a good candidate for use in data assimilation systems (DAS) and numerical weather prediction (NWP). In order to demonstrate its usefulness in a DAS or NWP system, a first step is to assess its impact oil the analysis. A one-dimensional variational off-line analysis (1DVAR), meaning the data are not assimilated 'In the 3D DAS, constitutes a starting approach to which further enhancements can be made. The chosen observable to be analyzed in this study is the refractivity. One way to extract temperature (humidity) from the refractivity, is to assume a humidity (temperature) profile. One variable may then be retrieved without any a priori information. The 1DVAR approach used here resolves the ambiguity problem raised in the interpretation of these data. It enables retrieving these two atmospheric variables at a reasonable computing cost.

Poli, Paul↗

Model of Atmospheric Links on Optical Communications from High Altitude

Optical communication links have the potential to solve many of the problems of current radio and microwave links to satellites and high-altitude aircraft. The higher frequency involved in optical systems allows for significantly greater signal bandwidth, and thus information transfer rate, in excess of 10 Gbps, and the highly directional nature of laser-based signals eliminates the need for frequency-division multiplexing seen in radio and microwave links today. The atmosphere, however, distorts an optical signal differently than a microwave signal. While the ionosphere is one of the most significant sources of noise and distortion in a microwave or radio signal, the lower atmosphere affects an optical signal more significantly. Refractive index fluctuations, primarily caused by changes in atmospheric temperature and density, distort the incoming signal in both deterministic and nondeterministic ways. Additionally, suspended particles, such as those in haze or rain, further corrupt the transmitted signal. To model many of the atmospheric effects on the propagating beam, we use simulations based on the beam-propagation method. This method, developed both for simulation of signals in waveguides and propagation in atmospheric turbulence, separates the propagation into a diffraction and refraction problem. The diffraction step is an exact solution, within the limits of numerical precision, to the problem of propagation in free space, and the refraction step models the refractive index variances over a segment of the propagation path. By applying refraction for a segment of the propagation path, then diffracting over that same segment, this method forms a good approximation to true propagation through the atmospheric medium. Iterating over small segments of the total propagation path gives a good approximation to the problem of propagation over the entire path. Parameters in this model, such as initial beam profile and atmospheric constants, are easily modified in a simulation such as this, which allows for the rapid analysis of different propagation scenarios. Therefore, this method allows the development of a near-optimal system design for a wide range of situations, typical of what would be seen in different atmospheric conditions over a receiving ground station. A simulation framework based upon this model was developed in FORTRAN, and for moderate grid sizes and propagation distances these simulations are computable in reasonable time on a standard workstation. This presentation will discuss results thus far.

Subich, Christopher↗

Numerical Study of the Turbulent Prandtl Number in Supersonic Plane-Channel Flow – the Effect of Thermal Boundary Conditions

Direct numerical simulation is used to investigate the turbulent Prandtl number Prt above cold (isothermal) and hot (adiabatic) walls in a family of supersonic channel flows. A range of Reynolds numbers and mean temperature/density variations is considered. The value of Prt away from the wall approaches 0.85 above both the isothermal and adiabatic walls, with the width of the Prt = 0.85 region increasing with wall-variable-based Reynolds number Reτw (the latter strongly affected by the thermal boundary condition). The variable, near-wall Prt profiles from both the present and previous, passive-scalar simulations collapse as a function of the semilocal y∗ wall scaling proposed by Huang et al., with the y∗-dependence weakly parameterized by Reτw.

Turbulence↗

Challenge Problem 1: Preliminary Results of the Direct Numerical Simulation of Transient Flows

This report presents the first direct numerical simulations (DNS) of transient mixed convection in an idealized downcomer-like channel (Challenge Problem 1, Phase II). Using the GPU-accelerated NekRS solver, we modeled a sudden decay in driving pressure, mimicking loss-of-flow events, and tracked the resulting evolution of Reynolds number, boundary-layer structure, turbulence statistics, and heat-transfer metrics. Key findings include the systematic thickening and eventual asymmetry of velocity and thermal boundary layers under buoyant deceleration; minimal “memory” lag in Reynolds shear stress and TKE profiles when sampled at matching Re, yet clear shifts of peak locations toward the cooled wall; overshoots in transient eddy-viscosity and eddy-diffusivity (and corresponding sub-unity turbulent Prandtl numbers) on the cooled side; and a pronounced transient Nusselt-number enhancement driven by wall-temperature inertia and residual eddy mixing. These effects combined to offer a temporary cooling margin above steady-state predictions during reactor LOF transients. Future work will extend this work to a more complex “Case II” geometry (90° turn + lower plenum) and generate multi-Re/Pr datasets for data-driven turbulence closures.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

The Crack-contact and the Free End Problem for a Strip Under Residual Stress

The plane problem for an infinite strip with two edge cracks under a given state of residual stress is considered. The residual stress is compressive near and at the surfaces and tensile in the interior of the strip. If the crack is deep enough to penetrate into the tensile zone, then the problem is one of crack-contact problem in which the depth of the contact area is an unknown which depends on the crack depth and the residual stress profile. The problem has applications to the static fatigue of glass plates and is solved for three typical residual stress profiles. In the limiting case of the crack crossing the entire plate thickness, the problem becomes a stress-free end problem for a semi-infinite strip under a given residual stress state away from the end. This is a typical stress diffusion problem in which decay behavior of the residual stress near and the nature of the normal displacement at the end of the semi-infinite strip are of special interest. For two typical residual stress states the solution is obtained, and some numerical results are given.

Bakioglu, M.↗

Parameterization of radiative flux profiles within layer clouds

This paper presents simple formulations for the vertical structure of solar and infrared radiative fluxes within a layer cloud overlying a partially reflective surface. The parameterized profile shapes are analytic, with governing parameters derived from bulk radiative properties and more physically based radiative transfer models. The bulk cloud, subcloud layer, and surface radiative properties are assumed to be known. The parameterizations are based on exponential functions of height, with decay scales related to cloud liquid water content. Although the results presented here are based on very simple assumptions about the cloud structure in the vertical, the method used is applicable to more general cases as well as to various other analytic and/or numerical radiative transfer calculations.

Hanson, Howard P.↗

Evolution of the solar nebula. I - Nonaxisymmetric structure during nebula formation

Numerical solutions of the equations of hydrodynamics, gravitation, and radiative transfer in three spatial dimensions are used to model the formation and time evolution of the early solar nebula in order to learn whether or not gravitational torques between nonaxisymmetric structures in the solar nebula can transport angular momentum rapidly enough to produce nebula clearing on astronomically indicated (10 to the 5 to 10 to the 7 yr) time scales. The models involve solutions for the collapse of spherical clouds with assumed initial density and rotation profiles onto protosuns of variable mass. Most of the models assume uniform initial density and rotation, and have variations in the initial parameters of cloud mass, cloud rotation rate, and protosun mass which are chosen to simulate a range of possible phases of early solar nebula evolution. The models show little tendency for directly forming small numbers of giant gaseous protoplanets through gaseous gravitational instability.

Boss, Alan P.↗

Computational Study of Chaotic and Ordered Solutions of the Kuramoto-Sivashinsky Equation

We report the results of extensive numerical experiments on the Kuramoto-Sivashinsky equation in the strongly chaotic regime as the viscosity parameter is decreased and increasingly more linearly unstable modes enter the dynamics. General initial conditions are used and evolving states do not assume odd-parity. A large number of numerical experiments are employed in order to obtain quantitative characteristics of the dynamics. We report on different routes to chaos and provide numerical evidence and construction of strange attractors with self-similar characteristics. As the 'viscosity' parameter decreases the dynamics becomes increasingly more complicated and chaotic. In particular it is found that regular behavior in the form of steady state or steady state traveling waves is supported amidst the time-dependent and irregular motions. We show that multimodal steady states emerge and are supported on decreasing windows in parameter space. In addition we invoke a self-similarity property of the equation, to show that these profiles are obtainable from global fixed point attractors of the Kuramoto-Sivashinsky equation at much larger values of the viscosity.

Smyrlis, Yiorgos S.↗

HST Studies of the Chromospheres, Wind, and Mass-Loss Rates of Cool Giant and Supergiant Stars

UV spectra of K-M giant and supergiant stars and of carbon stars have been acquired with the Goddard High Resolution Spectrograph (GHRS) on the Hubble Space Telescope (HST). These spectra have been used to measure chromospheric flow and turbulent velocities, study the acceleration of their stellar winds, acquire constraints on their outer atmospheric structure, and enable estimates of their mass-loss rates. Results from our observations of the giant stars Gamma Dra (K5 III hybrid), Alpha Tau (K5 III), Gamma Cru (M3.4 III), Mu Gem (M3 IIIab), and 30 Her (MG III), the supergiants Alpha Ori (M2 Iab) and Lambda Vel (K5 Ib), and the carbon stars TX Psc (NO; C6,2) and TW Hor (NO; C7,2) will be summarized and compared. The high resolution and wavelength accuracy of these data have allowed the direct measurement of the acceleration of the stellar winds in the chromospheres of several of these stars (from initial velocities of 3-9 km/s to upper velocities of 15-25 km/s) and of the chromospheric macroturbulence (-25-35 km/s). The high signal-to-noise and large dynamic range of these spectra have allowed the detection and identification of numerous new emission features, including weak C IV emission indicative of hot transition-region plasma in the non-coronal giant Alpha Tau, many new fluorescent lines of Fe II, and the first detection of fluorescent molecular hydrogen emission and of Ca II recombination lines in the UV spectrum of a giant star. The UV spectrum of two carbon stars have been studied with unprecedented resolution and reveal extraordinarily complicated Mg II lines nearly smothered by circumstellar absorptions. Finally, comparison of synthetic UV emission line profiles computed with the Lamers et al. (1987) Sobolev with Exact Integration (SEI) code with observations of chromospheric emission lines overlain with wind absorption features provides estimates of the mass-loss rates for four of these stars.

Carpenter, Kenneth G.↗

Mars Climate Orbiter's Investigation of the Atmosphere and Polar Caps

The Mars Climate Orbiter (MCO) is now on its way to Mars. It carries an atmospheric sounder whose observations will provide a continuous, global data set on weather and climate for a full Martian year. This paper describes the observation strategy and anticipated results from the Pressure Modulator Infrared Radiometer (PMIRR). PMIRR will measure vertical profiles of atmospheric infrared radiance in the 7 to 50 micron wavelength region extending from the surface of Mars to 80-km altitude. The observations have a vertical resolution of 5 km, or one-half the atmospheric scale height. From these radiance profiles we will retrieve profiles of atmospheric temperature, pressure, and the amounts of dust, condensates and water vapor. In addition, PMIRR will measure the radiative balance of the polar regions of Mars in an effort to better understand the short-term climate variability of the planet. The information obtained with PMIRR on MCO will be complementary to data obtained by the Thermal Emission Spectrometer (TES) and Radio Science (RS) experiments on the Mars Global Surveyor. A major emphasis of our research will be the assimilation of PMIRR data into numerical models of the Martian atmosphere. Assimilation schemes, of which several are currently in development, will permit the extension of measurements to spatial and temporal scales and to phenomena (e.g. winds) not observed directly by PMIRR.

McCleese, D. J.↗

Formation Of the Giant Planets By Concurrent Accretion Of Solids And Gas

New numerical simulations of the formation of the giant planets are presented, in which for the first time both the gas and planetesimal accretion rates are calculated in a self-consistent, interactive fashion. The simulations combine three elements: 1) three-body accretion cross-sections of solids onto an isolated planetary embryo, 2) a stellar evolution code for the planet's gaseous envelope, and 3) a planetesimal dissolution code within the envelope, used to evaluate the planet's effective capture radius and the energy deposition profile of accreted material. Major assumptions include: The planet is embedded in a disk of gas and small planetesimals with locally uniform initial surface mass density, and planetesimals are not allowed to migrate into or out of the planet's feeding zone. All simulations are characterized by three major phases. During the first phase, the planet's mass consists primarily of solid material. The planetesimal accretion rate, which dominates that of gas, rapidly increases owing to runaway accretion, then decreases as the planet's feeding zone is depleted. During the second phase, both solid and gas accretion rates are small and nearly independent of time. The third phase, marked by runaway gas accretion, starts when the solid and gas masses are about equal. It is engendered by a strong positive feedback on the gas accretion rates, driven by the rapid contraction of the gaseous envelope and the rapid expansion of the outer boundary, which depends on the planet's total mass. The overall evolutionary time scale is generally determined by the length of the second phase. The actual rates at which the giant planets accreted small planetesimals is probably intermediate between the constant rates assumed in most previous studies and the highly variable rates that we have used. Within the context, of the adopted model of planetesimal accretion, the joint constraints of the time scale for dissipation of the solar nebula and the current high-Z masses of the giant planets lead to estimates of the initial surface density (sigma(sub init)) of planetesimals in the outer region of the solar nebula. The results show sigma(sub init) approx. = 10 g/sq cm near Jupiter's orbit and that sigma(sub init) proportional to alpha(sup -2), where alpha is the distance from the Sun. These values are a factor of 3 - 4 times as high as that of the "minimum mass" solar nebula at Jupiter's distance and a factor of 2 - 3 times as high it Saturn's distance. Our estimates for the formation time of Jupiter and Saturn are 1 - 10 million years while those for Uranus fall in the range of 2 - 16 million years. These estimates follow from the properties of our Solar System and do not necessarily apply to giant planets in other planetary systems.

Pollack, James B.↗

A surface temperature and moisture parameterization for use in mesoscale numerical models

A modified multi-level soil moisture and surface temperature model is presented for use as in defining lower boundary conditions in mesoscale weather models. Account is taken of the hydraulic and thermal diffusion properties of the soil, their variations with soil type, and the mixing ratio at the surface. Techniques are defined for integrating the surface input into the multi-level scheme. Sample simulation runs were performed with the modified model and the original model defined by Pielke, et al. (1977, 1981). The models were applied to regional weather forecasting over soils composed of sand and clay loam. The new form of the model avoided iterations necessary in the earlier version of the model and achieved convergence at reasonable profiles for surface temperature and moisture in regions where the earlier version of the model failed.

Tremback, C. J.↗

Numerical Study of Wave Propagation in a Non-Uniform Flow

The propagation of acoustic waves originating from cylindrical and spherical pulses, in a non-uniform mean flow, and in the presence of a reflecting wall is investigated by Hardin and Pope approach using compact approximation of spatial derivatives. The 2-D and 3-D stagnation flows and a flow around a cylinder are taken as prototypes of real world flows with strong gradients of mean pressure and velocity. The intensity and directivity of acoustic wave patterns appear to be quite different from the benchmark solutions obtained in a static environment for the same geometry. The physical reasons for amplification and weakening of sound are discussed in terms of dynamics of wave profile and redistribution of acoustic energy and its potential and kinetic components. For an acoustic wave in the flow around a cylinder, the observed mean acoustic pressure is approximately doubled (upstream pulse position) and halved (downstream pulse position) in comparison with the sound propagation in static ambient conditions.

Povitsky, Alex↗

Validation and Sensitivity Analysis of a 1-D Lake Model Across Global Lakes

Lakes have important influence on weather and climate from local to global scales. However, their prediction using numerical models is notoriously difficult because lakes are highly heterogeneous across the globe, but observations are sparse. In this study, we assessed the performance of a 1-D lake model in simulating the thermal structures of 58 lakes with diverse morphometric and geographic characteristics by following the phase 2a local lake protocol of the Inter-sectoral Impact Model Intercomparison Project (ISIMIP2a). After calibration, the root-mean-square errors (RMSE) were below 2 °C for 70% and 75% of the lakes for epilimnion and full-profile temperature simulations, with an average of 1.71 °C and 1.43 °C, respectively. The model performance mainly depended on lake shape rather than location, supporting the possibility of grouping model parameters by lake shape for global applications. Furthermore, through machine-learning based parameter sensitivity tests, we identified turbulent heat fluxes, wind-driven mixing and water transparency as the major processes controlling lake thermal and mixing regimes. Snow density was also important for modeling the ice phenology of high-latitude lakes. The relative influence of the key processes and the corresponding parameters mainly depended on lake latitude and depth. Turbulent heat fluxes showed a decreasing importance in affecting epilimnion temperature with increasing latitude. Wind-driven mixing was less influential to lake stratification for deeper lakes while the impact of light extinction, on the contrary, showed a positive correlation with depth. Our findings may guide improvements in 1-D lake model parameterizations to achieve higher fidelity in simulating global lake thermal dynamics.

54 ENVIRONMENTAL SCIENCES↗

Evidence for the 3D Radiative Effects of Boundary-Layer Clouds From Observations of Direct and Diffuse Surface Solar Fluxes

Numerical experiments have revealed the importance of horizontal transport of light in the presence of clouds (“3D effects”), with consequences for climate, weather, and solar resource availability predictions. Yet, analysis of 3D effects from observations remain sparse because of the difficulty to isolate the effect of horizontal transport in radiation measurements. In this study, we provide observational evidence for 3D effects based on the direct-diffuse partition of surface solar fluxes. It is compared to outputs from the ecRad radiative transfer scheme run on retrieved cloud profiles. Furthermore, the direct-beam calculation takes careful account of the field-of-view of the pyrheliometer to ensure consistency between observed and modeled direct fluxes. Only the solver that accounts for 3D effects is able to reproduce the observed mean direct-diffuse partition as a function of solar zenith angle and cloud cover, in particular at large solar zenith angles where cloud sides intercept most of the direct beam.

54 ENVIRONMENTAL SCIENCES↗

Retention of Long-period Gas Giant Planets: Type II Migration Revisited

During their formation, emerging protoplanets tidally interact with their natal disks. Proto–gas giant planets, with Hill radii larger than the disk thickness, open gaps and quench gas flow in the vicinity of their orbits. It is usually assumed that their type II migration is coupled to the viscous evolution of the disk. Although this hypothesis provides an explanation for the origin of close-in planets, it also encounters a predicament on the retention of long-period orbits for most gas giant planets. Moreover, numerical simulations indicate that the planets' migrations are not solely determined by the viscous diffusion of their natal disk. In this work, we carry out a series of hydrodynamic simulations combined with analytic studies to examine the transition between different paradigms of type II migration. We find a range of planetary mass for which gas continues to flow through a severely depleted gap so that the surface density distribution in the disk region beyond the gap is maintained in a quasi-steady state. The associated gap profile modifies the location of corotation and Lindblad resonances. In the proximity of the planet's orbit, high-order Lindblad and corotation torque are weakened by the gas depletion in the gap, while low-order Lindblad torques near the gap walls preserve their magnitude. Consequently, the intrinsic surface density distribution of the disk delicately determines both the pace and direction of the planets' type II migration. We show that this effect might stall the inward migration of giant planets and preserve them in disk regions where the surface density is steep.

79 ASTRONOMY AND ASTROPHYSICS↗