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 415 records · Page 23

Ocean color algorithm for remote sensing of chlorophyll

An algorithm for the remote detection of chlorophyll a in the ocean was tested during a Nantucket Shoals experiment conducted by NASA. A set of Multichannel Ocean Color Sensor (MOCS) data selected from one flight for each of the two altitudes flown was used to calibrate the algorithm for chlorophyll a concentration. The equations were then applied to all unsaturated MOCS data collected during the 8-day experiment to generate contour maps of chlorophyll a concentration over the shoals. One additional flight was conducted away from the shoals both on and off the Continental Shelf. Although no solar elevation or environmental corrections were made to the original conversions, the equations in these tests determined chlorophyll a concentrations to an accuracy better than 1.0 micron g/L despite the fact that the solar elevation varied between 20 deg and 56 deg during the data collection periods of the experiments.

Grew, G. W.↗

A two-equation turbulence model for compressible reacting flows

A two-equation turbulence model (k-epsilon) has been modified to model compressible reacting flows. A compressibility correction model based on modeling the dilatational terms in the Reynolds stress equations has been used. Computations are carried out using the Navier-Stokes solver SPARK which includes a finite rate chemistry model for hydrogen-air combustion. Initially, planar compressible mixing layers are simulated. Comparisons of the predictions with available experimental data and the predictions of a compressible Reynolds stress closure indicate that the model is well suited for the study of such flows. The decrease in the growth rate of the mixing layer with increasing convective Mach number is well predicted. The model is then tested on a three-dimensional supersonic flow problem where solid boundaries are present. An axisymmetric nonreacting mixing layer flow involving dissimilar gases is evaluated next using the model and the predictions are compared with available experimental data. Finally, two reacting turbulent flow configurations have been studied and comparison of predictions and experimental data has been done. The comparisons indicate that the model predicts the overall features of the test cases very well.

Narayan, J. R.↗

A finite difference scheme for three-dimensional steady laminar incompressible flow

A finite difference scheme for three-dimensional steady laminar incompressible flows is presented. The Navier-Stokes equations are expressed conservatively in terms of velocity and pressure increments (delta form). First order upwind differences are used for first order partial derivatives of velocity increments resulting in a diagonally dominant matrix system. Central differences are applied to all other terms for second order accuracy. The SIMPLE pressure correction algorithm is used to satisfy the continuity equation. Numerical results are presented for cubic cavity flow problems for Reynolds numbers up to 2000 and are in good agreement with other numerical results.

Hwang, Danny P.↗

A finite difference scheme for three-dimensional steady laminar incompressible flow

A finite difference scheme for three-dimensional steady laminar incompressible flows is presented. The Navier-Stokes equations are expressed conservatively in terms of velocity and pressure increments (delta form). First order upwind differences are used for first order partial derivatives of velocity increments resulting in a diagonally dominant matrix system. Central differences are applied to all other terms for second order accuracy. The SIMPLE pressure correction algorithm is used to satisfy the continuity equation. Numerical results are presented for cubic cavity flow problems for Reynolds numbers up to 2000 and are in good agreement with other numerical results.

Hwang, Danny P.↗

Investigation of computational and spectral analysis methods for aeroacoustic wave propagation

Most computational fluid dynamics (CFD) schemes are not adequately accurate for solving aeroacoustics problems, which have wave amplitudes several orders of magnitude smaller yet with frequencies larger than the flow field variations generating the sound. Hence, a computational aeroacoustics (CAA) algorithm should have minimal dispersion and dissipation features. A dispersion relation preserving (DRP) scheme is, therefore, applied to solve the linearized Euler equations in order to simulate the propagation of three types of waves, namely: acoustic, vorticity, and entropy waves. The scheme is derived using an optimization procedure to ensure that the numerical derivatives preserve the wave number and angular frequency of the partial differential equations being discretized. Consequently, simulated waves propagate with the correct wave speeds and exhibit their appropriate properties. A set of radiation and outflow boundary conditions, compatible with the DRP scheme and derived from the asymptotic solutions of the governing equations, are also implemented. Numerical simulations are performed to test the effectiveness of the DRP scheme and its boundary conditions. The computed solutions are shown to agree favorably with the exact solutions. The major restriction appears to be that the dispersion relations can be preserved only for waves with wave lengths longer than four or five spacings. The boundary conditions are found to be transparent to the outgoing disturbances. However, when the disturbance source is placed closer to a boundary, small acoustic reflections start appearing. CAA generates enormous amounts of temporal data which needs to be reduced to understand the physical problem being simulated. Spectral analysis is one approach that helps us in extracting information which often can not be easily interpreted in the time domain. Thus, three different methods for the spectral analysis of numerically generated aeroacoustic data are studied. First, the capabilities of two traditional methods for spectral analysis, namely, the Blackman-Tukey method and periodogram method, are compared in estimating the spectra of a simple-periodic process. The periodogram is then applied to analyze transitory-deterministic processes. Finally, these two methods are compared with a more recent method, referred as the Weighted-Overlapped-Segment-Averaging (WOSA) method, in estimating the spectra of a chaotic (random-like) process. From the demonstrative case for the spectral analyses of data generated by simple-periodic process, the periodogram method is found to give a better estimate of the steep-sloped spectra than the Blackman-Tukey method. Also, for this problem, the Hanning window is found to perform better with the periodogram method than with the Blackman-Tukey method. Finally, for the spectral analysis of data generated by the chaotic process, the periodogram method does not perform well, whereas, the WOSA and Blackman-Tukey methods give equivalently good results.

