Search NASA⌕ Search

SEARCH · Search NASA

Results for “Relaxation method”

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 145 records · Page 8

Results of the integration of a transonic full-potential analysis program with a free-wake lifting-line program for hovering rotors

The hovering performance predictions of the TFAR1 and OPLIN codes and the experimental data are discussed. The TFAR1 program solves the full-potential equation in a rotor-fixed coordinate system by use of the line relaxation method. The OPLIN program calculates the positions of wake vortices and rotor performance using the influence-coefficient-and-lifing-line method. The two programs are combined by adding the induced velocities from the OPLIN code to the near flow-field of the rotor from the TFAR1 code. Results show that the TFAR1 program converges better with the downwash-coupling method than the twist correction method to include the wake downwash.

Chung, Song-Young↗

Crack growth in single-crystal silicon

Crack growth in single-crystal silicon at room temperature in air was evaluated by double torsion (DT) load-relaxation method and monitored by acoustic emission (AE) technique. Both DT and AE methods indicated lack of subcritical crack growth in silicon. At the critical stress intensity factor, the crack front was found to be jumping several times in a 'mirror' region and then followed by fast crack growth in a 'hackle' region. Hackle marks were found to be associated with plastic deformation at the tip of the fast moving crack. No dislocation etch pits were found in the 'mirror' region, in which crack growth may result from interatomic bonds broken at the crack tip under stress without any plastic deformation. Acoustic emission appears to be spontaneously generated from both interatomic bonds broken and dislocation generation at the moving crack tip during the crack growth in single-crystal silicon.

Chen, C. P.↗

Retrieval of constituent mixing ratios from limb thermal emission spectra

An onion-peeling iterative, least-squares relaxation method to retrieve mixing ratio profiles from limb thermal emission spectra is presented. The method has been tested on synthetic data, containing various amounts of added random noise for O3, HNO3, and N2O. The retrieval method is used to obtain O3 and HNO3 mixing ratio profiles from high-resolution thermal emission spectra. Results of the retrievals compare favorably with those obtained previously.

Shaffer, William A.↗

A review of reduced Navier-Stokes computations for compressible viscous flows

A reduced form of the Navier-Stokes equations, defined by a single composite of the Euler, boundary layer, and triple deck approximations, is considered for the computation of viscous interacting flows. Global pressure or pseudopotential relaxation methods are integrated with coupled sparse matrix direct solvers or coupled strongly implicit ILU algorithms to efficiently capture sharp shocks and regions of recirculation. Solutions are obtained for a variety of 2D and 3D geometries and for Mach numbers (M) spanning the range from incompressible (M = 0) to supersonic (M = 6) flow.

Rubin, S. G.↗

Crustal deformation measurements in central Japan determined by a Global Positioning System fixed-point network

Results are presented from temporally dense measurements of crustal deformation associated with the convergence of the Eurasian (EUR), Pacific, North American, and Philippine Sea (PHS) plates, carried out in April 1988 by a 10-station GPS fixed-point network established in central Japan. Using regional orbit relaxation methods, the analysis of the first 17-month data revealed significant horizontal deformation across the Suruga trough. Namely, it was found that a site in the northern tip of PHS plate moved nearly westward with a velocity of 28 +/-5 mm per year, and a site at the southeastern tip of EUR plate moved south-southwestward with a velocity of 18 +/-5 mm per year. A significant vertical uplift with a velocity of 20 mm/yr was detected at a site inland of the Tokai district located in the Akaishi uplift zone and at a site on the Hatsushima Island in Sagami Bay.

Shimada, Seiichi↗

Comprehensive Analyses of Data Collected from TEREK (Solar EUV Telescope) RES-C (Solar X-Ray Spectrometer) and SORS (Solar Radio Spectrometer) on board CORONAS-1 Using Magnetohydrodynamic Models

By using the observed magnetic field data obtained from the Wilcox Solar Observatory at Stanford University as the inputs to a two-dimensional plane-of-sky magnetohydrodynamic model, via numerical relaxation method, we have deduced the plasma and magnetic field parameters for the observed coronal hole by CORONAS-1. The method for this self-consistent MHD analysis will be discussed in detail. Numerical results for the magnetic field configuration, velocity distribution, density and temperature distributions will be presented. We have converted the computed density to polarization brightness in order to directly compare the MHD outputs with observations. Also included is a summary of achievements made during the grant period. This section is summarized into three categories: 1) Visit of Co-Investigators; 2) Presentations; and 3) Papers published, accepted and submitted for publication in journals.

