Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “weak grids”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
NASA NTRS1981

A synoptic analysis of the first AVE-SESAME '79 period

Key features of a severe convection observed during April 10-11, 1979 as part of the Atmospheric Variability Experiment-Severe Environmental Storms and Mesoscale Experiment (AVE-SESAME) are examined. Three-hourly rawinsonde readings from 23 stations were taken, and vertical motion and divergence parameters are considered. The data were converted into a 127 km grid at the surface, and at 50 mb intervals from 900 mb to 100 mb by an objective analysis scheme, while a kinematic method was used to compute vertical motion. A weak upper tropospheric short wave trough was found to propagate from New Mexico into the Texas panhandle, while a jet maximum propagated eastward. The development of a strong wind maximum over Oklahoma and Kansas was associated with a rapid increase in upper-level divergence and the development of a small-scale pressure perturbation in the Texas panhandle, as well as a low-level jet and convergence, which led to rapid changes over the Red River Valley, where stability was decreased.

Moore, J. T.↗
NASA NTRS2006

Use, Assessment, and Improvement of the Loci-CHEM CFD Code for Simulation of Combustion in a Single Element GO2/GH2 Injector and Chamber

This document is a viewgraph presentation of a paper that documents a continuing effort at Marshall Space Flight Center (MSFC) to use, assess, and continually improve CFD codes to the point of material utility in the design of rocket engine combustion devices. This paper describes how the code is presently being used to simulate combustion in a single element combustion chamber with shear coaxial injectors using gaseous oxygen and gaseous hydrogen propellants. The ultimate purpose of the efforts documented is to assess and further improve the Loci-CHEM code and the implementation of it. Single element shear coaxial injectors were tested as part of the Staged Combustion Injector Technology (SCIT) program, where detailed chamber wall heat fluxes were measured. Data was taken over a range of chamber pressures for propellants injected at both ambient and elevated temperatures. Several test cases are simulated as part of the effort to demonstrate use of the Loci-CHEM CFD code and to enable us to make improvements in the code as needed. The simulations presented also include a grid independence study on hybrid grids. Several two-equation eddy viscosity low Reynolds number turbulence models are also evaluated as part of the study. All calculations are presented with a comparison to the experimental data. Weaknesses of the code relative to test data are discussed and continuing efforts to improve the code are presented.

Westra, Douglas G.↗
NASA NTRS1996

Parallel Newton-Krylov-Schwarz algorithms for the transonic full potential equation

We study parallel two-level overlapping Schwarz algorithms for solving nonlinear finite element problems, in particular, for the full potential equation of aerodynamics discretized in two dimensions with bilinear elements. The overall algorithm, Newton-Krylov-Schwarz (NKS), employs an inexact finite-difference Newton method and a Krylov space iterative method, with a two-level overlapping Schwarz method as a preconditioner. We demonstrate that NKS, combined with a density upwinding continuation strategy for problems with weak shocks, is robust and, economical for this class of mixed elliptic-hyperbolic nonlinear partial differential equations, with proper specification of several parameters. We study upwinding parameters, inner convergence tolerance, coarse grid density, subdomain overlap, and the level of fill-in in the incomplete factorization, and report their effect on numerical convergence rate, overall execution time, and parallel efficiency on a distributed-memory parallel computer.

Cai, Xiao-Chuan↗
DOE OSTI2021

Holes in Optical Lightning Flashes: Identifying Poorly Transmissive Clouds in Lightning Imager Data

Space-based optical lightning sensors including the lightning imaging sensor (LIS) and geostationary lightning mapper (GLM) are pixelated imagers that detect lightning as transient increases in cloud top illumination. Detection requires optical emissions to escape the cloud top to space with sufficient energy to trigger a pixel on the imaging array. Through scattering and absorption, certain clouds are able to block most light from reaching the instrument, causing a reduction in detection efficiency (DE) and possibly location accuracy (LA). Radiant lightning emissions that illuminate large cloud top areas are used to examine scenarios where clouds block light from reaching orbit. In some cases, these anomalies in the spatial radiance distribution from the lightning pulse lead to “holes” in the optical lightning flash where certain pixels fail to trigger. Such holes are identified algorithmically in the Tropical Rainfall Measuring Mission satellite LIS record and the microphysical properties of the coincident storm region are queried. We find that holes primarily occur in tall (IR T b < 235 K) convection (87%) and overhanging anvil clouds (10%). The remaining 3% of holes occur in moderate-to-weak convection or in clear air breaks between stormclouds. We further demonstrate how an algorithm that assesses the spatial radiance patterns from energetic lightning pulses might be used to construct an optical transmission gridded stoplight product for GLM that could help operators identify clouds with a potentially reduced DE and LA.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Subaru Hyper Suprime-Cam excavates colossal over- and underdense structures over 360 deg 2 out to z = 1

