SEARCH · Search NASA

Results for “Numerical performance”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 109 records · Page 6

All sources NASA NTRS DOE OSTI

NASA NTRS2004

Microscopic and Metallurgical Aspects of the Space Shuttle Columbia Accident Investigation and Reconstruction

The Space Shuttle Columbia was descending for a landing at the Kennedy Space Center (KSC) on February 1, 2003. Approximately 20 minutes prior to touchdown, the Columbia began disintegrating over the western United States; the majority of debris eventually impacted in eastern Texas and western Louisiana. A monumental effort eventually recovered approximately 84,000 pieces of debris, approximately 38% of the Orbiter's original dry weight. The debris was transported to KSC, where the items were catalogued and evaluated. Critical areas of interest, such as the left and right leading edge surfaces and the underside of the ship, were placed upon a grid to aid in the reconstruction. Items of interest included metallic structures, reinforced carbon-carbon composites, and ceramic heat insulation tiles. Many of the leading edge elements had re-solidified metallic deposits spattered on them. These deposits became known as slag and were one of the main focuses of the investigation. In order to help determine the sequence of events inside the left wing during the accident, the slag's composition, layering order, and directionality of deposition were studied. A myriad of analytical tests were performed in an attempt to ascertain the compositional and depositional characteristics of selected slag deposits, including the ordering of deposited layers within each individual slag deposit harvested. Initially, Scanning Electron Microscopy and Energy Dispersive X-Ray Spectroscopy (SEM/EDX) were performed to quickly characterize the overall composition of individual slag deposits: SEM utilizes a narrowlyfocused high-energy electron beam impinging upon a specimen. The incident beam excites and liberates lower energy secondary electrons, which are detected and analyzed, providing a visual representation of the sample's surface topography. EDX also relies on an incident electron beam, except an EDX unit measures X-ray energies generated by the impinging beam. Each element generates a unique X-ray signature; the EDX detector measures these discreet energies. EDX actually penetrates approximately 2 microns into the bulk of the sample. However, random examination of various portions of slag, coupled with the semiquantitative nature of the SEM/EDX analysis, did not yield convincingly pertinent data. Therefore, X-ray dot mapping was conducted, which provided more understandable data, both in terms of slag layering and composition. An X-ray dot map is generated by performing numerous EDX scans for individual elements, then compiling the scans in a visual representation. Eventually, specimens consisting of not only the slag, but of the adjacent RCC substrate as well were cross-sectioned. X-Ray dot mapping of the materialographicallymounted and -polished cross- sections provided a visual representation of both the layering sequence and compositional characteristics of the slag. Contemporaneously, Electron Spectroscopy for Chemical Analysis/X-Ray Photoelectron Spectroscopy (ESCA/XPS) and powdered X-Ray Diffraction (XRD) were performed to further characterize the deposits and to attempt to identify what, if any, compounds were present. The ESCA/XPS analysis allowed the analyst to "sputter" into the sample with an electron gun, aiding in the identification of the layering sequence. XPS uses photons, rather than electrons, which impinge upon the surface of the sample. XPS measures the electrons emitted from within the first 5 nm of the sample's surface. The XRD measures the scatter angles of incident X-rays; the angle and intensity of scatter depend upon the crystalline structure of the pulverized sample. XRD is considered a qualitative rather than quantitative technique. ESCA/XPS revealed that the final layer to deposit was predominantly carbonaceous. XRD was successful in identifying specific compounds, such as Al 2O3, Al and/or Al3 21SiO47, mullite (3(Al2)O3 -SiO2), and nickel-aluminides. Eventually, Electron MicroProbe Analysis (EMPA) was conducted on the marialographically-prepared cross- sections of selected slag deposits. Microprobe combines SEM and Wavelength Dispersive X-Ray Spectroscopy (WDS), and, like EDX, uses a narrowly-focused high-energy electron beam impinging upon a specimen to elicit, in the case of EPMA, characteristic X-rays with specific wavelengths. This quantitative, analytical tool proved the most useful in determining depositional layering and composition of the slag deposits. This information was utilized in verifying the location of the breach in the left leading edge of the wing of the Columbia.

McDaniels, Steven J.↗

NASA NTRS1995

Chemical and Thermal Analysis

