Search NASA⌕ Search

SEARCH · Search NASA

Results for “EQUATIONS - MOTION”

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
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

A solution method for the unsteady incompressible Navier-Stokes equations in generalized coordinate systems

A solution method based on a fractional step approach is developed for obtaining time-dependent solutions of the three-dimensional, incompressible Navier-Stokes equations in generalized coordinate systems. The governing equations are discretized conservatively by finite volumes using a staggered mesh system. The primitive variable formulation uses the volume fluxes across the faces of each computational cell as dependent variables. This procedure, combined with accurate and consistent approximations of geometric parameters, is done to satisfy the discretized mass conservation equation to machine accuracy as well as to gain favorable convergence properties of the Poisson solver. The discretized equations are second-order-accurate in time and space and no smoothing terms are added. An approximate-factorization scheme is implemented in solving the momentum equations. A novel ZEBRA scheme with four-color ordering is devised for the efficient solution of the Poisson equation. Several two and three-dimensional solutions are compared with other numerical and experimental results to validate the present method.

Rosenfeld, Moshe↗
NASA NTRS1988

On Hilbert-Schmidt norm convergence of Galerkin approximation for operator Riccati equations

An abstract approximation framework for the solution of operator algebraic Riccati equations is developed. The approach taken is based on a formulation of the Riccati equation as an abstract nonlinear operator equation on the space of Hilbert-Schmidt operators. Hilbert-Schmidt norm convergence of solutions to generic finite dimensional Galerkin approximations to the Riccati equation to the solution of the original infinite dimensional problem is argued. The application of the general theory is illustrated via an operator Riccati equation arising in the linear-quadratic design of an optimal feedback control law for a 1-D heat/diffusion equation. Numerical results demonstrating the convergence of the associated Hilbert-Schmidt kernels are included.

Rosen, I. G.↗
NASA NTRS1989

A genuinely multi-dimensional upwind cell-vertex scheme for the Euler equations

The solution of the two-dimensional Euler equations is based on the two-dimensional linear convection equation and the Euler-equation decomposition developed by Hirsch et al. The scheme is genuinely two-dimensional. At each iteration, the data are locally decomposed into four variables, allowing convection in appropriate directions. This is done via a cell-vertex scheme with a downwind-weighted distribution step. The scheme is conservative, and third-order accurate in space. The derivation and stability analysis of the scheme for the convection equation, and the derivation of the extension to the Euler equations are given. Preconditioning techniques based on local values of the convection speeds are discussed. The scheme for the Euler equations is applied to two channel-flow problems. It is shown to converge rapidly to a solution that agrees well with that of a third-order upwind solver.

Powell, Kenneth G.↗
NASA NTRS1989

Unsteady density-current equations for highly curved terrain

New nonlinear partial differential equations containing terrain curvature and its rate of change are derived that describe the flow of an atmospheric density current. Unlike the classical hydraulic-type equations for density currents, the new equations are valid for two-dimensional, gradually varied flow over highly curved terrain, hence suitable for computing unsteady (or steady) flows over arbitrary mountain/valley profiles. The model assumes the atmosphere above the density current exerts a known arbitrary variable pressure upon the unknown interface. Later this is specialized to the varying hydrostatic pressure of the atmosphere above. The new equations yield the variable velocity distribution, the interface position, and the pressure distribution that contains a centrifugal component, often significantly larger than its hydrostatic component. These partial differential equations are hyperbolic, and the characteristic equations and characteristic directions are derived. Using these to form a characteristic mesh, a hypothetical unsteady curved-flow problem is calculated, not based upon observed data, merely as an example to illustrate the simplicity of their application to unsteady flows over mountains.

Sivakumaran, N. S.↗
NASA NTRS1991

On the wall-normal velocity of the compressible boundary-layer equations

Numerical methods for the compressible boundary-layer equations are facilitated by transformation from the physical (x,y) plane to a computational (xi,eta) plane in which the evolution of the flow is 'slow' in the time-like xi direction. The commonly used Levy-Lees transformation results in a computationally well-behaved problem for a wide class of non-similar boundary-layer flows, but it complicates interpretation of the solution in physical space. Specifically, the transformation is inherently nonlinear, and the physical wall-normal velocity is transformed out of the problem and is not readily recovered. In light of recent research which shows mean-flow non-parallelism to significantly influence the stability of high-speed compressible flows, the contribution of the wall-normal velocity in the analysis of stability should not be routinely neglected. Conventional methods extract the wall-normal velocity in physical space from the continuity equation, using finite-difference techniques and interpolation procedures. The present spectrally-accurate method extracts the wall-normal velocity directly from the transformation itself, without interpolation, leaving the continuity equation free as a check on the quality of the solution. The present method for recovering wall-normal velocity, when used in conjunction with a highly-accurate spectral collocation method for solving the compressible boundary-layer equations, results in a discrete solution which is extraordinarily smooth and accurate, and which satisfies the continuity equation nearly to machine precision. These qualities make the method well suited to the computation of the non-parallel mean flows needed by spatial direct numerical simulations (DNS) and parabolized stability equation (PSE) approaches to the analysis of stability.

