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 901 records · Page 50

SAGE III/ISS Rapid Data Analysis Through Dashboarding with Jupyter Notebooks

Spaceborne remote sensing observations of Earth’s atmosphere produce significant quantities of data over the life of each mission. In the case of the Stratospheric Aerosol and Gas Experiment III on the International Space Station (SAGE III/ISS) nearly four years of vertical profiles of atmospheric ozone, water vapor, and nitrogen dioxide concentrations as well as aerosol extinction coefficients have been released. The dichotomy of the desire for both long-term trends in the atmospheric state alongside the assessment of short-term impacts of major disruptive events such as volcanic eruptions and pyrocumulus injections requires agile tools to handle these cases in near real-time as new data are produced. The analysis landscape is further complicated by the desire to compare results between the numerous contemporary observations available for a given dataset. The SAGE III/ISS team has developed a suite of tools leveraging modern web-based frameworks allowing members to interact with a dashboard-style interface to load the data record, assess new profiles as they are generated and in ensemble, compare between species, and additionally add in measurements observed by other platforms as necessary. Leveraging a commonly packaged data format of NetCDF alongside the Python Jupyter Notebook framework, the data can be served to interested parties from an analysis server while still runnable on personal systems if required. This presentation illustrates the ecosystem developed by the SAGE III/ISS team, the applicability to measurements made by any limb-observing platform, and the benefit to transforming routine analyses into readily accessible dynamic plots. Frameworks currently exist at larger scales with projects such as GIOVANNI, and this illustration seeks to show that similar frameworks are accessible and possible within the local research environment while simultaneously unloading human processing cycles for more specialized analysis tasks.

Dashboarding↗

A cost-effective strategy for nonoscillatory convection without clipping

Clipping of narrow extrema and distortion of smooth profiles is a well known problem associated with so-called high resolution nonoscillatory convection schemes. A strategy is presented for accurately simulating highly convective flows containing discontinuities such as density fronts or shock waves, without distorting smooth profiles or clipping narrow local extrema. The convection algorithm is based on non-artificially diffusive third-order upwinding in smooth regions, with automatic adaptive stencil expansion to (in principle, arbitrarily) higher order upwinding locally, in regions of rapidly changing gradients. This is highly cost effective because the wider stencil is used only where needed-in isolated narrow regions. A recently developed universal limiter assures sharp monotonic resolution of discontinuities without introducing artificial diffusion or numerical compression. An adaptive discriminator is constructed to distinguish between spurious overshoots and physical peaks; this automatically relaxes the limiter near local turning points, thereby avoiding loss of resolution in narrow extrema. Examples are given for one-dimensional pure convection of scalar profiles at constant velocity.

Leonard, B. P.↗

Sharp monotonic resolution of discontinuities without clipping of narrow extrema

A strategy is presented for accurately simulating highly convective flows containing discontinuities such as density fronts or shock waves, without distorting smooth profiles or clipping narrow local extrema. The convection algorithm is based on non-artificially-diffusive third-order upwinding in smooth regions, with automatic adaptive stencil expansion to (in principle, arbitrarily) higher order upwinding locally, in regions of rapidly changing gradients. This is highly cost-effective because the wider stencil is used only where needed - in isolated narrow regions. A recently developed universal limiter assures sharp monotonic resolution of discontinuities without introducing artificial diffusion or numerical compression. An adaptive discriminator is constructed to distinguish between spurious overshoots and physical peaks; this automatically relaxes the limiter near local turning points, thereby avoiding loss of resolution in narrow extrema. Examples are given for one-dimensional pure convection of scalar profiles at constant velocity.

Leonard, B. P.↗

Computational Flow Field in Energy Efficient Engine (EEE)

In this paper, preliminary results for the recently-updated Open National Combustion Code (Open NCC) as applied to the EEE are presented. The comparison between two different numerical schemes, the standard Jameson-Schmidt-Turkel (JST) scheme and the advection upstream splitting method (AUSM), is performed for the cold flow and the reacting flow calculations using the RANS. In the cold flow calculation, the AUSM scheme predicts a much stronger reverse flow in the central recirculation zone. In the reacting flow calculation, we test two cases: gaseous fuel injection and liquid spray injection. In the gaseous fuel injection case, the overall flame structures of the two schemes are similar to one another, in the sense that the flame is attached to the main nozzle, but is detached from the pilot nozzle. However, in the exit temperature profile, the AUSM scheme shows a more uniform profile than that of the JST scheme, which is close to the experimental data. In the liquid spray injection case, we expect different flame structures in this scenario. We will give a brief discussion on how two numerical schemes predict the flame structures inside the EEE using different ways to introduce the fuel injection.

