Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dependability modeling”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
NASA NTRS2020

Development of a Daily GRACE Mascon Solution for Terrestrial Water Storage

The Gravity Recovery and Climate Experiment (GRACE) and GRACE-Follow On missions have provided a global history of terrestrial water storage changes since 2002. Traditional GRACE products resolve monthly time-variable gravity at spatial resolutions of 300-500 km, with many recent efforts focusing on regularized mass concentration (mascon) solutions to better resolve signals spatially. However, monthly resolution inhibits the applicability of GRACE to investigations of sub-monthly signals. This study presents a new daily mass change solution estimated as deviations from the NASA Goddard Space Flight Center (GSFC) converged monthly mascon product and quantifies the fundamental trade-off between temporal and spatial resolution in GRACE-only solutions. As an iteration of the GSFC product, this daily solution represents a mixing of monthly information at higher relative spatial resolution with daily information at lower relative spatial resolution, resolving each temporal scale at the best resolution achievable without requiring hydrological model dependencies in the estimation. The resolution of the daily mascons is 300-400 km at high latitudes and 600-1,000 km lower, depending on each daily orbital track and the proximity of each mascon to constraint region boundaries (e.g., coastlines). Through simulations and model comparisons, we demonstrate that daily signals over areas larger than 400,000 km(2) are recoverable, with basins larger than 800,000 km(2) exhibiting strong signal recovery relative to leakage. This analysis establishes baselines for daily signal recovery from GRACE in the context of longer time scale solutions, characterizes leakage inherent in daily GRACE information, and creates new opportunities for applying GRACE to investigations of sub-monthly signals.

M. J. Croteau↗
NASA NTRS1994

Time-dependent simulation of cosmic-ray shocks, including Alfven transport

Time evolution of plane, cosmic-ray modified shocks has been simulated numerically for the case with parallel magnetic fields. Computations were done in a 'three-fluid' dynamical model incorporating cosmic-ray and Alfven-wave energy transport equations. Nonlinear feedback from the cosmic rays and Alfven waves is included in the equation of motion for the underlying plasma, as is the finite propagation speed and energy dissipation of the Alfven waves. Exploratory results confirm earlier, steady state analyses that found these Alfven transport effects to be potentially important when the upstream Alfven speed and gas sound speeds are comparable. As noted earlier, Alfven transport effects tend to reduce the transfer of energy through a shock from gas to energetic particles. These studies show as well that the timescale for modification of the shock is altered in nonlinear ways. It is clear, however, that the consequences of Alfven transport are strongly model dependent and that both advection of cosmic rays by the waves and dissipation of wave energy in the plasma will be important to model correctly when quantitative results are needed. Comparison is made between simulations based on a constant diffusion coefficient and more realistic diffusion models allowing the diffusion coefficient to vary in response to changes in Alfven wave intensity. No really substantive differences were found between them.

Jones, T. W.↗
NASA NTRS1982

Clouds and wet removal as causes of variability in the trace-gas composition of the marine troposphere

A modeling study of the effects of clouds and wet removal on the chemistry of the remote marine troposphere is described. Using a time-dependent model with parameterized vertical transport to calculate trace-gas concentrations, it is found that large variations in key species (e.g., HNO3, H2CO, and H2O2) result from simulations of sporadic rainfall, changes in cloud cover, and external inputs such as surface NO sources. Depending on the frequency and intensity of an event, the effects of these perturbations may persist for several days, thereby invalidating assumptions of photochemical equilibrium in the interpretation of measurements. Long-term integrations with fixed boundary conditions and regularly occurring cloud and rain episodes demonstrate a strong sensitivity of the mean concentration of longer-lived soluble gases to precipitation frequency but also confirm the validity of using properly chosen parameterizations of wet removal in steady state calculations.

Thompson, A. M.↗
DOE OSTI2023

Search for leptoquarks decaying into the bτ final state in pp collisions at $\sqrt{\textrm{s}}$ = 13 TeV with the ATLAS detector

A search for leptoquarks decaying into the $bτ$ final state is performed using Run 2 proton-proton collision data from the Large Hadron Collider, corresponding to an integrated luminosity of 139 fb —1 at $\sqrt{\textrm{s}}$ = 13 TeV recorded by the ATLAS detector. The benchmark models considered in this search are vector leptoquarks with electric charge of 2/3e and scalar leptoquarks with an electric charge of 4/3e. No significant excess above the Standard Model prediction is observed, and 95% confidence level upper limits are set on the cross-section times branching fraction of leptoquarks decaying into bτ. For the vector leptoquark production two models are considered: the Yang-Mills and Minimal coupling models. In the Yang-Mills (Minimal coupling) scenario, vector leptoquarks with a mass below 1.58 (1.35) TeV are excluded for a gauge coupling of 1.0 and below 2.05 (1.99) TeV for a gauge coupling of 2.5. In the case of scalar leptoquarks, masses below 1.28 (1.53) TeV are excluded for a Yukawa coupling of 1.0 (2.5). Finally, an interpretation of the results with minimal model dependence is performed for each of the signal region categories, and limits on the visible cross-section for beyond the Standard Model processes are provided.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1995

