Search NASA⌕ Search

SEARCH · Search NASA

Results for “continuous relaxation”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Predicting long-term stress relaxation on Alloy 709 using the crystal plasticity finite element method

This report summarizes the results of physics-based, crystal plasticity simulations for the long-term stress relaxation behavior of Alloy 709. The purpose of the study was to provide insight into five key questions related to long-term behavior in high temperatures materials which are difficult or impossible to answer experimentally: (1) is there a threshold stress for long-term relaxation? (2) is there strain threshold for relaxation damage, below which significant damage does not accumulate? (3) does damage continue to accumulate as the material relaxes or will damage accumulation plateau under some loading conditions? (4) does stress relaxation loading inevitably lead to failure? and (5) which, if any, engineering models for relaxation damage accumulation reasonably match the simulation results? The report summarizes the numerical simulations used to address these five questions and provides at least partial answers to each question.

36 MATERIALS SCIENCE↗
NASA NTRS1998

Improved Modeling of Finite-Rate Turbulent Combustion Processes in Research Combustors

The objective of this thesis is to further develop and test a stochastic model of turbulent combustion in recirculating flows. There is a requirement to increase the accuracy of multi-dimensional combustion predictions. As turbulence affects reaction rates, this interaction must be more accurately evaluated. In this work a more physically correct way of handling the interaction of turbulence on combustion is further developed and tested. As turbulence involves randomness, stochastic modeling is used. Averaged values such as temperature and species concentration are found by integrating the probability density function (pdf) over the range of the scalar. The model in this work does not assume the pdf type, but solves for the evolution of the pdf using the Monte Carlo solution technique. The model is further developed by including a more robust reaction solver, by using accurate thermodynamics and by more accurate transport elements. The stochastic method is used with Semi-Implicit Method for Pressure-Linked Equations. The SIMPLE method is used to solve for velocity, pressure, turbulent kinetic energy and dissipation. The pdf solver solves for temperature and species concentration. Thus, the method is partially familiar to combustor engineers. The method is compared to benchmark experimental data and baseline calculations. The baseline method was tested on isothermal flows, evaporating sprays and combusting sprays. Pdf and baseline predictions were performed for three diffusion flames and one premixed flame. The pdf method predicted lower combustion rates than the baseline method in agreement with the data, except for the premixed flame. The baseline and stochastic predictions bounded the experimental data for the premixed flame. The use of a continuous mixing model or relax to mean mixing model had little effect on the prediction of average temperature. Two grids were used in a hydrogen diffusion flame simulation. Grid density did not effect the predictions except for peak temperature and tangential velocity. The hybrid pdf method did take longer and required more memory, but has a theoretical basis to extend to many reaction steps which cannot be said of current turbulent combustion models.

VanOverbeke, Thomas J.↗
NASA NTRS1993

Computation of axisymmetric and ionized flows using particle and continuum methods

Comparisons between particle and continuum simulations of hypersonic near-continuum flows are presented. The particle approach employs the direct simulation Monte Carlo method (DSMC), and the continuum approach solves the Euler equations. Both simulations have thermochemistry models for air implemented including ionization. A new axisymmetric DSMC code which is efficiently vectorized is developed for this study. In this DSMC code, particular attention is paid to matching the relaxation rates employed in the continuum approach. This investigation represents a continuation of a previous study which considered thermochemical relaxation in one-dimensional shock waves of nitrogen. Comparison of the particle and continuum methods is first made for an axisymmetric blunt-body flow of air at 7 km/s. Very good agreement is obtained for the two solutions. The two techniques also compare well for a one-dimensional shock wave in air at 10 km/s. In both applications, the results are found to be sensitive to various aspects of the chemistry models employed.

Boyd, Iain D.↗
DOE OSTI2024

Multi-parametric analysis for mixed integer linear programming: An application to transmission upgrade and congestion management

