Search NASA⌕ Search

SEARCH · Search NASA

Results for “RELAXATION METHOD (MATHEMATICS)”

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 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Automatic Abstraction in Planning

Traditionally, abstraction in planning has been accomplished by either state abstraction or operator abstraction, neither of which has been fully automatic. We present a new method, predicate relaxation, for automatically performing state abstraction. PABLO, a nonlinear hierarchical planner, implements predicate relaxation. Theoretical, as well as empirical results are presented which demonstrate the potential advantages of using predicate relaxation in planning. We also present a new definition of hierarchical operators that allows us to guarantee a limited form of completeness. This new definition is shown to be, in some ways, more flexible than previous definitions of hierarchical operators. Finally, a Classical Truth Criterion is presented that is proven to be sound and complete for a planning formalism that is general enough to include most classical planning formalisms that are based on the STRIPS assumption.

LOGIC PROGRAMMING↗
NASA NTRS2016

Cryovolcanic Emplacement of Domes on Europa

Here we explore the hypothesis that certain domes on Europa may have been produced by the extrusion of viscous cryolavas. A new mathematical method for the emplacement and relaxation of viscous lava domes is presented and applied to putative cryovolcanic domes on Europa. A similarity solution approach is applied to the governing equation for fluid flow in a cylindrical geometry, and dome relaxation is explored assuming a volume of cryolava has been rapidly emplaced onto the surface. Nonphysical sin- gularities inherent in previous models for dome relaxation have been eliminated, and cryolava cooling is represented by a time-variable viscosity. We find that at the onset of relaxation, bulk kinematic viscosities may lie in the range between 10(exp 3) and 10(exp 6) sq m/s, while the actual fluid lava viscosity may be much lower. Plausible relaxation times to form the domes, which are linked to bulk cryolava rheology, are found to range from 3.6 days to 7.5 years. We find that cooling of the cryolava, while dominated by conduction through an icy skin, should not prevent fluids from advancing and relaxing to form domes within the timescales considered. Determining the range of emplacement conditions for putative cryolava domes will shed light on Europa's resurfacing history. In addition, the rheologies and compositions of erupted cryolavas have implications for subsurface cryomagma ascent and local surface stress conditions on Europa.

volcanism↗
NASA NTRS1992

Uncertainty management by relaxation of conflicting constraints in production process scheduling

Mathematical-analytical methods as used in Operations Research approaches are often insufficient for scheduling problems. This is due to three reasons: the combinatorial complexity of the search space, conflicting objectives for production optimization, and the uncertainty in the production process. Knowledge-based techniques, especially approximate reasoning and constraint relaxation, are promising ways to overcome these problems. A case study from an industrial CIM environment, namely high-grade steel production, is presented to demonstrate how knowledge-based scheduling with the desired capabilities could work. By using fuzzy set theory, the applied knowledge representation technique covers the uncertainty inherent in the problem domain. Based on this knowledge representation, a classification of jobs according to their importance is defined which is then used for the straightforward generation of a schedule. A control strategy which comprises organizational, spatial, temporal, and chemical constraints is introduced. The strategy supports the dynamic relaxation of conflicting constraints in order to improve tentative schedules.

Dorn, Juergen↗
NASA NTRS1999

Phenomenological Theory of the Translational Relaxation Times in Gases

The exact solution to the classical equations governing the translational dispersion and absorption of sound in a gas obscures its relaxational character because of its mathematical complexity. The approach taken here is to solve the secular equation by the method of Pade approximants, which even to the relatively low order R(sub 11) yields a remarkably close approximation to the exact solution over a wide range of frequency/pressure (f/P) ratios. As a result, translational relaxation can be formulated in terms of a conventional relaxation process with well-defined relaxation times, relaxation strength, collision numbers, additivity relations, etc. To extend the theory to high values of f/P ratio, a model is proposed to account for the noncontinuum behavior of the transport coefficients (viscosity and thermal conductivity) as the molecular mean free path approaches the acoustical enclosure dimensions. The theoretical dispersion and absorption show good agreement with measurements in argon over the classical and transition regions of f/P, but a discrepancy appears at higher values of f/P, where collective propagating modes, assumed in the theory, give way to single-particle modes, prevailing in the experiments.