Pruett, C. David↗
NASA NTRS1991

Whipple shield sizing equations

Wipple shield sizing equations are given for use in current Space Station Freedom (SSF) Work Package 2 (WP-2) trade study activities. These equations are modifications of the 1969 Cour-Palais predictor equations which were used by McDonnell Douglas to size WP-2 meteoroid and debris shielding for the SSF preliminary design review. Recent hypervelocity impact (HVI) test results have shown that the original 1969 Cour-Palais predictor improperly scales to the particle sizes that the WP-2 shields must be designed to protect against. The original equations have been redesigned to correct this scaling deficiency. Substantial increases to WP-2 shielding weight estimates are indicated by the modified equations. Several possibilities exist, however, to reduce the weight of WP-2 shielding. These equations will be updated in the future as warranted by the results of ultra-high speed (greater than 19 km/sec) HVI tests and further analysis.

Christiansen, Eric L.↗
NASA NTRS1993

Comparison of total body water estimates from O-18 and bioelectrical response prediction equations

Identification of an indirect, rapid means to measure total body water (TBW) during space flight may aid in quantifying hydration status and assist in countermeasure development. Bioelectrical response testing and hydrostatic weighing were performed on 27 subjects who ingested O-18, a naturally occurring isotope of oxygen, to measure true TBW. TBW estimates from three bioelectrical response prediction equations and fat-free mass (FFM) were compared to TBW measured from O-18. A repeated measures MANOVA with post-hoc Dunnett's Test indicated a significant (p less than 0.05) difference between TBW estimates from two of the three bioelectrical response prediction equations and O-18. TBW estimates from FFM and the Kushner & Schoeller (1986) equation yielded results that were similar to those given by O-18. Strong correlations existed between each prediction method and O-18; however, standard errors, identified through regression analyses, were higher for the bioelectrical response prediction equations compared to those derived from FFM. These findings suggest (1) the Kushner & Schoeller (1986) equation may provide a valid measure of TBW, (2) other TBW prediction equations need to be identified that have variability similar to that of FFM, and (3) bioelectrical estimates of TBW may prove valuable in quantifying hydration status during space flight.

Barrows, Linda H.↗
NASA NTRS1993

Nonlinear modes of the tensor Dirac equation and CPT violation

Recently, it has been shown that Dirac's bispinor equation can be expressed, in an equivalent tensor form, as a constrained Yang-Mills equation in the limit of an infinitely large coupling constant. It was also shown that the free tensor Dirac equation is a completely integrable Hamiltonian system with Lie algebra type Poisson brackets, from which Fermi quantization can be derived directly without using bispinors. The Yang-Mills equation for a finite coupling constant is investigated. It is shown that the nonlinear Yang-Mills equation has exact plane wave solutions in one-to-one correspondence with the plane wave solutions of Dirac's bispinor equation. The theory of nonlinear dispersive waves is applied to establish the existence of wave packets. The CPT violation of these nonlinear wave packets, which could lead to new observable effects consistent with current experimental bounds, is investigated.

Reifler, Frank J.↗
NASA NTRS1992

A near-wall two-equation model for turbulent heat fluxes

A near-wall two-equation model for turbulent heat fluxes is derived from the temperature variance and its dissipation-rate equations and the assumption of gradient transport. Only incompressible flows with non-buoyant heat transfer are considered. The near-wall asymptotics of each term in the exact equations are examined and used to derive near-wall correction functions that render the modeled equations consistent with these behavior. Thus modeled, the equations are used to calculate fully-developed pipe and channel flows with heat transfer. It is found that the proposed two-equation model yields asymptotically correct near-wall behavior for the normal heat flux, the temperature variance and its near-wall budget and correct limiting wall values for these properties compared to direct simulation data and measurements obtained under different wall boundary conditions.

Sommer, T. P.↗
NASA NTRS1994

Theoretical study of the incompressible Navier-Stokes equations by the least-squares method

Usually the theoretical analysis of the Navier-Stokes equations is conducted via the Galerkin method which leads to difficult saddle-point problems. This paper demonstrates that the least-squares method is a useful alternative tool for the theoretical study of partial differential equations since it leads to minimization problems which can often be treated by an elementary technique. The principal part of the Navier-Stokes equations in the first-order velocity-pressure-vorticity formulation consists of two div-curl systems, so the three-dimensional div-curl system is thoroughly studied at first. By introducing a dummy variable and by using the least-squares method, this paper shows that the div-curl system is properly determined and elliptic, and has a unique solution. The same technique then is employed to prove that the Stokes equations are properly determined and elliptic, and that four boundary conditions on a fixed boundary are required for three-dimensional problems. This paper also shows that under four combinations of non-standard boundary conditions the solution of the Stokes equations is unique. This paper emphasizes the application of the least-squares method and the div-curl method to derive a high-order version of differential equations and additional boundary conditions. In this paper, an elementary method (integration by parts) is used to prove Friedrichs' inequalities related to the div and curl operators which play an essential role in the analysis.