During the past six months we have conducted significant research in several domains in order to clarify and understanding the aging and chemical failure mechanism of thermoplastics (PVDF or Tefzel) for pipes. We organized numerous analytical studies with methods including Fourier Transform Infrared Spectroscopy, Dynamic Mechanical Analysis, Differential Scanning Calorimetry, and Stress Relaxation experiments. In addition we have reanalyzed previous thermogravimetric data concerning the rate of deplasticization of Coflon pipe. We investigated a number of aged samples of both Tefzel and Coflon that were forwarded from MERL. We conducted stress relaxation experiments of Coflon pipe at several temperatures and determined an activation energy. We also examined the dynamic mechanical response PVDF during deplasticization and during methanol plasticization. We performed numerous DSC analyses to research the changing crystalline morphology. We have noted significant changes in crystallinity upon aging for both PVDF and Tefzel. Little variation in elemental composition was noted for many of the aged Coflon and Tefzel samples tested.

Bulluck, J. W.↗

NASA NTRS2006

NASA Langley Research Center Systems Analysis & Concepts Directorate Participation in the Exploration Systems Architecture Study

The NASA Langley Research Center (LaRC) Systems Analysis & Concepts Directorate (SACD) began studying human exploration missions beyond low Earth orbit (LEO) in the year 1999. This included participation in NASA s Decadal Planning Team (DPT), the NASA Exploration Team (NExT), Space Architect studies and Revolutionary Aerospace Systems Concepts (RASC) architecture studies that were used in formulating the new Vision for Space Exploration. In May of 2005, NASA initiated the Exploration Systems Architecture Study (ESAS). The primary outputs of the ESAS activity were concepts and functional requirements for the Crewed Exploration Vehicle (CEV), its supporting launch vehicle infrastructure and identification of supporting technology requirements and investments. An exploration systems analysis capability has evolved to support these functions in the past and continues to evolve to support anticipated future needs. SACD had significant roles in supporting the ESAS study team. SACD personnel performed the liaison function between the ESAS team and the Shuttle/Station Configuration Options Team (S/SCOT), an agency-wide team charged with using the Space Shuttle to complete the International Space Station (ISS) by the end of Fiscal Year (FY) 2010. The most significant of the identified issues involved the ability of the Space Shuttle system to achieve the desired number of flights in the proposed time frame. SACD with support from the Kennedy Space Center performed analysis showing that, without significant investments in improving the shuttle processing flow, that there was almost no possibility of completing the 28-flight sequence by the end of 2010. SACD performed numerous Lunar Surface Access Module (LSAM) trades to define top level element requirements and establish architecture propellant needs. Configuration trades were conducted to determine the impact of varying degrees of segmentation of the living capabilities of the combined descent stage, ascent stage, and other elements. The technology assessment process was developed and implemented by SACD as the ESAS architecture was refined. SACD implemented a rigorous and objective process which included (a) establishing architectural functional needs, (b) collection, synthesis and mapping of technology data, and (c) performing an objective decision analysis resulting in technology development investment recommendations. The investment recommendation provided budget, schedule, and center/program allocations to develop required technologies for the exploration architecture, as well as the identification of other investment opportunities to maximize performance and flexibility while minimizing cost and risk. A summary of the trades performed and methods utilized by SACD for the Exploration Systems Mission Directorate (ESAS) activity is presented along with how SACD is currently supporting the implementation of the Vision for Space Exploration.

Keyes, Jennifer↗

DOE OSTI2021

Evaluation of Energy Storage Potential of Unconventional Shale Reservoirs Using Numerical Simulation of Cyclic Gas Injection

Compressed air energy storage (CAES) stores energy as compressed air in underground formations, typically salt dome caverns. When electricity demand grows, the compressed air is released through a turbine to produce electricity. CAES in the US is limited to one plant built in 1991, due in part to the inherent risk and uncertainty of developing subsurface storage reservoirs. As an alternative to CAES, we propose using some of the hundreds of thousands of hydraulically fractured horizontal wells to store energy as compressed natural gas in unconventional shale reservoirs. To store energy, produced or “sales” natural gas is injected back into the formation using excess electricity and is later produced through an expander to generate electricity. To evaluate this concept, we performed numerical simulations of cyclic natural gas injection into unconventional shale reservoirs using cmg-gem commercial reservoir modeling software. We tested short-term (diurnal) and long-term (seasonal) energy storage potential by modeling well injection and production gas flowrates as a function of bottom-hole pressure. First, we developed a conceptual model of a single fracture stage in an unconventional shale reservoir to characterize reservoir behavior during cyclic injection and production. Next, we modeled cyclic injection in the Marcellus shale gas play using published data. Results indicate that Marcellus unconventional shale reservoirs could support both short- and long-term energy storage at capacities of 100–1000 kWe per well. The results indicate that energy storage in unconventional shale gas wells may be feasible and warrants further investigation.