Upgrading the capacity of existing transmission lines is essential for meeting the growing energy demands, facilitating the integration of renewable energy, and ensuring the security of the transmission system. This study focuses on the selection of lines whose capacities and by how much should be expanded from the perspective of the Independent System Operators (ISOs) to minimize the total system cost. We employ advanced multi-parametric programming and an enhanced branch-and-bound algorithm to address complex mixed-integer linear programming (MILP) problems, considering multi-period time constraints and physical limitations of generators and transmission lines. To characterize the various decisions in transmission expansion, we model the increased capacity of existing lines as parameters within a specified range. This study first relaxes the binary variables to continuous variables and applies the Lagrange method and Karush-Kuhn-Tucker (KKT) conditions to obtain optimal solutions and identify critical regions associated with active and inactive constraints. Moreover, we extend the traditional branch-and-bound (B&B) method by determining the problem’s upper and lower bounds at each node of the B&B decision tree, helping to manage computational challenges in large-scale MILP problems. Here, we compare the difference between the upper and lower bounds to obtain an approximate optimal solution within the decision-makers’ tolerable error range. In addition, the first derivative of the objective function on the parameters of each line is used to inform the selection of lines for easing congestion and maximizing social welfare. Finally, the capacity upgrades are selected by weighing the reductions in system costs against the expense of upgrading line capacities. The findings are supported by numerical simulations and provide transmission-line planners with decision-making guidance.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS1990

Mechanically induced orientation of adult rat cardiac myocytes in vitro

The present study describes the spatial orientation of a population of freshly isolated adult rat cardiac myocytes using a computerized mechanical cell stimulator device for tissue cultured cells. A continuous unidirectional stretch of the substratum at 60 to 400 microns/min for 120 to 30 min, respectively, during the cell attachment period in a serum-free medium was found to induce a significant threefold increase in the number of rod-shaped myocytes oriented parallel to the direction of movement. The myocytes orient less well with unidirectional substratum stretching after their adhesion to the substratum. Adult myocytes plated onto a substratum undergoing continuous 10-percent stretch-relaxation cycling show no significant change in the myocyte orientation or cytoskeletal organization. In addition to the type of mechanical activity, orientation of rod-shaped myocytes is dependent on the speed of the substratum, the final stretch amplitude, and the timing between initiation of substratum stretching and adhesion of myocytes to the substratum.

Samuel, J.-L.↗
NASA NTRS2023

Equilibrium Selection via Current Sheet Relaxation and Guide Field Amplification

Although there is a continuous spectrum of current sheet equilibria, how a particular equilibrium is selected by a given system remains a mystery. Yet, only a limited number of equilibrium solutions are used for analyses of magnetized plasma phenomena. Here we present the exact process of equilibrium selection, by analyzing the relaxation process of a disequilibrated current sheet under a finite guide field. It is shown via phase-space analyses and particle-in-cell simulations that the current sheet relaxes in such a way that the guide field is locally amplified, yielding a mixed equilibrium from the spectrum. Comparisons to spacecraft observations and solar wind current sheet statistics demonstrate that such mixed equilibria are ubiquitous and exist as underlying local structures in various physical environments.

Young-Dae Yoon↗
NASA NTRS2011

Lessons Learned from the Space Shuttle Engine Cutoff System (ECO) Anomalies

The Space Shuttle Orbiter's main engine cutoff (ECO) system first failed ground checkout in April, 2005 during a first tanking test prior to Return-to-Flight. Despite significant troubleshooting and investigative efforts that followed, the root cause could not be found and intermittent anomalies continued to plague the Program. By implementing hardware upgrades, enhancing monitoring capability, and relaxing the launch rules, the Shuttle fleet was allowed to continue flying in spite of these unexplained failures. Root cause was finally determined following the launch attempts of STS-122 in December, 2007 when the anomalies repeated, which allowed drag-on instrumentation to pinpoint the fault (the ET feedthrough connector). The suspect hardware was removed and provided additional evidence towards root cause determination. Corrective action was implemented and the system has performed successfully since then. This white paper presents the lessons learned from the entire experience, beginning with the anomalies since Return-to-Flight through discovery and correction of the problem. To put these lessons in better perspective for the reader, an overview of the ECO system is presented first. Next, a chronological account of the failures and associated investigation activities is discussed. Root cause and corrective action are summarized, followed by the lessons learned.

Martinez, Hugo E.↗
DOE OSTI2022

High-Temperature 2D Optical Relaxation Visualizes Enhanced Oxygen Exchange Kinetics at Metal-Mixed Conducting Oxide Interfaces

