Search NASA⌕ Search

SEARCH · Search NASA

Results for “Closure problem”

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 253 records · Page 14

Lunar and Planetary Science XXXV: Terrestrial Planets: Building Blocks and Differentiation

The session "Terrestrial Planets: Building Blocks and Differentiation: included the following topics:Magnesium Isotopes in the Earth, Moon, Mars, and Pallasite Parent Body: High-Precision Analysis of Olivine by Laser-Ablation Multi-Collector ICPMS; Meteoritic Constraints on Collision Rates in the Primordial Asteroid Belt and Its Origin; New Constraints on the Origin of the Highly Siderophile Elements in the Earth's Upper Mantle; Further Lu-Hf and Sm-Nd Isotopic Data on Planetary Materials and Consequences for Planetary Differentiation; A Deep Lunar Magma Ocean Based on Neodymium, Strontium and Hafnium Isotope Mass Balance Partial Resetting on Hf-W System by Giant Impacts; On the Problem of Metal-Silicate Equilibration During Planet Formation: Significance for Hf-W Chronometry ; Solid Metal-Liquid Metal Partitioning of Pt, Re, and Os: The Effect of Carbon; Siderophile Element Abundances in Fe-S-Ni-O Melts Segregated from Partially Molten Ordinary Chondrite Under Dynamic Conditions; Activity Coefficients of Silicon in Iron-Nickel Alloys: Experimental Determination and Relevance for Planetary Differentiation; Reinvestigation of the Ni and Co Metal-Silicate Partitioning; Metal/Silicate Paritioning of P, Ga, and W at High Pressures and Temperatures: Dependence on Silicate Melt Composition; and Closure of the Fe-S-Si Liquid Miscibility Gap at High Pressure and Its Implications for Planetary Core Formation.

Source record↗

Mixed Layer Sub-Mesoscale Parameterization - Part 1: Derivation and Assessment

Several studies have shown that sub-mesoscales (SM 1km horizontal scale) play an important role in mixed layer dynamics. In particular, high resolution simulations have shown that in the case of strong down-front wind, the re-stratification induced by the SM is of the same order of the de-stratification induced by small scale turbulence, as well as of that induced by the Ekman velocity. These studies have further concluded that it has become necessary to include SM in ocean global circulation models (OGCMs), especially those used in climate studies. The goal of our work is to derive and assess an analytic parameterization of the vertical tracer flux under baroclinic instabilities and wind of arbitrary directions and strength. To achieve this goal, we have divided the problem into two parts: first, in this work we derive and assess a parameterization of the SM vertical flux of an arbitrary tracer for ocean codes that resolve mesoscales, M, but not sub-mesoscales, SM. In Part 2, presented elsewhere, we have used the results of this work to derive a parameterization of SM fluxes for ocean codes that do not resolve either M or SM. To carry out the first part of our work, we solve the SM dynamic equations including the non-linear terms for which we employ a closure developed and assessed in previous work. We present a detailed analysis for down-front and up-front winds with the following results: (a) down-front wind (blowing in the direction of the surface geostrophic velocity) is the most favorable condition for generating vigorous SM eddies; the de-stratifying effect of the mean flow and re-stratifying effect of SM almost cancel each other out,

Canuto, V. M.↗

On the consistency of Reynolds stress turbulence closures with hydrodynamic stability theory

The consistency of second-order closure models with results from hydrodynamic stability theory is analyzed for the simplified case of homogeneous turbulence. In a recent study, Speziale, Gatski, and MacGiolla Mhuiris showed that second-order closures are capable of yielding results that are consistent with hydrodynamic stability theory for the case of homogeneous shear flow in a rotating frame. It is demonstrated in this paper that this success is due to the fact that the stability boundaries for rotating homogeneous shear flow are not dependent on the details of the spatial structure of the disturbances. For those instances where they are -- such as in the case of elliptical flows where the instability mechanism is more subtle -- the results are not so favorable. The origins and extent of this modeling problem are examined in detail along with a possible resolution based on rapid distortion theory (RDT) and its implications for turbulence modeling.

