Search NASA⌕ Search

SEARCH · Search NASA

Results for “Differentiable Programming”

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 1,225 records · Page 68

Monitoring Changes of Tropical Extreme Rainfall Events Using Differential Absorption Barometric Radar (DiBAR)

This work studies the potential of monitoring changes in tropical extreme rainfall events such as tropical storms from space using a Differential-absorption BArometric Radar (DiBAR) operating at 50-55 gigahertz O2 absorption band to remotely measure sea surface air pressure. Air pressure is among the most important variables that affect atmospheric dynamics, and currently can only be measured by limited in-situ observations over oceans. Analyses show that with the proposed radar the errors in instantaneous (averaged) pressure estimates can be as low as approximately 5 millibars (approximately 1 millibar) under all weather conditions. With these sea level pressure measurements, the forecasts, analyses and understanding of these extreme events in both short and long time scales can be improved. Severe weathers, especially hurricanes, are listed as one of core areas that need improved observations and predictions in WCRP (World Climate Research Program) and NASA Decadal Survey (DS) and have major impacts on public safety and national security through disaster mitigation. Since the development of the DiBAR concept about a decade ago, our team has made substantial progress in advancing the concept. Our feasibility assessment clearly shows the potential of sea surface barometry using existing radar technologies. We have developed a DiBAR system design, fabricated a Prototype-DiBAR (P-DiBAR) for proof-of-concept, conducted lab, ground and airborne P-DiBAR tests. The flight test results are consistent with our instrumentation goals. Observational system simulation experiments for space DiBAR performance show substantial improvements in tropical storm predictions, not only for the hurricane track and position but also for the hurricane intensity. DiBAR measurements will lead us to an unprecedented level of the prediction and knowledge on tropical extreme rainfall weather and climate conditions.

Lin, Bing↗

Landsat 8: The Plans, the Reality, and the Legacy

Landsat 8, originally known as the Landsat Data Continuity Mission (LDCM), is a National Aeronautics and Space Administration (NASA)-U.S. Geological Survey (USGS) partnership that continues the legacy of continuous moderate resolution observations started in 1972. The conception of LDCM to the reality of Landsat 8 followed an arduous path extending over nearly 13 years, but the successful launch on February 11, 2013 ensures the continuity of the unparalleled Landsat record. The USGS took over mission operations on May 30, 2013 and renamed LCDM to Landsat 8. Access to Landsat 8 data was opened to users worldwide. Three years following launch we evaluate the science and applications impact of Landsat 8. With a mission objective to enable the detection and characterization of global land changes at a scale where differentiation between natural and human-induced causes of change is possible, LDCM promised incremental technical improvements in capabilities needed for Landsat scientific and applications investigations. Results show that with Landsat 8, we are acquiring more data than ever before, the radiometric and geometric quality of data are generally technically superior to data acquired by past Landsat missions, and the new measurements, e.g., the coastal aerosol and cirrus bands, are opening new opportunities. Collectively, these improvements are sparking the growth of science and applications opportunities. Equally important, with Landsat 7 still operational, we have returned to global imaging on an 8-day cycle, a capability that ended when Landsat 5 ceased operational Earth imaging in November 2011. As a result, the Landsat program is on secure footings and planning is underway to extend the record for another 20 or more years.

Loveland, Thomas R.↗

First Measurements of Differential Cross Sections In Kinematic Imbalance Variables With The MicroBooNE Detector

Making high-precision measurements of neutrino oscillation parameters requires an unprecedented understanding of neutrino-nucleus scattering. In this presentation, we present the first muon neutrino charged current double-differential cross sections in kinematic imbalance variables. These variables characterize the imbalance in the plane transverse to an incoming neutrino. We use events with a single muon above 100 MeV/c, a single final state proton above 300 MeV/c, and no recorded final state pions. Thus, these variables act as a direct probe of nuclear effects such as final state interactions, Fermi motion, and multi-nucleon processes. We also present a complementary ongoing analysis using electron neutrinos. This channel is of the utmost importance for the extraction of neutrino oscillation parameters by making high-precision measurements. Our measurements allow us to constrain systematic uncertainties associated with neutrino oscillation results performed by near-future experiments of the Short Baseline Neutrino (SBN) program, as well as by future large-scale experiments like DUNE.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