Solid-state heterointerfaces are of interest for emergent local behavior that is distinct from either bulk parent compound. One technologically relevant example is the case of mixed ionic/electronic conductor (MIEC)–metal interfaces, which play an important role in electrochemistry. Metal–MIEC composite electrodes can demonstrate improved catalytic activity vs single-phase MIECs, improving fuel cell efficiency. Similarly, MIEC surface reaction kinetics are often evaluated using techniques that place metal current collectors in contact with the surface under evaluation, potentially altering the response vs the native surface. Here, techniques enabling direct and local in situ observation of the behavior at and around such heterointerfaces are needed. Here, we develop a spatially resolved optical transmission relaxation (2D-OTR) method providing continuous evaluation of local, high-temperature, controlled atmosphere defect kinetics across a ~1 cm2 sample area simultaneously in a contact-free manner. We apply it to observe the spatial variance of oxygen incorporation and evolution rates at ~525–620 °C, in response to step changes in oxygen partial pressure, on MIEC SrTi 0.65 Fe 0.35 O 3–x films as a function of distance from porous Pt and Au layers. Using this model geometry, we find significant enhancements in kinetics adjacent to the metals that decay over a few millimeter distance. To extract kinetic parameters, we fit the short-term optical data (initial portion of relaxations) with an exponential decay function appropriate for surface-exchange-limited kinetics, yielding apparent surface exchange coefficients (k chem ) with spatial resolution, decreasing with distance from the metal. To understand the kinetic processes governing the complete (long-term) optical relaxations, we performed COMSOL simulations, which demonstrated that a combination of laterally varying k chem and in-plane diffusion controls the observed kinetics over the full time range. Further support for spatially varying kchem comes from demonstrations of changing surface and bulk chemistry vs distance from the metal–MIEC interface, by X-ray photoelectron and optical absorption spectroscopies, respectively. Although microporous Pt and Au are not excellent electrodes in isolation, both metals exert a synergistic effect on the oxygen surface exchange rate in the presence of the mixed conducting film.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A direct transcription-based multiple shooting formulation for dynamic optimization

The growing need for fast and efficient solution techniques for solving dynamic optimization problems is driven by a broad spectrum of applications in scheduling and control. We suggest a novel framework for dynamic optimization that utilizes a multiple shooting “backbone” with discrete rather than continuous subproblems, thereby eliminating need for repeated time-integration. A Lagrangian relaxation (LR)-based decomposition scheme is proposed, which dualizes the state continuity requirements between subproblems and enables parallel solution of the problem. We demonstrate the applicability of our method on two case studies: the Van der Pol oscillator and a batch reactor.

42 ENGINEERING↗
DOE OSTI2022

AMM: Adaptive Multilinear Meshes

Adaptive representations are increasingly indispensable for reducing the in-memory and on-disk footprints of large-scale data. Usual solutions are designed broadly along two themes: reducing data precision, e.g., through compression, or adapting data resolution, e.g., using spatial hierarchies. Additionally, recent research suggests that combining the two approaches, i.e., adapting both resolution and precision simultaneously, can offer significant gains over using them individually. However, there currently exist no practical solutions to creating and evaluating such representations at scale. In this work, we present a new resolution-precision-adaptive representation to support hybrid data reduction schemes and offer an interface to existing tools and algorithms. Through novelties in spatial hierarchy, our representation, Adaptive Multilinear Meshes (AMM), provides considerable reduction in the mesh size. AMM creates a piecewise multilinear representation of uniformly sampled scalar data and can selectively relax or enforce constraints on conformity, continuity, and coverage, delivering a flexible adaptive representation. AMM also supports representing the function using mixed-precision values to further the achievable gains in data reduction. We describe a practical approach to creating AMM incrementally using arbitrary orderings of data and demonstrate AMM on six types of resolution and precision datastreams. By interfacing with state-of-the-art rendering tools through VTK, we demonstrate the practical and computational advantages of our representation for visualization techniques. With an open-source release of our tool to create AMM, we make such evaluation of data reduction accessible to the community, which we hope will foster new opportunities and future data reduction schemes.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1973

Vibrational relaxation of carbon monoxide studied by two-wavelength infrared emission

Experimental results are presented for the vibrational relaxation of pure carbon monoxide behind incident shock waves over the temperature range 4000 to 6300 K. The data were obtained as infrared emission from the fundamental and overtone vibrational band systems (in some of the experiments the two-band systems were recorded simultaneously). The data are consistent with present theories for the vibrational relaxation of diatomic molecules and can be interpreted in terms of an initial Boltzmann vibrational distribution relaxing toward final equilibrium via a continuous sequence of intermediate Boltzmann distributions.

Chackerian, C., Jr.↗
NASA NTRS2014

Forced Reconnection in the Near Magnetotail: Onset and Energy Conversion in PIC and MHD Simulations

