Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “distribution grid”

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 937 records · Page 52

Investigation of the effects of external current systems on the MAGSAT data utilizing grid cell modeling techniques

Progress made in reducing MAGSAT data and displaying magnetic field perturbations caused primarily by external currents is reported. A periodic and repeatable perturbation pattern is described that arises from external current effects but appears as unique signatures associated with upper middle latitudes on the Earth's surface. Initial testing of the modeling procedure that was developed to compute the magnetic fields at satellite orbit due to current distributions in the ionosphere and magnetosphere is also discussed. The modeling technique utilizes a linear current element representation of the large scale space current system.

Klumpar, D. M.↗

Characteristics procedure for supersonic flows including consideration of viscous contributions to flow rotationality

A characteristics procedure has been developed which includes viscous and conductive transport terms normal to streamlines. By introducing the transport terms as corrections, the equations retain their hyperbolic character. The solutions have been developed using an inverse grid scheme in a streamline-normal network. Inclusion of the transport terms enables tracing of shock waves through the supersonic portion of the boundary layer so that shock reflections are well-located. Excellent agreement is obtained with the results of other characteristics procedures and with the experimental pressure distributions for the waisted body of Winter, Smith and Rotta.

Anderson, B. H.↗

Calculation of steady and oscillating airfoil flow fields via the Navier Stokes equations

A Navier-Stokes calculation procedure is applied to high Reynolds number flow about steady and unsteady airfoils. The procedure solves the ensemble averaged governing equations via a linearized block implicit (LBI) technique which in general converges to a nominal steady flow within 120 time steps. The grid used in the computation is highly stretched thus resolving the turbulent boundary layers. Calculations have been compared with data for both NACA 4412 airfoil at high incidence and an NACA 0012 airfoil oscillating in dynamic stall. In both cases, good agreement is noted between measured and calculated surface pressure distributions.

Shamroth, S. J.↗

Viscous analysis of three-dimensional rotor flows using a multigrid method

A three-dimensional code for rotating blade-row flow analysis was developed. The space discretization uses a cell-centered scheme with eigenvalues scaling for the artificial dissipation. The computational efficiency of a four-stage Runge-Kutta scheme is enhanced by using variable coefficients, implicit residual smoothing, and a full-multigrid method. An application is presented for the NASA rotor 67 transonic fan. Due to the blade stagger and twist, a zonal, non-periodic H-type grid is used to minimize the mesh skewness. The calculation is validated by comparing it with experiments in the range from the maximum flow rate to a near-stall condition. A detailed study of the flow structure near peak efficiency and near stall is presented by means of pressure distribution and particle traces inside boundary layers.

Arnone, A.↗

Reduction of blade-vortex interaction noise through porous leading edge

The effect of the porous leading edge of an airfoil on the blade-vortex interaction noise, which dominates the far-field acoustic spectrum of the helicopter, is investigated. The thin-layer Navier-Stokes equations are solved with a high-order upwind-biased scheme and a multizonal grid system. The Baldwin-Lomax turbulence model is modified for considering transpiration on the surface. The amplitudes of the propagating acoustic wave in the near field are calculated directly from the computation. The porosity effect on the surface is modeled in two ways: (1) imposition of prescribed transpiration velocity distribution and (2) calculation of transpiration velocity distribution by Darcy's law. Results show leading-edge transpiration can suppress pressure fluctuations at the leading edge during blade-vortex interaction and consequently reduce the amplitude of propagating noise by 30% at a maximum in the near field.

Lee, Soogab↗

Effects of gravity on sheared and nonsheared turbulent nonpremixed flames

The present numerical study is concerned with the fundamental physics of the multiway interaction between turbulence, chemical reaction, and buoyancy in a nonpremixed flame. The method of direct numerical simulation (DNS) is used to solve the instantaneous, three-dimensional governing equations. Because of the present supercomputer limitations, we consider two simple flow geometries, namely an initially uniform flow without shear (equivalent to grid-generated turbulence) and an initially uniform shear flow. In each flow, the fuel and oxidant initially exist as two separate streams. As the reactants mix, chemical reaction takes place and exothermic energy is released causing variations in density. In the presence of a gravity field, the spatial and temporal distributions of the induced buoyancy forces depend on the local density gradients and the direction of the gravitational acceleration. The effects of buoyancy include the generation of local shear, baroclinic production or destruction of vorticity, and countergradient heat and mass transport. Increased vorticity and small-scale turbulence promote further mixing and reaction. However, if the strain-rates become too high, local flame extinction can occur. Our objective is to gain an understanding of the complex interactions between the physical phenomena involved, with particular attention to the effects of buoyancy on the turbulence structure, flame behavior, and factors influencing flame extinction.