Analysis of composite ablators using massively parallel computation

In this work, the feasibility of using massively parallel computation to study the response of ablative materials is investigated. Explicit and implicit finite difference methods are used on a massively parallel computer, the Thinking Machines CM-5. The governing equations are a set of nonlinear partial differential equations. The governing equations are developed for three sample problems: (1) transpiration cooling, (2) ablative composite plate, and (3) restrained thermal growth testing. The transpiration cooling problem is solved using a solution scheme based solely on the explicit finite difference method. The results are compared with available analytical steady-state through-thickness temperature and pressure distributions and good agreement between the numerical and analytical solutions is found. It is also found that a solution scheme based on the explicit finite difference method has the following advantages: incorporates complex physics easily, results in a simple algorithm, and is easily parallelizable. However, a solution scheme of this kind needs very small time steps to maintain stability. A solution scheme based on the implicit finite difference method has the advantage that it does not require very small times steps to maintain stability. However, this kind of solution scheme has the disadvantages that complex physics cannot be easily incorporated into the algorithm and that the solution scheme is difficult to parallelize. A hybrid solution scheme is then developed to combine the strengths of the explicit and implicit finite difference methods and minimize their weaknesses. This is achieved by identifying the critical time scale associated with the governing equations and applying the appropriate finite difference method according to this critical time scale. The hybrid solution scheme is then applied to the ablative composite plate and restrained thermal growth problems. The gas storage term is included in the explicit pressure calculation of both problems. Results from ablative composite plate problems are compared with previous numerical results which did not include the gas storage term. It is found that the through-thickness temperature distribution is not affected much by the gas storage term. However, the through-thickness pressure and stress distributions, and the extent of chemical reactions are different from the previous numerical results. Two types of chemical reaction models are used in the restrained thermal growth testing problem: (1) pressure-independent Arrhenius type rate equations and (2) pressure-dependent Arrhenius type rate equations. The numerical results are compared to experimental results and the pressure-dependent model is able to capture the trend better than the pressure-independent one. Finally, a performance study is done on the hybrid algorithm using the ablative composite plate problem. It is found that there is a good speedup of performance on the CM-5. For 32 CPU's, the speedup of performance is 20. The efficiency of the algorithm is found to be a function of the size and execution time of a given problem and the effective parallelization of the algorithm. It also seems that there is an optimum number of CPU's to use for a given problem.

Shia, David↗
NASA NTRS2018

Modeling the Proton Radiation Belt With Van Allen Probes Relativistic Electron-Proton Telescope Data

An empirical model of the proton radiation belt is constructed from data taken during 2013-2017 by the Relativistic Electron-Proton Telescopes on the Van Allen Probes satellites. The model intensity is a function of time, kinetic energy in the range 18-600 megaelectronvolts, equatorial pitch angle, and L shell of proton guiding centers. Data are selected, on the basis of energy deposits in each of the nine silicon detectors, to reduce background caused by hard proton energy spectra at low L. Instrument response functions are computed by Monte Carlo integration, using simulated proton paths through a simplified structural model, to account for energy loss in shielding material for protons outside the nominal field of view. Overlap of energy channels, their wide angular response, and changing satellite orientation require the model dependencies on all three independent variables be determined simultaneously. This is done by least squares minimization with a customized steepest descent algorithm. Model uncertainty accounts for statistical data error and systematic error in the simulated instrument response. A proton energy spectrum is also computed from data taken during the 8 January 2014 solar event, to illustrate methods for the simpler case of an isotropic and homogeneous model distribution. Radiation belt and solar proton results are compared to intensities computed with a simplified, on-axis response that can provide a good approximation under limited circumstances.

Selesnick, R. S.↗
NASA NTRS1978

Dynamic MHD modeling of solar wind corotating stream interaction regions observed by Pioneer 10 and 11

The use of the Pioneer 10 and 11 projects to test an MHD one-dimensional time-dependent model of corotating solar wind streams during the period from Sept. 30 to Nov. 25, 1973 is described. During this period, five or six corotating interaction regions streamed past the two spacecraft, and, as a result of multiple-spacecraft radial alignment and temporally varying conditions at the solar wind source, the pattern predicted by the Steinolfson et al. (1975) model could be compared with observations. The results, in general, support the validity of the model, although the neglect of thermal energy exchange leads to incorrect values for the proton temperature. A detailed analysis of a stream is discussed.