Wu, S. T.↗

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↗

A Pseubo-Temporal Multi-Grid Relaxation Scheme for Solving the Parabolized Navier-Stokes Equations

A multi-grid, flux-difference-split, finite-volume code, VULCAN, is presented for solving the elliptic and parabolized form of the equations governing three-dimensional, turbulent, calorically perfect and non-equilibrium chemically reacting flows. The space marching algorithms developed to improve convergence rate and or reduce computational cost are emphasized. The algorithms presented are extensions to the class of implicit pseudo-time iterative, upwind space-marching schemes. A full approximate storage, full multi-grid scheme is also described which is used to accelerate the convergence of a Gauss-Seidel relaxation method. The multi-grid algorithm is shown to significantly improve convergence on high aspect ratio grids.

Morrison, J. H.↗

Calculation of Supersonic Combustion Using Implicit Schemes

One of the technology goals of NASA for advanced space transportation is to develop highly efficient propulsion systems to reduce the cost of payload for space missions. Developments of rockets for the second generation Reusable Launch Vehicle (RLV) in the past several years have been focused on low-cost versions of conventional engines. However, recent changes in the Integrated Space Transportation Program to build a crew transportation vehicle to extend the life of the Space Shuttle fleet might suggest that air-breathing rockets could reemerge as a possible propulsion system for the third generation RLV to replace the Space Shuttle after 2015. The weight of the oxygen tank exceeds thirty percent of the total weight of the Space Shuttle at launch while the payload is only one percent of the total weight. The air-breathing rocket propulsion systems, which consume oxygen in the air, offer clear advantages by making vehicles lighter and more efficient. Experience in the National Aerospace Plane Program in the late 1980s indicates that scramjet engines can achieve high specific impulse for low hypersonic vehicle speeds. Whether taking a form of Rocket Based Combined Cycle (RBCC) or Turbine Based Combined Cycle (TBCC), the scramjet is an essential mode of operation for air-breathing rockets. It is well known that fuel-air mixing and rapid combustion are of crucial importance for the success of scramjet engines since the spreading rate of the supersonic mixing layer decreases as the Mach number increases. A factored form of the Gauss-Seidel relaxation method has been widely used in hypersonic flow research since its first application to non-equilibrium flows. However, difficulties in stability and convergence have been encountered when there is strong interaction between fluid motion and chemical reaction, such as multiple fuel injection problems. The present paper reports the results from investigation of the effect of modifications to the original algorithm on the performance for multiple injectors.

Yoon, Seokkwan↗

Cooling of Gas Turbines: Computed Temperature Distribution Through Cross Section of Water-Cooled Turbine Blade - 6

A theoretical analysis of the cross-sectional temperature distribution of a water-cooled turbine blade was made using the relaxation method to solve the differential equation derived from the analysis. The analysis was applied to specific turbine blade and the studies icluded investigations of the accuracy of simple methods to determine the temperature distribution along the mean line of the rear part of the blade, of the possible effect of varying the perimetric distribution of the hot gas-to -metal heat transfer coefficient, and of the effect of changing the thermal conductivity of the blade metal for a constant cross sectional area blade with two quarter inch diameter coolant passages.

Livingood, John N. B.↗

A Pseudo-Temporal Multi-Grid Relaxation Scheme for Solving the Parabolized Navier-Stokes Equations

A multi-grid, flux-difference-split, finite-volume code, VULCAN, is presented for solving the elliptic and parabolized form of the equations governing three-dimensional, turbulent, calorically perfect and non-equilibrium chemically reacting flows. The space marching algorithms developed to improve convergence rate and or reduce computational cost are emphasized. The algorithms presented are extensions to the class of implicit pseudo-time iterative, upwind space-marching schemes. A full approximate storage, full multi-grid scheme is also described which is used to accelerate the convergence of a Gauss-Seidel relaxation method. The multi-grid algorithm is shown to significantly improve convergence on high aspect ratio grids.

White, J. A.↗

Relaxation Revisited: A Fresh Look at Multigrid for Steady Flows