Using two-dimensional particle-in-cell (PIC) together with magnetohydrodynamic (MHD) Q1 simulations of magnetotail dynamics, we investigate the evolution toward onset of reconnection and the subsequent energy transfer and conversion. In either case, reconnection onset is preceded by a driven phase, during which magnetic flux is added to the tail at the high-latitude boundaries, followed by a relaxation phase, during which the configuration continues to respond to the driving. The boundary deformation leads to the formation of thin embedded current sheets, which are bifurcated in the near tail, converging to a single sheet farther out in the MHD simulations. The thin current sheets in the PIC simulation are carried by electrons and are associated with a strong perpendicular electrostatic field, which may provide a connection to parallel potentials and auroral arcs and an ionospheric signal even prior to the onset of reconnection. The PIC simulation very well satisfies integral entropy conservation (intrinsic to ideal MHD) during this phase, supporting ideal ballooning stability. Eventually, the current intensification leads to the onset of reconnection, the formation and ejection of a plasmoid, and a collapse of the inner tail. The earthward flow shows the characteristics of a dipolarization front: enhancement of Bz, associated with a thin vertical electron current sheet in the PIC simulation. Both MHD and PIC simulations show a dominance of energy conversion from incoming Poynting flux to outgoing enthalpy flux, resulting in heating of the inner tail. Localized Joule dissipation plays only a minor role.

two-dimensional↗
NASA NTRS2014

Onset of Reconnection in the near Magnetotail: PIC Simulations

Using 2.5-dimensional particle-in-cell (PIC) simulations of magnetotail dynamics, we investigate the onset of reconnection in two-dimensional tail configurations with finite Bz. Reconnection onset is preceded by a driven phase, during which magnetic flux is added to the tail at the high-latitude boundaries, followed by a relaxation phase, during which the configuration continues to respond to the driving. We found a clear distinction between stable and unstable cases, dependent on deformation amplitude and ion/electron mass ratio. The threshold appears consistent with electron tearing. The evolution prior to onset, as well as the evolution of stable cases, are largely independent of the mass ratio, governed by integral flux tube entropy conservation as imposed in MHD (magnetohydrodynamics). This suggests that ballooning instability in the tail should not be expected prior to the onset of tearing and reconnection. The onset time and other onset properties depend on the mass ratio, consistent with expectations for electron tearing. At onset,we found electron anisotropies T⊥⁄ T∥ (bottom tail divided by parallel tail) equals 1.1-1.3, raising growth rates and wavenumbers. Our simulations have provided a quantitative onset criterion that is easily evaluated in MHD simulations, provided the spatial resolution is sufficient. The evolution prior to onset and after the formation of a neutral line does not depend on the electron physics, which should permit an approximation by MHD simulations with appropriate dissipation terms.

PIC↗
NASA NTRS1995

Time-domain approach for the transient responses in stratified viscoelastic Earth models

We have developed the numerical algorithm for the computation of transient viscoelastic responses in the time domain for a radially stratified Earth model. Stratifications in both the elastic parameters and the viscosity profile have been considered. The particular viscosity profile employed has a viscosity maximum with a constrast of O(100) in the mid lower mantle. The distribution of relaxation times reveals the presence of a continuous spectrum situated between O(100) and O(exp 4) years. The principal mode is embedded within this continuous spectrum. From this initial-value approach we have found that for the low degree harmonics the non-modal contributions are comparable to the modal contributions. For this viscosity model the differences between the time-domain and normal-mode results are found to decrease strongly with increasing angular order. These calculations also show that a time-dependent effective relaxation time can be defined, which can be bounded by the relaxation times of the principal modes.

Hanyk, L.↗
NASA NTRS1987

Wall relaxation and the driving forces for cell expansive growth

When water uptake by growing cells is prevented, the turgor pressure and the tensile stress in the cell wall are reduced by continued wall loosening. This process, termed in vivo stress relaxation, provides a new way to study the dynamics of wall loosening and to measure the wall yield threshold and the physiological wall extensibility. Stress relaxation experiments indicate that wall stress supplies the mechanical driving force for wall yielding. Cell expansion also requires water absorption. The driving force for water uptake during growth is created by wall relaxation, which lowers the water potential of the expanding cells. New techniques for measuring this driving force show that it is smaller than believed previously; in elongating stems it is only 0.3 to 0.5 bar. This means that the hydraulic resistance of the water transport pathway is small and that rate of cell expansion is controlled primarily by wall loosening and yielding.

NASA Program Space Biology↗
DOE OSTI2020

