Search NASA⌕ Search

SEARCH · Search NASA

Results for “convergence diagnostics”

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 NTRS2001

Numerical Simulations and Diagnostic Studies of Meteorological Conditions During PEM-Tropics B

Provides a final report on the work accomplished by several meteorological scientists under a NASA grant in conjunction with the DC-8 component of Pacific Exploratory Mission (PEM)-Tropics B. The responsibilities of the principal investigator included collaboration with the Science Team on flight planning, presentation of forecasts, and the preparation of map discussions for each flight. In a published manuscript, the principal investigator summarized the meteorological conditions during PEM-TB which included mean flow patterns, subtropical anticyclones, the South Pacific Convergence Zone (SPCZ), and the Intertropical Convergence Zone (ITCZ). Methodologies used included streamlines, ten day backward trajectories, thermodynamic soundings, and satellite imagery. Other interests included air sampling for the purpose of determining pollution levels.

Fuelberg, Henry E.↗
NASA NTRS1984

Response of the annual and zonal mean winds and temperatures to variations in the heat and momentum sources

Comparisons are made among solutions to zonal-mean equations obtained with parameterized friction and radiative cooling and those forced from specified heat and momentum sources. Budget equations are defined for zonally averaged steady state responses in a thin spherically rotating atmosphere. The heat sources and sinks and mean meridional circulation that maintain observed annual and zonal mean temperatures are identified with a diagnostic calculation. Estimates are made of the surface sensible heating, atmospheric latent heating and vertical flux eddy divergences. The heat and moisture sources and sinks are varied to obtain the steady state responses. The Hadley circulation is fairly insensitive to changes in the strength of the eddy momentum flux when sufficient internal friction is present. Varying the width of the total precipitation of the intertropical convergence zone with fixed eddy fluxes and extratropical heat sources yields conditions similar to El Nino. Finally, a minimum speed is found for the jet stream after varying the horizontal eddy momentum fixing latent and eddy heat sources.

Schneider, E. K.↗
NASA NTRS1991

South Pacific convergence zone and global-scale circulations

The major goals of the research can be classified under three categories: (1) heat and moisture budget studies in the tropics; (2) the role of tropical forcing on subtropical wind maxima in the Southern Hemisphere; and (3) general circulation modeling of low latitude circulations. Research related to the first goal included budget estimates of precipitation, as well as a comparison to rainfall rates observed at numerous island stations. It also included a comparison between observed and satellite estimates of precipitable water. Results and methods used to achieve this goal are given. Included in the second research goal were diagnostic studies of the intraseasonal oscillation which is known to be an important component of organized tropical convection. Also, the role of extratropical forcing on subtropical wind maxima was examined through the use of Eliassen-Palm (EP) fluxes. Again, results and procedural details are discussed in section 3. To accomplish the third research goal, the GLA general circulation model (GCM) was ported to Purdue University and successfully used in a number of control and experimental runs to investigate the significance of tropical heating on subtropical and extratropical circulation features. As with the first two research goals, results and methodology are detailed. Documentation of publications, theses, conference preprint papers, and other scientific information relevant to the research results and significant accomplishments achieved under the contract are given. Finally, the conclusions of the sponsored research are given.

Vincent, Dayton G.↗
NASA NTRS1975

Numerical procedure for analyzing Langmuir probe data

A numerical procedure is proposed for analyzing Langmuir probe data in the presence of a two-group plasma containing both primary and Maxwellian electrons. The procedure is known as a least-squares differential-correction technique for determining the unknown coefficients of the governing equation. It is shown that for a given set of input data the results of the analysis are unique and independent of the initial estimate of electron temperature, that convergence is fastest when electron temperature is overestimated, that the results are sensitive to the region of the curve used as input and also to the voltage increment, and that plasma properties determined by the proposed numerical procedure are either consistent with those determined graphically or closer to expected values. With a suitable data acquisition system, the Fortran IV program worked out for this procedure could be used to provide real-time plasma diagnostic information for an operating ion thruster.

Beattie, J. R.↗
DOE OSTI2023

Numerical discreteness errors in multispecies cosmological N -body simulations