First Measurements of Differential Cross Sections In Kinematic Imbalance Variables With The MicroBooNE Detector

Making high-precision measurements of neutrino oscillation parameters requires an unprecedented understanding of neutrino-nucleus scattering. In this presentation, we present the first muon neutrino charged current double-differential cross sections in kinematic imbalance variables. These variables characterize the imbalance in the plane transverse to an incoming neutrino. We use events with a single muon above 100 MeV/c, a single final state proton above 300 MeV/c, and no recorded final state pions. Thus, these variables act as a direct probe of nuclear effects such as final state interactions, Fermi motion, and multi-nucleon processes. We also present a complementary ongoing analysis using electron neutrinos. This channel is of the utmost importance for the extraction of neutrino oscillation parameters by making high-precision measurements. Our measurements allow us to constrain systematic uncertainties associated with neutrino oscillation results performed by near-future experiments of the Short Baseline Neutrino (SBN) program, as well as by future large-scale experiments like DUNE.

43 PARTICLE ACCELERATORS↗

AN ACID BAKING APPROACH TO ENHANCE RARE EARTH ELEMENT RECOVERY FROM BITUMINOUS COAL SOURCES

Rare earth elements (REE) are a group of 17 elements typically classified as light and heavy rare earth elements, which play a crucial role in developing the latest technologies for energy, defense, and medical sectors. Even though REEs have been found in more than 200 minerals, only bastnaesite, monazite and xenotime are commercially exploited for REE extraction. However, the recent exponential increase in REE demand has spurred countries such as the United States into research for the extraction of REEs from secondary sources such as coal, acid mine drainage, and coal ash. Several coal sources (e.g., Fire Clay seam coal) across the United States have been identified to contain elevated concentrations of rare earth elements (>600 ppm), and various researchers have investigated the feasibility of both physical and hydrometallurgical extraction techniques for rare earth concentration and subsequent extraction. However, both physical and direct leaching were concluded to be inefficient for RE beneficiation and extraction due to low recoveries. Alternatively, thermal treatment provides promising means for RE recovery from bituminous coal sources. The positive impact of thermal treatment/calcination was established to be due to the decarbonization and dehydroxylation of the clays, which released entrapped rare earth elements within the dominant minerals and converted them into an acid-soluble form. Nonetheless, the improvement in recovery was limited to the light REEs (LREEs) with an insignificant increase in the heavy REEs (HREE). It was demonstrated that the light and heavy REEs in the material were associated with difficult-to-leach minerals such as monazite, xenotime, and zircon, which were not decomposed by simple calcination due to their high thermal stability. Hence, roasting in the presence of chemicals was necessary to ensure the decomposition of those REE-containing minerals. As such, this study was focused on the acid baking treatment of bituminous coals with an aim to enhance REE recovery, especially HREEs. Based on the presence of REE minerals like monazite and xenotime, three pre-leach treatment methods, i.e., 1) roasting, 2) direct acid baking, and 3) acid baking after roasting were investigated. Roasting tests at 600 ⁰C revealed that the recovery of light REEs (LREEs) was enhanced while the recovery of HREEs remained relatively unaffected. LREE and HREE recovery values of 38.3% and 21.3%, respectively, were achieved using a 50 g/L sulfuric acid solution at 5% solid concentration and a solution temperature of 75 ⁰C for 2 hours. Comparatively, direct acid baking at 250 °C provided substantial increased LREE and HREE recovery values to approximately 49.4% and 53.0%, respectively, using an equivalent acid dosage. Recoveries were maximized to 77.0% and 79.6% for LREE and HREE, respectively, by roasting followed by acid baking. Similar results were obtained from the treatment of a second bituminous coal source. Due to strong correlations between REE and Al recovery values, tests were performed on kaolinite and illite, which were prominent clay minerals within the source coals. These experiments revealed that the REE recovery improvements were likely a result of dehydroxylation of clays and subsequent release and decomposition of REE-bearing minerals such as monazite, xenotime and zircon. Subsequently, a parametric study was conducted to identify the impact of acid baking parameters on rare earth element recovery. The factors investigated using a three-level statistical experimental program were acid baking time, acid solution concentration, baking temperature, and acid solution-to-solids ratio which were found to significantly impact REE and contaminant element (Al, Fe, and Ca) recovery. An increase in baking temperature up to around 250 ⁰C improved the light and heavy REE recovery values by more than 50 absolute percentage points relative to performances achieved when direct leaching. As aforementioned, acid baking was needed to both decompose the clay minerals and liberate the REE minerals, which allowed access for the acid to solubilize the REEs. Acid concentration of the solution used for acid baking was studied as a means of minimizing the amount of acid needed to achieve a target REE recovery. However, thermo-gravimetric and differential scanning calorimetry analysis (TGA-DSC) of sulfuric acid under oxidizing atmosphere revealed that the addition of water decreased the evaporation temperature, which explains the lower REE recovery values obtained when using acid concentrations less than 100%. Using pure sulfuric acid at an acid-to-solid ratio of 0.8:1 resulted in recovery values of around 70% for both LREEs and HREEs. The decomposition reaction time was relatively quick with 65% of the TREEs recovered within the first 10 minutes. Following the identification of optimum operating conditions through the parametric study, a systematic leaching study was carried out to examine the impact of leaching parameters, such as solid-to-liquid (S/L) ratio, temperature, and time, on REE recovery using acid baking conditions of 1:1 acid to solids ratio at 250 °C for 30 minutes. The solid-to-liquid ratio was varied from 1-20% by weight at 25 °C, 50 °C, and 75 °C solution temperatures. The results indicated that reaction time and solution temperature considerably impacted the recovery of heavy and light REEs. Interestingly, LREE recovery reduced from 68% at 5% S/L to 58% at 20%S/L, whereas HREE recovery of 78% remained unaffected. The decrease in the LREE recovery was determined to be due to La and Ce precipitation, likely through isomorphic substitution with calcium in gypsum. The kinetic data indicated that 67% LREE and 77% HREE recovery could be obtained within the first 15 minutes of the reaction, suggesting fast reaction kinetics. Furthermore, raising the solution temperature from 25 °C to 75 °C increased the LREE and HREE recovery from 60% and 32% to 67% and 79%, respectively. The kinetic modeling results demonstrated that the rate-limiting step in the LREE dissolution was diffusion and chemical reaction, whereas the HREE extraction was controlled by only chemical reaction. The leaching study concluded that using 20% S/L at 75 °C for 15 minutes maximized LREE and HREE recovery. The lab-scale precipitation study showed that Fe and Al in solution could be removed at pH 4.5 followed by REE precipitation at pH 6.0 using 6 mol/L NaOH. Finally, the bench-scale data was used to develop a process flowsheet for REE recovery from low-grade bituminous coal sources using acid baking. Finally, based on the proposed flowsheet, a concentrated RE-cake obtained through selective precipitation at pH 6.5 was re-leached using HCl at pH 1.5. The resultant leachate was used to identify the impact of various operating parameters on REE recovery and purity with an aim to maximize REE precipitation efficiency while minimizing the oxalic acid dosage. The operating parameters for this investigation were oxalic acid dosage, iron (III) contamination, solution pH and temperature. The resultant model suggested that oxalic acid dosage and reaction pH are the most significant factors for the REE precipitation efficiency, followed by the interaction of oxalic dosage and Fe concentration. Test results indicated that increasing the oxalic acid concentration from 0g/L to 80g/L improved the REE precipitation efficiency from approximately 4.2% to 95.0%. Furthermore, raising the solution pH from 0.5 to 2.5 considerably enhanced the precipitation efficiency from 0.0% to 98.9%. A solution temperature elevation decreased REE recovery, which indicated an exothermic reaction between REEs and oxalate anions. Finally, a high level of Fe contamination adversely impacted REE precipitation efficiency. The speciation analysis revealed that the dominant iron species in the solution system were Fe-(C₂O₄)₃³⁻, Fe-(C₂O₄)²⁻, and Fe-(C₂O₄)⁺, which consumed the majority of the oxalate anions