CFD↗

Aeroelastic oscillations of a pitching flexible wing with structural geometric nonlinearities: Theory and numerical simulation

In this paper, we focus on the derivation of an analytical model of the aeroelastic dynamics of an elastically mounted flexible wing. This equations of motion obtained serve to help understand the behaviour of the aeroelastic wind tunnel setup in question, which consists of a rectangular wing with a uniform NACA 0012 airfoil profile, whose base is free to rotate rigidly about a longitudinal axis. Of particular interest are the structural geometric nonlinearities primarily introduced by the coupling between the rigid body pitch degree-of-freedom and the continuous system. A coupled system of partial differential equations (PDEs) coupled with an ordinary differential equation (ODE) describing axial-bending-bending-torsion-pitch motion is derived using Hamilton's principle. A finite dimensional approximation of the system of coupled differential equations is obtained using the Galerkin method, leading to a system of coupled nonlinear ODEs. Subsequently, these nonlinear ODEs are solved numerically using Houbolt's method. The results that are obtained are verified by comparison with the results obtained by direct integration of the equations of motion using a finite difference scheme. Adopting a linear unsteady aerodynamic model, it is observed that the system undergoes coalescence flutter due to coupling between the rigid body pitch rotation dominated mode and the first flapwise bending dominated mode. Finally, the behaviour of the limit cycle oscillations is primarily influenced by the structural geometric nonlinear terms in the coupled system of PDEs and ODE.

42 ENGINEERING↗

Effect of Inflow Boundary Conditions on Hovering Tilt-Rotor Flows

A study on the effect of inflow turbulence boundary conditions on the local flow on and around a body in flight has been carried out. The study has been carried out using OVERFLOW2 flow solver using the default Spalart-Allmaras turbulence model in OVERFLOW. Many OVERFLOW turbulent flow simulations have been reported using the SA-fv3 model1 over the years. The present study demonstrates that the turbulence levels imposed as a boundary condition (b.c.) in the far-field as implemented in OVERFLOW are not correct. In fact, very low level of turbulence at the far-field boundary as implemented in OVERFLOW results in the SA-fv3 model predicting transition-like profiles on a given body. By choosing sufficiently high levels of the Reynolds number of turbulence, Rt, as a boundary condition, this anomalous behavior of SA-fv3 model is eliminated. Since numerous papers using OVERFLOW have been presented in the literature including that by the author2 using low level of inflow turbulence, it will be beneficial to the CFD community at large and in particular to the OVERFLOW community to understand this effect of the inflow b.c. in the SA-fv3 model. Various results reported over the years in domains such as high-lift applications, drag prediction applications and rotorcraft flow applications can be revisited using the new turbulence inflow boundary condition prescription as suggested in this study. It has been demonstrated in the literature that the effect of inflow turbulence level on the downstream development of flow is too important to be ignored. Demonstrative results are shown and compared with experiment for the case of hovering XV-15 rotor flow at a tip Reynolds number of 4.9 x 10(exp 6) and Mach number of 0.69 in Figs. 1 and 2 below. Fig. 1(a,b) shows skin friction predictions at radial station of r/R = 0.28 corresponding to 5 different inflow boundary conditions obtained with the SA-fv3 model. As the inflow turbulence is increased monotonically, from a value of Rt = 0.1 to 100.0 as shown in Fig. 1, a fully developed turbulent skin friction profile is realized progressively as shown in Fig. 1(b). Similarly skin friction results shown in Fig. 2(a,b) corresponding to radial stations of r/R = 0.72 reinforce the postulate made here that as the inflow turbulence level is increased gradually, fully turbulent on-the-blade simulation is realized uniformly. In the full paper extensive results on skin friction, pressure and velocity profiles as well as the Figure of Merit will be presented and compared with experiment that will illustrate the substantial effect of the inflow boundary conditions on the development of hovering rotor flow downstream.

Kaul, Upender K.↗

The Advanced Technology Microwave Sounder (ATMS): First Year On-Orbit

