Search NASA⌕ Search

SEARCH · Search NASA

Results for “Criticality Calculations”

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

Calculating Shocks In Flows At Chemical Equilibrium

Boundary conditions prove critical. Conference paper describes algorithm for calculation of shocks in hypersonic flows of gases at chemical equilibrium. Although algorithm represents intermediate stage in development of reliable, accurate computer code for two-dimensional flow, research leading up to it contributes to understanding of what is needed to complete task.

Eberhardt, Scott↗

Calculation of solar attenuation coefficient using ACCOS V along a critical scattering path

An appraisal of the optical properties of the AVHRR is made for a proposed orbital trajectory which results in an unexpected solar stray-light path; i.e., bypassing the external scan mirror and Cassegrain telescope and proceeding directly to the inner conical baffle. This stray-light path is directed onto an internal beamsplitter which, in turn, is directly observed by the detectors. Stray-light analysis, as a goal, seeks to remove or minimize the influence of such critical scattering paths. The AVHRR path is evaluated using a standard optical ray-tracing program, ACCOS V. To determine the expected sensor degradation, a calculation of the expected attenuation coefficient of scattered sunlight in the AVHRR sensor is estimated based on this important critical scattering path.

Ames, Alan J.↗

Critical exponents and scaling relations for self-organized critical phenomena

Critical indices beta, gamma delta, nv, etc. are defined and calculated for self-organized critical phenomena. Scaling relations are derived and checked numerically. The order-parameter exponent beta describes the spontaneous current and the relaxation to the criticl point. The power spectrum has 'l/f' behavior with the exponent phi = nv x z, where z is the dynamical critical exponent.

Tang, Chao↗

Utilization of Airborne Observations to Assess Model Parametrizations of Critical RH Profiles in the Arctic Ocean

Observational paucity is a reality in the Arctic Ocean. This is especially true for near-surface variables such as temperature, moisture, heat fluxes and BL clouds. As a result, modeling has become one of the major avenues for understanding current and future Arctic trends. Reanalyses are frequently used to force global ocean circulation and sea-ice models. But in northern high latitudes, model integrations and reanalyses are known to have large uncertainties in temperature and humidity profiles, and in boundary layer cloudiness. These are common sources of error in the surface radiative budget terms.An important way to diagnose these biases spatially and temporally is by using satellite remote sensing data. However, remotely-sensed observations also have large uncertainties, especially in near-surface temperature and relative humidity profiles.In situ observational studies are important in bridging our knowledge gap in regions such as the Arctic Ocean. Here, we utilize airborne and ship observations from the ARISE, ACME-IV, and ASCOS campaigns to construct critical relative humidity (RH) profiles over the Beaufort Sea. Such profiles are used as parameterization inputs in the NASA GOES global model to derive the total water condensate in a model grid-box, which determines the cloud fraction. Currently, the critical RH profiles are derived by global AIRS data, relaying mostly on mid-latitude regions, which are not necessarily relevant to the Arctic.We derive campaign-wide mean, standard deviation, and critical RH values, for grid size of 50x50 km and altitude bins between 50 to 400 m, covering both open-ocean and sea-ice covered regions. We compare profiles over open ocean and sea-ice, and look at correlations between the observed critical RH values and water condensate (by cloud number concentration) from observations versus the modeled ones. We then input our calculated values of minimal critical RH values into a set of GEOS single column model (SCM) simulations over the ARISE and ASCOS regions and compare the differences between the predicted values of cloud liquid water path (LWP), ice water path (IWP) and surface fluxes with the observed ones under the range of input parameterizations. Finally, we discuss the implications on surface radiative budget predictions in this region.

Segal-Rozenhaimer, Michal↗

Critical speeds and forced response solutions for active magnetic bearing turbomachinery, part 1

The prediction of critical speeds and forced response of active magnetic bearing turbomachinery is of great interest due to the increased use of this new and promising technology. Calculating the system undamped critical speeds and forced response is important to all those who are involved in the design of the active magnetic bearing system. An extended Jeffcott model which was used as an approximate solution to a more accurate transfer matrix procedure is presented. Theory behind a two-degree-of freedom extended Jeffcoat model is presented. Results of the natural frequency calculation are shown followed by the results of the forced response calculation. The system response was predicted for two types of forcing. A constant magnitude excitation with a wide frequency variation was applied at the bearings as one forcing function. The normal unbalance force at the midspan was the second source of excitation. The results of this extended Jeffcott solution gives useful design guidance for the influence of the first and third modes of a symmetric rotor system.