Dryer, M.↗
NASA NTRS1990

Distributed deformation and block rotation in 3D

The authors address how block rotation and complex distributed deformation in the Earth's shallow crust may be explained within a stationary regional stress field. Distributed deformation is characterized by domains of sub-parallel fault-bounded blocks. In response to the contemporaneous activity of neighboring domains some domains rotate, as suggested by both structural and paleomagnetic evidence. Rotations within domains are achieved through the contemporaneous slip and rotation of the faults and of the blocks they bound. Thus, in regions of distributed deformation, faults must remain active in spite of their poor orientation in the stress field. The authors developed a model that tracks the orientation of blocks and their bounding faults during rotation in a 3D stress field. In the model, the effective stress magnitudes of the principal stresses (sigma sub 1, sigma sub 2, and sigma sub 3) are controlled by the orientation of fault sets in each domain. Therefore, adjacent fault sets with differing orientations may be active and may display differing faulting styles, and a given set of faults may change its style of motion as it rotates within a stationary stress regime. The style of faulting predicted by the model depends on a dimensionless parameter phi = (sigma sub 2 - sigma sub 3)/(sigma sub 1 - sigma sub 3). Thus, the authors present a model for complex distributed deformation and complex offset history requiring neither geographical nor temporal changes in the stress regime. They apply the model to the Western Transverse Range domain of southern California. There, it is mechanically feasible for blocks and faults to have experienced up to 75 degrees of clockwise rotation in a phi = 0.1 strike-slip stress regime. The results of the model suggest that this domain may first have accommodated deformation along preexisting NNE-SSW faults, reactivated as normal faults. After rotation, these same faults became strike-slip in nature.

Scotti, Oona↗
NASA NTRS2000

Towards LES Models of Jets and Plumes

As pointed out by Rodi standard integral solutions for jets and plumes developed for discharge into infinite, quiescent ambient are difficult to extend to complex situations, particularly in the presence of boundaries such as the sea floor or ocean surface. In such cases the assumption of similarity breaks down and it is impossible to find a suitable entrainment coefficient. The models are also incapable of describing any but the most slowly varying unsteady motions. There is therefore a need for full time-dependent modeling of the flow field for which there are three main approaches: (1) Reynolds averaged numerical simulation (RANS), (2) large eddy simulation (LES), and (3) direct numerical simulation (DNS). Rodi applied RANS modeling to both jets and plumes with considerable success, the test being a match with experimental data for time-averaged velocity and temperature profiles as well as turbulent kinetic energy and rms axial turbulent velocity fluctuations. This model still relies on empirical constants, some eleven in the case of the buoyant jet, and so would not be applicable to a partly laminar plume, may have limited use in the presence of boundaries, and would also be unsuitable if one is after details of the unsteady component of the flow (the turbulent eddies). At the other end of the scale DNS modeling includes all motions down to the viscous scales. Boersma et al. have built such a model for the non-buoyant case which also compares well with measured data for mean and turbulent velocity components. The model demonstrates its versatility by application to a laminar flow case. As its name implies, DNS directly models the Navier-Stokes equations without recourse to subgrid modeling so for flows with a broad spectrum of motions (high Re) the cost can be prohibitive - the number of required grid points scaling with Re(exp 9/4) and the number of time steps with Re(exp 3/4). The middle road is provided by LES whereby the Navier-Stokes equations are formally filtered with the filter chosen to only exclude the smallest turbulent motions. If successful, LES should provide much of the detail available to DNS but at more bearable cost. Fatica et al. in comparing LES with DNS for a low Reynolds number jet showed that the LES could simulate the temporally evolving behavior including growth of the jet thickness. It is the intention of this report to explore the application of an LES model to jets and plumes. As always, before tackling complex situations, the model must be tested for the simplest of cases and so we address only two, a non-buoyant axisymmetric jet issuing steadily from an orifice into a semi-infinite stationary environment and a buoyant jet in the same environment. The work is a continuation of Basu and Mansour.

Webb, A. T.↗
NASA NTRS2003

Development of Erosive Burning Models for CFD Predictions of Solid Rocket Motor Internal Environments