Vanel, Florence O.↗

Solving Upwind-Biased Discretizations: Defect-Correction Iterations

This paper considers defect-correction solvers for a second order upwind-biased discretization of the 2D convection equation. The following important features are reported: (1) The asymptotic convergence rate is about 0.5 per defect-correction iteration. (2) If the operators involved in defect-correction iterations have different approximation order, then the initial convergence rates may be very slow. The number of iterations required to get into the asymptotic convergence regime might grow on fine grids as a negative power of h. In the case of a second order target operator and a first order driver operator, this number of iterations is roughly proportional to h-1/3. (3) If both the operators have the second approximation order, the defect-correction solver demonstrates the asymptotic convergence rate after three iterations at most. The same three iterations are required to converge algebraic error below the truncation error level. A novel comprehensive half-space Fourier mode analysis (which, by the way, can take into account the influence of discretized outflow boundary conditions as well) for the defect-correction method is developed. This analysis explains many phenomena observed in solving non-elliptic equations and provides a close prediction of the actual solution behavior. It predicts the convergence rate for each iteration and the asymptotic convergence rate. As a result of this analysis, a new very efficient adaptive multigrid algorithm solving the discrete problem to within a given accuracy is proposed. Numerical simulations confirm the accuracy of the analysis and the efficiency of the proposed algorithm. The results of the numerical tests are reported.

Diskin, Boris↗

On uniformly valid high-frequency far-field asymptotic solutions of the Helmholtz equation

An asymptotic, large wave number approximation for the Helmholtz equation is derived. The theory is an extension of the geometric acoustic theory, and provides corrections to that theory in the form of multiplicative functions which satisfy parabolic equations. A simple example is used both to illustrate failure of the geometric theory for large propagation distances, and to show the improvement obtained by use of the new theory.

Mcaninch, G. L.↗

Equation-of-Motion Coupled-Cluster Cumulant Green’s Function for Excited States and X-ray Spectra

Green’s function methods within many-body perturbation theory provide a general framework for treating electronic correlations in excited states and spectra. Conventional methods using the Dyson equation or the cumulant expansion are typically based on the GW self-energy approximation. In order to extend this approximation in molecular systems, a non-perturbative real-time coupled-cluster cumulant Green’s function approach has been introduced, where the cumulant is obtained as the solution to a set of coupled first order, non-linear differential equations. This approach naturally includes non-linear corrections to conventional cumulant Green’s function techniques where the cumulant is linear in the GW self-energy. The method yields the spectral function associated with the core Green’s function, which is directly related to the x-ray photoemission spectra (XPS) of molecular systems. The approach also yields very good results for binding energies and satellite excitations. The x-ray absorption spectrum (XAS) is then obtained as a convolution of the core spectral function and an effective one-body absorption spectrum. Here this approach is extended to include the full coupled-cluster-singles (CCS) core Green’s function by including the complete form of the non-linear contributions to the cumulant as well as all single, double, and triple cluster excitations in the CC amplitude equations. This approach naturally builds in orthogonality and shake-up effects analogous to those in the Mahan-Noizeres-de Dominicis edge singularity corrections that enhance the XAS near the edge. Themethod is illustrated for the XPS and XAS of NH 3 .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

D 3 -creatine dilution for skeletal muscle mass measurement: historical development and current status