rare earth elements, acid baking, high-temperature↗

Phenotypically anchored transcriptomics across diverse agrichemicals reveals conserved pathways and unique gene expression signatures in zebrafish

Agrichemicals such as herbicides, fungicides, insecticides, and biocides are widely used in agriculture, yet some are associated with adverse effects in humans and the environment. While many of these chemicals have been extensively studied in vitro and are included in the EPA’s ToxCast program, comprehensive in vivo comparisons using RNA sequencing across structurally diverse agrichemicals, in a single screening platform, are lacking. In this study, we examined structurally diverse agrichemicals found in the U.S. Environmental Protection Agency’s (EPA) Toxcast Phase I and II library by statically exposing early life stage zebrafish at 6 h post fertilization (hpf) until 120 hpf at concentrations ranging from 0.25 to 100 µM. Morphological outcomes were assessed at 120 hpf across 10 endpoints, including yolk sac edema, craniofacial malformations, and axis abnormalities. Chemicals that produced robust concentration-response relationships were selected for transcriptomic profiling. For transcriptomic analysis, zebrafish were statically exposed to each chemical and sampled at 48 hpf, prior to the onset of morphological effects observed at 120 hpf. Differential expression analysis identified between 0 and 4,538 differentially expressed genes (DEGs) per chemical, with no clear correlation to morphological severity. Both DEG and co-expression network analyses revealed chemical-specific expression patterns that converged on shared biological pathways, including neurodevelopment and cytoskeletal organization. Key regulatory genes such as mylpfa and krt4 were identified within co-expression modules, suggesting their potential role in conserved toxicity mechanisms. Semantic similarity analysis of enriched gene ontology (GO) terms, when compared to existing datasets, highlighted gaps in the annotation of neurodevelopmental processes, indicating that some in vivo effects may not be fully captured by current curated resources. The results provide new insights into the modes of action of diverse agrichemicals and establish a framework for understanding how agrichemical structure relates to biological function in a vertebrate model.