25 ENERGY STORAGE↗

DOE OSTI2021

Preconditioners for multiphase poromechanics with strong capillarity

This paper aims to enhance the performance of Newton–Krylov solvers for coupled poromechanical problems with two-phase flow. In particular, we investigate the impact of capillary pressure on preconditioning strategies. Capillarity complicates the coupling between the solid deformation and fluid pressure degrees of freedom, as well as increases the nonlinearity of the system. Depending on the capillary pressure relation used in the constitutive formulation, the flow equations may exhibit a spectrum of advection-dominated to diffusion-dominated behavior. We propose preconditioning approaches that account for this behavior and lead to robust numerical performance within a broad range of regimes.

42 ENGINEERING↗

DOE OSTI2020

Non-local model for surface tension in fluid-fluid simulations

We propose a non-local model for surface tension obtained in the form of an integral of a molecular-force-like function with support 3.5e added to the Navier-Stokes momentum conservation equation. We demonstrate analytically and numerically that with the non-local model interfaces with a radius of curvature larger than the support length behave macroscopically and microscopically, otherwise. For static droplets, the pressure difference Pe,in - Pe,out satisfies the Young-Laplace law for droplet radius greater than 3.5e and other- wise deviates from the Young-Laplace law. The latter indicates that the surface tension in the proposed model decreases with decreasing radius of curvature, which agrees with molecular dynamics and experimental studies of nanodroplets. Using the non-local model we perform numerical simulations of droplets under dynamic conditions, including a rising droplet, a droplet in shear flow, and two colliding droplets in shear flow, and compare results with a standard Navier-Stokes model subject to the Young-Laplace boundary condition at the fluid-fluid interface implemented via the Conservative Level Set (CLS) method. We find good agreement with existing numerical methods and analytical results for a rising macroscopic droplet and a droplet in a shear flow. For colliding droplets in shear flow, the non- local model converges (with respect to the grid size) to the correct behavior, including sliding, coalescence, and merging and breaking of two droplets depending on the capillary number. In contrast, we find that the results of the CLS model are highly grid-size dependent.

Howard, Amanda A.↗

DOE OSTI2025

Enabling Low-Temperature (LTP) Ignition Technologies for Multi-Mode Engines through the Development of a Validated High-Fidelity LTP Model for Predicative Simulations Tools

The goal of multi-mode engine architectures is to extend current lean-burn dilution limits with renewable fuels, which requires spark plugs to deposit high energies (hundreds of mJ) in order to initiate ignition and complete combustion. At elevated energy deposition rates, spark plugs experience increased electrode erosion and thermal losses, which ultimately shortens the spark-plug lifetime and lowers ignition efficiency. As such, in order to safeguard the efficiency gains of multi-mode concepts, new and improved ignition technologies are required. Recently, non-equilibrium low-temperature plasmas (LTP) have been shown to promote energy-efficient ignition via quenching and transport of electronically excited atoms and molecules, selective radical production and fast heating of hydrocarbon/air mixtures [1-2]. Thus, LTP is seen as a technology that can potentially improve the energy extraction efficiency of fuels, while enabling kinetically controlled combustion modes towards fuel leaner conditions to realize current DOE VTO goals of improving the sustainability of future mobility [3]. Although many previous studies have demonstrated the efficacy of plasma-assisted ignition to enhance combustion, the detailed enhancement mechanisms remain largely unknown, especially for oxygenated fuels and at elevated pressures that are most relevant to practical engine conditions. These barriers hinder the development of accurate and comprehensive numerical models that seek to describe LTP-based ignition in existing engine design software tools and methods. Current state-of-the-art simulation capabilities for LTP ignition systems are in need of improvements since they deliver qualitative results only due to important limitations of existing approaches. Firstly, validated kinetic models with elementary steps for plasma discharges in oxygenated fuel/air mixtures of relevance to the transportation sector are required. Such kinetic models do not exist at present and will be developed and validated within this project. Secondly, plasma discharges and reactive mixture ignition are multi-scale, unsteady processes requiring high-performance numerical methods and software that execute efficiently on DOE supercomputers. Such software does not exist at present and will be developed and applied to practical LTP ignition scenarios as part of this project. Thirdly, experimental databases that are tailored to serve as benchmark in support of the development of predictive computational models of LTP ignition do not exist and will be part of this project.