Speziale, Charles G.↗

The Behavior of Filters and Smoothers for Strongly Nonlinear Dynamics

The Kalman filter is the optimal filter in the presence of known Gaussian error statistics and linear dynamics. Filter extension to nonlinear dynamics is non trivial in the sense of appropriately representing high order moments of the statistics. Monte Carlo, ensemble-based, methods have been advocated as the methodology for representing high order moments without any questionable closure assumptions (e.g., Miller 1994). Investigation along these lines has been conducted for highly idealized dynamics such as the strongly nonlinear Lorenz (1963) model as well as more realistic models of the oceans (Evensen and van Leeuwen 1996) and atmosphere (Houtekamer and Mitchell 1998). A few relevant issues in this context are related to the necessary number of ensemble members to properly represent the error statistics and, the necessary modifications in the usual filter equations to allow for correct update of the ensemble members (Burgers 1998). The ensemble technique has also been applied to the problem of smoothing for which similar questions apply. Ensemble smoother examples, however, seem to quite puzzling in that results of state estimate are worse than for their filter analogue (Evensen 1997). In this study, we use concepts in probability theory to revisit the ensemble methodology for filtering and smoothing in data assimilation. We use Lorenz (1963) model to test and compare the behavior of a variety implementations of ensemble filters. We also implement ensemble smoothers that are able to perform better than their filter counterparts. A discussion of feasibility of these techniques to large data assimilation problems will be given at the time of the conference.

Zhu, Yanqiu↗

Integrated Crewed Mars Mission Analysis, Part II: The Trajectory Strikes Back

Crewed missions to Mars present a challenging problem involving strong interactions between a variety of spacecraft systems. Design and architecture decisions made by one system can have powerful implications for another. Siloed subsystem development and optimization efforts are unlikely to generate optimal results for the integrated system and may lead to counterproductive investments at the mission and system level. Mismatched assumptions between different subsystems can also negatively impact the final result. Integrated systems analysis, trade studies, and sensitivity analysis are a critical element of successful Mars mission design. Previous work examined the linkages between Environmental Control and Life Support System (ECLSS) loop closure, food water content, spares Probability of Sufficiency (POS), and other factors by performing an integrated mass and cost assessment of a variety of Mars habitat and transportation system designs. This paper expands upon that work by evaluating the impact of mission duration, trajectory, and propulsion technology alongside key habitat variables. Transportation system and trajectory design decisions have enormous impacts on overall system mass, and the impacts of those decisions are themselves sensitive to the mass of the habitat that is being transported. This paper examines a variety of mission cases and evaluates total departure mass for each case. An estimate of relative costs for specific cases is also provided. The results are discussed, with an emphasis on the sensitivity of mass, cost, and other figures of merit to key parameters. Interactions and dependencies between various decisions, particularly those that relate to technology investments, are also highlighted. Overall, the various elements, systems, and subsystems required to support human spaceflight should not be optimized and evaluated in isolation, but instead should be analyzed as part of an integrated system. This paper demonstrates how integrated analysis can improve system understanding and lead to more optimal and effective crewed exploration systems.

Human Spaceflight↗

Integrated Crewed Mars Mission Analysis, Part II: The Trajectory Strikes Back