agrichemical↗

Natural and Engineered Features Supporting Environmental Performance of Idaho National Laboratory's Remote-Handled Low-Level Waste Disposal Facility - 19436

The Idaho National Laboratory's (INL's) Remote-Handled Low-Level Waste (RH LLW) Disposal Facility has been designed and constructed to receive waste generated at the INL site in support of the US DOE Office of Nuclear Energy (NE) and Naval Reactors (NR) missions. The facility has been designed to receive legacy RH LLW currently in storage at INL facilities and new waste generated by nuclear research programs. RH LLW to be received at the facility is solid waste consisting of surface contaminated materials; ion exchange resins used to purify cooling and canal water; and activated metal components generated during reactor core change outs and Naval spent fuel management; and nuclear R and D. Disposal of INL's RH LLW poses unique disposal challenges due to high radiation levels (contact rates up to 600 Sv/hr [60,000 Rem/hr]) that require remote handling and adequate shielding for transportation and disposal. The facility has been designed to meet the requirements of DOE Order 435.1, Radioactive Waste Management, including the all-pathways dose, air pathway dose, inadvertent intruder doses, and radon emissions limits. This paper presents an overview of natural and engineered design features that support the safety case and enhance the environmental performance of the facility. These features include location selection, use of reinforced pre-cast concrete vaults; a robust hydraulic drainage system; use of steel waste canisters for waste disposal; the facility monitoring system; facility layout; and crediting waste forms. The facility performance assessment (PA) developed to meet the requirements of DOE Order 435.1, credits each of these protective features. Releases from waste forms disposed of at this facility are limited by corrosion of the steel waste canisters which will ultimately allow contact with infiltrating water followed by transport of radionuclides from the vault system. The steel waste canisters were credited for their ability to limit water contact with the waste based on the demonstrated ability of the reinforced pre-caste concrete vaults to provide long-term structural protection of the steel waste canisters. Verification of hydraulic and concrete performance was conducted through laboratory data, field infiltration tests supported by a robust monitoring system, and numerical modeling. Releases into the vault system from activated metals are based on site-specific corrosion data. Radionuclide transport from the vault system from all waste forms is based on site-specific climatological, hydrologic, and stratigraphic data. At INL, the groundwater pathway historically and appropriately has received the most attention. By crediting aspects of the natural system, the engineered vault system, and the waste, estimated radiologic doses via the groundwater pathway are minimal. Air pathway doses were reduced by differentially crediting corrosion of aluminum and stainless steel activated metals resulting in the pathway being screened out of the PA. Groundwater and inadvertent intruder doses were reduced through selective facility layout. This paper provides an overview of the system-wide protective features incorporated into the facility design along with an extensive reference list to supporting technical papers and analyses that provide data and interpretation. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Bit threads and holographic entanglement of purification