The French chemist Michel Eugène Chevreul discovered creatine in meat two centuries ago. Extensive biochemical and physiological studies of this organic molecule followed with confirmation that creatine is found within the cytoplasm and mitochondria of human skeletal muscles. Two groups of investigators exploited these relationships five decades ago by first estimating the creatine pool size in vivo with 14 C and 15 N labelled isotopes. Skeletal muscle mass (kg) was then calculated by dividing the creatine pool size (g) by muscle creatine concentration (g/kg) measured on a single muscle biopsy or estimated from the literature. This approach for quantifying skeletal muscle mass is generating renewed interest with the recent introduction of a practical stable isotope (creatine-(methyl-d 3 )) dilution method for estimating the creatine pool size across the full human lifespan. The need for a muscle biopsy has been eliminated by assuming a constant value for whole-body skeletal muscle creatine concentration of 4.3 g/kg wet weight. The current single compartment model of estimating creatine pool size and skeletal muscle mass rests on four main assumptions: tracer absorption is complete; tracer is all retained; tracer is distributed solely in skeletal muscle; and skeletal muscle creatine concentration is known and constant. Three of these assumptions are false to varying degrees. Not all tracer is retained with urinary isotope losses ranging from 0% to 9%; an empirical equation requiring further validation is used to correct for spillage. Not all tracer is distributed in skeletal muscle with non-muscle creatine sources ranging from 2% to 10% with a definitive value lacking. Lastly, skeletal muscle creatine concentration is not constant and varies between muscles (e.g. 3.89–4.62 g/kg), with diets (e.g. vegetarian and omnivore), across age groups (e.g. middle-age, ~4.5 g/kg; old-age, 4.0 g/kg), activity levels (e.g. athletes, ~5 g/kg) and in disease states (e.g. muscular dystrophies, <3 g/kg). Some of the variability in skeletal muscle creatine concentrations can be attributed to heterogeneity in the proportions of wet skeletal muscle as myofibres, connective tissues, and fat. These observations raise serious concerns regarding the accuracy of the deuterated-creatine dilution method for estimating total body skeletal muscle mass as now defined by cadaver analyses of whole wet tissues and in vivo approaches such as magnetic resonance imaging. A new framework is needed in thinking about how this potentially valuable method for measuring the creatine pool size in vivo can be used in the future to study skeletal muscle biology in health and disease.

body composition↗

Code-verification techniques for the method-of-moments implementation of the combined-field integral equation

Code verification plays an important role in establishing the credibility of computational simulations by assessing the correctness of the implementation of the underlying numerical methods. In computational electromagnetics, the numerical solution to integral equations incurs multiple interacting sources of numerical error, as well as other challenges, which render traditional code-verification approaches ineffective. In this paper, we provide approaches to separately measure the numerical errors arising from these different error sources for the method-of-moments implementation of the combined-field integral equation. Here, we demonstrate the effectiveness of these approaches for cases with and without coding errors.

97 MATHEMATICS AND COMPUTING↗

Particle-number fluctuations near the critical point of nuclear matter

Equation of state with the quantum statistics corrections is used for particle-number fluctuations ω of the isotopically symmetric nuclear matter with interparticle van der Waals and Skyrme local density interactions. Here, the fluctuations, ω ∝ 1/K, are analytically derived through the isothermal incompressibility K at first order over a small quantum-statistics parameter. Our approximate analytical results appear to be in good agreement with the results of accurate numerical calculations. These results are also close to those obtained by using more accurate Tolman and Rowlinson expansions of the incompressibility K near the critical point. More general formula for fluctuations ω, improved at the critical point, was obtained for a finite particle-number average $\langle$N$\rangle$ by neglecting, for simplicity, small quantum statistics effects. It is shown that for a large dimensionless parameter, α ∝ K 2 $\langle$N$\rangle$/K", where K" is the second derivative of the incompressibility K as function of the average particle density n, far from the critical point (α >> 1), one finds the traditional asymptote, ω ∝ 1/K, for the fluctuations ω. For a small parameter, α << 1, near the critical point, where K = 0 and α = 0, one obtains another asymptote of ω. These fluctuations, having a maximum near the critical point as function of the average density n, for finite values of $\langle$N$\rangle$ are finite and relatively small, in contrast to the results of the traditional calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Nonperturbative resolution of strong coupling singularities in 4D 𝒩 = 1 heterotic M-theory

