Search NASA⌕ Search

SEARCH · Search NASA

Results for “Continuity 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 541 records · Page 30

Combustion-acoustic stability analysis for premixed gas turbine combustors

Lean, prevaporized, premixed combustors are susceptible to combustion-acoustic instabilities. A model was developed to predict eigenvalues of axial modes for combustion-acoustic interactions in a premixed combustor. This work extends previous work by including variable area and detailed chemical kinetics mechanisms, using the code LSENS. Thus the acoustic equations could be integrated through the flame zone. Linear perturbations were made of the continuity, momentum, energy, chemical species, and state equations. The qualitative accuracy of our approach was checked by examining its predictions for various unsteady heat release rate models. Perturbations in fuel flow rate are currently being added to the model.

Darling, Douglas↗

Lumped Parameter Models of the Central Nervous System for VIIP Research

INTRODUCTION: Current long-duration missions to the International Space Station and future exploration-class missions beyond low-Earth orbit, such as to Mars and asteroids, expose astronauts to increased risk of Visual Impairment and Intracranial Pressure (VIIP) syndrome [1]. It has been hypothesized that the headward shift of cerebral spinal fluid (CSF) and blood in microgravity may cause significant elevation of intracranial pressure (ICP), which in turn induces VIIP syndrome through biomechanical pathways [1, 2]. However, there is insufficient evidence to confirm this hypothesis. In this light, we are developing lumped-parameter models of fluid transport in the central nervous system (CNS) as a means to simulate the influence of microgravity on ICP. The CNS models will also be used in concert with the lumped parameter and finite element models of the eye described in the realted IWS abstracts submitted by Nelson et al., Feola et al. and Ethier et al. METHODS: We have developed a nine compartment CNS model (Figure 1) capable of both time-dependent and steady state fluid transport simulations, based on the works of Stevens et al. [3]. The breakdown of compartments within the model includes: vascular (3), CSF (2), brain (1) and extracranial (3). The boundary pressure in the Central Arteries [A] node is prescribed using an oscillating pressure function PA(t) simulating the carotid pulsatile pressure wave as developed by Linninger et al. [4]. For each time step, pressures are integrated through time using an adaptive-timestep 4th and 5th order Runga-Kutta solver. Once pressures are found, constitutive equations are used to solve for flowrates (Q) between each compartment. In addition to fluid flow between the different compartments, compliance (C) interactions between neighboring compartments are represented. We are also developing a second CNS model based on the works of Linninger et al. [4] which takes a more granular approach to represent the interactions of the intracranial and spinal compartments with the inclusion of arteries, arterioles, capillaries, venules, veins, venous sinus, and ventricles. The flow through the arteries, veins and CSF compartments are governed by continuity, momentum and distensibility balance equations. Furthermore, unlike the Stevens et al. approach, the Monro-Kellie doctrine of constant cranial volume and the bi-phasic nature of the brain parenchyma are implemented. These features appear to be more consistent with the physiologic and anatomical behavior of the CNS, and follow a modeling philosophy similar to the lumped parameter eye model that is intended to be integrated with the CNS model. However, Linninger’s approach has never been implemented to include hydrostatic gradient and microgravity simulation capabilities. Therefore, we aim at implement this modeling approach for spaceflight simulations and assess its overall applicability to VIIP research. OBJECTIVES: We will present verification and validation test results for both models, as well as head-to-head comparison to explore their strengths and limitations with respect to mathematical implementation and physiological significance for VIIP research. In doing so, we hope to provide some guidance to the HRP research community on how to appropriately leverage lumped parameter models for space biomedical research.

Vera, J.↗

Full Multigrid Flow Solver

FMG3D (full multigrid 3 dimensions) is a pilot computer program that solves equations of fluid flow using a finite difference representation on a structured grid. Infrastructure exists for three dimensions but the current implementation treats only two dimensions. Written in Fortran 90, FMG3D takes advantage of the recursive subroutine feature, dynamic memory allocation, and structured-programming constructs of that language. FMG3D supports multi-block grids with three types of block-to-block interfaces: periodic, C-zero, and C-infinity. For all three types, grid points must match at interfaces. For periodic and C-infinity types, derivatives of grid metrics must be continuous at interfaces. The available equation sets are as follows: scalar elliptic equations, scalar convection equations, and the pressure-Poisson formulation of the Navier-Stokes equations for an incompressible fluid. All the equation sets are implemented with nonzero forcing functions to enable the use of user-specified solutions to assist in verification and validation. The equations are solved with a full multigrid scheme using a full approximation scheme to converge the solution on each succeeding grid level. Restriction to the next coarser mesh uses direct injection for variables and full weighting for residual quantities; prolongation of the coarse grid correction from the coarse mesh to the fine mesh uses bilinear interpolation; and prolongation of the coarse grid solution uses bicubic interpolation.

Mineck, Raymond E.↗

Burgers approximation for two-dimensional flow past an ellipse