33 ADVANCED PROPULSION SYSTEMS↗

NASA NTRS1998

Incompressible Turbulent Wing-Body Junction Flow

The overall objective of this study is to contribute to the optimized design of fan bypass systems in advanced turbofan engines. Increasing the engine bypass ratios have provided a major boost in engine performance improvement over the last fifty years. An engine with high bypass ratio (11-16:1) such as the Advanced Ducted Propulsion (ADP) is being developed and is expected to provide an additional 25% improvement in overall efficiency over the early turbofans. Such significant improvements in overall efficiency would reduce the cost per seat mile, which is a major government and Industry challenge for the 21th century. The research is part of the Advanced Subsonic Technology (AST) program that involves a NASA, U.S. Industry and FAA partnership with the goal of a safe and highly productive global air transportation system. The immediate objective of the study is to perform numerical simulation of duct-strut interactions to elucidate the loss mechanisms associated with this configuration that is typical of advanced turbofan engines such as ADP. However, at present experimental data for a duct-strut configuration are not available. Thus, as a first step a wing-body junction flow would be studied and is the specific objective of the present study. At the outset it is to be recognized that while duct-strut interaction flow is similar to that of wing-body junction flows, there are some differences owing to the presence of a wall at both ends of the strut. Likewise, some differences are due to the sheared inflow (as opposed to a uniform inflow) velocity profile. It is however expected that some features of a wing-body junction flow would persist. Next, some of the salient aspects of the complex flow near a wing-body junction, as revealed by various studies reported in the literature will be reviewed. One of the principle characteristics of the juncture flow, is the presence of the mean flow components in a plane perpendicular to the direction of the oncoming free-stream flow. The lateral curvature of the wing/strat causes the oncoming turbulent layer to skew about am axis (x-axis) parallel to the plane (xz-plane) of the mean shear. This is the principle mechanism for the generation of secondary flow. Such skew-induced secondary flows are slow to be attenuated by Reynolds stresses. Additional contribution to the generation of secondary flow comes from anisotropies in Reynolds stresses. Upstream of the strut, the mean-vorticity is directed span wise (along the y-direction). The presence of secondary flow in the vicinity of the strut causes the vorticity to stretch around the obstacle in a horse-shoe shape, with each leg having a vorticity of the opposite sense. The blockage effect of the strut imposes a severe adverse pressure gradient on the oncoming turbulent shear layer, causing boundary layer separation ahead of the leading edge, resulting in a vortex that rolls up and flows downstream into the juncture region. The separation vortices trailing in the wake of the wing can alter the lift or drag characteristics of the surfaces downstream of the wing-body juncture. Likewise, on submarines, the wake flow behind the appendage can degrade the performance of the propeller located downstream. The complex nature of this flow is caused by the presence of all six components of Reynolds stresses. Devenport and Simpson report that in the vicinity of the horse-shoe vortex there is intense recirculation with turbulent stresses being much larger than those normally observed in turbulent flows. These features contribute to making this flow a challenge to predict numerically. Some of the past studies provide useful insights into this flow that would guide our numerical efforts. In measurements reported by Shabaka and Bradshaw, the eddy viscosity tensor is seen to be non-isotropic and has negative components in certain regions. In an effort to evaluate the closure assumptions of various turbulence models, Devenport and Simpson used their own extensive measurements in juncture flows around the nose of a wing-body junction. Measured values of mean-velocity and/or turbulence kinetic energy was used to predict the magnitude of the shear stress vector. Algebraic stress models performed the best followed by Cebeci-Smith eddy viscosity model. The flow is reported to be dominated by a pressure field produced by the wing and the velocity field generated by the horseshoe vortex that is wrapped around the junction between the wing and wall. Kubendran et al. conclude from an experimental study that the shape of leading edge of the wing as characterized by its slenderness ratio is a major factor in determining the flow fields in the juncture region. The more thinner the leading edge of the juncture, the weaker the horseshoe vortex is. Also, with a slender leading edge, the secondary flow in the juncture would be due mainly to the cross-stream gradients of Reynolds stresses rather than due to a lateral skewing of the shear layer.

Krishnamurthy, R.↗

DOE OSTI2024

Mortar-based Cohesive Zone Model with Application to TRISO Particle Debonding