ABSTRACT We present a detailed analysis of numerical discreteness errors in two-species, gravity-only, cosmological simulations using the density power spectrum as a diagnostic probe. In a simple set-up where both species are initialized with the same total matter transfer function, biased growth of power forms on small scales when the solver force resolution is finer than the mean interparticle separation. The artificial bias is more severe when individual density and velocity transfer functions are applied. In particular, significant large-scale offsets in power are measured between simulations with conventional offset grid initial conditions when compared against converged high-resolution results where the force resolution scale is matched to the interparticle separation. These offsets persist even when the cosmology is chosen so that the two particle species have the same mass, indicating that the error is sourced from discreteness in the total matter field as opposed to unequal particle mass. We further investigate two mitigation strategies to address discreteness errors: the frozen potential method and softened interspecies short-range forces. The former evolves particles under the approximately ‘frozen’ total matter potential in linear theory at early times, while the latter filters cross-species gravitational interactions on small scales in low-density regions. By modelling closer to the continuum limit, both mitigation strategies demonstrate considerable reductions in large-scale power spectrum offsets.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1991

Retrieval of wind temperature and pressure from single Doppler radar and a numerical model

A 4D data assimilation algorithm to obtain 3D wind and thermodynamic fields from radial velocities and a numerical model is derived. It is a blend of the continuous updating technique of Charney et al. (1969) and a diagnostic pressure and temperature retrieval technique suggested by Gal-Chen (1978). This model is tested only against model-generated data which are viewed as 'data observed from a real atmosphere'. A thermal bubble developing in a dry, neutral environment is simulated. The thermal recovery technique is found to be a necessary procedure to obtain successful data assimilation results. The present assimilation method is found to be capable of reducing the 'observational error' and making the assimilation run converge toward the control run. Without data insertion, errors in the nonobserved wind component are large.

Liou, Yu-Chieng↗
NASA NTRS1995

The mesoscale forcing of a midlatitude upper-tropospheric jet streak by a simulated convective system. 1: Mass circulation and ageostrophic processes

The mutual forcing of a midlatitude upper-tropospheric jet streak by organized mesoscale adiabatic and diabatic processes within a simulated convective system (SCS) is investigated. Using isentropic diagnostics, results from a three-dimensional numerical simulation of an SCS are examined to study the isallobaric flow field, modes of dominant ageostrophic motion, and stability changes in relation to the mutual interdependence of adiabatic processes and latent heat release. Isentropic analysis affords an explicit isolation of a component of isallobaric flow associated with diabatic processes within the SCS. Prior to convective development within the simulations, atmospheric destabilization occurs through adiabatic ageostrophic mass adjustment and low-level convergence in association with the preexisting synoptic-scale upper-tropospheric jet streak. The SCS develops in a baroclinic zone and quickly initiates a vigorous mass circulation. By the mature stage, a pronounced vertical couplet of low-level convergence and upper-level mass divergence is established, linked by intense midtropospoheric diabatic heating. Significant divergence persists aloft for several hours subsequent to SCS decay. The dominant role of ageostrophic motion within which the low-level mass convergence develops is the adiabatic isallobaric component, while the mass divergence aloft develops principally through the diabatic isallobaric component. Both compnents are intrinsically linked to the convectively forced vertical mass transport. The inertial diabatic ageostrophic component is largest near the level of maximum heating and is responsible for the development of inertial instability to the north of SCS, resulting in this quadrant being preferred for outflow. The inertial advective component, the dominant term that produces the new downstream wind maximum, rapidly develops north of the SCS and through mutual adjustment creates the baroclinic support for the new jet streak.

Wolf, Bart J.↗
DOE OSTI2021

The Trinity High Explosive Implosion System: The Foundation for Precision Explosive Applications

This article is set during the 1944 and 1945 final push to complete Project Y—the Manhattan Project at Los Alamos—and focuses primarily on overcoming the challenge of creating and demonstrating a successful convergent explosive implosion to turn a subcritical quantity of plutonium into a critical mass. The critical mass would then efficiently yield kilotons of trinitrotoluene (TNT)-equivalent energy in about a microsecond, demonstrating the implosion atomic bomb concept. This work culminated in the Trinity atomic test near Alamogordo on July 16, 1945. This implosion effect demarcated the approach to explosive science and technology the Laboratory has followed ever since, including development of high-explosive synthesis and formulation, small and large test and diagnostic facilities, shock dynamics theory, high-explosive system design engineering, and three-dimensional implosion modeling and simulation using some of the fastest computers in the world. This work also ushered in a period of broader application of precision high explosives in conventional munitions, demolition, mining and oil exploration, and space travel.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
NASA NTRS1994

Photonics: From target recognition to lesion detection

Since 1989, Martin Marietta has invested in the development of an innovative concept for robust real-time pattern recognition for any two-dimensioanal sensor. This concept has been tested in simulation, and in laboratory and field hardware, for a number of DOD and commercial uses from automatic target recognition to manufacturing inspection. We have now joined Rose Health Care Systems in developing its use for medical diagnostics. The concept is based on determining regions of interest by using optical Fourier bandpassing as a scene segmentation technique, enhancing those regions using wavelet filters, passing the enhanced regions to a neural network for analysis and initial pattern identification, and following this initial identification with confirmation by optical correlation. The optical scene segmentation and pattern confirmation are performed by the same optical module. The neural network is a recursive error minimization network with a small number of connections and nodes that rapidly converges to a global minimum.

