Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “consistent scheme”

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 667 records · Page 37

Parallel processing of real-time dynamic systems simulation on OSCAR (Optimally SCheduled Advanced multiprocessoR)

Parallel processing of real-time dynamic systems simulation on a multiprocessor system named OSCAR is presented. In the simulation of dynamic systems, generally, the same calculation are repeated every time step. However, we cannot apply to Do-all or the Do-across techniques for parallel processing of the simulation since there exist data dependencies from the end of an iteration to the beginning of the next iteration and furthermore data-input and data-output are required every sampling time period. Therefore, parallelism inside the calculation required for a single time step, or a large basic block which consists of arithmetic assignment statements, must be used. In the proposed method, near fine grain tasks, each of which consists of one or more floating point operations, are generated to extract the parallelism from the calculation and assigned to processors by using optimal static scheduling at compile time in order to reduce large run time overhead caused by the use of near fine grain tasks. The practicality of the scheme is demonstrated on OSCAR (Optimally SCheduled Advanced multiprocessoR) which has been developed to extract advantageous features of static scheduling algorithms to the maximum extent.

Kasahara, Hironori↗

Winter QPF Sensitivities to Snow Parameterizations and Comparisons to NASA CloudSat Observations

Steady increases in computing power have allowed for numerical weather prediction models to be initialized and run at high spatial resolution, permitting a transition from larger scale parameterizations of the effects of clouds and precipitation to the simulation of specific microphysical processes and hydrometeor size distributions. Although still relatively coarse in comparison to true cloud resolving models, these high resolution forecasts (on the order of 4 km or less) have demonstrated value in the prediction of severe storm mode and evolution and are being explored for use in winter weather events . Several single-moment bulk water microphysics schemes are available within the latest release of the Weather Research and Forecast (WRF) model suite, including the NASA Goddard Cumulus Ensemble, which incorporate some assumptions in the size distribution of a small number of hydrometeor classes in order to predict their evolution, advection and precipitation within the forecast domain. Although many of these schemes produce similar forecasts of events on the synoptic scale, there are often significant details regarding precipitation and cloud cover, as well as the distribution of water mass among the constituent hydrometeor classes. Unfortunately, validating data for cloud resolving model simulations are sparse. Field campaigns require in-cloud measurements of hydrometeors from aircraft in coordination with extensive and coincident ground based measurements. Radar remote sensing is utilized to detect the spatial coverage and structure of precipitation. Here, two radar systems characterize the structure of winter precipitation for comparison to equivalent features within a forecast model: a 3 GHz, Weather Surveillance Radar-1988 Doppler (WSR-88D) based in Omaha, Nebraska, and the 94 GHz NASA CloudSat Cloud Profiling Radar, a spaceborne instrument and member of the afternoon or "A-Train" of polar orbiting satellites tasked with cataloguing global cloud characteristics. Each system provides a unique perspective. The WSR-88D operates in a surveillance mode, sampling cloud volumes of Rayleigh scatterers where reflectivity is proportional to the sixth moment of the size distribution of equivalent spheres. The CloudSat radar provides enhanced sensitivity to smaller cloud ice crystals aloft, as well as consistent vertical profiles along each orbit. However, CloudSat reflectivity signatures are complicated somewhat by resonant Mie scattering effects and significant attenuation in the presence of cloud or rain water. Here, both radar systems are applied to a case of light to moderate snowfall within the warm frontal zone of a cold season, synoptic scale storm. Radars allow for an evaluation of the accuracy of a single-moment scheme in replicating precipitation structures, based on the bulk statistical properties of precipitation as suggested by reflectivity signatures.

Molthan, Andrew↗

Electromagnetic particle simulation codes

Electromagnetic particle simulations solve the full set of Maxwell's equations. They thus include the effects of self-consistent electric and magnetic fields, magnetic induction, and electromagnetic radiation. The algorithms for an electromagnetic code which works directly with the electric and magnetic fields are described. The fields and current are separated into transverse and longitudinal components. The transverse E and B fields are integrated in time using a leapfrog scheme applied to the Fourier components. The particle pushing is performed via the relativistic Lorentz force equation for the particle momentum. As an example, simulation results are presented for the electron cyclotron maser instability which illustrate the importance of relativistic effects on the wave-particle resonance condition and on wave dispersion.