Generalizing the bit thread formalism, we use convex duality to derive dual flow programs to the bipartite and multipartite holographic entanglement of purification proposals and then prove several inequalities using these constructions. In the multipartite case we find the flows exhibit novel behavior which allows for a constrained flux on the boundary of the homology region. We show this flux can be made distinct from bi-partite terms and reflects the truly multipartite portion of the holographic entanglement of purification.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA Robotic Neurosurgery Testbed

The detection of tissue interface (e.g., normal tissue, cancer, tumor) has been limited clinically to tactile feedback, temperature monitoring, and the use of a miniature ultrasound probe for tissue differentiation during surgical operations. In neurosurgery, the needle used in the standard stereotactic CT (Computational Tomography) or MRI (Magnetic Resonance Imaging) guided brain biopsy provides no information about the tissue being sampled. The tissue sampled depends entirely upon the accuracy with which the localization provided by the preoperative CT or MRI scan is translated to the intracranial biopsy site. In addition, no information about the tissue being traversed by the needle (e.g., a blood vessel) is provided. Hemorrhage due to the biopsy needle tearing a blood vessel within the brain is the most devastating complication of stereotactic CT/MRI guided brain biopsy. A robotic neurosurgery testbed has been developed at NASA Ames Research Center as a spin-off of technologies from space, aeronautics and medical programs. The invention entitled 'Robotic Neurosurgery Leading to Multimodality Devices for Tissue Identification' is nearing a state ready for commercialization. The devices will: 1) improve diagnostic accuracy and precision of general surgery, with near term emphasis on stereotactic brain biopsy, 2) automate tissue identification, with near term emphasis on stereotactic brain biopsy, to permit remote control of the procedure, and 3) reduce morbidity for stereotactic brain biopsy. The commercial impact from this work is the potential development of a whole new generation of smart surgical tools to increase the safety, accuracy and efficiency of surgical procedures. Other potential markets include smart surgical tools for tumor ablation in neurosurgery, general exploratory surgery, prostate cancer surgery, and breast cancer surgery.

Mah, Robert↗

Attitude dynamics simulation subroutines for systems of hinge-connected rigid bodies with nonrigid appendages

Three computer subroutines designed to solve the vector-dyadic differential equations of rotational motion for systems that may be idealized as a collection of hinge-connected rigid bodies assembled in a tree topology, with an optional flexible appendage attached to each body are reported. Deformations of the appendages are mathematically represented by modal coordinates and are assumed small. Within these constraints, the subroutines provide equation solutions for (1) the most general case of unrestricted hinge rotations, with appendage base bodies nominally rotating at a constant speed, (2) the case of unrestricted hinge rotations between rigid bodies, with the restriction that those rigid bodies carrying appendages are nominally nonspinning, and (3) the case of small hinge rotations and nominally nonrotating appendages. Sample problems and their solutions are presented to illustrate the utility of the computer programs.

Fleischer, G. E.↗

Validation of Test Methods for Air Leak Rate Verification of Spaceflight Hardware