The year 1971 saw the publication of one of the landmark papers in computational aerodynamics, that of Murman and Cole. As with many seminal works, its significance lies not so much in the specific problem that it addressed| small disturbance, plane transonic flow - but in the identification of a general approach to the solution of a technically important and theoretically difficult problem. The key features of Murman and Cole's work were the use of type- dependent differencing to correctly account for the proper domain of dependence of a mixed elliptic/hyperbolic equation, and the introduction of line relaxation to solve the steady flow equation. All subsequent work in transonic potential flows was based on these concepts. Jameson extended Murman and Cole's ideas to the full potential equation with two important contributions. First, he introduced the rotated difference stencil, which generalized the Murman and Cole type-dependent difference operator to general coordinates. Second, he used the interpretation, introduced by Garabedian, of relaxation as an iteration in artificial time to construct stable relaxation schemes, generalizing the original line relaxation method of Reference. The decade of the 1970s saw an explosion of activity in the solution of transonic potential flows, which has been summarized in the review article of Caughey.

Roberts, Thomas W.↗

MHD Modeling of the Solar Wind with Turbulence Transport and Heating

We have developed a magnetohydrodynamic model that describes the global axisymmetric steady-state structure of the solar wind near solar minimum with account for transport of small-scale turbulence associated heating. The Reynolds-averaged mass, momentum, induction, and energy equations for the large-scale solar wind flow are solved simultaneously with the turbulence transport equations in the region from 0.3 to 100 AU. The large-scale equations include subgrid-scale terms due to turbulence and the turbulence (small-scale) equations describe the effects of transport and (phenomenologically) dissipation of the MHD turbulence based on a few statistical parameters (turbulence energy, normalized cross-helicity, and correlation scale). The coupled set of equations is integrated numerically for a source dipole field on the Sun by a time-relaxation method in the corotating frame of reference. We present results on the plasma, magnetic field, and turbulence distributions throughout the heliosphere and on the role of the turbulence in the large-scale structure and temperature distribution in the solar wind.

Goldstein, M. L.↗

Three-Fluid Magnetohydrodynamic Modeling of the Solar Wind in the Outer Heliosphere

We have developed a three-fluid, fully three-dimensional magnetohydrodynamic model of the solar wind plasma in the outer heliosphere as a co-moving system of solar wind protons, electrons, and interstellar pickup protons, with separate energy equations for each species. Our approach takes into account the effects of electron heat conduction and dissipation of Alfvenic turbulence on the spatial evolution of the solar wind plasma and interplanetary magnetic fields. The turbulence transport model is based on the Reynolds decomposition of physical variables into mean and fluctuating components and uses the turbulent phenomenologies that describe the conversion of fluctuation energy into heat due to a turbulent cascade. We solve the coupled set of the three-fluid equations for the mean-field solar wind and the turbulence equations for the turbulence energy, cross helicity, and correlation length. The equations are written in the rotating frame of reference and include heating by turbulent dissipation, energy transfer from interstellar pickup protons to solar wind protons, and solar wind deceleration due to the interaction with the interstellar hydrogen. The numerical solution is constructed by the time relaxation method in the region from 0.3 to 100 AU. Initial results from the novel model are presented.

Usmanov, Arcadi V.↗

Numerical Solution of the Radiative Transfer Equation: X-Ray Spectral Formation from Cylindrical Accretion onto a Magnetized Neutron Star

Predicting the emerging X-ray spectra in several astrophysical objects is of great importance, in particular when the observational data are compared with theoretical models. This requires developing numerical routines for the solution of the radiative transfer equation according to the expected physical conditions of the systems under study. Aims. We have developed an algorithm solving the radiative transfer equation in the Fokker-Planck approximation when both thermal and bulk Comptonization take place. The algorithm is essentially a relaxation method, where stable solutions are obtained when the system has reached its steady-state equilibrium. Methods. We obtained the solution of the radiative transfer equation in the two-dimensional domain defined by the photon energy E and optical depth of the system pi using finite-differences for the partial derivatives, and imposing specific boundary conditions for the solutions. We treated the case of cylindrical accretion onto a magnetized neutron star. Results. We considered a blackbody seed spectrum of photons with exponential distribution across the accretion column and for an accretion where the velocity reaches its maximum at the stellar surface and at the top of the accretion column, respectively. In both cases higher values of the electron temperature and of the optical depth pi produce flatter and harder spectra. Other parameters contributing to the spectral formation are the steepness of the vertical velocity profile, the albedo at the star surface, and the radius of the accretion column. The latter parameter modifies the emerging spectra in a specular way for the two assumed accretion profiles. Conclusions. The algorithm has been implemented in the XPEC package for X-ray fitting and is specifically dedicated to the physical framework of accretion at the polar cap of a neutron star with a high magnetic field (approx > 10(exp 12) G). This latter case is expected to be of typical accreting systems such as X-ray pulsars and supergiant fast X ray transients.