Jiang, Bo-Nan↗
NASA NTRS1994

A modified restricted Euler equation for turbulent flows with mean velocity gradients

The restricted Euler equation captures many important features of the behavior of the velocity gradient tensor observed in direct numerical simulations (DNS) of isotropic turbulence. However, in slightly more complex flows the agreement is not good, especially in regions of low dissipation. In this paper, it is demonstrated that the Reynolds-averaged restricted Euler equation violates the balance of mean momentum for virtually all homogeneous turbulent flows with only two major exceptions: isotropic and homogeneously-sheared turbulence. A new model equation which overcomes this shortcoming and is more widely applicable is suggested. This modele is derived from the Navier-Stokes equation with a restricted Euler type approximation made on the fluctuating velocity gradient field. Analytical solutions of the proposed modified restricted Euler equation appear to be difficult to obtain. Hence, a strategy for numerically calculating the velocity gradient tensor is developed. Preliminary calculations tend to indicate that the modified restricted Euler equation captures many important aspects of the behavior of the fluctuating velocity gradients in anisotropic homogeneous turbulence.

Girimaji, Sharath S.↗
NASA NTRS1994

Radiation boundary condition and anisotropy correction for finite difference solutions of the Helmholtz equation

In this paper finite-difference solutions of the Helmholtz equation in an open domain are considered. By using a second-order central difference scheme and the Bayliss-Turkel radiation boundary condition, reasonably accurate solutions can be obtained when the number of grid points per acoustic wavelength used is large. However, when a smaller number of grid points per wavelength is used excessive reflections occur which tend to overwhelm the computed solutions. Excessive reflections are due to the incompability between the governing finite difference equation and the Bayliss-Turkel radiation boundary condition. The Bayliss-Turkel radiation boundary condition was developed from the asymptotic solution of the partial differential equation. To obtain compatibility, the radiation boundary condition should be constructed from the asymptotic solution of the finite difference equation instead. Examples are provided using the improved radiation boundary condition based on the asymptotic solution of the governing finite difference equation. The computed results are free of reflections even when only five grid points per wavelength are used. The improved radiation boundary condition has also been tested for problems with complex acoustic sources and sources embedded in a uniform mean flow. The present method of developing a radiation boundary condition is also applicable to higher order finite difference schemes. In all these cases no reflected waves could be detected. The use of finite difference approximation inevita bly introduces anisotropy into the governing field equation. The effect of anisotropy is to distort the directional distribution of the amplitude and phase of the computed solution. It can be quite large when the number of grid points per wavelength used in the computation is small. A way to correct this effect is proposed. The correction factor developed from the asymptotic solutions is source independent and, hence, can be determined once and for all. The effectiveness of the correction factor in providing improvements to the computed solution is demonstrated in this paper.

Tam, Christopher K. W.↗
NASA NTRS1995

The limitation and applicability of Musher-Sturman equation to two dimensional lower hybrid wave collapse

The existence of localized regions of intense lower hybrid waves in the auroral ionosphere recently observed by rocket and satellite experiments can be understood by the study of a non-linear two-timescale coupling process. In this Letter, we demonstrate that the leading non-linear term in the standard Musher-Sturman equation vanishes identically in strict two-dimensions (normal to the magnetic field). Instead, the new two-dimensional equation is characterized by a much weaker non-linear term which arises from the ponderomotive force perpendicular to the magnetic field, particularly that due to the ions. The old and new equations are compared by means of time-evolution calculations of wave fields. The results exhibit a remarkable difference in the evolution of the waves as governed by the two equations. Such dissimilar outcomes motivate our investigation of the limitation of Musher-Sturman equation in quasi-two-dimensions. Only within all these limits can Musher-Sturman equation adequately describe the collapse of lower hybrid waves.

Tam, Sunny W. Y.↗
NASA NTRS1995

On absorbing boundary conditions for linearized Euler equations by a perfectly matched layer