Four erosive burning models, equations (11) to (14). are developed in this work by using a power law relationship to correlate (1) the erosive burning ratio and the local velocity gradient at propellant surfaces; (2) the erosive burning ratio and the velocity gradient divided by centerline velocity; (3) the erosive burning difference and the local velocity gradient at propellant surfaces; and (4) the erosive burning difference and the velocity gradient divided by centerline velocity. These models depend on the local velocity gradient at the propellant surface (or the velocity gradient divided by centerline velocity) only and, unlike other empirical models, are independent of the motor size. It was argued that, since the erosive burning is a local phenomenon occurring near the surface of the solid propellant, the erosive burning ratio should be independent of the bore diameter if it is correlated with some local flow parameters such as the velocity gradient at the propellant surface. This seems to be true considering the good results obtained by applying these models, which are developed from the small size 5 inch CP tandem motor testing, to CFD simulations of much bigger motors.

Wang, Qun-Zhen↗
DOE OSTI2025

Dark Energy Survey Year 6 Results: improved mitigation of spatially varying observational systematics with masking

As photometric surveys reach unprecedented statistical precision, systematic uncertainties increasingly dominate large-scale structure probes relying on galaxy number density. Defining the final survey footprint is critical, as it excludes regions affected by artefacts or suboptimal observing conditions. For galaxy clustering, spatially varying observational systematics, such as seeing, are a leading source of bias. Template maps of contaminants are used to derive spatially dependent corrections, but extreme values may fall outside the applicability range of mitigation methods, compromising correction reliability. The complexity and accuracy of systematics modelling depend on footprint conservativeness, with aggressive masking enabling simpler, robust mitigation. We present a unified approach to define the DES Year 6 joint footprint, integrating observational systematics templates and artefact indicators that degrade mitigation performance. This removes extreme values from an initial seed footprint, leading to the final joint footprint. By evaluating the DES Year 6 lens sample MagLim++ plus plus on this footprint, we enhance the Iterative Systematics Decontamination (ISD) method, detecting non-linear systematic contamination and improving correction accuracy. While the mask's impact on clustering is less significant than systematics decontamination, it remains non-negligible, comparable to statistical uncertainties in certain w(theta) scales and redshift bins. Supporting coherent analyses of galaxy clustering and cosmic shear, the final footprint spans 4031.04 deg2, setting the basis for DES Year 6 1x2pt, 2x2pt, and 3x2pt analyses. This work highlights how targeted masking strategies optimise the balance between statistical power and systematic control in Stage-III and -IV surveys.

Rodríguez-Monroy, M. [Madrid, IFT; IJCLab, Orsay]↗
DOE OSTI2026

Model-Agnostic Signal Discovery with Machine Learning: Bridging the Gap Between Theory and Practice

Searches for new phenomena in complex scientific data are predominantly model-dependent, optimized for specific hypotheses, and therefore limited in their coverage of the space of possible signals. Recently, new AI-based model-agnostic search strategies, many of which have been pioneered in high-energy physics, have been proposed which provide a complementary paradigm, prioritizing broad exploration over tailored analyses. These techniques offer an opportunity to enhance the overall discovery potential of modern experiments, especially in regimes where theoretical guidance is scarce. In this document, we review the conceptual framework behind the main classes of AI-based model-agnostic strategies. We discuss the potential pitfalls of these methods, and strategies for their validation and interpretation. We aim for this document to serve as a useful reference both for practitioners and for researchers interested in learning more about these model-agnostic search strategies.

Amram, Oz [Fermilab] (ORCID:0000000237653123)↗
NASA NTRS1978

Release of potential instability by mesoscale triggering - An objective model simulation

The effects of mesoscale triggering on organized nonsevere convective cloud systems in the High Plains are considered. Two experiments were conducted to determine if a one-dimensional quasi-time dependent model could (1) detect soundings which were sensitive to mesoscale triggering, and (2) discriminate between cases which had mesoscale organized convection and those with no organized convection. The MESOCU model was used to analyze the available potential instability and thermodynamic potential for cloud growth. It is noted that lifting is a key factor in the release of available potential instability on the High Plains.

Matthews, D. A.↗
DOE OSTI2021

Constraints on the nuclear symmetry energy from asymmetric-matter calculations with chiral 𝑁⁢𝑁 and 3⁢𝑁 interactions