Subaru Strategic Program with the Hyper-Suprime Cam (HSC-SSP) has proven to be successful with its extremely wide area coverage in past years. Taking advantages of this feature, we report initial results from exploration and research of expansive over- and underdense structures at z = 0.3–1 based on the second Public Data Release where optical 5-band photometric data for ~ eight million sources with i < 23 mag are available over ~360 deg 2 . We not only confirm known superclusters but also find candidates of titanic over- and underdense regions out to z = 1. The mock data analysis suggests that the density peaks would involve one or more massive dark matter haloes (>1014 M ⊙ ) of the redshift, and the density troughs tend to be empty of massive haloes over >10 comoving Mpc. Besides, the density peaks and troughs at z ≲ 0.6 are in part identified as positive and negative weak lensing signals respectively, in mean tangential shear profiles, showing a good agreement with those inferred from the full-sky weak lensing simulation. The coming extensive spectroscopic surveys will be able to resolve these colossal structures in 3D space. Finally, the number density information over the entire survey field is available as grid-point data on the website of the HSC-SSP data release (https://hsc.mtk.nao.ac.jp/ssp/data-release/).

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2003

A Cartesian Adaptive Level Set Method for Two-Phase Flows

In the present contribution we develop a level set method based on local anisotropic Cartesian adaptation as described in Ham et al. (2002). Such an approach should allow for the smallest possible Cartesian grid capable of resolving a given flow. The remainder of the paper is organized as follows. In section 2 the level set formulation for free surface calculations is presented and its strengths and weaknesses relative to the other free surface methods reviewed. In section 3 the collocated numerical method is described. In section 4 the method is validated by solving the 2D and 3D drop oscilation problem. In section 5 we present some results from more complex cases including the 3D drop breakup in an impulsively accelerated free stream, and the 3D immiscible Rayleigh-Taylor instability. Conclusions are given in section 6.

Ham, F.↗
NASA NTRS1992

Turbulence and deterministic chaos

Several turbulent and nonturbulent solutions of the Navier-Stokes equations are obtained. The unaveraged equations are used numerically in conjunction with tools and concepts from nonlinear dynamics, including time series, phase portraits, Poincare sections, largest Liapunov exponents, power spectra, and strange attractors. Initially neighboring solutions for a low Reynolds number fully developed turbulence are compared. Several flows are noted: fully chaotic, complex periodic, weakly chaotic, simple periodic, and fixed-point. Of these, only fully chaotic is classified as turbulent. Besides the sustained flows, a flow which decays as it becomes turbulent is examined. For the finest grid, 128(exp 3) points, the spatial resolution appears to be quite good. As a final note, the variation of the velocity derivatives skewness of a Navier-Stokes flow as the Reynolds number goes to zero is calculated numerically. The value of the skewness is shown to become small at low Reynolds numbers, in agreement with intuitive arguments that nonlinear terms should be negligible.

Deissler, Robert G.↗
NASA NTRSDate not supplied

Cloud-Type Mean Cloud Properties and Associated Cloud Radiative Effects in the Tropical “Chimney” Zones Using 20 Years High-Resolution CERES Satellite Data

It has been well known that convection tends to be more intense over land than over ocean and continental convection generally contains wider cores that are protected from entrainment than their oceanic counterparts. How does this difference in convective intensity impact the rest of cloud types and associated cloud radiative effects (CREs)? We select five regions in the tropical “chimney” zones, with two regions over land (Africa and Amazon) and three regions over ocean (eastern and western Pacific, and Atlantic) to understand the differences in the cloud-type mean cloud properties and associated CREs using 20-years high-resolution (2x2 km2) CERES satellite data. The cloud types are based upon the ISCCP classification using the joint cloud-top pressure and cloud optical depth distribution. These five regions are compared to the entire tropics (25°S to 25°N). First, we compare the frequencies of occurrence of cloud types among the regions. Compared to the entire tropics, the “chimney” zones are cloudier, particularly, for cloud types with moderate cloud optical depths in the middle and upper troposphere. Except for the eastern Pacific region, the low-level cloud types are less abundant over land (due to higher boundary layers) and ocean (due to weaker lower-tropospheric subsidence). Over Africa, upper-level clouds are scarcer, due perhaps to weak production of convective anvils related to the drier atmosphere. Second, liquid water path (LWP) and ice water path (IWP) are compared among the regions. Compared to the entire tropics, mean LWPs for moderate cloud optical depths are higher for oceanic regions but lower for Africa while they are slightly lower in the lower troposphere but slightly higher in the middle troposphere for Amazon. Not surprisingly, mean IWPs are higher over land than over ocean, which is directly related to the difference in convective intensity. That is, the explanations to this result are the relatively weak convective intensity over the eastern Pacific and the abundance of anvil clouds over the western Pacific. Both reduce the mean IWPs. The CRE results are still being analyzed because the totally-clear grids (1° x 1°) of the CERES data product should be removed before taking a regional mean of clear-sky fluxes, which cause the wrong signs of CREs for a few cloud types.

cloud radiative effects↗
NASA NTRS1981

Lanning 14, a new DO white dwarf - Spectrophotometry and atmospheric parameters

Optical spectrophotometric observations are presented and analyzed for Lanning 14, which was recently classified as a DO white dwarf and has one of the strongest He II 4868 absorption lines of any known star, with strong He Bracket lines and weak He I. An effective temperature of 55,000 plus or minus 2500 K, log g near 8, and He/H ratio equal to or greater than 100, have been derived by means of a new grid of hot, high-gravity, mixed-composition models. It is concluded that Lanning 14 is both hotter and more helium-rich than the prototype DO star HZ 21.

Liebert, J.↗
NASA NTRS1986

MR Cygni revisited

New analysis tools and additional unanalyzed observations justify a reanalysis of MR Cygni. The reanalysis applied successively more restrictive physical models, each with an optimization program. The final model assigned separate first and second order limb darkening coefficients, from model atmospheres, to individual grid points. Proper operation of the optimization procedure was tested on simulated observational data, produced by light synthesis with assigned system parameters, and modulated by simulated observational error. The iterative solution converged to a weakly-determined mass ratio of 0.75. Assuming the B3 primary component is on the main sequence, the HR diagram location of the secondary from the light ratio (ordinate) and adjusted T sub eff (abscissa) was calculated. The derived mass ratio, together with a main-sequence mass for the B3 component, implies a main-sequence secondary spectral type of B4. The photometrically-determined secondary radii agree with this spectral type, in marginal disagreement with the B7 type from the HR diagram analysis. The individual masses, derived from the radial velocity curve of the primary component, the photometrically-determined i, and alternative values of derived mass ratio are seriously discrepant with main sequence objects. The imputed physical status of the system is in disagreement with representations that have appeared in the literature.

Linnell, Albert P.↗
DOE OSTI2025

A moment-conserving discontinuous Galerkin representation of the relativistic Maxwellian distribution

Kinetic simulations of relativistic gases and plasmas are critical for understanding diverse astrophysical and terrestrial systems, but the accurate construction of the relativistic Maxwellian, the Maxwell–Jüttner distribution, on a discrete simulation grid is challenging. Difficulties arise from the finite velocity bounds of the domain, which may not capture the entire distribution function, as well as errors introduced by projecting the function onto a discrete grid. Here, we present a novel scheme for iteratively correcting the moments of the projected distribution applicable to all grid-based discretizations of the relativistic kinetic equation. In addition, we describe how to compute the needed nonlinear quantities, such as Lorentz boost factors, in a discontinuous Galerkin scheme through a combination of numerical quadrature and weak operations. The resulting method accurately captures the distribution function and ensures that the moments match the desired values to machine precision.

astrophysical plasmas↗
DOE OSTI2024

DCMIP2016: the tropical cyclone test case

This paper describes and analyzes the Reed–Jablonowski (RJ) tropical cyclone (TC) test case used in the 2016 Dynamical Core Model Intercomparison Project (DCMIP2016). This intermediate-complexity test case analyzes the evolution of a weak vortex into a TC in an idealized tropical environment. Reference solutions from nine general circulation models (GCMs) with identical simplified physics parameterization packages that participated in DCMIP2016 are analyzed in this study at 50 km horizontal grid spacing, with five of these models also providing solutions at 25 km grid spacing. Evolution of minimum surface pressure (MSP) and maximum 1 km azimuthally averaged wind speed (MWS), the wind–pressure relationship, radial profiles of wind speed and surface pressure, and wind composites are presented for all participating GCMs at both horizontal grid spacings. While all TCs undergo a similar evolution process, some reach significantly higher intensities than others, ultimately impacting their horizontal and vertical structures. TCs simulated at 25 km grid spacings retain these differences but reach higher intensities and are more compact than their 50 km counterparts. These results indicate that dynamical core choice is an essential factor in GCM development, and future work should be conducted to explore how specific differences within the dynamical core affect TC behavior in GCMs.

Willson, Justin L.↗
NASA NTRS1997

Flow Quality and Boundary Layer Transition

The widely held view is that transition to turbulence in the Blasius boundary layer occurs via amplification and eventual nonlinear breakdown of initially small amplitude instabilities i.e. Tollmien-Schlichting (TS) waves. However this scenario is only observed for low amplitude free-stream turbulence levels, i.e. u/U < 0.1%. Bypass of linear TS instability mechanism occurs for higher EST levels, yet considerable differences exist between the few experiments carefully designed to assess the effect of EST on transition. The consensus is that EST leads to longitudinal streaks that form near the leading edge in the boundary layer . These streaks appeal to be regions of concentrated streamwise vorticity and they are often referred to as Klebanoff modes. The importance of mean flow free-stream nonuniformity (FSN) is not as widely appreciated as EST for characterizing wind tunnel flow quality. Here it is shown that, although the v like generated by a d=50micron wire located upstream of the contraction (Re(sub d)=6.6, x/d=45,000) is immeasurably small by the time it interacts with the leading edge in the test section, it is responsible for generation of a pair of weak streamwise vortices in the boundary layer downstream. The characteristics of these wake-induced vortices and their effect on TS waves are demonstrated. Small remnant FSN variations are also shown to exist downstream of a turbulence grid. The question arises Are the adverse effects introduced by the turbulence grid caused by FST or by small remnant FSN variations?

Watmuff, Jonathan H↗
DOE OSTI2025

Atmospheric River Frequency-Category Characteristics Shape U.S. West Coast Runoff

Abstractrunoff response to atmospheric rivers (ARs) over the U.S. West Coast. We focused on runoff time series variations impacted by AR characteristics (e.g., category and frequency) and land preconditions during Northern Hemisphere cool seasons in the period of 1940–2023. Results show that high-category ARs significantly increase local runoff with higher hourly precipitation rates leading to a greater incremental rate and peak runoff. Extreme runoff increases greatly with the AR category with an increase rate up to 12.5 times stronger than non-extreme runoff. Besides the AR category, land preconditions such as soil moisture and snowpack also play crucial roles in modulating runoff response. We found that runoff induced by weak-category ARs is more sensitive to land preconditions than high-category ARs, with high peak runoff occurring when soil is nearly saturated. Additionally, more than 50% of high-peak-runoff events in snow-covered grid cells are associated with rain-on-snow events particularly for the events associated with weaker ARs. Regression analysis reveals that AR precipitation and land preconditions jointly influence runoff, emphasizing the importance of including soil moisture and snowpack levels in AR impact assessments. The study also highlights the intensified runoff response to back-to-back ARs with short intervals, which may become more frequent with climate warming, posing increased flood risks via facilitating wet soil conditions. Our findings have significant implications for AR risk predictions and the development of prediction models for AR-induced runoff.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Updraft and Downdraft Core Size and Intensity as Revealed by Radar Wind Profilers: MCS Observations and Idealized Model Comparisons

This study explores the updraft and downdraft properties of mature stage mesoscale convective systems (MCSs) in terms of draft core width, shape, intensity, and mass flux characteristics. The observations use extended radar wind profiler (RWP) and surveillance radar datasets from the U.S. Department of Energy Atmospheric Radiation Measurement program for midlatitude (Oklahoma, USA) and tropical (Amazon, Brazil) sites. MCS drafts behave qualitatively similar to previous aircraft and RWP cloud summaries. The Oklahoma MCSs indicate larger and more intense convective up- and downdraft cores, and greater mass flux than Amazon MCS counterparts. However, similar size-intensity relationships and draft vertical profile behaviors are observed for both regions. Additional similarities include weak positive correlations between core intensity and core width (correlation coefficient r ~ 0.5), and increases in draft intensity with altitude. A model-observational intercomparison for draft properties (core width, intensity, mass flux) is also performed to illustrate the potential usefulness of statistical observed draft characterizations. Idealized simulations with the Weather Research and Forecasting model aligned with midlatitude MCS conditions are performed at model grid spacings (?x) that range from 4 km to 250 m. It is shown that the simulations performed at ?x = 250 m at similar mature MCS lifecycle stages are those that exhibit draft intensity, width, mass flux, and shape parameter performances best matching with observed properties.

radar, arm, wind profiler, goamazon↗
NASA NTRS1996

The Determination of the Large-Scale Circulation of the Pacific Ocean from Satellite Altimetry using Model Green's Functions

A Green's function method for obtaining an estimate of the ocean circulation using both a general circulation model and altimetric data is demonstrated. The fundamental assumption is that the model is so accurate that the differences between the observations and the model-estimated fields obey a linear dynamics. In the present case, the calculations are demonstrated for model/data differences occurring on very a large scale, where the linearization hypothesis appears to be a good one. A semi-automatic linearization of the Bryan/Cox general circulation model is effected by calculating the model response to a series of isolated (in both space and time) geostrophically balanced vortices. These resulting impulse responses or 'Green's functions' then provide the kernels for a linear inverse problem. The method is first demonstrated with a set of 'twin experiments' and then with real data spanning the entire model domain and a year of TOPEX/POSEIDON observations. Our present focus is on the estimate of the time-mean and annual cycle of the model. Residuals of the inversion/assimilation are largest in the western tropical Pacific, and are believed to reflect primarily geoid error. Vertical resolution diminishes with depth with 1 year of data. The model mean is modified such that the subtropical gyre is weakened by about 1 cm/s and the center of the gyre shifted southward by about 10 deg. Corrections to the flow field at the annual cycle suggest that the dynamical response is weak except in the tropics, where the estimated seasonal cycle of the low-latitude current system is of the order of 2 cm/s. The underestimation of observed fluctuations can be related to the inversion on the coarse spatial grid, which does not permit full resolution of the tropical physics. The methodology is easily extended to higher resolution, to use of spatially correlated errors, and to other data types.

Stammer, Detlef↗
NASA NTRSDate not supplied

Cloud-Precipitation Hybrid Regimes and their Projection onto IMERG Precipitation Data

We extend and enhance the concept of the Cloud Regimes (CRs) developed from two-dimensional joint histograms of cloud optical thickness and cloud top pressure from the Moderate Resolution Imaging Spectroradiometer (MODIS), by adding precipitation information in order to better understand cloud-precipitation relationships. Taking advantage of the high-resolution Integrated Multi-satellitE Retrievals for GPM (IMERG) precipitation dataset, cloud-precipitation “hybrid” regimes are derived by implementing the k-means clustering algorithm with advanced initialization and objective measures to determine the most optimal clusters. By expressing precipitation rates within 1-degree grid cell as histograms and making choices on the relative weight of cloud and precipitation, we could obtain several editions of hybrid cloud-precipitation regimes (CPRs), and examine their characteristics. In the deep tropics, when precipitation is weighted weakly, the cloud part of the hybrid entroids resembles the centroid of cloud-only regimes, but still tightens the cloud-precipitation relationship by decreasing the precipitation variability of each regime. As precipitation weight progressively increases, the shape of the cloudy part of the hybrid centroids becomes blunter, while the precipitation part of the centroids sharpens. In the case where cloud and precipitation are weighted equally, the CPRs representing high clouds with intermediate to heavy precipitation exhibit distinct features in the precipitation parts of the centroids, which allows us to project them onto the 30-minly IMERG domain. Such a projection can be used to overcome the temporal sparseness of MODIS cloud observations, which leads to great application potential for various convection-focused studies, including diurnal cycle analysis.

cloud-precipitation↗
DOE OSTI2020

Analysis of the SBP-SAT Stabilization for Finite Element Methods Part I: Linear Problems

In the hyperbolic community, discontinuous Galerkin (DG) approaches are mainly applied when finite element methods are considered. As the name suggested, the DG framework allows a discontinuity at the element interfaces, which seems for many researchers a favorable property in case of hyperbolic balance laws. On the contrary, continuous Galerkin methods appear to be unsuitable for hyperbolic problems and there exists still the perception that continuous Galerkin methods are notoriously unstable. To remedy this issue, stabilization terms are usually added and various formulations can be found in the literature. However, this perception is not true and the stabilization terms are unnecessary, in general. In this paper, we deal with this problem, but present a different approach. We use the boundary conditions to stabilize the scheme following a procedure that are frequently used in the finite difference community. Here, the main idea is to impose the boundary conditions weakly and specific boundary operators are constructed such that they guarantee stability. This approach has already been used in the discontinuous Galerkin framework, but here we apply it with a continuous Galerkin scheme. No internal dissipation is needed even if unstructured grids are used. Further, we point out that we do not need exact integration, it suffices if the quadrature rule and the norm in the differential operator are the same, such that the summation-by-parts property is fulfilled meaning that a discrete Gauss Theorem is valid. This contradicts the perception in the hyperbolic community that stability issues for pure Galerkin scheme exist. In numerical simulations, we verify our theoretical analysis.

97 MATHEMATICS AND COMPUTING↗