UNRAVEL: A Robust Modular Velocity Dealiasing Technique for Doppler Radar

Unfold Radar Velocity (UNRAVEL) is an open-source modular Doppler velocity dealiasing algorithm for weather radars. UNRAVEL is an algorithm that does not need external reference velocity data, making it easily applicable. The proposed algorithm includes 11 core modules and 2 dealiasing strategies. UNRAVEL is an iterative algorithm. The goal is to build the dealiasing results starting with the strictest possible continuity tests in azimuth and range and, after each step, relaxing the parameters to include more results from a progressively growing number of reference points. UNRAVEL also has modules that perform 3D continuity checks. Thanks to this modular design, the number of dealiasing strategies can be expanded in order to optimize the dealiasing results. While the first driver dealiases Doppler velocity from each tilt independently from one another, the second driver also performs a three-dimensional continuity check of the velocity using successive elevations. The proposed dealiasing algorithm is tested using severe weather data from an S-band Doppler radar that have been aliased to mimic aliased radial velocity patterns that would be observed by a C-band Doppler radar. Artificially aliasing S-band data permits creation of a reference to which the performance of various dealiasing techniques can be compared. Comparisons show that UNRAVEL consistently outperforms other established dealiasing algorithms for the test period selected in this work.

47 OTHER INSTRUMENTATION↗
NASA NTRS1993

Next generation initiation techniques

Four-dimensional data assimilation strategies can generally be classified as either current or next generation, depending upon whether they are used operationally or not. Current-generation data-assimilation techniques are those that are presently used routinely in operational-forecasting or research applications. They can be classified into the following categories: intermittent assimilation, Newtonian relaxation, and physical initialization. It should be noted that these techniques are the subject of continued research, and their improvement will parallel the development of next generation techniques described by the other speakers. Next generation assimilation techniques are those that are under development but are not yet used operationally. Most of these procedures are derived from control theory or variational methods and primarily represent continuous assimilation approaches, in which the data and model dynamics are 'fitted' to each other in an optimal way. Another 'next generation' category is the initialization of convective-scale models. Intermittent assimilation systems use an objective analysis to combine all observations within a time window that is centered on the analysis time. Continuous first-generation assimilation systems are usually based on the Newtonian-relaxation or 'nudging' techniques. Physical initialization procedures generally involve the use of standard or nonstandard data to force some physical process in the model during an assimilation period. Under the topic of next-generation assimilation techniques, variational approaches are currently being actively developed. Variational approaches seek to minimize a cost or penalty function which measures a model's fit to observations, background fields and other imposed constraints. Alternatively, the Kalman filter technique, which is also under investigation as a data assimilation procedure for numerical weather prediction, can yield acceptable initial conditions for mesoscale models. The third kind of next-generation technique involves strategies to initialize convective scale (non-hydrostatic) models.

Warner, Tom↗
DOE OSTI2021

Antiferromagnetic ordering and possible lattice response to dynamic uranium valence in U 3 O 8

Determining the correct electronic structure of U 3 O 8 remains a formidable experimental and theoretical challenge. In the low-temperature phase, two crystallographic U sites are separated into a distinct 2U(V)+1U(VI) oxidation configuration. At low temperatures, the U(V) sites form a distorted honeycomb lattice, but the U(VI) sit on a triangular sublattice, suggesting potential for magnetic frustration effects. The spin configuration of the unpaired f electrons on the U(V) sites is likely antiferromagnetic (AFM) from susceptibility measurements, but this has not been confirmed. Here in this paper, we present a neutron scattering investigation of the structure and dynamics of U 3 O 8 from 1.7 to 600 K. We confirm static AFM ordering onset at between 22 and 25 K, which is present down to at least 1.7 K with AFM peaks corresponding to [0.5 1 1] and [0.5 2 2] in the orthorhombic phase. These measurements rule out static AFM order along the a axis of the Amm2 phase, a configuration previously suggested by theory. Above 100 K a quasielastic scattering channel opens that we speculate arises from a lattice relaxation response to thermally activated electron hopping. This term does not conform to a magnetic form factor, so it is not related to spin relaxations. If correct, this mechanism stabilizes a continuous valence transition from 2U(V)+1U(VI) in the low-temperature (T<600 K) orthorhombic phase to the hexagonal phase that contains only one degenerate U site, wherein the U valence can be dynamically stabilized between U(V)↔U(VI) by phonon-assisted electron hopping.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