The Advanced Technology Microwave Sounder (ATMS) is a new satellite microwave sounding sensor designed to provide operational weather agencies with atmospheric temperature and moisture profile information for global weather forecasting and climate applications. A TMS will continue the microwave sounding capabilities first provided by its predecessors, the Microwave Sounding Unit (MSU) and Advanced Microwave Sounding Unit (AMSU). The first flight unit was launched a year ago in October, 2011 aboard the Suomi-National Polar-Orbiting Partnership (S-NPP) satellite, part of the new Joint Polar-Orbiting Satellite System (JPSS). Microwave soundings by themselves are the highest-impact input data used by Numerical Weather Prediction models; and A TMS, when combined with the Cross-track Infrared Sounder (CrIS), forms the Cross-track Infrared and Microwave Sounding Suite (CrIMSS). The microwave soundings help meet sounding requirements under cloudy sky conditions and provide key profile information near the surface. ATMS was designed & built by Aerojet Corporation in Azusa, California, (now Northrop Grumman Electronic Systems). It has 22 channels spanning 23-183 GHz, closely following the channel set of the MSU, AMSU-AI/2, AMSU-B, Microwave Humidity Sounder (MHS), and Humidity Sounder for Brazil (HSB). It continues their cross-track scanning geometry, but for the first time, provides Nyquist sample spacing. All this is accomplished with approximately V. the volume, Y, the mass, and Y, the power of the three AMSUs. A description will be given of its performance from its first year of operation as determined by post-launch calibration activities. These activities include radiometric calibration using the on-board warm targets and cold space views, and geolocation determination. Example imagery and zooms of specific weather events will be shown. The second ATMS flight model is currently under construction and planned for launch on the "Jl" satellite of the JPSS program in approximately 2016. Additional units are expected on the J2 and 13 satellites, as well as potentially on future European METOP satellites.

Kim, Edward J.↗

Numerical Analysis of Liquid Metal MHD Flow and Heat Transfer for Open-Surface Li Divertor in FNSF

Within the ongoing U.S.-based program on the development of liquid metal plasma-facing components, numerical simulations and analyses are performed to address the feasibility of the open-surface Li divertor. In the previous scoping studies (Smolentsev, 2021), heat-removal capabilities of the divertor were assessed using a simplified flow model for a slug-type velocity profile and constant flow thickness. Here, new analyses take into account forces acting on the flowing Li layer. Three reduced-order mathematical models are applied under the conditions of the U.S. Fusion Nuclear Science Facility (FNSF) to access magnetohydrodynamic (MHD) flow development effects, velocity distribution, and surface waves: 1) fully developed MHD flow; 2) quasi-2-D developing MHD flow; and 3) multiphase MHD flow. The obtained results for MHD flows and the surface heat flux computed with the plasma code scrape-off layer plasma simulation for ITER (SOLPS-ITER) are then used as input data to compute the temperature distribution in the divertor by solving the convection–diffusion energy equation.

Smolentsev, Sergey↗

The effect of ram-pressure stripping on dwarf galaxies

ABSTRACT Ram-pressure stripping (RPS) is a well observed phenomenon of massive spiral galaxies passing through the hot intracluster medium (ICM) of galaxy clusters. For dwarf galaxies (DGs) within a cluster, the transformation from gaseous to gas-poor systems by RPS is not easily observed and must happen in the outskirts of clusters. In a few objects in close by galaxy clusters and the field, RPS has been observed. Since cluster early-type DGs also show a large variety of internal structures (unexpected central gas reservoirs, blue stellar cores, composite radial stellar profiles), we aim in this study to investigate how ram pressure (RP) affects the interstellar gas content and therefore the star formation (SF) activity. Using a series of numerical simulations, we quantify the dependence of the stripped-off gas on the velocity of the infalling DGs and on the ambient ICM density. We demonstrated that SF can be either suppressed or triggered by RP depending on the ICM density and the DGs mass. Under some conditions, RP can compress the gas, so that it is unexpectedly retained in the central DG region and forms stars. When gas clouds are still bound against stripping but lifted from a thin disc and fall back, their new stars form an ellipsoidal (young) stellar population already with a larger velocity dispersion without the necessity of harassment. Most spectacularly, star clusters can form downstream in stripped-off massive gas clouds in the case of strong RP. We compare our results to observations.

Steyrleithner, P.↗

Development of a molecularly informed biogeochemical framework for reactive transport modeling of subsurface carbon inventories, transformations and fluxes (Final Report)