The U.S. Department of Energy’s Nuclear Energy Advanced Modeling and Simulation (NEAMS) program aims to develop predictive capabilities by applying computational methods to the analysis and design of advanced reactor and fuel cycle systems. This program has been providing engineering-scale support for the development of BISON, a high-fidelity and high-resolution fuel performance tool. This report outlines work towards the development of an integrated, mortar-based cohesive zone framework with mechanical contact. We model debonding of tri-structural isotropic (TRISO) layers using a mortar approach, extending the framework used for mechanical contact, using weighted quantities. A bilinear mixed mode traction approach, already used in Jiang et al. (2021), is applied analogously herein to model debonding between TRISO particle layers. The examples include the debonding of the inner pyrolytic carbon (IPyC) from the silicon carbide (SiC) layers, as well as partial and full debonding of the buffer from the IPyC layers. Both 2D and 3D examples are utilized to demonstrate the superior numerical performance of the new mortar approach, which allows for modeling arbitrarily meshed interfaces (i.e., non-matching discretizations).

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2021

Notes on the Variety of Multipactor Patterns

The experimental results on multipactor often are not in a good agreement with theoretical predictions and numerical simulations. The experimental results can noticeably vary from test to test even in the same geometry of electrodes, depending on the material they made of and condition of their surfaces, i.e., secondary emission properties. On the other hand, an actual secondary emission yield (SEY) of electrode material is never known a priory in the experiments. In practice, the SEY data for given material used for multipactor evaluation are obtained in specialized measurements. A difference between actual and tabulated SEY of material may result in a disagreement between pre-dictions and measurements. In this work the impact of the basic SEY characteristic variations on the multipactor dynamics in a variety of coaxial and rectangular waveguides is studied. The study was performed numerically with the use of CST Particle Studio.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRS2011

Venting of a Water/Inhibited Propylene Glycol Mixture in a Vacuum Environment-Characterization and Representative Test Results

A planned use of the Orion space vehicle involves its residence at the International Space Station for six months at a time. One concept of operations involves temporarily venting portions of the idle Orion active thermal control system (ATCS) during the docked phase, preventing freezing. The venting would have to be reasonably complete with few, if any, completely filled pockets of frozen liquid. Even if pockets of frozen liquid did not damage the hardware during the freezing process, they could prevent the system from filling completely prior to its reactivation. The venting of single component systems in a space environment has been performed numerous times and is well understood. Local nucleation occurs at warm, relatively massive parts of the system, which creates vapor and forces the bulk liquid out of the system. The remnants of the liquid will freeze, then evaporate over time through local heating. Because the Orion ATCS working fluid is a 50/50 mixture of water and inhibited propylene glycol, its boiling behavior was expected to differ from that of a pure fluid. It was thought that the relatively high vapor pressure water might evaporate preferentially, leaving behind a mixture enriched with the low vapor pressure propylene glycol, which would be vaporization ]resistant. Owing to this concern, a test was developed to compare the evaporation behavior of pure water, a 50/50 mixture of water and inhibited propylene glycol, and inhibited propylene glycol. The test was performed using room temperature fluids in an instrumented thin walled stainless steel vertical tube. The 1 in x 0.035 in wall tube was instrumented with surface thermocouples and encased in closed cell polyurethane foam. Reticulated polyurethane foam was placed inside the tube to reduce the convection currents. A vacuum system connected to the top of the tube set the pressure boundary condition. Tests were run for the three fluids at back pressures ranging from 1 to 18 torr. During each test, the mass of the test article was measured as it changed over time, as was its temperature and backpressure. The tests were successful. Somewhat surprisingly, the results showed that the evaporation behavior of the three fluids had more similarities than differences. The 50/50 mixture evaporated similarly to the pure water - albeit at a slower rate. The test results indicate that our extensive space - based experience with venting of single component fluids can be applied to the problem of Orion ATCS venting as long as the appropriate puts, takes, and caveats are applied.

Ungar, Eugene K.↗

DOE OSTI2022

Towards sharp error analysis of extended Lagrangian molecular dynamics