Crewed missions to Mars present a challenging problem involving strong interactions between a variety of spacecraft systems. Design and architecture decisions made by one system can have powerful implications for another. Siloed subsystem development and optimization efforts are unlikely to generate optimal results for the integrated system and may lead to counterproductive investments at the mission and system level. Mismatched assumptions between different subsystems can also negatively impact the final result. Integrated systems analysis, trade studies, and sensitivity analysis are a critical element of successful Mars mission design. Previous work examined the linkages between Environmental Control and Life Support System (ECLSS) loop closure, food water content, spares Probability of Sufficiency (POS), and other factors by performing an integrated mass and cost assessment of a variety of Mars habitat and transportation system designs. This paper expands upon that work by evaluating the impact of mission duration, trajectory, and propulsion technology alongside key habitat variables. Transportation system and trajectory design decisions have enormous impacts on overall system mass, and the impacts of those decisions are themselves sensitive to the mass of the habitat that is being transported. This paper examines a variety of mission cases and evaluates total departure mass for each case. An estimate of relative costs for specific cases is also provided. The results are discussed, with an emphasis on the sensitivity of mass, cost, and other figures of merit to key parameters. Interactions and dependencies between various decisions, particularly those that relate to technology investments, are also highlighted. Overall, the various elements, systems, and subsystems required to support human spaceflight should not be optimized and evaluated in isolation, but instead should be analyzed as part of an integrated system. This paper demonstrates how integrated analysis can improve system understanding and lead to more optimal and effective crewed exploration systems.

Human Spaceflight↗

Parallel expansion of a fuel pellet plasmoid

The problem of the assimilation of a cryogenic fuel pellet injected into a hot plasma is considered. Due to the transparency to ambient particles of the plasmoid, the localised region of high-density plasma created by ionisation of the ablated pellet material, electrons reach a ‘quasiequilibrium’ (QE) state which is characterised by a steady-state on the fastest collisional time scale. The simplified electron kinetic equation of the QE state is solved. Taking a velocity moment of the higher-order electron kinetic equation, which is valid on the expansion time scale, permits a fluid closure, yielding an evolution equation for the macroscopic parameters describing the QE distribution function. In contrast to the Braginskii equations, the closure does not require that electrons have a short mean free path compared with the size of density perturbations, and permits an anisotropic and highly non-Maxwellian distribution function. As the QE distribution function accounts for both trapped and passing electrons, the self-consistent electric potential that causes the expansion can be properly described, in contrast to earlier models of pellet plasmoid expansion with an unbounded potential. The plasmoid expansion is simulated using both a Vlasov model and a cold-fluid model for the ions. During the expansion plasmoid ions and electrons obtain nearly equal amounts of energy; as hot ambient electrons provide this energy in the form of collisional heating of plasmoid electrons, the expansion of a pellet plasmoid is expected to be a potent mechanism for the transfer of energy from electrons to ions on a time scale shorter than that of ion–electron thermalisation.

Physics↗

Alternative Representations of Convective Processes in the NASA GEOS-5 AGCM

The gap in explicit resolution of phenomena between global climate models and cloud resolving models is shrinking at a steady pace with global integrations of several km in resolution now practical for at least some time scales. In moving toward finer resolution the long standing problem of convective parameterization is being examined along with the assumption of convective quasi-equilibrium and how this can be reconciled with the stochastic and intermittent nature of convection. In this context we examine the nature of parameterized convection in the NASA Goddard Earth Observing System (GEOS-5) Atmospheric General Circulation Model. Our analysis uses both coarse (2.5 degree) and fine scale 0.25 degree spatial resolution integrations. Two basic formulations are compared: the default option is the Relaxed Arakawa-Schubert (RAS) scheme which invokes a sequence of linearly entraining plumes and quasi-equilibrium closure. An optional modification of this method (the "Stochastic Tokioka" constraint) places a random lower limit on plume entrainment. An alternative representation is the Kain-Fritsch parameterization which was originally developed for mesoscale numerical modeling strategies. Here entrainment is determined by a crude buoyancy sorting approach that allows the plume spectrum to be more responsive to ambient vertical stratification of moisture. Diagnostics of the model behavior are referenced to recent observational evidence of continuous phase transition behavior. In particular we examine the relationship between column water vapor and probablility of convective presence and intensity. Sensitivity of the statistics of convective behavior to parcel mixing/entrainment formulations and parcel initial thermodynamics are considered. Observational statistics from A-Train and TRMM sensors provide validation of the model integrations.

