Search NASA⌕ Search

SEARCH · Search NASA

Results for “Corrected equation”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

A Gauss-Radau-Laguerre Discrete Variable Representation for Use in Continuum Electron Dynamics

In this work, we detail an implementation, suitable for calculations on highly correlated ionizing systems, of a modified finite element discrete variable representation (FE-DVR) appended with a Gauss-Radau-Laguerre element. The appended element includes exterior complex scaling (ECS) to impose outgoing wave boundary conditions on treatments of processes involving continuum electrons. In this “infinite range” ECS (irECS), the complications that introduce reflections from the end of the grid when the last ECS finite element has finite range are avoided by the use of the Laguerre-weighted exponentially decaying tails, while outgoing wave boundary conditions are still imposed via the ECS transformation. For highly correlated systems in the absence of strong external fields we find that accurate two-electron integrals are essential in this modified FE-DVR. To accurately compute the two-electron integrals over the entire ECS contour, we present a detailed examination of the implications from the boundary terms that arise in a solution of Poisson’s equation with the Radau-Laguerre basis. A boundary term correction is necessary, and when included, the Radau-Laguerre DVR can accurately describe highly correlated states such as the doubly excited states of helium over the entire ECS contour.

elements↗
DOE OSTI2025

Nonresonant two-photon x-ray absorption in Cu

We present a real-space Green's function theory and calculations of two-photon x-ray absorption (TPA). Our focus is on nonresonant 𝐾-shell TPA in metallic Cu, which has been observed experimentally at intense x-ray free electron laser (XFEL) sources. The theory is based on an independent particle Green's function treatment of the Kramers-Heisenberg equation and an approximation for the sum over nonresonant intermediate states in terms of a static quadrupole transition operator. XFEL effects are modeled by a partially depleted 𝑑 band. This approach is shown to give results for 𝐾-shell TPA in quantitative agreement with XFEL experiment and with a Bethe-Salpeter equation approach. Furthermore, we also briefly discuss many-body corrections and TPA sum rules.

Approximation methods for many-body systems↗
NASA NTRS1975

Numerical procedure for analyzing Langmuir probe data

A numerical procedure is proposed for analyzing Langmuir probe data in the presence of a two-group plasma containing both primary and Maxwellian electrons. The procedure is known as a least-squares differential-correction technique for determining the unknown coefficients of the governing equation. It is shown that for a given set of input data the results of the analysis are unique and independent of the initial estimate of electron temperature, that convergence is fastest when electron temperature is overestimated, that the results are sensitive to the region of the curve used as input and also to the voltage increment, and that plasma properties determined by the proposed numerical procedure are either consistent with those determined graphically or closer to expected values. With a suitable data acquisition system, the Fortran IV program worked out for this procedure could be used to provide real-time plasma diagnostic information for an operating ion thruster.

Beattie, J. R.↗
NASA NTRS1980

A study of the effects of state transition matrix approximations

The effects of using an approximate state transition matrix in orbit estimation are investigated. The approximate state transition matrix results when higher order geopotential terms in the equations of motion are ignored in the formation of the variational equations. Two methods of orbit estimation are considered: the differential correction procedure (DC) and the extended Kalman filter (EKF). The system used for the study is the Research & Development version of the Goddard Trajectory Determination System. The effects of the approximation are analyzed on a number of orbits. These orbits of various inclinations and semimajor axes. Other parameters studied include geopotential models and DC arc length.

May, J. A.↗
NASA NTRS1981

The aerodynamics of circulation control

Two dimensional subsonic wind tunnel tests were conducted on a 20% thickness: chord ratio circulation controlled elliptic aerofoil section equipped with forward and reverse blowing slots. Overall performance measurements were made over a range of trailing edge blowing momentum coefficients from 0 to 0.04; some included the effect of leading edge blowing. A detailed investigation of the trailing edge wall jet, using split film probes, hot wire probes and total head tubes, provided measurements of mean velocity components, Reynolds normal and shear stresses, and radial static pressure. The closure of the two dimensional angular momentum and continuity equations was examined using the measured data, with and without correction, and the difficulty of obtaining a satisfactory solution illustrated. Suggestions regarding the nature of the flow field which should aid the understanding of Coanda effect and the theoretical solution of highly curved wall jet flows are presented.