As deep space exploration continues to be the goal of NASAs human spaceflight program, verification of the performance of spaceflight hardware becomes increasingly critical. Suitable test methods for verifying the leak rate of sealing systems are identified in program qualification testing requirements. One acceptable method for verifying the air leak rate of gas pressure seals is the tracer gas leak detector method. In this method, a tracer gas (commonly helium) leaks past the test seal and is transported to the leak detector where the leak rate is quantified. To predict the air leak rate, a conversion factor of helium-to-air is applied depending on the magnitude of the helium flow rate. The conversion factor is based on either the molecular mass ratio or the ratio of the dynamic viscosities. The current work was aimed at validating this approach for permeation-level leak rates using a series of tests with a silicone elastomer O-ring. An established pressure decay method with constant differential pressure was used to evaluate both the air and helium leak rates of the O-ring under similar temperature and pressure conditions. The results from the pressure decay tests showed, for the elastomer O-ring, that neither the molecular flow nor the viscous flow helium-to-air conversion factors were applicable. Leak rate tests were also performed using nitrogen and argon as the test gas. Molecular mass and viscosity based helium-to-test gas conversion factors were applied, but did not correctly predict the measured leak rates of either gas. To further this study, the effect of pressure boundary conditions was investigated. Often, pressure decay leak rate tests are performed at a differential pressure of 101.3 kPa with atmospheric pressure on the downstream side of the test seal. In space applications, the differential pressure is similar, but with vacuum as the downstream pressure. The same O-ring was tested at four unique differential pressures ranging from 34.5 to 137.9 kPa. Up to six combinations of upstream and downstream pressures for each differential pressure were compared. For a given differential pressure, the various combinations of upstream and downstream dry air pressures did not significantly affect the leak rate. As expected, the leak rate of the O-ring increased with increasing differential pressure. The results suggested that the current leak test pressure conditions, used to verify spacecraft sealing systems with elastomer seals, produce accurate values even though the boundary conditions do not model the space application.

helium leak detector↗

Nuclear Data: From Measurement, Evaluation, Processing to Practical Applications [Slides]

Any analysis, design of a nuclear system will be as good as the knowledge of the nuclear data used whether differential and integral. ORNL is engaged in generating high level nuclear data for reactor, shielding, and safety applications. Nuclear data remains a major source of data uncertainties. Uncertainties on nuclear data rely upon measurements and evaluations. Tight connections with experimental facilities around the world have helped ORNL to address issues on ND and uncertainty. It is extremely important for ORNL to retain expertise on nuclear data in support of the research and expertise activities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Intracellular stress tomography reveals stress focusing and structural anisotropy in cytoskeleton of living cells

We describe a novel synchronous detection approach to map the transmission of mechanical stresses within the cytoplasm of an adherent cell. Using fluorescent protein-labeled mitochondria or cytoskeletal components as fiducial markers, we measured displacements and computed stresses in the cytoskeleton of a living cell plated on extracellular matrix molecules that arise in response to a small, external localized oscillatory load applied to transmembrane receptors on the apical cell surface. Induced synchronous displacements, stresses, and phase lags were found to be concentrated at sites quite remote from the localized load and were modulated by the preexisting tensile stress (prestress) in the cytoskeleton. Stresses applied at the apical surface also resulted in displacements of focal adhesion sites at the cell base. Cytoskeletal anisotropy was revealed by differential phase lags in X vs. Y directions. Displacements and stresses in the cytoskeleton of a cell plated on poly-L-lysine decayed quickly and were not concentrated at remote sites. These data indicate that mechanical forces are transferred across discrete cytoskeletal elements over long distances through the cytoplasm in the living adherent cell.

NASA Discipline Cell Biology↗

Airborne gravity and other geophysical techniques for understanding the lithosphere beneath the West Antarctic Ice Sheet