Robertson, Franklin↗

Investigation of Navier-Stokes Code Verification and Design Optimization

With rapid progress made in employing computational techniques for various complex Navier-Stokes fluid flow problems, design optimization problems traditionally based on empirical formulations and experiments are now being addressed with the aid of computational fluid dynamics (CFD). To be able to carry out an effective CFD-based optimization study, it is essential that the uncertainty and appropriate confidence limits of the CFD solutions be quantified over the chosen design space. The present dissertation investigates the issues related to code verification, surrogate model-based optimization and sensitivity evaluation. For Navier-Stokes (NS) CFD code verification a least square extrapolation (LSE) method is assessed. This method projects numerically computed NS solutions from multiple, coarser base grids onto a freer grid and improves solution accuracy by minimizing the residual of the discretized NS equations over the projected grid. In this dissertation, the finite volume (FV) formulation is focused on. The interplay between the xi concepts and the outcome of LSE, and the effects of solution gradients and singularities, nonlinear physics, and coupling of flow variables on the effectiveness of LSE are investigated. A CFD-based design optimization of a single element liquid rocket injector is conducted with surrogate models developed using response surface methodology (RSM) based on CFD solutions. The computational model consists of the NS equations, finite rate chemistry, and the k-6 turbulence closure. With the aid of these surrogate models, sensitivity and trade-off analyses are carried out for the injector design whose geometry (hydrogen flow angle, hydrogen and oxygen flow areas and oxygen post tip thickness) is optimized to attain desirable goals in performance (combustion length) and life/survivability (the maximum temperatures on the oxidizer post tip and injector face and a combustion chamber wall temperature). A preliminary multi-objective optimization study is carried out using a geometric mean approach. Following this, sensitivity analyses with the aid of variance-based non-parametric approach and partial correlation coefficients are conducted using data available from surrogate models of the objectives and the multi-objective optima to identify the contribution of the design variables to the objective variability and to analyze the variability of the design variables and the objectives. In summary the present dissertation offers insight into an improved coarse to fine grid extrapolation technique for Navier-Stokes computations and also suggests tools for a designer to conduct design optimization study and related sensitivity analyses for a given design problem.

Vaidyanathan, Rajkumar↗

A case study of view-factor rectification procedures for diffuse-gray radiation enclosure computations

The view factors which are used in diffuse-gray radiation enclosure calculations are often computed by approximate numerical integrations. These approximately calculated view factors will usually not satisfy the important physical constraints of reciprocity and closure. In this paper several view-factor rectification algorithms are reviewed and a rectification algorithm based on a least-squares numerical filtering scheme is proposed with both weighted and unweighted classes. A Monte-Carlo investigation is undertaken to study the propagation of view-factor and surface-area uncertainties into the heat transfer results of the diffuse-gray enclosure calculations. It is found that the weighted least-squares algorithm is vastly superior to the other rectification schemes for the reduction of the heat-flux sensitivities to view-factor uncertainties. In a sample problem, which has proven to be very sensitive to uncertainties in view factor, the heat transfer calculations with weighted least-squares rectified view factors are very good with an original view-factor matrix computed to only one-digit accuracy. All of the algorithms had roughly equivalent effects on the reduction in sensitivity to area uncertainty in this case study.

Taylor, Robert P.↗

Quantifying the Impacts of PM2.5 Constituents and Relative Humidity on Visibility Impairment in a Suburban Area of Eastern Asia Using Long-Term In-Situ Measurements.