Henry, E. Michael↗
DOE OSTI2024

A fast, matrix-based method to perform omnidirectional pressure integration

Abstract Experimentally-measured pressure fields play an important role in understanding many fluid dynamics problems. Unfortunately, pressure fields are difficult to measure directly with non-invasive, spatially resolved diagnostics, and calculations of pressure from velocity have proven sensitive to error in the data. Omnidirectional line integration methods are usually more accurate and robust to these effects as compared to implicit Poisson equations, but have seen slower uptake due to the higher computational and memory costs, particularly in 3D domains. This paper demonstrates how omnidirectional line integration approaches can be converted to a matrix inversion problem. This novel formulation uses an iterative approach so that the boundary conditions are updated each step, preserving the convergence behavior of omnidirectional schemes while also keeping the computational efficiency of Poisson solvers. This method is implemented in Matlab and also as a GPU-accelerated code in CUDA-C++. The behavior of the new method is demonstrated on 2D and 3D synthetic and experimental data. Three-dimensional grid sizes of up to 125 million grid points are tractable with this method, opening exciting opportunities to perform volumetric pressure field estimation from 3D PIV measurements.

42 ENGINEERING↗
DOE OSTI2024

Automated 3D cytoplasm segmentation in soft X-ray tomography

Cells’ structure is key to understanding cellular function, diagnostics, and therapy development. Soft X-ray tomography (SXT) is a unique tool to image cellular structure without fixation or labeling at high spatial resolution and throughput. Fast acquisition times increase demand for accelerated image analysis, like segmentation. Currently, segmenting cellular structures is done manually and is a major bottleneck in the SXT data analysis. This paper introduces ACSeg, an automated 3D cytoplasm segmentation model. ACSeg is generated using semi-automated labels and 3D U-Net and is trained on 43 SXT tomograms of immune T cells, rapidly converging to high-accuracy segmentation, therefore reducing time and labor. Furthermore, adding only 6 SXT tomograms of other cell types diversifies the model, showing potential for optimal experimental design. ACSeg successfully segmented unseen tomograms and is published on Biomedisa, enabling high-throughput analysis of cell volume and structure of cytoplasm in diverse cell types.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2008

The ALI-ARMS Code for Modeling Atmospheric non-LTE Molecular Band Emissions: Current Status and Applications

The Accelerated Lambda Iteration (ALI) technique was developed in stellar astrophysics at the beginning of 1990s for solving the non-LTE radiative transfer problem in atomic lines and multiplets in stellar atmospheres. It was later successfully applied to modeling the non-LTE emissions and radiative cooling/heating in the vibrational-rotational bands of molecules in planetary atmospheres. Similar to the standard lambda iterations ALI operates with the matrices of minimal dimension. However, it provides higher convergence rate and stability due to removing from the iterating process the photons trapped in the optically thick line cores. In the current ALI-ARMS (ALI for Atmospheric Radiation and Molecular Spectra) code version additional acceleration of calculations is provided by utilizing the opacity distribution function (ODF) approach and "decoupling". The former allows replacing the band branches by single lines of special shape, whereas the latter treats non-linearity caused by strong near-resonant vibration-vibrational level coupling without additional linearizing the statistical equilibrium equations. Latest code application for the non-LTE diagnostics of the molecular band emissions of Earth's and Martian atmospheres as well as for the non-LTE IR cooling/heating calculations are discussed.

Kutepov, A. A.↗
DOE OSTI2024

Algorithm 1049: The Delaunay Density Diagnostic

Accurate approximation of a real-valued function depends on two aspects of the available data: the density of inputs within the domain of interest and the variation of the outputs over that domain. There are few methods for assessing whether the density of inputs is sufficient to identify the relevant variations in outputs—i.e., the “geometric scale” of the function—despite the fact that sampling density is closely tied to the success or failure of an approximation method. In this article, we introduce a general purpose, computational approach to detecting the geometric scale of real-valued functions over a fixed domain using a deterministic interpolation technique from computational geometry. The algorithm is intended to work on scalar data in moderate dimensions (2–10). Our algorithm is based on the observation that a sequence of piecewise linear interpolants will converge to a continuous function at a quadratic rate (in L 2 norm) if and only if the data are sampled densely enough to distinguish the feature from noise (assuming sufficiently regular sampling). We present numerical experiments demonstrating how our method can identify feature scale, estimate uncertainty in feature scale, and assess the sampling density for fixed (i.e., static) datasets of input–output pairs. Finally, we include analytical results in support of our numerical findings and have released lightweight code that can be adapted for use in a variety of data science settings.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Kinetic equilibrium reconstruction and the impact on stability analysis of KSTAR plasmas