As part of a program entitled Corridor Aerogeophysics of the Southeastern Ross Transect Zone (CASERTZ), an aerogeophysical platform was developed to study the interaction of geological and glaciological processes in West Antarctica. A de Havilland Twin Otter was equipped with an ice-penetrating radar, a proton precession magnetometer, an airborne gravity system, and a laser altimeter. The 60-MHz ice-penetrating radar can recover sub-ice topography with an accuracy of about 10 m through 3 km of comparatively warm West Antarctic ice, while the laser altimeter profiling of the ice surface is accurate to approximately 1 m. The magnetic field observations are accurate to several nT, and the gravity measurements are accurate to better than 3 mGal. The aircraft is navigated by a local radio transponder network, while differential positioning techniques based on the Global Positioning System (GPS) satellites are used for recovering high-resolution horizontal and vertical positions. Attitude information from an inertial navigation system is used to correct the laser altimetry and a digital pressure transducer is used to recover vertical positions and accelerations in the absence of satellite positioning. Continuous base-station observations are made for the differential GPS positioning and the removal of ionospheric noise from the airborne magnetometer measurements.

Bell, Robin E.↗

Development of Multiobjective Optimization Techniques for Sonic Boom Minimization

A discrete, semi-analytical sensitivity analysis procedure has been developed for calculating aerodynamic design sensitivities. The sensitivities of the flow variables and the grid coordinates are numerically calculated using direct differentiation of the respective discretized governing equations. The sensitivity analysis techniques are adapted within a parabolized Navier Stokes equations solver. Aerodynamic design sensitivities for high speed wing-body configurations are calculated using the semi-analytical sensitivity analysis procedures. Representative results obtained compare well with those obtained using the finite difference approach and establish the computational efficiency and accuracy of the semi-analytical procedures. Multidisciplinary design optimization procedures have been developed for aerospace applications namely, gas turbine blades and high speed wing-body configurations. In complex applications, the coupled optimization problems are decomposed into sublevels using multilevel decomposition techniques. In cases with multiple objective functions, formal multiobjective formulation such as the Kreisselmeier-Steinhauser function approach and the modified global criteria approach have been used. Nonlinear programming techniques for continuous design variables and a hybrid optimization technique, based on a simulated annealing algorithm, for discrete design variables have been used for solving the optimization problems. The optimization procedure for gas turbine blades improves the aerodynamic and heat transfer characteristics of the blades. The two-dimensional, blade-to-blade aerodynamic analysis is performed using a panel code. The blade heat transfer analysis is performed using an in-house developed finite element procedure. The optimization procedure yields blade shapes with significantly improved velocity and temperature distributions. The multidisciplinary design optimization procedures for high speed wing-body configurations simultaneously improve the aerodynamic, the sonic boom and the structural characteristics of the aircraft. The flow solution is obtained using a comprehensive parabolized Navier Stokes solver. Sonic boom analysis is performed using an extrapolation procedure. The aircraft wing load carrying member is modeled as either an isotropic or a composite box beam. The isotropic box beam is analyzed using thin wall theory. The composite box beam is analyzed using a finite element procedure. The developed optimization procedures yield significant improvements in all the performance criteria and provide interesting design trade-offs. The semi-analytical sensitivity analysis techniques offer significant computational savings and allow the use of comprehensive analysis procedures within design optimization studies.

Chattopadhyay, Aditi↗

PETSc/TAO Users Manual (Rev. 3.19)