The deterioration of visibility due to air pollutants and relative humidity has been a serious environmental problem in eastern Asia. In most previous studies, chemical compositions of atmospheric particles were provided using filter-based offline analyses, which were unable to provide long-term and in-situ measurements that resolve sufficient temporal variations of air pollution and meteorology, hindering the resolution of the relationship between air quality and visibility. Here, we present a year-long continuously measured data from a comprehensive suite of online instruments to investigate diurnal and seasonal impacts of the aerosol chemical compositions in PM2.5 on visibility seasonally and diurnally. The measured dry aerosol extinction at λ = 550 nm reached a closure with that predicted by aerosol compositions within 12%. However, the hygroscopic growth of particles under ambient RH could enhance the aerosol extinction by a factor of 2 – 6, matching the perceptive visibility of the public. Particulate ammonium nitrate was most sensitive to reducing visibility, while ammonium sulfate contributed the most to the light extinction. In spring and winter, the monsoon and stagnant air masses reduced the visibility and increased PM2.5 (> 35 μg m-3).The moisture was found to substantially enhance the light extinction under RH = 60 – 90%,reducing visibility by approximately 15 km, largely attributed to hygroscopic inorganic salts.This study serves as a metric to highlight the need to consider the influence of RH, and aqueous reactions in producing secondary inorganic aerosols on atmospheric visibility, underpinning the more accurate mitigation strategies of air pollution.

visibility↗

A robust rotorcraft flight control system design methodology utilizing quantitative feedback theory

Rotorcraft flight control systems present design challenges which often exceed those associated with fixed-wing aircraft. First, large variations in the response characteristics of the rotorcraft result from the wide range of airspeeds of typical operation (hover to over 100 kts). Second, the assumption of vehicle rigidity often employed in the design of fixed-wing flight control systems is rarely justified in rotorcraft where rotor degrees of freedom can have a significant impact on the system performance and stability. This research was intended to develop a methodology for the design of robust rotorcraft flight control systems. Quantitative Feedback Theory (QFT) was chosen as the basis for the investigation. Quantitative Feedback Theory is a technique which accounts for variability in the dynamic response of the controlled element in the design robust control systems. It was developed to address a Multiple-Input Single-Output (MISO) design problem, and utilizes two degrees of freedom to satisfy the design criteria. Two techniques were examined for extending the QFT MISO technique to the design of a Multiple-Input-Multiple-Output (MIMO) flight control system (FCS) for a UH-60 Black Hawk Helicopter. In the first, a set of MISO systems, mathematically equivalent to the MIMO system, was determined. QFT was applied to each member of the set simultaneously. In the second, the same set of equivalent MISO systems were analyzed sequentially, with closed loop response information from each loop utilized in subsequent MISO designs. The results of each technique were compared, and the advantages of the second, termed Sequential Loop Closure, were clearly evident.

Gorder, Peter James↗

Next-Generation MKIII Lightweight HUT/Hatch Assembly

The MK III (H-1) carbon-graphite/ epoxy Hard Upper Torso (HUT)/Hatch assembly was designed, fabricated, and tested in the early 1990s. The spacesuit represented an 8.3 psi (≈58 kPa) technology demonstrator model of a zero prebreathe suit. The basic torso shell, brief, and hip areas of the suit were composed of a carbon-graphite/epoxy composite lay-up. In its current configuration, the suit weighs approximately 120 lb (≈54 kg). However, since future planetary suits will be designed to operate at 0.26 bar (≈26 kPa), it was felt that the suit's re-designed weight could be reduced to 79 lb (≈35 kg) with the incorporation of lightweight structural materials. Many robust, lightweight structures based on the technologies of advanced honeycomb materials, revolutionary new composite laminates, metal matrix composites, and recent breakthroughs in fullerene fillers and nanotechnology lend themselves well to applications requiring materials that are both light and strong. The major problem involves the reduction in weight of the HUT/ Hatch assembly for use in lunar and/or planetary applications, while at the same time maintaining a robust structural design. The technical objective is to research, design, and develop manufacturing methods that support fa b rica - tion of a lightweight HUT/Hatch assembly using advanced material and geometric redesign as necessary. Additionally, the lightweight HUT/Hatch assembly will interface directly with current MK III hardware. Using the new operating pressure and current MK III (H-1) interfaces as a starting block, it is planned to maximize HUT/Hatch assembly weight reduction through material selection and geometric redesign. A hard upper torso shell structure with rear-entry closure and corresponding hatch will be fabricated. The lightweight HUT/Hatch assembly will retrofit and interface with existing MK III (H-1) hardware elements, providing NASA with immediate "plug-andplay" capability. NASA crewmembers will have a lightweight, robust, life-support system that will minimize fatigue during extraterrestrial surface sojourns. Its unique feature is the utilization of a new and innovative family of materials used by the aerospace industry, which at the time of this reporting has not been used for the proposed application.