Fairnelli, R.↗

Implicit Coupling Approach for Simulation of Charring Carbon Ablators

This study demonstrates that coupling of a material thermal response code and a flow solver with nonequilibrium gas/surface interaction for simulation of charring carbon ablators can be performed using an implicit approach. The material thermal response code used in this study is the three-dimensional version of Fully Implicit Ablation and Thermal response program, which predicts charring material thermal response and shape change on hypersonic space vehicles. The flow code solves the reacting Navier-Stokes equations using Data Parallel Line Relaxation method. Coupling between the material response and flow codes is performed by solving the surface mass balance in flow solver and the surface energy balance in material response code. Thus, the material surface recession is predicted in flow code, and the surface temperature and pyrolysis gas injection rate are computed in material response code. It is demonstrated that the time-lagged explicit approach is sufficient for simulations at low surface heating conditions, in which the surface ablation rate is not a strong function of the surface temperature. At elevated surface heating conditions, the implicit approach has to be taken, because the carbon ablation rate becomes a stiff function of the surface temperature, and thus the explicit approach appears to be inappropriate resulting in severe numerical oscillations of predicted surface temperature. Implicit coupling for simulation of arc-jet models is performed, and the predictions are compared with measured data. Implicit coupling for trajectory based simulation of Stardust fore-body heat shield is also conducted. The predicted stagnation point total recession is compared with that predicted using the chemical equilibrium surface assumption

Ablation↗

The Force-Free Magnetosphere of a Rotating Black Hole

We revisit the Blandford-Znajek process and solve the fundamental equation that governs the structure of the steady-state force-free magnetosphere around a Kerr black hole. The solution depends on the distributions of the magnetic field angular velocity and the poloidal electric current. These are not arbitrary. They are determined self-consistently by requiring that magnetic field lines cross smoothly the two singular surfaces of the problem: the inner "light surface" located inside the ergosphere and the outer "light surface" which is the generalization of the pulsar light cylinder.We find the solution for the simplest possible magnetic field configuration, the split monopole, through a numerical iterative relaxation method analogous to the one that yields the structure of the steady-state axisymmetric force-free pulsar magnetosphere. We obtain the rate of electromagnetic extraction of energy and confirm the results of Blandford and Znajek and of previous time-dependent simulations. Furthermore, we discuss the physical applicability of magnetic field configurations that do not cross both "light surfaces."

Contopoulos, Ioannis↗

Evaluation of Finite-Rate Gas/Surface Interaction Models for a Carbon Based Ablator

Two sets of finite-rate gas-surface interaction model between air and the carbon surface are studied. The first set is an engineering model with one-way chemical reactions, and the second set is a more detailed model with two-way chemical reactions. These two proposed models intend to cover the carbon surface ablation conditions including the low temperature rate-controlled oxidation, the mid-temperature diffusion-controlled oxidation, and the high temperature sublimation. The prediction of carbon surface recession is achieved by coupling a material thermal response code and a Navier-Stokes flow code. The material thermal response code used in this study is the Two-dimensional Implicit Thermal-response and Ablation Program, which predicts charring material thermal response and shape change on hypersonic space vehicles. The flow code solves the reacting full Navier-Stokes equations using Data Parallel Line Relaxation method. Recession analyses of stagnation tests conducted in NASA Ames Research Center arc-jet facilities with heat fluxes ranging from 45 to 1100 wcm2 are performed and compared with data for model validation. The ablating material used in these arc-jet tests is Phenolic Impregnated Carbon Ablator. Additionally, computational predictions of surface recession and shape change are in good agreement with measurement for arc-jet conditions of Small Probe Reentry Investigation for Thermal Protection System Engineering.

Ablation↗