The extended Lagrangian molecular dynamics (XLMD) method provides a useful framework for reducing the computational cost of a class of molecular dynamics simulations with constrained latent variables. The XLMD method relaxes the constraints by introducing a fictitious mass ε for the latent variables and solving a set of singularly perturbed ordinary differential equations. While favorable numerical performance of XLMD has been demonstrated in several different contexts in the past decade, mathematical analysis of the method remains scarce. Here, we propose the first error analysis of the XLMD method in the context of a classical polarizable force field model. While the dynamics with respect to the atomic degrees of freedom are general and nonlinear, the key mathematical simplification of the polarizable force field model is that the constraints on the latent variables are given by a linear system of equations. We prove that when the initial value of the latent variables is compatible in a sense that we define, XLMD converges as the fictitious mass ε is made small with $\mathscr{O}$(ε) error for the atomic degrees of freedom and with $\mathscr{O}$($\sqrt{ε}$) error for the latent variables, when the dimension of the latent variable d' is 1. Furthermore, when the initial value of the latent variables is improved to be optimally compatible in a certain sense, we prove that the convergence rate can be improved to $\mathscr{O}$(ε) for the latent variables as well. Numerical results verify that both estimates are sharp not only for d'=1, but also for arbitrary d'. In the setting of general d', we do obtain convergence, but with the non-sharp rate of $\mathscr{O}$($\sqrt{ε}$) for both the atomic and latent variables.

74 ATOMIC AND MOLECULAR PHYSICS↗

NASA NTRS1993

Stanford Aerospace Research Laboratory research overview

Over the last ten years, the Stanford Aerospace Robotics Laboratory (ARL) has developed a hardware facility in which a number of space robotics issues have been, and continue to be, addressed. This paper reviews two of the current ARL research areas: navigation and control of free flying space robots, and modelling and control of extremely flexible space structures. The ARL has designed and built several semi-autonomous free-flying robots that perform numerous tasks in a zero-gravity, drag-free, two-dimensional environment. It is envisioned that future generations of these robots will be part of a human-robot team, in which the robots will operate under the task-level commands of astronauts. To make this possible, the ARL has developed a graphical user interface (GUI) with an intuitive object-level motion-direction capability. Using this interface, the ARL has demonstrated autonomous navigation, intercept and capture of moving and spinning objects, object transport, multiple-robot cooperative manipulation, and simple assemblies from both free-flying and fixed bases. The ARL has also built a number of experimental test beds on which the modelling and control of flexible manipulators has been studied. Early ARL experiments in this arena demonstrated for the first time the capability to control the end-point position of both single-link and multi-link flexible manipulators using end-point sensing. Building on these accomplishments, the ARL has been able to control payloads with unknown dynamics at the end of a flexible manipulator, and to achieve high-performance control of a multi-link flexible manipulator.

Ballhaus, W. L.↗

NASA NTRS1995

Robust integration schemes for generalized viscoplasticity with internal-state variables. Part 2: Algorithmic developments and implementation

This two-part report is concerned with the development of a general framework for the implicit time-stepping integrators for the flow and evolution equations in generalized viscoplastic models. The primary goal is to present a complete theoretical formulation, and to address in detail the algorithmic and numerical analysis aspects involved in its finite element implementation, as well as to critically assess the numerical performance of the developed schemes in a comprehensive set of test cases. On the theoretical side, the general framework is developed on the basis of the unconditionally-stable, backward-Euler difference scheme as a starting point. Its mathematical structure is of sufficient generality to allow a unified treatment of different classes of viscoplastic models with internal variables. In particular, two specific models of this type, which are representative of the present start-of-art in metal viscoplasticity, are considered in applications reported here; i.e., fully associative (GVIPS) and non-associative (NAV) models. The matrix forms developed for both these models are directly applicable for both initially isotropic and anisotropic materials, in general (three-dimensional) situations as well as subspace applications (i.e., plane stress/strain, axisymmetric, generalized plane stress in shells). On the computational side, issues related to efficiency and robustness are emphasized in developing the (local) interative algorithm. In particular, closed-form expressions for residual vectors and (consistent) material tangent stiffness arrays are given explicitly for both GVIPS and NAV models, with their maximum sizes 'optimized' to depend only on the number of independent stress components (but independent of the number of viscoplastic internal state parameters). Significant robustness of the local iterative solution is provided by complementing the basic Newton-Raphson scheme with a line-search strategy for convergence. In the present second part of the report, we focus on the specific details of the numerical schemes, and associated computer algorithms, for the finite-element implementation of GVIPS and NAV models.

Li, Wei↗

NASA NTRS2003

The Filtered Abel Transform and Its Application in Combustion Diagnostics