Recently, Berenger introduced a Perfectly Matched Layer (PML) technique for absorbing electromagnetic waves. In the present paper, a perfectly matched layer is proposed for absorbing out-going two-dimensional waves in a uniform mean flow, generated by linearized Euler equations. It is well known that the linearized Euler equations support acoustic waves, which travel with the speed of sound relative to the mean flow, and vorticity and entropy waves, which travel with the mean flow. The PML equations to be used at a region adjacent to the artificial boundary for absorbing these linear waves are defined. Plane waves solutions to the PML equations are developed and wave propagation and absorption properties are given. It is shown that the theoretical reflection coefficients at an interface between the Euler and PML domains are zero, independent of the angle of incidence and frequency of the waves. As such, the present study points out a possible alternative approach for absorbing out-going waves of the Euler equations with little or no reflection in computation. Numerical examples that demonstrate the validity of the proposed PML equations are also presented.

Hu, Fang Q.↗
NASA NTRS1998

Fresnel Lens Solar Concentrator Design Based on Geometric Optics and Blackbody Radiation Equations

Fresnel lenses have been used for years as solar concentrators in a variety of applications. Several variables effect the final design of these lenses including: lens diameter, image spot distance from the lens, and bandwidth focused in the image spot. Defining the image spot as the geometrical optics circle of least confusion, a set of design equations has been derived to define the groove angles for each groove on the lens. These equations allow the distribution of light by wavelength within the image spot to be calculated. Combining these equations with the blackbody radiation equations, energy distribution, power, and flux within the image spot can be calculated. In addition, equations have been derived to design a lens to produce maximum flux in a given spot size. Using these equations, a lens may be designed to optimize the spot energy concentration for given energy source.

Watson, Michael D.↗
NASA NTRS1994

Progress in Computational Aeroelasticity Using High Fidelity Flow and Structural Equations on Parallel Computers

Aeroelasticity which involves strong coupling of fluids, structures and controls is an important element in designing an aircraft. Computational aeroelasticity using low fidelity methods such as the linear aerodynamic flow equations coupled with the modal structural equations are well advanced. Though these low fidelity approaches are computationally less intensive, they are not adequate for the analysis of modern aircraft such as High Speed Civil Transport (HSCT) and Advanced Subsonic Transport (AST) which can experience complex flow/structure interactions. HSCT can experience vortex induced aeroelastic oscillations whereas AST can experience transonic buffet associated structural oscillations. Both aircraft may experience a dip in the flutter speed at the transonic regime. For accurate aeroelastic computations at these complex fluid/structure interaction situations, high fidelity equations such as the Navier-Stokes for fluids and the finite-elements for structures are needed. Computations using these high fidelity equations require large computational resources both in memory and speed. Current conventional supercomputers have reached their limitations both in memory and speed. As a result, parallel computers have evolved to overcome the limitations of conventional computers. This paper will address the transition that is taking place in computational aeroelasticity from conventional computers to parallel computers. The paper will address special techniques needed to take advantage of the architecture of new parallel computers. Results will be illustrated from computations made on iPSC/860 and IBM SP2 computer by using ENASERO code that directly couples the Euler/Navier-Stokes flow equations with high resolution finite-element structural equations.

Guruswamy, Guru P.↗
NASA NTRS1994

On the Connection Between One-and Two-Equation Models of Turbulence

A formalism will be presented that allows the transformation of two-equation eddy viscosity turbulence models into one-equation models. The transformation is based on an assumption that is widely accepted over a large range of boundary layer flows and that has been shown to actually improve predictions when incorporated into two-equation models of turbulence. Based on that assumption, a new one-equation turbulence model will be derived. The new model will be tested in great detail against a previously introduced one-equation model and against its parent two-equation model.

Menter, F. R.↗
NASA NTRS1994

Fluid/Structure Interaction Studies of Aircraft Using High Fidelity Equations on Parallel Computers

Abstract Aeroelasticity which involves strong coupling of fluids, structures and controls is an important element in designing an aircraft. Computational aeroelasticity using low fidelity methods such as the linear aerodynamic flow equations coupled with the modal structural equations are well advanced. Though these low fidelity approaches are computationally less intensive, they are not adequate for the analysis of modern aircraft such as High Speed Civil Transport (HSCT) and Advanced Subsonic Transport (AST) which can experience complex flow/structure interactions. HSCT can experience vortex induced aeroelastic oscillations whereas AST can experience transonic buffet associated structural oscillations. Both aircraft may experience a dip in the flutter speed at the transonic regime. For accurate aeroelastic computations at these complex fluid/structure interaction situations, high fidelity equations such as the Navier-Stokes for fluids and the finite-elements for structures are needed. Computations using these high fidelity equations require large computational resources both in memory and speed. Current conventional super computers have reached their limitations both in memory and speed. As a result, parallel computers have evolved to overcome the limitations of conventional computers. This paper will address the transition that is taking place in computational aeroelasticity from conventional computers to parallel computers. The paper will address special techniques needed to take advantage of the architecture of new parallel computers. Results will be illustrated from computations made on iPSC/860 and IBM SP2 computer by using ENSAERO code that directly couples the Euler/Navier-Stokes flow equations with high resolution finite-element structural equations.

Guruswamy, Guru↗