McCarthy, Mike↗

Computational Investigation of Fluidic Counterflow Thrust Vectoring

A computational study of fluidic counterflow thrust vectoring has been conducted. Two-dimensional numerical simulations were run using the computational fluid dynamics code PAB3D with two-equation turbulence closure and linear Reynolds stress modeling. For validation, computational results were compared to experimental data obtained at the NASA Langley Jet Exit Test Facility. In general, computational results were in good agreement with experimental performance data, indicating that efficient thrust vectoring can be obtained with low secondary flow requirements (less than 1% of the primary flow). An examination of the computational flowfield has revealed new details about the generation of a countercurrent shear layer, its relation to secondary suction, and its role in thrust vectoring. In addition to providing new information about the physics of counterflow thrust vectoring, this work appears to be the first documented attempt to simulate the counterflow thrust vectoring problem using computational fluid dynamics.

Hunter, Craig A.↗

New approaches to virtual environment surgery

This research focused on two main problems: 1) low cost, high fidelity stereoscopic imaging of complex tissues and organs; and 2) virtual cutting of tissue. A further objective was to develop these images and virtual tissue cutting methods for use in a telemedicine project that would connect remote sites using the Next Generation Internet. For goal one we used a CT scan of a human heart, a desktop PC with an OpenGL graphics accelerator card, and LCD stereoscopic glasses. Use of multiresolution meshes ranging from approximately 1,000,000 to 20,000 polygons speeded interactive rendering rates enormously while retaining general topography of the dataset. For goal two, we used a CT scan of an infant skull with premature closure of the right coronal suture, a Silicon Graphics Onyx workstation, a Fakespace Immersive WorkBench and CrystalEyes LCD glasses. The high fidelity mesh of the skull was reduced from one million to 50,000 polygons. The cut path was automatically calculated as the shortest distance along the mesh between a small number of hand selected vertices. The region outlined by the cut path was then separated from the skull and translated/rotated to assume a new position. The results indicate that widespread high fidelity imaging in virtual environment is possible using ordinary PC capabilities if appropriate mesh reduction methods are employed. The software cutting tool is applicable to heart and other organs for surgery planning, for training surgeons in a virtual environment, and for telemedicine purposes.

NASA Discipline Neuroscience↗

Heterogeneous Mixtures of Dictionary Functions to Approximate Subspace Invariance in Koopman Operators: Why Deep Koopman Operators Work