A linearization of the Navier-Stokes equation due to Burgers in which vorticity is transported by the velocity field corresponding to continuous potential flow is examined. The governing equations are solved exactly for the two dimensional steady flow past an ellipse of arbitrary aspect ratio. The requirement of no slip along the surface of the ellipse results in an infinite algebraic system of linear equations for coefficients appearing in the solution. The system is truncated at a point which gives reliable results for Reynolds numbers R in the range 0 R 5. Predictions of the Burgers approximation regarding separation, drag and boundary layer behavior are investigated. In particular, Burgers linearization gives drag coefficients which are closer to observed experimental values than those obtained from Oseen's approximation. In the special case of flow past a circular cylinder, Burgers approximation predicts a boundary layer whose thickness is roughly proportional to R-1/2.

Dorrepaal, J. M.↗

Maneuver simulations of flexible spacecraft by solving TPBVP

The optimal control of large angle rapid maneuvers and vibrations of a Shuttle mast reflector system is considered. The nonlinear equations of motion are formulated by using Lagrange's formula, with the mast modeled as a continuous beam. The nonlinear terms in the equations come from the coupling between the angular velocities, the modal coordinates, and the modal rates. Pontryagin's Maximum Principle is applied to the slewing problem, to derive the necessary conditions for the optimal controls, which are bounded by given saturation levels. The resulting two point boundary value problem (TPBVP) is then solved by using the quasilinearization algorithm and the method of particular solutions. In the numerical simulations, the structural parameters and the control limits from the Spacecraft Control Lab Experiment (SCOLE) are used. In the 2-D case, only the motion in the plane of an Earth orbit or the single axis slewing motion is discussed. In the 3-D slewing, the mast is modeled as a continuous beam subjected to 3-D deformations. The numerical results for both the linearized system and the nonlinear system are presented to compare the differences in their time response.

Bainum, Peter M.↗

Local Discontinuous Galerkin Methods for Partial Differential Equations with Higher Order Derivatives

In this paper we review the existing and develop new continuous Galerkin methods for solving time dependent partial differential equations with higher order derivatives in one and multiple space dimensions. We review local discontinuous Galerkin methods for convection diffusion equations involving second derivatives and for KdV type equations involving third derivatives. We then develop new local discontinuous Galerkin methods for the time dependent bi-harmonic type equations involving fourth derivatives, and partial differential equations involving fifth derivatives. For these new methods we present correct interface numerical fluxes and prove L(exp 2) stability for general nonlinear problems. Preliminary numerical examples are shown to illustrate these methods. Finally, we present new results on a post-processing technique, originally designed for methods with good negative-order error estimates, on the local discontinuous Galerkin methods applied to equations with higher derivatives. Numerical experiments show that this technique works as well for the new higher derivative cases, in effectively doubling the rate of convergence with negligible additional computational cost, for linear as well as some nonlinear problems, with a local uniform mesh.

Yan, Jue↗

On the continuous dependence with respect to sampling of the linear quadratic regulator problem for distributed parameter systems

The convergence of solutions to the discrete or sampled time linear quadratic regulator problem and associated Riccati equation for infinite dimensional systems to the solutions to the corresponding continuous time problem and equation, as the length of the sampling interval (the sampling rate) tends toward zero (infinity) is established. Both the finite and infinite time horizon problems are studied. In the finite time horizon case, strong continuity of the operators which define the control system and performance index together with a stability and consistency condition on the sampling scheme are required. For the infinite time horizon problem, in addition, the sampled systems must be stabilizable and detectable, uniformly with respect to the sampling rate. Classes of systems for which this condition can be verified are discussed. Results of numerical studies involving the control of a heat/diffusion equation, a hereditary of delay system, and a flexible beam are presented and discussed.

Rosen, I. G.↗

Turbulence models from the point of view of an industrial user

Since there appears to be little prospect of the development of methods for the solution of the full, time dependent Navier-Stokes equations in the near future, the need continues for turbulence models to approximate the Reynolds stress terms that appear in the time averaged Navier-Stokes equations. Even if methods were available for solving the full equations, this would not necessarily be the optimum choice in all cases. As the cost of numerical computations decreases, the trend toward the use of more complex methods is likely to continue, but there will always be a need for a range of methods, depending on the accuracy and detail required from the calculation.

Birch, S. F.↗

Theory of Electrostatic Fields in the Ionosphere at Equatorial Latitudes