High fidelity kinetic plasma equilibrium reconstructions are an essential requirement for accurate stability and disruption prediction analyses to support continuous operation of high beta tokamak plasmas. The present kinetic equilibrium reconstructions of plasmas in the KSTAR device include plasma density and temperature profiles from Thomson scattering and ion temperature from charge exchange spectroscopy diagnostics, and allowance for fast particle pressure. In addition, up to 25 channels of motional Stark effect diagnostic data are used to constrain the magnetic field pitch angle profile in the plasma to produce a reliable computation of the safety factor, q, profile. In this work, H-mode plasmas exhibit clear pedestal characteristics in the reconstructed pressure profile compared to internal transport barrier or L-mode plasmas. The plasma configuration and vertical position of inner strike points are validated by CCD and infrared camera images. Ideal and resistive magnetohydrodynamic (MHD) stability analyses using the DCON and resistive DCON codes utilize these kinetic equilibrium reconstructions to compare to experimental plasma stability. Equilibria with sufficiently low convergence error can provide reliable computation of ideal and resistive MHD stability analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Process-oriented evaluation of quasi-stationary Rossby waves and their impact on surface air temperature extremes in dynamical downscaling over North America

Quasi-stationary Rossby waves are a crucial component of the general circulation and play a significant role in regional water and energy cycles, as well as in extreme events. However, process-oriented evaluation for Rossby waves is rarely performed for dynamical downscaling simulations. To close this gap, we evaluate three classes of dynamical downscaling approaches, with a focus on quasi-stationary Rossby waves and their impact on surface air temperature over North America during Northern Hemisphere summer. The three classes of models differ in the way large-scale forcing is provided: a limited-area model (LAM) constrained only by lateral boundary conditions, represented by RegCM4 from the North American branch of the Coordinated Regional Downscaling Experiment (NA-CORDEX), a LAM with spectral nudging to maintain consistency in large-scale dynamics with the forcing data, represented by the Weather Research and Forecasting (WRF) model simulation in NA-CORDEX, and a global variable-resolution model with smoothly varying grid spacings, represented by the Community Atmosphere Model version 5.4, with the Model for Prediction Across Scales (MPAS) as its dynamical core (CAM-MPAS). With no constraints on the atmospheric dynamics, CAM-MPAS exhibits several mean biases in the upper-level circulations over the Pacific Coast region: a weaker subtropical jet, a northward-shifted mid-latitude jet, and an overestimated southerly flow. With the lateral boundary constraint alone, RegCM4 also exhibits weaker jets and overestimated southerly winds off the West Coast. Rossby ray theory reveals that those wind biases direct incoming Rossby waves northward. The erroneously routed Rossby waves distort the relationship between the accumulation of wave activity over the US West Coast and surface temperature anomalies over the Southern Great Plains, which emerges approximately 4 d after the convergence of wave-activity flux in the ERA-Interim reanalysis. Furthermore, the response of heatwaves to the extreme wave activity flux is not reproduced by the two models, a serious drawback as a dynamical downscaling framework is expected to connect large-scale forcing to local-scale phenomena. The WRF model employing spectral nudging is largely free from the aforementioned problems. A pair of sensitivity simulations suggests that spectral nudging is the key to improving the dynamics of quasi-stationary Rossby waves and their impact on surface air temperature. Our results also demonstrate the effectiveness of Rossby wave diagnostics that allow for realistic background flows for assessing the credibility of dynamical downscaling over North America, where incoming Rossby waves propagate through complex circulation patterns before traveling across the continent.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Assessment of DFT functionals for a minimal nitrogenase [Fe(SH) 4 H]- model employing state-of-the-art ab initio methods