The nuclear symmetry energy is a key quantity in nuclear (astro)physics. It describes the isospin dependence of the nuclear equation of state, which is commonly assumed to be almost quadratic. Here, in this work, we confront this standard quadratic expansion of the equation of state with explicit asymmetric nuclear-matter calculations based on a set of commonly used Hamiltonians including two- and three-nucleon forces derived from chiral effective-field theory. We study, in particular, the importance of nonquadratic contributions to the symmetry energy, including the nonanalytic logarithmic term introduced by Kaiser [Phys. Rev. C 91, 065201 (2015)]. Our results suggest that the nonquadratic contribution to the symmetry energy can be systematically determined from the various Hamiltonians employed, and we obtain 0.74$^{+0.11}_{−0.08}$ MeV (or −1.02$^{+0.11}_{−0.08}$ MeV for the potential term with the effective-mass contribution) at nuclear saturation density, while the logarithmic contribution to the symmetry energy is relatively small and model-dependent. We also employ the meta-model approach to study the impact of the higher-order contributions on the neutron-star crust-core transition density, and find a 5% correction.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Strong gravity extruding peaks in speed of sound profiles of massive neutron stars

The speed of sound squared (SSS) s 2 in massive neutron stars (NSs) characterizes not only the stiffness of supradense neutron-rich matter within but also equivalently properties of the curved geometry due to the strong-field gravity and matter-geometry coupling. A peaked density or radius profile of s 2 has been predicted for massive NSs using various NS equation of state (EOS) models. However, the nature, cause, location and size of the peak in s 2 profiles are still very EOS model dependent. In this work, we investigate systematically s 2 profiles in massive NSs in a new approach that is independent of the nuclear EOS model and without any presumption about the NS structure and/or composition. In terms of the small quantities (reduced radius, the energy density and pressure scaled by their central values), we perform double-element perturbative expansions in solving perturbatively the scaled Tolman-Oppenheimer-Volkoff (TOV) equations and analyzing s 2 profiles from the Newtonian limit to the general relativistic (GR) case. The GR term in the TOV equations plays a twofold role: it compresses NS matter and modifies the pressure/energy density ratio from small values in Newtonian stars showing no s 2 peak to large ones for massive NSs possessing a peak in their s 2 profiles, and eventually takes away the peak in extremely compact/massive NSs approaching the causality limit. In particular, the peaked behavior in s 2 is expected to emerge near the center of massive NSs like PSR J⁢0740+6620, while a sharp phase transition is unlikely to occur there. Furthermore, these features revealed from our analyses are universal as they are intrinsic properties of the GR stellar structure equations independent of the still very uncertain EOS of supradense neutron-rich matter in NSs.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1987

A pulsar emission model - Observational tests

The usual model for coherent radio emission, curvature radiation from bunches, is applied here to a modified 'hollow cone' model in which only one or a few geometrically fixed flux tubes are radiating into a fan beam. The way that the model handles the basic problems of sufficient radio power output, reasonable spectrum, and narrow pulses as well as the secondary issues of orthogonal mode changing, complex pulse profiles is addressed. Issues that seem more model-dependent, including flux tube scales, nulling, self-absorption, drifting subpulses, and radius-to-frequency mapping, are examined by considering a specific global model.

Michel, F. Curtis↗
NASA NTRS1989

Numerical study of the current sheet and PSBL in a magnetotail model

The current sheet and plasma sheet boundary layer (PSBL) in a magnetotail model are discussed. A test particle code is used to study the response of ensembles of particles to a two-dimensional, time-dependent model of the geomagnetic tail, and test the proposition (Coroniti, 1985a, b; Buchner and Zelenyi, 1986; Chen and Palmadesso, 1986; Martin, 1986) that the stochasticity of the particle orbits in these fields is an important part of the physical mechanism for magnetospheric substorms. The realistic results obtained for the fluid moments of the particle distribution with this simple model, and their insensitivity to initial conditions, is consistent with this hypothesis.

Doxas, I.↗
DOE OSTI2022

Universal trend of charge radii of even-even Ca–Zn nuclei

Radii of nuclear charge distributions carry information about the strong and electromagnetic forces acting inside the atomic nucleus. Whereas the global behavior of nuclear charge radii is governed by the bulk properties of nuclear matter, their local trends are affected by quantum motion of proton and neutron nuclear constituents. The measured differential charge radii δ$\langle$$r^{2}_{c}$$\rangle$ between neutron numbers N = 28 and N = 40 exhibit a universal pattern as a function of n = N – 28 that is independent of the atomic number. In this work, we analyze this remarkable behavior in even-even nuclei from calcium to zinc using two state-of-the-art theories based on quantified nuclear interactions: the ab initio coupled cluster theory and nuclear density functional theory. Both theories reproduce the smooth rise of differential charge radii and their weak dependence on the atomic number. By considering a large set of isotopic chains, we show that this trend can be captured by just two parameters: the slope and curvature of δ$\langle$$r^{2}_{c}$$\rangle$(n). We demonstrate that these parameters show appreciable model dependence, and the statistical analysis indicates that they are not correlated with any single model property, i.e., they are impacted by both bulk nuclear properties as well as shell structure.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