The properties of the elongated electrostatic fields that are required to provide the coupling mechanism in the dynamo-motor concept of the E and F regions of the ionosphere are examined theoretically for the conditions that prevail in equatorial latitudes. The analysis is developed for an electrostatic field of arbitrary horizontal scale in a horizontally stratified partly ionized gas subject to an imposed magnetic field having the form of a parabolic arch over the equator. The anisotropic character and continuous variation with height of the conductivity are retained throughout , and numerical solutions are determined for the attenuation of the electric field with distance along the field line. The results are similar qualitatively to those found previously upon analysis of the corresponding problem for middle latitudes, but the attenuation of the electrostatic field with height is considerably greater. It develops, in particular, that the coupling between E and F regions is very small for fields having horizontal wavelengths of a few kilometers. It follows that the dynamo-motor concept could not be used to account for the presence of irregularities of this scale. On the other hand, it is found that almost all the attenuation occurs at heights near that of the E region. If an electrostatic field having a horizontal wavelength of a few kilometers could be produced at heights of, say 200 km or greater, it follows that it would be very elongated and extend with little change in amplitude from hemisphere to hemisphere. The results also indicate that effective coupling could be achieved for fields having horizontal wavelengths of the order of tens of kilometers, but the assumptions introduced to simplify the analysis may impair the quantitative reliability of the results for fields of this scale.

Briggs, Benjamin R.↗

FORTRAN 4 computer program for calculating critical speeds of rotating shafts

A FORTRAN 4 computer program, written for the IBM DCS 7094/7044 computer, that calculates the critical speeds of rotating shafts is described. The shaft may include bearings, couplings, extra masses (nonshaft mass), and disks for the gyroscopic effect. Shear deflection is also taken into account, and provision is made in the program for sections of the shaft that are tapered. The boundary conditions at the ends of the shaft can be fixed (deflection and slope equal to zero) or free (shear and moment equal to zero). The fixed end condition enables the program to calculate the natural frequencies of cantilever beams. Instead of using the lumped-parameter method, the program uses continuous integration of the differential equations of beam flexure across different shaft sections. The advantages of this method over the usual lumped-parameter method are less data preparation and better approximation of the distribution of the mass of the shaft. A main feature of the program is the nature of the output. The Calcomp plotter is used to produce a drawing of the shaft with superimposed deflection curves at the critical speeds, together with all pertinent information related to the shaft.

Trivisonno, R. J.↗

Seasonal and magnetic storm related changes in the thermosphere induced by eddy mixing

The possibility of explaining the seasonal and magnetic storm related changes in the thermosphere and the ionosphere through variation in eddy diffusion coefficient in the lower thermosphere is investigated theoretically by obtaining simultaneous numerical solutions of the relevant continuity, momentum, and energy balance equations in a self-consistent manner. It is shown that various important features of thermospheric seasonal behavior, including the winter helium bulge phenomenon, can be explained by assuming effective reduction in the winter-time eddy diffusion coefficient by about one order of magnitude. Storm related changes in the thermal structure and composition are seen to arise by assuming a downward shift of turbulence. Physical implications and relative significance of variation in thermospheric turbulence and large scale inter-hemispheric circulation are briefly discussed.

Sinha, A. K.↗

The influence of convection electric fields on thermal proton outflow from the ionosphere

The continuity, momentum and energy hydrodynamic equations for an O(+)-H(+) ionosphere have been solved self-consistently for steady state conditions when a perpendicular (convection) electric field is present. Comparison of the H(+) temperature profiles obtained with and without the electric field show that the effect of the electric field is to enhance the H(+) temperature at high altitudes from about 3600 to 6400 K. Due to ion heating by the electric field, there is a net reduction of O(+) in the F2-region as compared with the case of a nonconvecting ionosphere. When the reduction of O(+) is neglected, the electric field acts to increase the H(+) outward flux. However, when the reduction of O(+) is included, there is a net reduction in the outward H(+) flux. Nevertheless, the convection electric field still results in an increase in the rate of depletion of the F-region ionization due to H(+) outflow, by a factor of 2.2 for a 100 mV/m electric field.

Raitt, W. J.↗

A comprehensive model of the Venus ionosphere

Coupled time-dependent continuity-momentum and energy balance equations for the Venus ionosphere were simultaneously solved for CO2(+), O2(+), O(+), He(+), and H(+) densities and electron and ion temperatures for an altitude range of 120-500 km. Values of the solar zenith angle varied from 0 deg (subsolar point) to 90 deg (terminator). The calculations include the horizontal bulk transport of ions by neutral winds but not the horizontal diffusion. The two-stream photoelectron transport method was used to find the heating rates for the ambient electrons. Different boundary conditions were considered, and a nightside ionosphere was calculated. The results of these model calculations are in good agreement with measurements in the region of maximum electron density. Characteristics of topside and nightside densities and temperatures are discussed.

Chen, R. H.↗

Numerical analysis of turbulent wakes of axial flow turbomachinery rotor blades

The wakes of rotor blades are turbulent, three-dimensional, and are subjected to curvature and rotation effects. The objective of this study is to predict the development of such wakes and compare the predictions with the existing data. A finite difference procedure is employed in the numerical analysis of the wake utilizing the continuity, momentum, and turbulence closure equations in the rotating curvilinear and non-orthogonal coordinate system. The turbulence closure is affected by the two equations model and the simplified Reynolds stress model. The predictions from both the models agree well with the mean velocity and the turbulence data.

Hah, C.↗