We investigate the interior of the moduli space of four-dimensional 𝒩 = 1 theories of gravity arising from compactifications of the 𝐸 8 × 𝐸 8 heterotic string on Calabi-Yau threefolds. By studying the threshold corrections to the coupling of the heterotic gauge groups, we infer the existence of a strong coupling singularity for one of the perturbative heterotic gauge groups, which effectively yields an additional finite distance boundary of the classical scalar field space. In heterotic M-theory, this boundary maps to a domain wall solution for which the gauge coupling and the warp factor on one of the Horava-Witten 9-branes diverge, thus highlighting the gravitational origin of the classical strong coupling singularity. The divergence of the warp factor is, however, regulated once nonperturbative effects are taken into account, as we demonstrate by studying the instanton corrections to the 5D BPS domain wall equations. This regularization implies that the classical strong coupling boundary of the scalar field 4D 𝒩 = 1 heterotic M-theory is resolved, indicating that, at the quantum level, the field space can be extended beyond this classical boundary.

Cvetič, Mirjam↗

Some properties of unstable particle gases

Elastic scattering formation and decay of unstable particle gases, discussing density and corrections to chemical and nuclear statistical equilibrium equations

Delano, M. D.↗

Fokker-Planck equations for charged-particle transport with a discussion of some higher-order effects

A derivation of the Fokker-Planck equation, based on the central limit theorem, is presented which clearly illustrates the conditions for its validity. It is reiterated that previous use of the Fokker-Planck equation in cosmic-ray transport is correct. Higher-order effects associated with magnetic mirroring and field line random walk at low energies are discussed heuristically.

Jokipii, J. R.↗

Error in temperature measurements due to conduction along the sensors

When a sensor is embedded in a solid body to measure its internal temperature, any conduction to or from its sensing element may cause the indicated temperature to be different from the true temperature. This paper describes an analysis of the error caused by conduction when there is an arbitrary temperature distribution in the body along the sensor. The sensor is modeled as a cylindrical fin and the appropriate conduction equation is solved. The solution gives a correction for the error which depends on such parameters as, depth of immersion, thermocouple wire and insulation properties, contact between the sensor and the body, and temperature distribution in the body. The latter may not be known, but the measured temperature distribution can be used as a first approximation. The corrected value can then be used to obtain a better estimate of the error. The results show good agreement with experimental observations.

Singh, B. S.↗

MAGSAT satellite magnetic anomaly map over South America

A scalar magnetic anomaly map was prepared for South America and adjacent marine areas directly from original MAGSAT orbits. Special problems associated with the separation of external field and crustal anomalies, and the reduction of data to a common altitude are addressed. External fields are manifested in a long-wavelength ring current effect, a medium-wavelength equatorial electrojet, and short-wavelength noise. The noise is reduced by selecting profiles from quiet periods (Kp or = 3), and the effect of the electrojet is minimized by drawing the data set from dawn profiles only. The ring current is corrected through the use of a standard equation, augmented by further digital band-pass filtering. Profiles thus filtered differ primarily in amplitude due solely to satellite altitude differences. These differences are normalized by an inversion of the profile data onto a grid of equivalent point dipoles, and recalculated at an altitude of 350 km. The low altitudes in the study area cause instability in the inversion, necessitating separate inversions of several sub-areas which are subsequently merged. Crustal anomalies reduced-to-the-pole exhibit marked correlations to known tectonic features.

Ridgway, J. R.↗

An entropy correction method for unsteady full potential flows with strong shocks

An entropy correction method for the unsteady full potential equation is presented. The unsteady potential equation is modified to account for entropy jumps across shock waves. The conservative form of the modified equation is solved in generalized coordinates using an implicit, approximate factorization method. A flux-biasing differencing method, which generates the proper amounts of artificial viscosity in supersonic regions, is used to discretize the flow equations in space. Comparisons between the present method and solutions of the Euler equations and between the present method and experimental data are presented. The comparisons show that the present method more accurately models solutions of the Euler equations and experiment than does the isentropic potential formulation.

Whitlow, W., Jr.↗

An entropy correction method for unsteady full potential flows with strong shocks

An entropy correction method for the unsteady full potential equation is presented. The unsteady potential equation is modified to account for entropy jumps across shock waves. The conservative form of the modified equation is solved in generalized coordinates using an implicit, approximate factorization method. A flux-biasing differencing method, which generates the proper amounts of artificial viscosity in supersonic regions, is used to discretize the flow equations in space. Comparisons between the present method and solutions of the Euler equations and between the present method and experimental data are presented. The comparisons show that the present method more accurately models solutions of the Euler equations and experiment than does the isentropic potential formulation.

Whitlow, W., Jr.↗