Wood, N. J.↗
NASA NTRS1985

Development and application of a unified balancing approach with multiple constraints

The development of a general analytic approach to constrained balancing that is consistent with past influence coefficient methods is described. The approach uses Lagrange multipliers to impose orbit and/or weight constraints; these constraints are combined with the least squares minimization process to provide a set of coupled equations that result in a single solution form for determining correction weights. Proper selection of constraints results in the capability to: (1) balance higher speeds without disturbing previously balanced modes, thru the use of modal trial weight sets; (2) balance off-critical speeds; and (3) balance decoupled modes by use of a single balance plane. If no constraints are imposed, this solution form reduces to the general weighted least squares influence coefficient method. A test facility used to examine the use of the general constrained balancing procedure and application of modal trial weight ratios is also described.

Zorzi, E. S.↗
NASA NTRS1986

A computer code for three-dimensional incompressible flows using nonorthogonal body-fitted coordinate systems

In this report, a numerical method for solving the equations of motion of three-dimensional incompressible flows in nonorthogonal body-fitted coordinate (BFC) systems has been developed. The equations of motion are transformed to a generalized curvilinear coordinate system from which the transformed equations are discretized using finite difference approximations in the transformed domain. The hybrid scheme is used to approximate the convection terms in the governing equations. Solutions of the finite difference equations are obtained iteratively by using a pressure-velocity correction algorithm (SIMPLE-C). Numerical examples of two- and three-dimensional, laminar and turbulent flow problems are employed to evaluate the accuracy and efficiency of the present computer code. The user's guide and computer program listing of the present code are also included.

Chen, Y. S.↗
NASA NTRS1986

A numerical method for three-dimensional incompressible flows using nonorthogonal body-fitted coordinate systems

In the present paper, a numerical method for solving the equations of motion of three-dimensional incompressible flows in nonorthogonal body-fitted coordinate (BEC) systems has been developed and evaluated. The equations of motion are transformed to a generalized curvilinear coordinate system from which the transformed equations are discretized using finite difference approximations in the transformed domain. The hybrid scheme and a central differencing plus artificial dissipation scheme are used to approximate the convection terms in the governing equations. Effects of these two schemes on the accuracy of numerical predictions are studied. Solutions of the finite difference equations are obtained iteratively by using a pressure-velocity correction algorithm, SIMPLE-C. Numerical examples of two- and three-dimensional, laminar and turbulent flow problems are employed to evaluate the accuracy and efficiency of the present numerical method.

Chen, Y. S.↗
NASA NTRS1987

A linearized Euler analysis of unsteady flows in turbomachinery

A method for calculating unsteady flows in cascades is presented. The model, which is based on the linearized unsteady Euler equations, accounts for blade loading shock motion, wake motion, and blade geometry. The mean flow through the cascade is determined by solving the full nonlinear Euler equations. Assuming the unsteadiness in the flow is small, then the Euler equations are linearized about the mean flow to obtain a set of linear variable coefficient equations which describe the small amplitude, harmonic motion of the flow. These equations are discretized on a computational grid via a finite volume operator and solved directly subject to an appropriate set of linearized boundary conditions. The steady flow, which is calculated prior to the unsteady flow, is found via a Newton iteration procedure. An important feature of the analysis is the use of shock fitting to model steady and unsteady shocks. Use of the Euler equations with the unsteady Rankine-Hugoniot shock jump conditions correctly models the generation of steady and unsteady entropy and vorticity at shocks. In particular, the low frequency shock displacement is correctly predicted. Results of this method are presented for a variety of test cases. Predicted unsteady transonic flows in channels are compared to full nonlinear Euler solutions obtained using time-accurate, time-marching methods. The agreement between the two methods is excellent for small to moderate levels of flow unsteadiness. The method is also used to predict unsteady flows in cascades due to blade motion (flutter problem) and incoming disturbances (gust response problem).

Hall, Kenneth C.↗
NASA NTRS1987