Zuckerwar, Allan J.↗
NASA NTRS1973

Analysis of embedded shock waves calculated by relaxation methods.

The requirements for uniqueness of the calculated jump conditions across embedded shock waves are investigated for type-dependent difference systems used in transonic flow studies. A mathematical analysis shows that sufficient conditions are (1) the equations should be differenced in conservative form and (2) a special difference operator should be used when switching from a hyperbolic to an elliptic operator. The latter results in a consistency condition on the integral equations, rather than the differential, at these points. Calculated jump conditions for several embedded and detached shock waves are analyzed in the physical and hodograph planes. Comparisons are made with previous results, a time-dependent calculation, and data.

Murman, E. M.↗
NASA NTRS1991

Transient modeling/analysis of hyperbolic heat conduction problems employing mixed implicit-explicit alpha method

This paper describes the evaluation of mixed implicit-explicit finite element formulations for hyperbolic heat conduction problems involving non-Fourier effects. In particular, mixed implicit-explicit formulations employing the alpha method proposed by Hughes et al. (1987, 1990) are described for the numerical simulation of hyperbolic heat conduction models, which involves time-dependent relaxation effects. Existing analytical approaches for modeling/analysis of such models involve complex mathematical formulations for obtaining closed-form solutions, while in certain numerical formulations the difficulties include severe oscillatory solution behavior (which often disguises the true response) in the vicinity of the thermal disturbances, which propagate with finite velocities. In view of these factors, the alpha method is evaluated to assess the control of the amount of numerical dissipation for predicting the transient propagating thermal disturbances. Numerical test models are presented, and pertinent conclusions are drawn for the mixed-time integration simulation of hyperbolic heat conduction models involving non-Fourier effects.

Tamma, Kumar K.↗
DOE OSTI2026

Asymptotic-preserving semi-implicit finite volume scheme for extended magnetohydrodynamics

A Finite Volume (FV) scheme is developed for solving the extended magnetohydrodynamic (XMHD) equations, yielding accurate results in the ideal, resistive, and Hall MHD limits. This is accomplished by first re-writing the XMHD equations such that it allows the algorithm to retain the use of ideal MHD Riemann solvers and the constrained transport method to preserve divergence-free magnetic fields. Incorporation of electron inertia and displacement current introduces additional numerical stiffness which motivates a semi-implicit FV scheme that re-formulates the XMHD model as a relaxation system. The equations are then advanced in time using an explicit 2nd-order Runge–Kutta scheme with operator splitting applied to the implicit source term updates at each sub-stage. For additional numerical stability, a density-dependent slope limiter is implemented to increase flux diffusivity at low density regions where non-ideal effects become significant. The algorithm is subsequently implemented in a scalable adaptive mesh refinement (AMR) framework. As the new algorithm retains many aspects of the ideal MHD formulations, it asymptotes naturally to the ideal MHD limit. Moreover, it shows promising results at the resistive and Hall MHD limits. This is verified against reference test problems for ideal, resistive and Hall MHD.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Multigrid Reduction in Time for Chaotic and Hyperbolic Problems (Final Report)

The coming massive parallelism of exascale computing presents a pressing challenge for the many DOE simulations of time-dependent partial differential equations (PDEs), which typically use traditional sequential time stepping methods. Since this traditional approach is inherently serial, it presents a sequential bottleneck when moving to exascale computing, because future performance gains will come through greater concurrency, not faster clock speeds. Thus, the goal of this work is to research parallelism in time, i.e., methods that compute multiple time values simultaneously, not sequentially. The focus will be on hyperbolic and chaotic problems of interest to DOE, with the goal of enabling scalable simulations of time-dependent hyperbolic and chaotic problems on future architectures. The chosen methodology for solving these problems parallel-in-time is multigrid, because multigrid (when it works) is a powerful, optimal, and scalable solver for discretized PDEs. Multigrid is already commonly used in many DOE simulations for scalably and optimally solving space-only PDE problems. The areas of hyperbolic and chaotic problems are chosen because of their relevance to problems of programmatic interest to DOE. However, these problems are also well-known to be difficult for parallel-in-time methods, with the most common method, parareal, diverging in many cases. The current state of-the-art for parallel-in-time at LLNL is the multigrid reduction in time (MGRIT) XBraid package, which also struggles for such problems, while still showing some improvement over parareal. In summary, new methods are needed for an efficient parallel-in-time scheme for hyperbolic and chaotic problems, and this work shall research promising new multigrid methods in this area. In particular, this work shall continue researching the directions from the current collaboration with Dr. Falgout, which are laid out in the work Toward Parallel in Time for Chaotic Dynamical Systems and showed the first known results of a parallel-in-time speedup for a chaotic problem. This work outlines two key improvements to XBraid for chaotic problems, the so-called “theta” and “delta-correction” methods. Here, these two improvements will be further researched and improved (including with a new relaxation method inspired by on Least Squares Shadowing (LSS)) and explored for more complicated problems.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1976