Keesee, J.↗

A numerical study of capillary stability in a circular cylindrical container with a concave spheroidal bottom

The paper describes a theoretical investigation of the stability under gravitational and surface forces of a liquid in a circular cylindrical container with a concave spheroidal bottom for the case in which the volume of the liquid is sufficiently small so that the bottom is not covered completely. The gravitational field is assumed to be directed along the symmetry axis of the container, and for a specific container shape the critical Bond number is calculated as a function of liquid volume for contact angles of gamma = 0, 1, 2, and 4 deg. For gamma = 0 deg, some critical equilibrium configurations and corresponding perturbation modes are presented.

Concus, P.↗

Adapting CLUTCH methodology to multigroup TSUNAMI-3D for eigenvalue sensitivity calculations

The sensitivity of the eigenvalue to uncertainties in nuclear data and its evaluation are important for nuclear criticality safety. TSUNAMI-3D sequences within the SCALE code system offer several options to the user community for calculating eigenvalue sensitivity coefficients with multigroup (MG) and continuous energy (CE) 3D transport capabilities. TSUNAMI-3D sequences implement the adjoint-based perturbation theory with MG KENO code, the Contributon Linked eigenvalue sensitivity/Uncertainty estimation via Track length importance CHaracterization (CLUTCH) method with CE KENO code, and the Iterated Fission Probability (IFP) method with CE KENO and Shift codes. Each method has benefits and limitations depending on the problem that is run. The work presented here aims to adapt the CLUTCH method, which enables the Contributon method's mesh-free, memory-efficient approach for calculating adjoint-weighted tallies for sensitivity calculations, to the MG TSUNAMI-3D sequence. This application would eliminate the explicit adjoint KENO calculation, as well as the memory-consuming mesh flux moment tallies required by the conventional MG TSUNAMI-3D. Smaller memory footprints in the CLUTCH methodology and relatively shorter runtimes in MG KENO transport can make MG TSUNAMI-3D a viable method for some complex problems. Moreover, this adaptation allows MG sensitivity calculations with Shift, ORNL's next-generation high-performance Monte Carlo transport code, which currently does not offer any sensitivity capabilities with MG particle transport simulations. Initial implementation of the new MG TSUNAMI-3D sequence and its preliminary results with a selected critical benchmark experiment in the Verified, Archived Library of Inputs and Data (VALID) are presented in this study.

KENO↗

Verification of the Uniformly-Ordered Binary Decision Algorithm in Correlated-Benchmark Whisper Calculations

Whisper is a nuclear criticality safety code package that aids analysts in validation exercises by computing upper subcritical limits (USL) for applications of interest. To obtain statistically meaningful, significant, and conservative USLs, the analyst must ensure that Whisper selects a sufficient number of benchmarks that are neutronically similar to the application. Many of the available benchmarks are correlated but are currently treated as independent, leading to an artificially small sample size, as their individual information contributions will be overestimated. To aid the analyst in obtaining a sufficient sample size, prior work [2] demonstrated application of the Uniformly-Ordered Binary Decision (UOBD) algorithm in adjusting benchmark weights to account for benchmark correlations. This work provides verification of the Whisper implementation and considers the impact of updated benchmark correlations compared to those available previously. We demonstrate that the UOBD algorithm performs as expected with an analytic example. With HEU-SOL-THERM-001 cases 1 through 10 as the applications, we compare the USLs computed with benchmark correlations available in the Whisper 1.1 release only to those computed with additional benchmark correlations from DICE 2023 and demonstrate substantive differences.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Modal analysis of turborotors using planar modes - Theory