Re-analysis procedure based on the mixed formulation

An analysis procedure is presented for large-scale structural systems, with large numbers of degrees of freedom and design variables. The procedure uses a mixed formulation with the fundamental unknowns consisting of both stress and displacement parameters. Other elements of the procedure include: (1) lumping the design variables into a single tracing parameter; (2) operator splitting or additive decomposition of different arrays in the finite element equations into the corresponding arrays of the original structure plus correction terms; and (3) application of a reduction method through the use of the finite element method and the classical Bubnov-Galerkin technique. The re-analysis procedure is applied to the linear static and free vibration problems of plate and shell structures.

Noor, Ahmed K.↗
NASA NTRS1989

Calculation of unsteady flows in turbomachinery using the linearized Euler equations

A method for calculating unsteady flows in cascades is presented. The model, which is based on the linearized unsteady Euler equations, accounts for blade loading shock motion, wake motion, and blade geometry. The mean flow through the cascade is determined by solving the full nonlinear Euler equations. Assuming the unsteadiness in the flow is small, then the Euler equations are linearized about the mean flow to obtain a set of linear variable coefficient equations which describe the small amplitude, harmonic motion of the flow. These equations are discretized on a computational grid via a finite volume operator and solved directly subject to an appropriate set of linearized boundary conditions. The steady flow, which is calculated prior to the unsteady flow, is found via a Newton iteration procedure. An important feature of the analysis is the use of shock fitting to model steady and unsteady shocks. Use of the Euler equations with the unsteady Rankine-Hugoniot shock jump conditions correctly models the generation of steady and unsteady entropy and vorticity at shocks. In particular, the low frequency shock displacement is correctly predicted. Results of this method are presented for a variety of test cases. Predicted unsteady transonic flows in channels are compared to full nonlinear Euler solutions obtained using time-accurate, time-marching methods. The agreement between the two methods is excellent for small to moderate levels of flow unsteadiness. The method is also used to predict unsteady flows in cascades due to blade motion (flutter problem) and incoming disturbances (gust response problem).

Hall, Kenneth C.↗
NASA NTRS1989

3-D composite velocity solutions for subsonic/transonic flow over afterbodies

A composite velocity procedure for the three-dimensional reduced Navier-Stokes equations is developed. The velocity components are written as a combined multiplicative and additive composite of viscous like velocities and pseudo-potential or inviscid velocities. The solution procedure is then consistent with both asymptotic inviscid flow and boundary layer theory. For transonic flow cases, the Enquist-Osher flux biasing scheme developed for the full potential equation is used. A quasi-conservation form of the governing equations is used in the shock region to capture the correct rotational behavior. The composite velocity procedure is applied for the solution of three-dimensional afterbody problems.

Gornier, Raymond E.↗
NASA NTRS1990

Flat plate at incidence as a waverider in rarefied hypersonic flow

The physical validity of continuum models and their ability to predict the critical aerothermodynamic properties of a waverider at high altitudes are examined using a flat plate at angle of attack as a generic hypersonic lifting vehicle. For a shock layer far from local translational equilibrium, a theoretical study based on Grad's thirteen-moment equations shows that the Navier-Stokes based solutions can correctly predict the drag, lift, and surface heat transfer rate, with the prediction error comparable to that of the standard shock-layer theory. The conclusion is supported by a comparison with direct simulation Monte Carlo calculations.

Cheng, H. K.↗
NASA NTRS1992

Diffusive shock acceleration - Acceleration rate, magnetic-field direction and the diffusion limit

This paper reviews the concept of diffusive shock acceleration, showing that the acceleration of charged particles at a collisionless shock is a straightforward consequence of the standard cosmic-ray transport equation, provided that one treats the discontinuity at the shock correctly. This is true for arbitrary direction of the upstream magnetic field. Within this framework, it is shown that acceleration at perpendicular or quasi-perpendicular shocks is generally much faster than for parallel shocks. Paradoxically, it follows also that, for a simple scattering law, the acceleration is faster for less scattering or larger mean free path. Obviously, the mean free path can not become too large or the diffusion limit becomes inapplicable. Gradient and curvature drifts caused by the magnetic-field change at the shock play a major role in the acceleration process in most cases. Recent observations of the charge state of the anomalous component are shown to require the faster acceleration at the quasi-perpendicular solar-wind termination shock.