Solutions to divergence form equations using the method of partial implicitization

Partial implicitization is well suited for use on vector-processing computers since the method is an explicit unconditionally stable numerical technique. So far equations cast in nondivergence form have been used, though equations cast in the divergence form would have inherently greater accuracy. Since the method of partial implicitization is applicable only to relaxing problems to their steady state and is not applicable to true transient problems, a simple mathematical manipulation can be carried out which will produce the desired partial implicitization solution to divergence form equations.

Graves, R. A., Jr.↗
NASA NTRS1985

A survey of unified constitutive theories

The state of the art of time temperature dependent elastic viscoplastic constitutive theories which are based on the unified approach werre assessed. This class of constitutive theories is characterized by the use of kinetic equations and internal variables with appropriate evolutionary equations for treating all aspects of inelastic deformation including plasticity, creep, and stress relaxation. More than 10 such unified theories which are shown to satisfy the uniqueness and stability criteria imposed by Drucker's postulate and Ponter's inequalities are identified. The theories are compared for the types of flow law, kinetic equation, evolutionary equation of the internal variables, and treatment of temperature dependence. The similarities and differences of these theories are outlined in terms of mathematical formulations and illustrated by comparisons of theoretical calculations with experimental results which include monotonic stress-strain curves, cyclic hysteresis loops, creep and stress relaxation rates, and thermomechanical loops. Numerical methods used for integrating these stiff time temperature dependent constitutive equations are reviewed.

Chan, K. S.↗
DOE OSTI2025

Near-ideal relaxed MHD in slab geometry

We investigate the solutions of the relaxed magnetohydrodynamic (MHD) model (RxMHD) of R. Dewar and Z. Qu. This model generalizes Taylor relaxation by including the ideal Ohm's law constraint using an augmented Lagrangian method, providing a pathway to extend the multi-region relaxed MHD (MRxMHD) model. We present the first numerical solution of the RxMHD model by Dewar and Qu, demonstrating that it is mathematically well-defined and computationally feasible for constructing MHD equilibria in slab geometry. We also show that a cross-field flow can exist without enforcing an arbitrary constraint on the angular momentum, as is done in the case of MRxMHD with flow. Our results also demonstrate the self-organization of fully relaxed regions during the optimization, which was an important motivation behind developing this model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1991

Characterization of frequency standard instability by estimation of their covariance matrix

The popular 3-cornered hat method used for evaluating the noise contributions of individual frequency standards is revisited. This method is used in several cases, but sometimes the results are not consistent because one or more estimated clock variances turn out to be negative. Different causes of this unacceptable result have been conjectured: among them one regards the hypothesis of uncorrelated clocks, essential in this method. Since recently realistic cases of correlation between clocks, mainly due to the environmental conditions, have been observed, this paper proposes an entirely revisited version of the 3-cornered hat method which permits to evaluate the individual variances and also the possible covariances between clocks, by relaxing the hypothesis of uncorrelation. The uncertainty and the lack of contemporaneity of the measurement series are assumed to be negligible. The lack of the uncorrelation hypothesis calls for a more general mathematical model leading to an underdetermined linear system. The estimates of the (co)variances of the measurement series us well us those of the individual clocks are introduced by means of the scalar product of the related time series and arranged in the respective covariance matrices S and R. Since covariance matrix is positive definite by definition, the problem consists in estimating the unknown R, subject to the constraint of positive definiteness, from the known S. Unfortunately, this constraint is not sufficient to estimate R. Therefore a suitable optimization criterion is proposed, which assures the positive definiteness of R and, at the same time, minimizes the global correlation among clocks. Examples of frequency instability measurements processed by the "classical" 3-cornered hat method and the here-revisited method are presented showing that the solutions are identical only when the uncorrelation hypothesis doesn't violate the positive definiteness of R.