Many non-intrusive combustion diagnosis methods generate line-of-sight projections of a flame field. To reconstruct the spatial field of the measured properties, these projections need to be deconvoluted. When the spatial field is axisymmetric, commonly used deconvolution method include the Abel transforms, the onion peeling method and the two-dimensional Fourier transform method and its derivatives such as the filtered back projection methods. This paper proposes a new approach for performing the Abel transform method is developed, which possesses the exactness of the Abel transform and the flexibility of incorporating various filters in the reconstruction process. The Abel transform is an exact method and the simplest among these commonly used methods. It is evinced in this paper that all the exact reconstruction methods for axisymmetric distributions must be equivalent to the Abel transform because of its uniqueness and exactness. Detailed proof is presented to show that the two dimensional Fourier methods when applied to axisymmetric cases is identical to the Abel transform. Discrepancies among various reconstruction method stem from the different approximations made to perform numerical calculations. An equation relating the spectrum of a set of projection date to that of the corresponding spatial distribution is obtained, which shows that the spectrum of the projection is equal to the Abel transform of the spectrum of the corresponding spatial distribution. From the equation, if either the projection or the distribution is bandwidth limited, the other is also bandwidth limited, and both have the same bandwidth. If the two are not bandwidth limited, the Abel transform has a bias against low wave number components in most practical cases. This explains why the Abel transform and all exact deconvolution methods are sensitive to high wave number noises. The filtered Abel transform is based on the fact that the Abel transform of filtered projection data is equal to an integral transform of the original projection data with the kernel function being the Abel transform of the filtering function. The kernel function is independent of the projection data and can be obtained separately when the filtering function is selected. Users can select the best filtering function for a particular set of experimental data. When the kernal function is obtained, it can be used repeatedly to a number of projection data sets (rovs) from the same experiment. When an entire flame image that contains a large number of projection lines needs to be processed, the new approach significantly reduces computational effort in comparison with the conventional approach in which each projection data set is deconvoluted separately. Computer codes have been developed to perform the filter Abel transform for an entire flame field. Measured soot volume fraction data of a jet diffusion flame are processed as an example.

Simons, Stephen N.↗

NASA NTRS1998

Tone Noise and Nearfield Pressure Produced by Jet-Cavity Interaction

Cavity flow resonance can cause numerous problems in aerospace applications. While our long-term goal is to understand cavity flows well enough to devise effective cavity resonance suppression techniques, this paper describes a fundamental study of resonant tones produced by jet-cavity interaction at subsonic and supersonic speeds. Our specific jet-cavity configuration can also be used as a test bed for evaluating active and passive flow resonance control concepts. Two significant findings emerge from this study. 1) Originally, we expected that tones produced by jet-cavity interaction would resemble cavity tones or jet tones or would involve some simple combinations of each. The experimental data do not support these expectations: instead, the jet cavity interaction produce a unique set of tones. We propose simple yet and physically insightful correlations for these tones. Although the pressure patterns on the cavity floor display very complex variations with the Mach number for a length/depth = 8 cavity, the tones correspond to the acoustic modes of the cavity-independent of flow. For a length/ depth = 3 cavity, however, a surprise emerges: the pressure patterns on the cavity floor are not so complex but the tones depend significantly on the flow. Additionally, we examine the role of external feedback unique to jet-cavity interaction. 2) Previous research led us to expect that traditional classifications (open, transitional, or closed) for cavities in an infinite flight stream would be insensitive to small changes in Mach number and would depend primarily on cavity length/depth ratios. Use of the novel high resolution photoluminescent pressure sensitive paint shows that the classifications are actually quite sensitive to jet Mach number for a length/depth = 8 cavity. However, these classifications provide no guidance whatsoever for tone amplitude or frequency. Detailed experimental data and insights presented here will assist researchers who are performing numerical simulations of jet-cavity flows as a first step toward devising resonance suppression methods.

Raman, Ganesh↗

DOE OSTI2026

Axisymmetric hydrodynamics in numerical relativity: treating coordinate singularity, artificial heating and modeling MHD instabilities

Two-dimensional axisymmetric simulations of binary neutron star (BNS) merger remnant are a cheap alternative to 3D simulations. To maintain realism for secular timescales, simulations must avoid accumulated errors from drifts in conserved quantities and artificial heating, and they must model turbulent transport in a way that remains plausible throughout the evolution. It is also crucial to avoid numerical artifacts due to the polar coordinate axis singularity. Methods that behave well near the axis often break flux-conservative form of the hydrodynamic equations, resulting in significant drifts in conserved quantities. We present a flux-conservative scheme that maintains smoothness near the axis without sacrificing conservative formulation of the equations or incurring drifts in conserved global quantities. We compare the numerical performance of different treatments of the hydrodynamic equations when evolving a hypermassive neutron star resembling the remnant of a BNS merger. These simulations demonstrate that the new scheme combines the axis smoothness of non-conservative methods with the mass and angular momentum conservation of other conservative methods on $\sim 10^2$ ms timescales of viscous and neutrino-driven evolution. Because fluid profiles remain smooth in the remnant interior, it is possible to remove artificial heating by evolving the entropy density. We show how physical heating and cooling terms can be easily calculated from source terms of the conservative evolution variables and demonstrate our implementation. Finally, we discuss and implement improvements to the effective viscosity scheme to better model the effect of magnetohydrodynamic instabilities as the remnant evolves.