Jokipii, J. R.↗
NASA NTRS1996

An Algebraic Multigrid Solver for Navier-Stokes Problems in the Discrete Second-Order Approximation

An algebraic multigrid scheme is presented for solving the discrete Navier-Stokes equations to second-order accuracy using the defect-correction method. Solutions have been obtained for problems involving both structured and unstructured meshes, with the resolution and resolution grading controlled by global and local mesh refinements. The solver is efficient and robust to the extent that no underrelaxation of variables has been required to ensure convergence, but rates of convergence can be improved with small amounts of underrelaxation of the velocity-pressure coupling. Provided that the computational mesh can resolve the flow field, convergence characteristics are almost mesh independent. Rates of convergence actually improve with refinement, asymptotically approaching mesh independent values. For extremely coarse meshes where dispersive truncation errors would be expected to prevent convergence (or even induce divergence), solutions can still be obtained by using explicit underrelaxation in the iterative cycle.

Webster, R.↗
NASA NTRS1996

Heat Capacity Anomaly Near the Lower Critical Consolute Point of Triethylamine-Water

The heat capacity of the binary liquid mixture triethylamine-water has been measured near its lower critical consolute point using a scanning, adiabatic calorimeter. Two data runs are analyzed to provide heat capacity and enthalpy data that are fitted by equations with background terms and a critical term that includes correction to scaling. The critical exponent a was determined to be 0.107 +/- 0.006, consistent with theoretical predictions. When alpha was fixed at 0.11 to determine various amplitudes consistently, our values of A(+) and A(-) agreed with a previous heat capacity measurement, but the value of A(-) was inconsistent with values determined by density or refractive index measurements. While our value for the amplitude ratio A(+)/ A(-) = 0.56 +/- 0.02 was consistent with other recent experimental determinations in binary liquid mixtures, it was slightly larger than either theoretical predictions or recent experimental values in liquid-vapor systems. The correction to scaling amplitude ratio D(+)/D(-) = 0.5 +/- 0.1 was half of that predicted. As a result of several more precise theoretical calculations and experimental determinations, the two-scale-factor universality ratio X, which we found to be 0.019 +/- 0.003, now is consistent among experiments and theories. A new 'universal' amplitude ratio R(sup +/-)(sub Bcr) involving the amplitudes for the specific heat was tested. Our determination of R(sup +/-)(sub Bcr) = -0.5 +/- 0.1 and R(sup -)(sub Bcr) = 1.1 +/- 0.1 is smaller in magnitude than predicted and is the first such determination in a binary fluid mixture.

Flewelling, Anne C.↗
NASA NTRS2003

Lunar Rotation and the Lunar Interior

Variations in rotation and orientation of the Moon are sensitive to solid-body tidal dissipation, dissipation due to relative motion at the fluid-core/ solid-mantle boundary, and tidal Love number k2. There is weaker sensitivity to flattening of the core-mantle boundary (CMB) and fluid core moment of inertia. Accurate Lunar Laser Ranging (LLR) measurements of the distance from observatories on the Earth to four retroreflector arrays on the Moon are sensitive to lunar rotation and orientation variations and tidal displacements. Past solutions using the LLR data have given results for dissipation due to solid-body tides and fluid core plus Love number. Past detection of CMB flattening has been marginal but is improving, while direct detection of the core moment has not yet been achieved. Three decades of Lunar Laser Ranging (LLR) data are analyzed using a weighted least-squares approach. The lunar solution parameters include dissipation at the fluid-core/solid-mantle boundary, tidal dissipation, dissipation-related coefficients for rotation and orientation terms, potential Love number k2, a correction to the constant term in the tilt of the equator to the ecliptic which is meant to approximate the influence of core-mantle boundary flattening, and displacement Love numbers h2 and l2. Several solutions, with different combinations of solution parameters and constraints, are considered.

Williams, J. G.↗