The overall objective of the project was to combine new molecular-level characterization strategies with soil carbon flux measurements to develop and evaluate model representations of subsurface carbon cycling. We expanded on existing studies in the East River watershed, Colorado, in collaboration with the Berkeley Lab Watershed Function Scientific Focus Area (SFA), SLAC Groundwater Quality SFA, and Rocky Mountain Biological Laboratory (RMBL) to develop an elevation and vegetation gradient that is now the subject of long-term monitoring by the USGS. To achieve the overall objective, we combined field studies of soil respiration with laboratory analyses, ranging from spectroscopy to incubation studies. This combination of techniques enabled us to develop a new understanding of the drivers of high-elevation soil respiration. In the process, we developed three new modeling approaches that improve our ability to conceptualize and ultimately to represent soil respiration in numerical models. The first approach is a plot-scale transient inverse model that can be used to determine in situ CO 2 production rates from measured concentration profiles and surface fluxes. Application of this method revealed the importance of plant phenology and deep CO 2 production in moderating CO 2 fluxes to the atmosphere. The second modeling approach is a molecular-scale tool that enables spectroscopic and elemental data for carbon speciation to be transformed into functional group abundances, or the ‘SOC-fga model’. This method uniquely enables carbon speciation to be tracked within a reactive transport framework to partition carbon among different pathways and storage zones within the soil. The third modeling tool builds strongly on the previous approaches and captures the microbial processes driving heterotrophic respiration. This ‘dormancy model’ allows the native soil microbial population to respond transiently to the presence or absence of water in order to catalyze carbon respiration. This approach was also compared to the simpler and more widely used first-order model using two experimental datasets with different temporal and spatial resolutions. Our results illustrate that the simpler first order model provides a robust and efficient representation of deep (>1 meter) soil respiration, but that shallow soils, where most respiration occurs, require explicit representation of moisture-dependent activation and dormancy rates. By assessing soil organic carbon turnover at multiple scales, we see a complex array of controls emerge. Atte scale of a hillslope, spatial heterogeneity in soil respiration rates dominates and is uncorrelated with instantaneous soil moisture and plant community. At the profile scale, the balance between plant inputs and water availability is the dominant control. At the microbial to molecular scale, physiological processes associated with carbon use and carbon speciation are important controls. These scale-dependent controls emphasize the need for new modeling approaches that examine their interactions and hierarchies.

54 ENVIRONMENTAL SCIENCES↗

Quantifying Variability and Controls of Riverine Dissolved Organic Carbon Exported to Arctic Coastal Margins of North America (Final Report)

This project involved implementation and application of a coupled permafrost hydrology and dissolved organic carbon process models to investigate how spatial and seasonal variations in terrestrial hydrology and soil freeze/thaw dynamics influence the mobilization, loading, and export of organic carbon to the stream network for selected arctic basins across northern Alaska and northwest Canada. The model simulations were constrained by detailed observations of in- stream chemistry, soil active layer profile moisture and temperature dynamics, streamflow, soil carbon inventories and satellite microwave remote sensing based assessments of surface soil freeze-thaw dynamics. We developed and applied the numerical modeling and data analysis, incorporating observed data for calibration and validation, in order to investigate the terrestrial hydrology, permafrost dynamics, and associated DOC production and loading to rivers across a region encompassing watersheds draining to the coast. The project produced six publications and three datasets archived in public repositories.

54 ENVIRONMENTAL SCIENCES↗

A proposed origin of the Olympus Mons escarpment

Olympus Mons (Nix Olympica) on Mars is delimited by a unique steep, nearly circular scarp. A pyroclastic model is proposed for the construct's origin. It is postulated that the Olympus Mons plateau is constructed predominantly of numerous ash-flow tuffs which were erupted from central sources over an extended period of time. Lava flows may be intercalated with the tuffs. A schematic radial profile incorporating the inferred compaction zones for an ash sheet is proposed. Following emplacement, eolian (and possibly fluvial) erosion and abrasion during dust storms would act on the ash sheets. Interior portions of the sheets would spall and slump following eolian erosion, generating steep, relatively smooth boundary scarps. The scarp would be circular due to symmetrical distribution of compaction zones. The model implies further that the Olympus Mons plateau rests on a more resistant rock substrate.

King, J. S.↗

The interpretations of ultraviolet observations of comets