Pritchett, P. L.↗

Computation of viscous blast wave solutions with an upwind finite volume method

A fully conservative, viscous, implicit, upwind, finite-volume scheme for the thin-layer Navier-Stokes equations is described with application to blast wave flow fields. In this scheme, shocks are captured without the oscillations typical of central differencing techniques and wave speeds are accurately predicted. The finite volume philosophy ensures conservation and since boundary conditions are also treated conservatively, accurate reflections of waves from surfaces are assured. Viscous terms in the governing equations are treated in a manner consistent with the finite volume philosophy, resulting in very accurate prediction of boundary layer quantities. Numerical results are presented for four viscous problems: a steady boundary layer, a shock-induced boundary layer, a blast wave/cylinder interaction and a blast wave/supersonic missile interaction. Comparisons of the results with an established boundary layer code, similarity solution, and experimental data show excellent agreement.

Molvik, Gregory A.↗

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↗

A representational framework and user-interface for an image understanding workstation

Problems in image understanding involve a wide variety of data (e.g., image arrays, edge maps, 3-D shape models) and processes or algorithms (e.g., convolution, feature extraction, rendering). The underlying structure of an Image Understanding Workstation designed to support mulitple levels and types of representations for both data and processes is described, also the user-interface. The Image Understanding Workstation consists of two parts: the Image Understanding (IU) Framework, and the user-interface. The IU Framework is the set of data and process representations. It includes multiple levels of representation for data such as images (2-D), sketches (2-D), surfaces (2 1/2 D), and models (3-D). The representation scheme for processes characterizes their inputs, outputs, and parameters. Data and processes may reside on different classes of machines. The user-interface to the IU Workstation gives the user convenient access for creating, manipulating, transforming, and displaying image data. The user-interface follows the structure of the IU Framework and gives the user control over multiple types of data and processes. Both the IU Framework and user-interface are implemented on a LISP machine.

Schroeder, Joyce D.↗

Numerical simulation and analysis of a prefrontal squall line. I - Observations and basic simulation results. II - Propagation of the squall line as an internal gravity wave

Results of an observational and numerical study of the squall line that occurred on June 17-18, 1978, are presented. This squall line was initially triggered by the strong surface convergence along a cold front and stretched from Illinois to the Texas Panhandle. The squall line was aligned with the surface front during its initial development, but then propagated faster than the front, resulting in a separation of approximately 200 km by 0300 UTC and 300-400 km by 0600 UTC. The Colorado State University Regional Atmospheric Modeling System is used to model the squall-line development and propagation. Results are described from several experiments that tested the sensitivity to the use of the Kuo-type cumulus parameterization scheme and grid-scale microphysical processes. The movement of the squall line in the model is shown to be due to the propagation of a deep tropospheric internal gravity wave in a CISK-like process. The thermal and dynamic perturbations associated with the hypothesized wave are shown to be consistent with internal gravity wave theory.

Cram, Jennifer M.↗

Unsteady transonic Euler solutions using finite elements

A finite element solution of the unsteady Euler equations is presented and demonstrated for 2D airfoil configurations oscillating in transonic flows. Computations are performed by spatially discretizing the conservation equations using the Galerkin weighted residual method and then employing a multistage Runge-Kutta scheme to march forward in time. A mesh deformation scheme has been developed to efficiently move interior points in a smooth fashion as the airfoil undergoes rigid body pitch and plunge motion. Both steady and unsteady results are presented, and a comparison is made with solutions obtained using finite-volume techniques. The effects of using either a lumped or consistent mass matrix are presented; the finite element method provides an accurate solution for unsteady transonic flows about isolated airfoils.

Davis, Gary A.↗

Vortex ring instability and its sound

This work carries earlier finite-difference calculations of the Widnall instability of vortex rings into the late non-linear stage. Plots of energy in azimuthal Fourier modes indicate that low-order modes dominate at large times; their structure and dynamics remain unexplored, however. An attempt was made to calculate the acoustic signal using the theory of Mohring (1978), valid for unbounded flow. This theory shows that only low-order azimuthal modes contribute to the sound. As a check on the effects of axial periodicity and a slip wall at large radius imposed by the numerical scheme, the acoustic integrals were also computed in a truncated region. Half of the terms contributing to the sound have large differences between the two regions, and the results are therefore unreliable. The error is less severe for a contribution involving only the m = 2 mode, and its low frequency is consistent with a free elliptic bending wave on a thin ring.