PatriziaTravella↗
NASA NTRS1967

Functions of relaxed controls

Mathematical control theory problems involving solutions of certain partial differential equations, nonadditive set functions, or other functionals - approximation and existence theorems

FUNCTIONAL ANALYSIS↗
NASA NTRS1975

Theoretical effect of modifications to the upper surface of two NACA airfoils using smooth polynomial additional thickness distributions which emphasize leading edge profile and which vary linearly at the trailing edge

An investigation was conducted on a CDC 7600 digital computer to determine the effects of additional thickness distributions to the upper surface of airfoils. The additional thickness distribution had the form of a continuous mathematical function which disappears at both the leading edge and the trailing edge. Results were obtained at a Mach number of 0.2 with an angle of attack of 6 deg. All calculations employed the full potential flow equations for two dimensional flow. The relaxation method of Jameson was used for solution of the potential flow equations. It is shown that increasing the thickness and variations in shape increases the lift and the adverse pitching moment coefficients.

Hague, D. S.↗
DOE OSTI2024

Semi-Lagrangian nodal discontinuous Galerkin method for the BGK model

In this work, we propose a semi-Lagrangian (SL) nodal discontinuous Galerkin (DG) solver for the BGK equation. The BGK model was introduced by Bhatnagar, Gross, and Krook [1] as a relaxation model for the fundamental Boltzmann equation [5], which describes the kinetic dynamic of rarefied gases with a probability distribution function. The challenges of designing efficient numerical schemes for the Boltzmann equation mainly come from its high dimensionality and complicated nonlinear collision operator. The BGK model gains interests since it has much lower computational cost, due to the relatively simple structure of the relaxation operator in replacement of the collision operator, while simultaneously preserving several important physical properties, such as macroscopic quantities and dissipation of entropy.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1976

Transonic flow calculations

The development of relaxation methods to calculate transonic flows is discussed. Rather accurate predictions can be made for a number of flows of interest using the transonic potential flow equation, which may be derived from the Euler equations for inviscid compressible flow by introducing the assumption that the flow is irrotational. Once the choice of a mathematical model has been settled, the numerical procedure for actually computing a solution contains two main elements: the construction of a discrete approximation which converges to the solution of the continuous problem in the limit as the mesh width is reduced to zero, and the solution of the resulting set of nonlinear difference equations by a convergent iterative scheme. The choice of an appropriate coordinate system and its influence on the accuracy of the discrete approximation is discussed. Several applications of the general method are described.

Jameson, A.↗
NASA NTRS1982

A finite element formulation of Euler equations for the solution of steady transonic flows

The main objective of the considered investigation is related to the development of a relaxation scheme for the analysis of inviscid, rotational, transonic flow problems. To formulate the equations of motion for inviscid flows in a fixed coordinate system, an Eulerian type variational principle is required. The derivation of an Eulerian variational principle which is employed in the finite element formulation is discussed. The presented numerical method describes the mathematical formulation and the application of a numerical process for the direct solution of steady Euler equations. The development of the procedure as an extension of existing potential flow formulations provides the applicability of previous procedures, e.g., proper application of the artificial viscosity for supersonic elements, and the accurate modeling of the shock.

Ecer, A.↗
NASA NTRS1982

Investigation of rotational transonic flows through ducts using a finite element scheme

An application of the finite element method is presented in order to study the flow through a two-dimensional channel including the choked and nearly-choked flow conditions. The mathematical procedure provides a combined treatment of potential and Euler equations, where the steady Euler equations are solved through the integration of a pseudo-time system which is equivalent to a relaxation scheme; the isentropic, potential flow equations are embedded in this formulation. To analyze the transonic flow through a channel with near-choked flow conditions, the flow through a parallel channel with a 10 percent circular arc bump is calculated as a sample problem.

Akay, H. U.↗