The generalized dynamic equations of motion have been obtained by the direct stiffness method for multimass flexible rotor bearing systems including the effects of gyroscopic moments, disk skew, and rotor acceleration. A set of undamped critical speed mode shapes calculated from the average horizontal and vertical bearing stiffness is used to transform the equations of motion into a set of coupled modal equations of motion. The modal equations are coupled by the generalized bearing coefficients and the gyroscopic moments. An analysis using only undamped critical speeds or decoupled modal analysis assuming proportional damping may lead to erroneous results. This paper presents a rapid method of calculating rotor resonance speeds with their corresponding amplification factors, stability and unbalance response of turborotors. Examples of the application of this modal approach are presented and results are compared to those of other methods such as matrix transfer analysis.

Gunter, E. J.↗

Variation of the pressure limits of flame propagation with tube diameter for propane-air mixtures

An investigation was made of the variation of the pressure limits of flame propagation with tube diameter for quiescent propane with tube diameter for quiescent propane-air mixtures. Pressure limits were measured in glass tubes of six different inside diameters, with a precise apparatus. Critical diameters for flame propagation were calculated and the effect of pressure was determined. The critical diameters depended on the pressure to the -0.97 power for stoichiometric mixtures. The pressure dependence decreased with decreasing propane concentration. Critical diameters were related to quenching distance, flame speeds, and minimum ignition energy.

TEST APPARATUS - FLAME PROPOGATION - LFPL↗

Single event upset vulnerability of selected 4K and 16K CMOS static RAM's

Upset thresholds for bulk CMOS and CMOS/SOS RAMS were deduced after bombardment of the devices with 140 MeV Kr, 160 MeV Ar, and 33 MeV O beams in a cyclotron. The trials were performed to test prototype devices intended for space applications, to relate feature size to the critical upset charge, and to check the validity of computer simulation models. The tests were run on 4 and 1 K memory cells with 6 transistors, in either hardened or unhardened configurations. The upset cross sections were calculated to determine the critical charge for upset from the soft errors observed in the irradiated cells. Computer simulations of the critical charge were found to deviate from the experimentally observed variation of the critical charge as the square of the feature size. Modeled values of series resistors decoupling the inverter pairs of memory cells showed that above some minimum resistance value a small increase in resistance produces a large increase in the critical charge, which the experimental data showed to be of questionable validity unless the value is made dependent on the maximum allowed read-write time.

Kolasinski, W. A.↗

A measurement of $H_0$ from DESI DR1 using energy densities

We present a new measurement of the Hubble constant, independent of standard rulers and robust to pre-recombination modifications such as Early Dark Energy (EDE), obtained by calibrating the total energy density of the Universe. We start using the present-day photon density as an anchor, and use the baryon-to-photon ratio from Big Bang Nucleosynthesis based measurements and the baryon-to-matter ratio from the baryons' imprint on galaxy clustering to translate to a physical matter density at present day. We then compare this to measurements of the ratio of the matter density to the critical density ($Ω_{\mathrm{m}}$), calculated using the relative positions of the baryon acoustic oscillations, to measure the critical density of the universe and hence $H_0$. The important measurements of the evolution of the energy density all happen at low redshift, so we consider this a low-redshift measurement. We validate our method both on a suite of $N$-body mocks and on noiseless theory vectors generated across a wide range of Hubble parameters in both $Λ$CDM and EDE cosmologies. Using DESI DR1 data combined with the angular CMB acoustic scale and the latest BBN constraints, we find $H_0 = 69.0 \pm 2.5$ km s$^{-1}$ Mpc$^{-1}$, consistent with existing early and late-time determinations of the Hubble constant. We consider the impact of non-standard dark energy evolution on our measurement. Future data, including that from further iterations of DESI and from Euclid, will add to these results providing a powerful test of the Hubble tension.