This manual describes the use of the Portable, Extensible Toolkit for Scientific Computation (PETSc) and the Toolkit for Advanced Optimization (TAO) for the numerical solution of partial differential equations and related problems on high-performance computers. PETSc/TAO is a suite of data structures and routines that provide the building blocks for the implementation of large-scale application codes on parallel (and serial) computers. PETSc uses the MPI standard for all distributed memory communication. PETSc/TAO includes a large suite of parallel linear solvers, nonlinear solvers, time integrators, and opti mization that may be used in application codes written in Fortran, C, C++, and Python (via petsc4py; see Getting Started). PETSc provides many of the mechanisms needed within parallel application codes, such as parallel matrix and vector assembly routines. The library is organized hierarchically, enabling users to employ the level of abstraction that is most appropriate for a particular problem. By using techniques of object-oriented programming, PETSc provides enormous flexibility for users. PETSc is a sophisticated set of software tools; as such, for some users it initially has a much steeper learning curve than packages such as MATLAB or a simple subroutine library. In particular, for individuals without some computer science background, experience programming in C, C++, python, or Fortran and experience using a debugger such as gdb or lldb, it may require a significant amount of time to take full advantage of the features that enable efficient software use. However, the power of the PETSc design and the algorithms it incorporates may make the efficient implementation of many application codes simpler than “rolling them” yourself. For many tasks a package such as MATLAB is often the best tool; PETSc is not intended for the classes of problems for which effective MATLAB code can be written. There are several packages, built on PETSc, that may satisfy your needs without requiring directly using PETSc. We recommend reviewing these packages functionality before starting to code directly with PETSc. PETSc can be used to provide a “MPI parallel linear solver” in an otherwise sequential, or OpenMP parallel code. This approach cannot provide extremely large improvements in the application time by utilizing large numbers of MPI processes but can still improve the performance. Certainly all parts of a previously sequential code need not be parallelized but the matrix generation portion must be parallelized to expect true scalability to large numbers of MPI processes. See PCMPI for details on how to utilize the PETSc MPI linear solver server. Since PETSc is under continued development, small changes in usage and calling sequences of routines will occur. PETSc has been supported for twenty-five years; see mailing list information on our website for information on contacting support.

97 MATHEMATICS AND COMPUTING↗

PETSc Users Manual (Rev. 3.13)

This manual describes the use of PETSc for the numerical solution of partial differential equations and related problems on high-performance computers. The Portable, Extensible Toolkit for Scientific Computation (PETSc) is a suite of data structures and routines that provide the building blocks for the implementation of large-scale application codes on parallel (and serial) computers. PETSc uses the MPI standard for all message-passing communication. PETSc includes an expanding suite of parallel linear solvers, nonlinear solvers, and time integrators that may be used in application codes written in Fortran, C, C++, and Python. PETSc provides many of the mechanisms needed within parallel application codes, such as parallel matrix and vector assembly routines. The library is organized hierarchically, enabling users to employ the level of abstraction that is most appropriate for a particular problem. By using techniques of object-oriented programming, PETSc provides enormous flexibility for users. PETSc is a sophisticated set of software tools; as such, for some users it initially has a much steeper learning curve than a simple subroutine library. In particular, for individuals without some computer science background, experience programming in C, C++, python, or Fortran and experience using a debugger such as gdb or dbx, it may require a significant amount of time to take full advantage of the features that enable efficient software use. However, the power of the PETSc design and the algorithms it incorporates may make the efficient implementation of many application codes simpler than \rolling them" yourself.

97 MATHEMATICS AND COMPUTING↗

Development of the Functional Flow Block Diagram for the J-2X Rocket Engine System

The J-2X program calls for the upgrade of the Apollo-era Rocketdyne J-2 engine to higher power levels, using new materials and manufacturing techniques, and with more restrictive safety and reliability requirements than prior human-rated engines in NASA history. Such requirements demand a comprehensive systems engineering effort to ensure success. Pratt & Whitney Rocketdyne system engineers performed a functional analysis of the engine to establish the functional architecture. J-2X functions were captured in six major operational blocks. Each block was divided into sub-blocks or states. In each sub-block, functions necessary to perform each state were determined. A functional engine schematic consistent with the fidelity of the system model was defined for this analysis. The blocks, sub-blocks, and functions were sequentially numbered to differentiate the states in which the function were performed and to indicate the sequence of events. The Engine System was functionally partitioned, to provide separate and unique functional operators. Establishing unique functional operators as work output of the System Architecture process is novel in Liquid Propulsion Engine design. Each functional operator was described such that its unique functionality was identified. The decomposed functions were then allocated to the functional operators both of which were the inputs to the subsystem or component performance specifications. PWR also used a novel approach to identify and map the engine functional requirements to customer-specified functions. The final result was a comprehensive Functional Flow Block Diagram (FFBD) for the J-2X Engine System, decomposed to the component level and mapped to all functional requirements. This FFBD greatly facilitates component specification development, providing a well-defined trade space for functional trades at the subsystem and component level. It also provides a framework for function-based failure modes and effects analysis (FMEA), and a rigorous baseline for the functional architecture.

White, Thomas↗