The paper summarizes recent cometary UV observations, most of which were made in Ly-alpha light with instruments aboard earth-orbiting satellites. These include OAO-2 observations of comets Bennett and Tago-Sato-Kosaka, OGO-5 observations of comets Bennett and Encke, and numerous observations of comet Kohoutek. Models for the production of cometary hydrogen atoms are described, including the fountain, syndyname, and parent-daughter models. Calculations of emission line profiles and multiple-scattering effects are also discussed. Results of observations and interpretations are reviewed for each cited comet, far-UV observations in other emission lines are noted, and the use of comets as solar-wind probes is considered. It is concluded that the results of the present cometary Ly-alpha observations strongly support the concept of an icy conglomerate solid cometary nucleus and suggest water to be one of the most abundant molecules in comets.

Keller, H. U.↗

The use of computer-generated color graphic images for transient thermal analysis

Color computer graphics techniques were investigated as a means of rapidly scanning and interpreting large sets of transient heating data. The data presented were generated to support the conceptual design of a heat-sink thermal protection system (TPS) for a hypersonic research airplane. Color-coded vector and raster displays of the numerical geometry used in the heating calculations were employed to analyze skin thicknesses and surface temperatures of the heat-sink TPS under a variety of trajectory flight profiles. Both vector and raster displays proved to be effective means for rapidly identifying heat-sink mass concentrations, regions of high heating, and potentially adverse thermal gradients. The color-coded (raster) surface displays are a very efficient means for displaying surface-temperature and heating histories, and thereby the more stringent design requirements can quickly be identified. The related hardware and software developments required to implement both the vector and the raster displays for this application are also discussed.

Edwards, C. L. W.↗

Large Eddy Simulation of turbulent shear flows

The conceptual foundation underlying Large Eddy Simulation (LES) is summarized, and the numerical methods developed for simulation of the time-developing turbulent mixing layer and turbulent plane Poiseuille flow are discussed. Computational results show that the average Reynolds stress profile nearly attains the equilibrium shape which balances the downstream mean pressure gradient in the regions away from the walls. In the vicinity of the walls, viscous stresses are shown to be significant; together with the Reynolds stresses, these stresses balance the mean pressure gradient. It is stressed that the subgrid scale contribution to the total Reynolds stress is significant only in the vicinity of the walls. The continued development of LES is urged.

Moin, P.↗

Higher modes of the Orr-Sommerfeld problem for boundary layer flows

The discrete spectrum of the Orr-Sommerfeld problem of hydrodynamic stability for boundary layer flows in semi-infinite regions is examined. Related questions concerning the continuous spectrum are also addressed. Emphasis is placed on the stability problem for the Blasius boundary layer profile. A general theoretical result is given which proves that the discrete spectrum of the Orr-Sommerfeld problem for boundary layer profiles (U(y), 0,0) has only a finite number of discrete modes when U(y) has derivatives of all orders. Details are given of a highly accurate numerical technique based on collocation with splines for the calculation of stability characteristics. The technique includes replacement of 'outer' boundary conditions by asymptotic forms based on the proper large parameter in the stability problem. Implementation of the asymptotic boundary conditions is such that there is no need to make apriori distinctions between subcases of the discrete spectrum or between the discrete and continuous spectrums. Typical calculations for the usual Blasius problem are presented.

Lakin, W. D.↗

Thermal Analysis of the X-band 34-meter Antenna Feedcone

A thermal analysis was done to describe feedcone shell temperature profiles and induced thermal stresses in waveguide components. IMP-SS, a steady-state version of the thermal analysis computer program employing the implicit alterating-direction numerical technique documented in a previous report, was used to analyze the two dimensional regions of the feedcone. This portion of the analysis determined the free relative displacement of the feedcone for a worst-case thermal condition.

Hughes, R. D.↗

On the Kelvin-Helmholtz instabilities of supersonic shear layers

Analytical consideration is given to a Kelvin-Helmholtz instability in a finite width supersonic shear layer with linear and sinusoidal velocity profiles in the transition region. An analogy is explored between the pressure perturbation and the bound and scattering wave states of the one-dimensional Schroedinger equation. A numerical model is obtained for generating contour maps of the real and imaginary parts of the perturbation frequency in the (A,B)-plane. A comparison is made between the numerical results and analytical expressions derived for very small Mach numbers. Two unstable regimes are characterized in the (A,B)-plane. The method is concluded useful for both laboratory and astrophysics problems.

Choudhury, S. R.↗