Elghobashi, Said↗

Effects of Gravity on Sheared and Nonsheared Turbulent Nonpremixed Flames

The present numerical study is concerned with the fundamental physics of the multi-way interaction between turbulence, chemical reaction and buoyancy in a nonpremixed flame. The method of direct numerical simulation (DNS) is used to solve the instantaneous three-dimensional governing equations (continuity, Navier-Stokes, species mass fractions and energy) under the zero-Mach number assumption. Because of the present supercomputer limitations, we consider two simple flow geometries, namely an initially uniform flow without shear (equivalent to grid-generated turbulence) and an initially uniform shear flow. In each flow, the fuel and oxidant initially exist as two separate streams. As the reactants mix, chemical reaction takes place and exothermic energy is released causing variations in density. In the presence of a gravity field, the spatial and temporal distributions of the induced buoyancy forces depend on the local density gradients and the direction of the gravitational acceleration.

Elghobashi, Said↗

Evaluation of Modeling NO2 Concentrations Driven by Satellite-Derived and Bottom-Up Emission Inventories Using In-Situ Measurements Over China

Chemical transport models together with emission inventories are widely used to simulate NO2 concentrations over China, but validation of the simulations with in situ measurements has been extremely limited. Here we use ground measurements obtained from the air quality monitoring network recently developed by the Ministry of Environmental Protection of China to validate modeling surface NO2 concentrations from the CHIMERE regional chemical transport model driven by the satellite-derived DECSO and the bottom-up MIX emission inventories. We applied a correction factor to the observations to account for the interferences of other oxidized nitrogen compounds (NOz), based on the modeled ratio of NO2 to NOz. The model accurately reproduces the spatial variability in NO2 from in situ measurements, with a spatial correlation coefficient of over 0.7 for simulations based on both inventories. A negative and positive bias is found for the simulation with the DECSO (slopeD0.74 and 0.64 for the daily mean and daytime only) and the MIX (slopeD1.3 and 1.1) inventories, respectively, suggesting an underestimation and overestimation of NOx emissions from corresponding inventories. The bias between observed and modeled concentrations is reduced, with the slope dropping from 1.3 to 1.0 when the spatial distribution of NOx emissions in the DECSO inventory is applied as the spatial proxy for the MIX inventory, which suggests an improvement of the distribution of emissions between urban and suburban or rural areas in the DECSO inventory compared to that used in the bottom-up inventory. A rough estimate indicates that the observed concentrations, from sites predominantly placed in the populated urban areas, may be 10-40% higher than the corresponding model grid cell mean. This reduces the estimate of the negative bias of the DECSO-based simulation to the range of -30 to 0% on average and more firmly establishes that the MIX inventory is biased high over major cities. The performance of the model is comparable over seasons, with a slightly worse spatial correlation in summer due to the difficulties in resolving the more active NOx photochemistry and larger concentration gradients in summer by the model. In addition, the model well captures the daytime diurnal cycle but shows more significant disagreement between simulations and measurements during nighttime, which likely produces a positive model bias of about 15% in the daily mean concentrations. This is most likely related to the uncertainty in vertical mixing in the model at night.

NO2 modelling↗

Mapping implicit spectral methods to distributed memory architectures

Spectral methods were proven invaluable in numerical simulation of PDEs (Partial Differential Equations), but the frequent global communication required raises a fundamental barrier to their use on highly parallel architectures. To explore this issue, a 3-D implicit spectral method was implemented on an Intel hypercube. Utilization of about 50 percent was achieved on a 32 node iPSC/860 hypercube, for a 64 x 64 x 64 Fourier-spectral grid; finer grids yield higher utilizations. Chebyshev-spectral grids are more problematic, since plane-relaxation based multigrid is required. However, by using a semicoarsening multigrid algorithm, and by relaxing all multigrid levels concurrently, relatively high utilizations were also achieved in this harder case.

Overman, Andrea L.↗