Verzicco, R.↗

Multigrid solution of compressible turbulent flow on unstructured meshes using a two-equation model

The steady state solution of the system of equations consisting of the full Navier-Stokes equations and two turbulence equations has been obtained using a multigrid strategy of unstructured meshes. The flow equations and turbulence equations are solved in a loosely coupled manner. The flow equations are advanced in time using a multistage Runge-Kutta time-stepping scheme with a stability-bound local time step, while turbulence equations are advanced in a point-implicit scheme with a time step which guarantees stability and positivity. Low-Reynolds-number modifications to the original two-equation model are incorporated in a manner which results in well-behaved equations for arbitrarily small wall distances. A variety of aerodynamic flows are solved, initializing all quantities with uniform freestream values. Rapid and uniform convergence rates for the flow and turbulence equations are observed.

MULTIGRID SOLUTION↗

A Novel Four-Node Quadrilateral Smoothing Element for Stress Enhancement and Error Estimation

A four-node, quadrilateral smoothing element is developed based upon a penalized-discrete-least-squares variational formulation. The smoothing methodology recovers C1-continuous stresses, thus enabling effective a posteriori error estimation and automatic adaptive mesh refinement. The element formulation is originated with a five-node macro-element configuration consisting of four triangular anisoparametric smoothing elements in a cross-diagonal pattern. This element pattern enables a convenient closed-form solution for the degrees of freedom of the interior node, resulting from enforcing explicitly a set of natural edge-wise penalty constraints. The degree-of-freedom reduction scheme leads to a very efficient formulation of a four-node quadrilateral smoothing element without any compromise in robustness and accuracy of the smoothing analysis. The application examples include stress recovery and error estimation in adaptive mesh refinement solutions for an elasticity problem and an aerospace structural component.

Tessler, A.↗

Unsupervised classification of remote multispectral sensing data

The new unsupervised classification technique for classifying multispectral remote sensing data which can be either from the multispectral scanner or digitized color-separation aerial photographs consists of two parts: (a) a sequential statistical clustering which is a one-pass sequential variance analysis and (b) a generalized K-means clustering. In this composite clustering technique, the output of (a) is a set of initial clusters which are input to (b) for further improvement by an iterative scheme. Applications of the technique using an IBM-7094 computer on multispectral data sets over Purdue's Flight Line C-1 and the Yellowstone National Park test site have been accomplished. Comparisons between the classification maps by the unsupervised technique and the supervised maximum liklihood technique indicate that the classification accuracies are in agreement.

Su, M. Y.↗

Digital data reformatter/deserializer

A method and apparatus is presented for reformatting and de-serializing a serially-received sequence of data words, each consisting of a fixed number of binary data bits. A block of nm bits is serially fed into a shift register or serially-connected group of shift registers. In lieu of the(nm-1)th shifts, the bits are rearranged within the shift register in parallel fashion, according to a prescribed scheme. Shifting then continues, until the first bit of each data word appears in the last bit position in the shift register, at which time that data word is shifted in parallel into an output buffer stage, from which it is outputted in parallel, after a fixed delay.

Engel, A.↗

Electrostatic emissions between electron gyroharmonics in the outer magnetosphere

Electrostatic emissions are often observed by spacecraft in the outer magnetosphere at frequencies between low electron gyroharmonics. The emissions include the well-known odd half-harmonic emissions as well as weaker emissions which occur near the electron plasma frequency. A scheme is constructed for classifying the emissions, and it is shown that a theoretical model which has been previously used to explain '3/2' emissions can be extended to the other classes of emissions. All of the emissions therefore appear to be generated by the same basic mechanism: an unstable electron plasma distribution consisting of cold electrons (below 100 eV) and hot loss cone electrons (about 1 keV). Each emission class is associated with a particular range of model parameters; the wide-band electric field data can thus be used to infer the density and temperature of the cold plasma component. The model predicts that gyroharmonic emissions near the plasma frequency require large cold plasma densities.

Hubbard, R. F.↗

A study of the response of nonlinear springs