axisymmetry↗

NASA NTRS2016

Exercise Effects on the Brain and Sensorimotor Function in Bed Rest

Long duration spaceflight microgravity results in cephalad fluid shifts and deficits in posture control and locomotion. Effects of microgravity on sensorimotor function have been investigated on Earth using head down tilt bed rest (HDBR). HDBR serves as a spaceflight analogue because it mimics microgravity in body unloading and bodily fluid shifts. Preliminary results from our prior 70 days HDBR studies showed that HDBR is associated with focal gray matter (GM) changes and gait and balance deficits, as well as changes in brain functional connectivity. In consideration of the health and performance of crewmembers we investigated whether exercise reduces the effects of HDBR on GM, functional connectivity, and motor performance. Numerous studies have shown beneficial effects of exercise on brain health. We therefore hypothesized that an exercise intervention during HDBR could potentially mitigate the effects of HDBR on the central nervous system. Eighteen subjects were assessed before (12 and 7 days), during (7, 30, and ~70 days) and after (8 and 12 days) 70 days of 6-degrees HDBR at the NASA HDBR facility in UTMB, Galveston, TX, US. Each subject was randomly assigned to a control group or one of two exercise groups. Exercise consisted of daily supine exercise which started 20 days before the start of HDBR. The exercise subjects participated either in regular aerobic and resistance exercise (e.g. squat, heel raise, leg press, cycling and treadmill running), or aerobic and resistance exercise using a flywheel apparatus (rowing). Aerobic and resistance exercise intensity in both groups was similar, which is why we collapsed the two exercise groups for the current experiment. During each time point T1-weighted MRI scans and resting state functional connectivity scans were obtained using a 3T Siemens scanner. Focal changes over time in GM density were assessed using voxel based morphometry (VBM8) under SPM. Changes in resting state functional connectivity was assessed using both a region of interest (ROI, or seed-to-voxel) approach as well as a whole brain intrinsic connectivity (i.e., voxel-to-voxel) analysis. For the ROI analysis we selected 11 ROIs of brain regions that are involved in sensorimotor function (i.e., L. Insular C., L. Putamen, R. Premotor C., L.+R. Primary Motor C., R. Vestibular C., L. Posterior Cingulate G., R. Cerebellum Lobule V + VIIIb + Crus I, and the R. Superior Parietal G.) and correlated their time course of brain activation during rest with all other voxels in the brain. The whole brain connectivity analysis tests changes in the strength of the global connectivity pattern between each voxel and the rest of the brain. Functional mobility was assessed using an obstacle course. Vestibular contribution to balance was measured using Neurocom Sensory Organization Test 5. Behavioral measures were assessed pre-HDBR, and 0, 8 and 12 days post-HDBR. Linear mixed models were used to test for effects of time, group, and group-by-time interactions. Family-wise error corrected VBM revealed significantly larger increases in GM volume in the right primary motor cortex in bed rest control subjects than in bed rest exercise subjects. No other significant group by time interactions in gray matter changes with bed rest were observed. Functional connectivity MRI revealed that the increase in connectivity during bed rest of the left putamen with the bilateral midsagittal precunes and the right cingulate gyrus was larger in bed rest control subjects than in bed rest exercise subjects. Furthermore, the increase in functional connectivity with bed rest of the right premotor cortex with the right inferior frontal gyrus and the right primary motor cortex with the bilateral premotor cortex was smaller in bed rest control subjects than in bed rest exercise subjects. Functional mobility performance was less affected by HDBR in exercise subjects than in control subjects and post HDBR exercise subjects recovered faster than control subjects. The group performance differences and GM changes were not related. Exercise in HDBR partially mitigates the adverse effect of HDBR on functional mobility, particularly during the post-bed rest recovery phase. In addition, exercise appears to result in differential brain structural and functional changes in motor regions such as the primary motor cortex, the premotor cortex and the putamen. Whether these central nervous system changes are related to motor behavioral changes including gait and balance warrants further research.

Koppelmans, V.↗