Numerical solutions of the Navier-Stokes equations for transonic afterbody flows

The time dependent Navier-Stokes equations in mass averaged variables are solved for transonic flow over axisymmetric boattail plume simulator configurations. Numerical solution of these equations is accomplished with the unsplit explict finite difference algorithm of MacCormack. A grid subcycling procedure and computer code vectorization are used to improve computational efficiency. The two layer algebraic turbulence models of Cebeci-Smith and Baldwin-Lomax are employed for investigating turbulence closure. Two relaxation models based on these baseline models are also considered. Results in the form of surface pressure distribution for three different circular arc boattails at two free stream Mach numbers are compared with experimental data. The pressures in the recirculating flow region for all separated cases are poorly predicted with the baseline turbulence models. Significant improvements in the predictions are usually obtained by using the relaxation models.

Swanson, R. C., Jr.↗

Flow aerodynamics modeling of an MHD swirl combustor - Calculations and experimental verification

The paper describes a computer code for calculating the flow dynamics of a constant-density flow in the second-stage trumpet shaped nozzle section of a two-stage MHD swirl combustor for application to a disk generator. The primitive pressure-velocity variable, finite-difference computer code has been developed for the computation of inert nonreacting turbulent swirling flows in an axisymmetric MHD model swirl combustor. The method and program involve a staggered grid system for axial and radial velocities, and a line relaxation technique for the efficient solution of the equations. The code produces as output the flow field map of the nondimensional stream function, axial and swirl velocity. It was found that the best location for seed injection to obtain a uniform distribution at the combustor exit is in the central location for seed injected at the entrance to the second stage combustor.

Gupta, A. K.↗

Procedures for analysis of spatial relationships among ship survey data and sea surface temperature

The establishment of a limited spatial data base for the U.S. eastern seaboard vicinity is discussed along with the demonstration of computer assisted analysis techniques for investigating spatial patterns and relationships among ship survey data and remotely sensed sea surface temperature. Ship survey variables included concentrations of two zooplankton, two icthyoplankton, and two fish species, in addition to physical data of depth to bottom and surface and bottom water temperatures. Continuous spatial distributions of these data were created by both weighted nearest neighbor and iterative smoothing interpolation techniques. Maps of surface water temperature were created by digitizing GOES satellite images. All mapped data were spatially registered by conversion of latitude and longitude coordinates to rotated Lambert conic conformal rectilinear coordinates and stored in grid format of approximately one hundred square kilometers per cell. The analysis of these data include the generation of statistical summaries and maps describing the joint occurrence among variables.

Berry, J. K.↗

A critical assessment of viscous models of trench topography and corner flow

Stresses for Newtonian viscous flow in a simple geometry (e.g., corner flow, bending flow) are obtained in order to study the effect of imposed velocity boundary conditions. Stress for a delta function velocity boundary condition decays as 1/R(2); for a step function velocity, stress goes as 1/R; for a discontinuity in curvature, the stress singularity is logarithmic. For corner flow, which has a discontinuity of velocity at a certain point, the corresponding stress has a 1/R singularity. However, for a more realistic circular-slab model, the stress singularity becomes logarithmic. Thus the stress distribution is very sensitive to the boundary conditions, and in evaluating the applicability of viscous models of trench topography it is essential to use realistic geometries. Topography and seismicity data from northern Hoshu, Japan, were used to construct a finite element model, with flow assumed tangent to the top of the grid, for both Newtonian and non-Newtonian flow (power law 3 rheology). Normal stresses at the top of the grid are compared to the observed trench topography and gravity anomalies. There is poor agreement. Purely viscous models of subducting slables with specified velocity boundary conditions do not predict normal stress patterns compatible with observed topography and gravity. Elasticity and plasticity appear to be important for the subduction process.

Zhang, J.↗

A critical assessment of viscous models of trench topography and corner flow