Krolewski, Alex [U. Waterloo (main); Waterloo U., ↗

Advanced Materials for the Lunar Surface: Multiscale Computational Design of Refractory Alloys and Carbides

Emerging operational environments, such as the lunar surface, present novel challenges for NASA and drive the need for advanced materials in applications like fission surface power systems. To address these demands, computational materials science is rapidly evolving to augment or replace costly and hazardous empirical testing. Although materials selection at NASA remains predominantly experimentally driven, advanced simulation methodologies are being steadily integrated into the engineering lifecycle. This work details the application of multiscale simulation techniques—including first-principles calculations, CALPHAD, dislocation dynamics, and molecular dynamics—at NASA's Ames Research Center to evaluate advanced materials for extreme environments. First, we present contributions to the Space Nuclear Propulsion Project. Be-cause propellant channel coatings in nuclear thermal rockets must withstand high-pressure, high-temperature hydro-gen, optimizing these materials is critical. First-principles calculations were employed to establish a rigorous quantitative and qualitative understanding of the behavior of the refractory carbides ZrC, NbC, and their mixtures in high-enthalpy hydrogen environments. This necessitated the generation of high-fidelity thermodynamic models for both stoichiometric and carbon-depleted carbides, both with and without the presence of hydrogen. Furthermore, we highlight efforts under the Refractory Alloy Additive Manufacturing Build Optimization (RAAMBO) project, where existing and novel alloy compositions were assessed for additive manufacturing printability and subsequent performance in applications such as heat pipes and rocket nozzle extensions. This was accomplished through a comprehensive multiscale simulation framework that bridged the gap from the nanometer to the millimeter scale. Across both initiatives, rigorous validation against empirical data was prioritized. By systematically employing a verified and validated computational frame-work, we demonstrate how simulation effectively supports multidisciplinary engineering efforts, builds project-wide confidence, and drives critical materials development.

computational materials↗

Mineralogical Studies of Experimentally Shocked Dolomite: Implications for the Outgassing of Carbonates

Common rock-forming rhombohedral carbonates - calcite and dolomite - constitute a considerable fraction of terrestrial sediments that may be shocked during hypervelocity impacts, such as during the Cretaceous-Tertiary (K/T) event or at the Haughton structure. The relatively modest temperatures needed to decompose carbonates and to release their CO2 are easily attained during such impacts. However, detailed and quantitative understanding of the CO2 release of carbonates as a function of shock stress is still the subject of controversy, as are a number of other reactions and phase transitions that were suggested for carbonates. The first devolatilization studies of carbonates in the early 1980s suggested that incipient devolatilization of calcite and/or dolomite commences at modest pressures, in the range < 10-18 GPa (amounting to 0.03-0.3% CO2 loss), with massive decarbonation (30-40% CO2 loss) occurring at shock pressure of 20 GPa and above. However, most later studies revealed that both calcite and dolomite are unexpectedly stable under shock conditions and no significant outgassing has been observed at pressures as high as 40 GPa for calcite and 60 GPa for dolomite. Despite these uncertainties, the calculations of O'Keefe and Ahrens (and others) serve to illustrate that the K/T bolide liberated sufficient CO2 that a substantial temperature increase of the global atmosphere is possible due to CO2-triggered greenhouse effects. The quantitative understanding of the devolatilization of carbonates as a function of shock stress is obviously critical to refine such calculations. Additional information is contained in the original extended abstract.

Skala, R.↗

A Wall-Distance Method for Turbulence Modeling

The distance from a grid point to the closest wall surface, wall distance, is a funda- mental quantity in turbulence modeling. Efficiency of wall-distance calculations has become more critical as the size of computational grids has significantly increased in recent years. This paper reports on an initial implementation of a new search-based wall-distance method that is suitable for general unstructured computational fluid dynamics (CFD) grids and tailored for requirements specific for turbulence modeling. The method represents a two-step approach to calculate the wall distance. In the first step, the wall distance is approximated for each grid point as the minimum distance from this point to a vertex of a triangular face at the wall. The point-to-vertex distance calculation is relatively inexpensive but may lead to a significant error in the wall-distance ap- proximation, especially for grid points near the wall. In the second step, for grid points located within a predefined distance ( threshold ) from the wall, the wall distance is computed as the minimum distance to wall faces. As a result, the wall distance is exact for all grid points within the threshold. This two-step approach reduces the computational cost yet achieves high and controllable accuracy in the evaluation of the wall distance. Algorithmic enhancements are presented to improve efficiency of wall-distance computations. Comprehensive assessment of the new method is reported for large-scale unstructured CFD grids generated for the Fifth AIAA CFD High-Lift Prediction Workshop. The performance of the new wall-distance method compares favorably with performance of two established methods implemented in high-performance CFD codes.

Wall Distance↗