The various phases to developing a methodology for studying the response of a spring-reinforced arch subjected to a point load are discussed. The arch is simply supported at its ends with both the spring and the point load assumed to be at midspan. The spring is present to off-set the typical snap through behavior normally associated with arches, and to provide a structure that responds with constant resistance over a finite displacement. The various phases discussed consist of the following: (1) development of the closed-form solution for the shallow arch case; (2) development of a finite difference analysis to study (shallow) arches; and (3) development of a finite element analysis for studying more general shallow and nonshallow arches. The two numerical analyses rely on a continuation scheme to move the solution past limit points, and to move onto bifurcated paths, both characteristics being common to the arch problem. An eigenvalue method is used for a continuation scheme. The finite difference analysis is based on a mixed formulation (force and displacement variables) of the governing equations. The governing equations for the mixed formulation are in first order form, making the finite difference implementation convenient. However, the mixed formulation is not well-suited for the eigenvalue continuation scheme. This provided the motivation for the displacement based finite element analysis. Both the finite difference and the finite element analyses are compared with the closed form shallow arch solution. Agreement is excellent, except for the potential problems with the finite difference analysis and the continuation scheme. Agreement between the finite element analysis and another investigator's numerical analysis for deep arches is also good.

Hyer, M. W.↗

Tropical Atlantic Dust and Smoke Aerosol Variations Related to the Madden-Julian Oscillation in MODIS and MISR Observations

In this study, Moderate Resolution Imaging Spectroradiometer (MODIS) fine mode fraction and Multi-angle Imaging SpectroRadiometer (MISR) nonspherical fraction data are used to derive dust and smoke aerosol optical thickness (T(sub dust) and T(sub smoke)) over the tropical Atlantic in a complementary way: due to its wider swath, MODIS has 3-4 times greater sampling than MISR, but MISR dust discrimination is based on particle shape retrievals, whereas an empirical scheme is used for MODIS. MODIS and MISR show very similar dust and smoke winter climatologies. T(sub dust) is the dominant aerosol component over the tropical Atlantic, accounting for 40-70 percent of the total aerosol optical thickness (AOT), whereas T(sub smoke) is significantly smaller than T(sub dust). The consistency and high correlation between these climatologies and their daily variations lends confidence to their use for investigating the relative dust and smoke contributions to the total AOT variation associated with the Madden-Julian Oscillation (MJO). The temporal evolution and spatial patterns of the τdus anomalies associated with the MJO are consistent between MODIS and MISR: the magnitude of MJO-realted T(sub dust) anomalies is comparable to or even larger than that of the total T, while the T(sub smoke) anomaly represents about 15 percent compared to the total, which is quite different from their relative magnitudes to the total T on the climatological time scale. This suggests that dust and smoke are not influenced by the MJO in the same way. Based on correlation analysis, dust is strongly influenced by the MJO-modulated trade wind and precipitation anomalies, and can last as long as one MJO phase, whereas smoke is less affected.

Madden-Julian Oscillation↗

Hubble Space Telescope electrical power system simulation breadboard

The components and operation of the Hubble Space Telescope (HST) electrical power system (EPS) breadboard are described. The HST EPS breadboard consists of: (1) the power distribution breadboard, (2) the batteries, (3) the solar array simulators, (4) the load banks, (5) the charge current control hardware, and (6) a control computer; the functions of each of these elements are discussed. The HST EPS has been used to simulate changes in charge level, charge control schemes, and beta angle, failed cells, batteries, solar panel assemblies (SPAs), beginning of life SPAs, and end of life SPAs, and temperature differences. The test data reveal that the HST EPS is capable of operating over a wide range of conditions.

Bush, John R., Jr.↗

Three dimensional thermal analysis of rocket thrust chambers

A numerical model for the three-dimensional thermal analysis of rocket thrust chambers and nozzles has been developed. The input to the model consists of the composition of the fuel/oxidant mixture and flow rates, chamber pressure, coolant entrance temperature and pressure, dimensions of the engine, materials and the number of nodes in different parts of the engine. The model allows for temperature variation in three dimensions: axial, radial and circumferential directions and by implementing an iterative scheme, it provides nodal temperature distribution, rates of heat transfer, hot gas and coolant thermal and transport properties.

Naraghi, M. H. N.↗