Stresses for Newtonian viscous flow in a simple geometry (e.g., corner flow, bending flow) are obtained in order to study the effect of imposed velocity boundary conditions. Stress for a delta function velocity boundary condition decays as 1/R(2); for a step function velocity, stress goes as 1/R; for a discontinuity in curvature, the stress singularity is logarithmic. For corner flow, which has a discontinuity of velocity at a certain point, the corresponding stress has a 1/R singularity. However, for a more realistic circular-slab model, the stress singularity becomes logarithmic. Thus the stress distribution is very sensitive to the boundary conditions, and in evaluating the applicability of viscous models of trench topography it is essential to use realistic geometries. Topography and seismicity data from northern Hoshu, Japan, were used to construct a finite element model, with flow assumed tangent to the top of the grid, for both Newtonian and non-Newtonian flow (power law 3 rheology). Normal stresses at the top of the grid are compared to the observed trench topography and gravity anomalies. There is poor agreement. Purely viscous models of subducting slables with specified velocity boundary conditions do not predict normal stress patterns compatible with observed topography and gravity. Elasticity and plasticity appear to be important for the subduction process.

Zhang, J.↗

Turbofan forced mixer lobe flow modeling. 2: Three-dimensional inviscid mixer analysis (FLOMIX)

A three-dimensional potential analysis (FLOMIX) was formulated and applied to the inviscid flow over a turbofan foced mixer. The method uses a small disturbance formulation to analytically uncouple the circumferential flow from the radial and axial flow problem, thereby reducing the analysis to the solution of a series of axisymmetric problems. These equations are discretized using a flux volume formulation along a Cartesian grid. The method extends earlier applications of the Cartesian method to complex cambered geometries. The effects of power addition are also included within the potential formulation. Good agreement is obtained with an alternate small disturbance analysis for a high penetration symmetric mixer in a planar duct. In addition, calculations showing pressure distributions and induced secondary vorticity fields are presented for practical trubofan mixer configurations, and where possible, comparison was made with available experimental data. A detailed description of the required data input and coordinate definition is presented along with a sample data set for a practical forced mixer configuration. A brief description of the program structure and subroutines is also provided.

Barber, T.↗

The correlated-k method for radiation calculations in nonhomogeneous atmospheres

The accuracy of the correlated-k method, which is a technique for radiation calculations with spectrally averaged data in nonhomogeneous atmospheres, is investigated. Comparisons are made for scattering and absorbing atmospheres containing CO2, H2O, and O3, and it is concluded that: (1) the errors in correlated-k are generally of order of magnitude 1 percent, (2) much larger errors occur only when a radiative quantity is very much smaller than its average value, (3) errors do not depend systematically on the gas molecule, the distributions of gases and aerosols, or on the aerosol optical properties, and (4) errors do not systematically increase with the order of differencing. It is shown that the multiplication property for transmission by overlapping bands can be incorporated into correlated-k, that temperature effects can be interpolated on a coarse grid, and that 10 quadrature points are often sufficient to average over complex spectral intervals containing thousands of lines.

Goody, Richard↗

Numerical simulation of aerothermal loads in hypersonic engine inlets due to shock impingement

The effect of shock impingement on an axial corner simulating the inlet of a hypersonic vehicle engine is modeled using a finite-difference procedure. A three-dimensional dynamic grid adaptation procedure is utilized to move the grids to regions with strong flow gradients. The adaptation procedure uses a grid relocation stencil that is valid at both the interior and boundary points of the finite-difference grid. A linear combination of spatial derivatives of specific flow variables, calculated with finite-element interpolation functions, are used as adaptation measures. This computational procedure is used to study laminar and turbulent Mach 6 flows in the axial corner. The description of flow physics and qualitative measures of heat transfer distributions on cowl and strut surfaces obtained from the analysis are compared with experimental observations. Conclusions are drawn regarding the capability of the numerical scheme for enhanced modeling of high-speed compressible flows.

Ramakrishnan, R.↗

Finite difference solution for transient radiative cooling of a conducting semitransparent square region

Transient solutions were obtained for a square region of heat conducting semitransparent material cooling by thermal radiation. The region is in a vacuum environment, so energy is dissipated only by radiation from within the medium leaving through its boundaries. The effect of heat conduction during the transient is to partially equalize the internal temperature distribution. As the optical thickness of the region is increased, the temperature gradients increase near the boundaries and corners, unless heat conduction is large. The solution procedure must provide accurate temperature distributions in these regions to prevent error in the calculated radiation losses. Two-dimensional numerical Gaussian integration is used to obtain the local radiative source term. A finite difference procedure with variable space and time increments is used to solve the transient energy equation. Variable spacing was used to concentrate grid points in regions with large temperature gradients.

Siegel, R.↗