We have designed a [Fe(SH)4H]- model with the fifth proton binding either to Fe or S. We show that the energy difference between these two isomers (ΔE) is hard to estimate with quantum-mechanical (QM) methods. For example, different density functional theory (DFT) methods give ΔE estimates that vary by almost 140 kJ/mol, mainly depending on the amount of exact Hartree–Fock included (0%–54%). The model is so small that it can be treated by many high-level QM methods, including coupled-cluster (CC) and multiconfigurational perturbation theory approaches. With extrapolated CC series (up to fully connected coupled-cluster calculations with singles, doubles, and triples) and semistochastic heat-bath configuration interaction methods, we obtain results that seem to be converged to full configuration interaction results within 5 kJ/mol. Our best result for ΔE is 101 kJ/mol. With this reference, we show that M06 and B3LYP-D3 give the best results among 35 DFT methods tested for this system. Brueckner doubles coupled cluster with perturbaitve triples seems to be the most accurate coupled-cluster approach with approximate triples. CCSD(T) with Kohn–Sham orbitals gives results within 4–11 kJ/mol of the extrapolated CC results, depending on the DFT method. Single-reference CC calculations seem to be reasonably accurate (giving an error of ~5 kJ/mol compared to multireference methods), even if the D1 diagnostic is quite high (0.25) for one of the two isomers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Simulation of a rapid compression machine for evaluation of ignition chemistry and soot formation using gasoline/ethanol blends

Due to the projected decline of demand for gasoline in light duty engines and the advent of ethanol as a green fuel, the use of gasoline/ethanol blend fuels in heavy duty applications are being investigated as they are projected to have lower cost and lower lifecycle green house gas (GHG) emissions. In heavy duty engines, the primary mode of combustion is mixing controlled combustion where wide range of mixture conditions (equivalence ratio) exist. Soot emissions of these fuels in richer conditions are not well understood. The goal of this research is to evaluate some commercially available soot modeling codes for the particulate matter emissions from gasoline/ethanol fuel blends, especially at fuel rich conditions. A Rapid Compression Machine (RCM) is modeled in a three-dimensional numerical simulation using CONVERGE computational software using a reduced chemical kinetic mechanism with SAGE chemistry solver and a RANS k-ϵturbulence model with a sector model including the creviced piston. The creviced piston is used in the experimental setup to reduce boundary layer effects and to maintain a homogeneous core in the reaction cylinder. Computational fluid dynamics simulations are conducted for different gasoline-ethanol fuel blends from E10 (10% ethanol v/v) to E100. The fuel blend is modeled as a surrogate mixture of toluene, iso-octane, n-heptane for gasoline content, and ethanol. The computational results were validated against experimental results using pressure measurements and laser extinction diagnostics. Different soot models are investigated to evaluate their capability of predicting the sooting tendencies of fuel blends, especially in richer conditions experienced during mixing-controlled combustion. The experimental combustion characteristics such as the ignition delay of different blends of fuel are reasonably well predicted. The Particulate Size Mimic (PSM) model accurately predicts the soot generation characteristics of the different fuels, but the Hiroyasu-NSC model falls short in this regard. For accurate prediction of soot with the PSM model, the thermodynamic conditions during combustion must be accurately modeled. While the current computational modeling tools can produce accurate results for the prediction of particulate matter emissions, there is much work to be done in improving our understanding of the underlying fundamental processes.

Energy & Fuels↗
NASA NTRS2013

Consistency of Estimated Global Water Cycle Variations Over the Satellite Era

Motivated by the question of whether recent indications of decadal climate variability and a possible "climate shift" may have affected the global water balance, we examine evaporation minus precipitation (E-P) variability integrated over the global oceans and global land from three points of view-remotely sensed retrievals / objective analyses over the oceans, reanalysis vertically-integrated moisture convergence (MFC) over land, and land surface models forced with observations-based precipitation, radiation and near-surface meteorology. Because monthly variations in area-averaged atmospheric moisture storage are small and the global integral of moisture convergence must approach zero, area-integrated E-P over ocean should essentially equal precipitation minus evapotranspiration (P-ET) over land (after adjusting for ocean and land areas). Our analysis reveals considerable uncertainty in the decadal variations of ocean evaporation when integrated to global scales. This is due to differences among datasets in 10m wind speed and near-surface atmospheric specific humidity (2m qa) used in bulk aerodynamic retrievals. Precipitation variations, all relying substantially on passive microwave retrievals over ocean, still have uncertainties in decadal variability, but not to the degree present with ocean evaporation estimates. Reanalysis MFC and P-ET over land from several observationally forced diagnostic and land surface models agree best on interannual variations. However, upward MFC (i.e. P-ET) reanalysis trends are likely related in part to observing system changes affecting atmospheric assimilation models. While some evidence for a low-frequency E-P maximum near 2000 is found, consistent with a recent apparent pause in sea-surface temperature (SST) rise, uncertainties in the datasets used here remain significant. Prospects for further reducing uncertainties are discussed. The results are interpreted in the context of recent climate variability (Pacific Decadal Oscillation, Atlantic Meridional Overturning), and efforts to distinguish these modes from longer-term trends.

Robertson, F. R.↗