Abstract Koopman operators model nonlinear dynamics as a linear dynamic system acting on a nonlinear function as the state. This nonstandard state is often called a Koopman observable and is usually approximated numerically by a superposition of functions drawn from a dictionary . In a widely used algorithm, extended dynamic mode decomposition (EDMD), the dictionary functions are drawn from a fixed class of functions. Deep learning combined with EDMD has been used to learn novel dictionary functions in an algorithm called deep dynamic mode decomposition (deepDMD). The learned representation both (1) accurately models and (2) scales well with the dimension of the original nonlinear system. In this paper, we analyze the learned dictionaries from deepDMD and explore the theoretical basis for their strong performance. We explore State-Inclusive Logistic Lifting (SILL) dictionary functions to approximate Koopman observables. Error analysis of these dictionary functions show they satisfy a property of subspace approximation, which we define as uniform finite approximate closure. Typically, a Koopman dictionary’s nonlinear functions are homogeneous. In this paper, we discover that structured mixing of heterogeneous dictionary functions drawn from different classes of nonlinear functions achieve the same accuracy and dimensional scaling as the deep-learning-based deepDMD algorithm Yeung et al. ( In: 2019 American Control Conference (ACC), 2019). We specifically show this by building a heterogeneous dictionary comprised of SILL functions and conjunctive radial basis functions (RBFs). This mixed dictionary achieves similar accuracy and dimensional scaling to deepDMD with an order of magnitude reduction in parameters, while maintaining geometric interpretability. These results strengthen the viability of dictionary-based Koopman models to solving high-dimensional nonlinear learning problems.

Johnson, Charles A.↗

thornado+FLASH-X: A Hybrid Discontinuous Galerkin–Implicit-explicit and Finite-volume Framework for Neutrino-radiation Hydrodynamics in Core-collapse Supernovae

We present neutrino-transport algorithms implemented in the toolkit for high-order neutrino-radiation hydrodynamics (thornado) and their coupling to self-gravitating hydrodynamics within the adaptive mesh refinement–based multiphysics simulation framework FLASH-X. thornado, developed primarily for simulations of core-collapse supernovae (CCSNe), employs a spectral, six-species two-moment formulation with algebraic closure and special-relativistic observer corrections accurate to $\mathcal{O}(v/c)$, and uses discontinuous Galerkin (DG) methods for phase-space discretization combined with implicit-explicit time stepping. A key development is a nonlinear neutrino–matter coupling algorithm based on nested fixed-point iteration with Anderson acceleration, enabling fully implicit treatment of collisional processes, including energy-coupling interactions such as neutrino–electron scattering and pair production. Coupling to finite-volume (FV) hydrodynamics is achieved through a hybrid DG-FV representation of the fluid variables and operator-split evolution within FLASH-X. The implementation is verified using basic transport tests with idealized opacities and relaxation and deleptonization problems with tabulated microphysics. Spherically symmetric CCSN simulations demonstrate accuracy and robustness of the coupled scheme, including close agreement with the CCSN simulation code Chimera. An axisymmetric CCSN simulation further demonstrates the viability of DG-based neutrino transport for multidimensional supernova modeling within FLASH-X. thornado’s neutrino-transport solver is GPU-enabled using OpenMP offloading or OpenACC, and all CCSN applications included in this work use the GPU implementation. Together, these results establish a foundation for future enhancements in physics fidelity, numerical algorithms, and computational performance, for increasingly realistic large-scale CCSN simulations.

Endeve, Eirik [Oak Ridge National Laboratory (ORNL↗

Software Fault Tolerance: A Tutorial

Because of our present inability to produce error-free software, software fault tolerance is and will continue to be an important consideration in software systems. The root cause of software design errors is the complexity of the systems. Compounding the problems in building correct software is the difficulty in assessing the correctness of software for highly complex systems. After a brief overview of the software development processes, we note how hard-to-detect design faults are likely to be introduced during development and how software faults tend to be state-dependent and activated by particular input sequences. Although component reliability is an important quality measure for system level analysis, software reliability is hard to characterize and the use of post-verification reliability estimates remains a controversial issue. For some applications software safety is more important than reliability, and fault tolerance techniques used in those applications are aimed at preventing catastrophes. Single version software fault tolerance techniques discussed include system structuring and closure, atomic actions, inline fault detection, exception handling, and others. Multiversion techniques are based on the assumption that software built differently should fail differently and thus, if one of the redundant versions fails, it is expected that at least one of the other versions will provide an acceptable output. Recovery blocks, N-version programming, and other multiversion techniques are reviewed.

Torres-Pomales, Wilfredo↗