Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “coordinate computing”

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 469 records · Page 26

On the character of turbulent-like flows in self-consistent models of core-collapse supernovae

Neutrino-driven convection plays a crucial role in the development of core-collapse supernova (CCSN) explosions. However, the complex mechanism that triggers the shock revival and the subsequent explosion has remained inscrutable for many decades. Multidimensional simulations suggest that the growth of fluid instabilities, especially development of convection, will determine the morphology of the explosion. We have performed 3D simulations using spherical-polar coordinates covering a reduced angular extent (90° computational domain), and with angular resolutions of 2°, 1°, 1/2°, and 1/4°, to study the development of turbulent-like flows in core-collapse supernova explosions on a time scale of order 100 ms. We have employed the multi-physics Chimera code that includes detailed nuclear physics and spectral neutrino transport. Coarse resolution models do not develop an inertial range, presumably due to the bottleneck effect, such that the energy is prevented from cascading down to small scales and tends to accumulate at large scales. High-resolution models instead, start to recover the k −5/3 scaling of Kolmogorov's theory. Stochasticity and few simulation samples limit our ability to predict the development of explosions. Over the simulated time period, our models show no clear trend in improving (or diminishing) conditions for explosion as the angular resolution is increased. However, we find that disordered flow provides an effective pressure, as characterized by the Reynolds stress, behind the shock. This contribution, ~40%–50% of the thermal pressure, aids shock revival and thus the development of the explosion. Finally, we show that the kinetic energy power spectrum of reduced angular extent and full 4π models are consistent, thus indicating that a 90° computational domain is an adequate configuration to study the character of turbulence in CCSNe.

Casanova Bustamante, Jordi↗

MONKES: a fast neoclassical code for the evaluation of monoenergetic transport coefficients in stellarator plasmas

Abstract MONKES is a new neoclassical code for the evaluation of monoenergetic transport coefficients in stellarators. By means of a convergence study and benchmarks with other codes, it is shown that MONKES is accurate and efficient. The combination of spectral discretization in spatial and velocity coordinates with block sparsity allows MONKES to compute monoenergetic coefficients at low collisionality, in a single core, in approximately one minute. MONKES is sufficiently fast to be integrated into stellarator optimization codes for direct optimization of the bootstrap current and to be included in predictive transport suites. The code and data from this paper are available at https://github.com/JavierEscoto/MONKES/ .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Nearly hyperuniform, nonhyperuniform, and antihyperuniform density fluctuations in two-dimensional transition metal dichalcogenides with defects

Hyperuniform many-body systems in d-dimensional Euclidean space R d are characterized by completely suppressed (normalized) infinite-wavelength density fluctuations, and appear to be endowed with novel exotic physical properties. Recently, hyperuniform systems of disordered varieties have been observed in the context of various atomic-scale two-dimensional (2D) materials. Here, we analyze the effects of localized defects on the density fluctuations across length scales and on the hyperuniformity property of experimental samples of two-dimensional transition metal dichalcogenides. In particular, we extract atomic coordinates from time series annular dark field-scanning transmission electron microscopy (ADF-STEM) imaging data of 2D tungsten chalcogenides with the 2H structure (Te-doped 2H-WSe 2 ) showing continuous development and evolution of electron-beam induced defects, and construct the corresponding chemical-bonding informed coordination networks between the atoms. We then compute a variety of pair statistics and bond-orientational statistics to characterize the samples. At low defect concentrations, the corresponding materials are nearly hyperuniform, characterized by significantly suppressed scattering at the zero wave-number limit (omitting forward/ballistic scattering). As more defects are introduced, the (approximate) hyperuniformity of the materials is gradually destroyed, and the system becomes non-hyperuniform even when the material still contains a significant amount of crystalline regions. At high defect concentrations, the structures become antihyperuniform with diverging (normalized) large-scale density fluctuations, mimicking those typically observed at the thermal critical points associated with phase transitions. Overall, the defected materials possess varying degrees of orientation order, and there is apparently no intermediate hexatic phase emerging. To understand the observed nearly hyperuniform density fluctuations in the slightly defected materials, we construct a minimalist structural model and demonstrate that the experimental samples can be essentially viewed as perturbed honeycomb crystals with small correlated displacements and double chalcogen vacancies. Moreover, the small correlated displacements alone can significantly degrade hyperuniformity of the perfect honeycomb structure. Therefore, even a small amount of vacancies, when coupled with correlated displacements, can completely destroy hyperuniformity of the system.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Computer program utilizes FORTRAN 4 subroutines for contour plotting

Computer program constructs lists of xy-coordinate pairs that define contour curves for an arbitrary given function of two variables and transmits these lists to plotting equipment to produce contour plots. The principal subroutine, CONTUR, is independent of any specific system of plotting subroutines and equipment.

Block, N.↗

Numerical solutions of the Navier-Stokes equations for the supersonic laminar flow over a two-dimensional compression corner

Numerical solutions have been obtained for the supersonic, laminar flow over a two-dimensional compression corner. These solutions were obtained as steady-state solutions to the unsteady Navier-Stokes equations using the finite difference method of Brailovskaya, which has second-order accuracy in the spatial coordinates. Good agreement was obtained between the computed results and wall pressure distributions measured experimentally for Mach numbers of 4 and 6.06, and respective Reynolds numbers, based on free-stream conditions and the distance from the leading edge to the corner. In those calculations, as well as in others, sufficient resolution was obtained to show the streamline pattern in the separation bubble. Upstream boundary conditions to the compression corner flow were provided by numerically solving the unsteady Navier-Stokes equations for the flat plate flow field, beginning at the leading edge. The compression corner flow field was enclosed by a computational boundary with the unknown boundary conditions supplied by extrapolation from internally computed points.

Carter, J. E.↗

Medical Automation System at the Marshall Space Flight Center

A computer system is reported for scheduling and coordinating physical examinations for groups of people who work with chemicals, lasers, X-rays, isotopic sources, toxic fuels, adhesives, and exotic metals. Complete medical data on the examined population are continuously updated for a broad medical master file. Statistical methods are employed to project progressive changes in the health status of these employees for possible clinical interventions.

Spraul, J. H.↗

Finite-volume solution of the Euler equations for steady three-dimensional transonic flow

A numerical finite-volume method is proposed for solving the Euler equations of a steady three-dimensional flow. The subsonic regions embedded in steady supersonic flow about aircraft are computed by introducing a time-dependent coordinate system that maps the fuselage and wing into one coordinate surface. The suitability of the method to a delta-winged aircraft is demonstrated, and its applicability to rectangular-winged aircraft and problems in internal transonic flow is pointed out.

Rizzi, A.↗

Development of a three-dimensional Navier-Stokes code on CDC star-100 computer

A three-dimensional code in body-fitted coordinates was developed using MacCormack's algorithm. The code is structured to be compatible with any general configuration, provided that the metric coefficients for the transformation are available. The governing equations are developed in primitive variables in order to facilitate the incorporation of physical boundary conditions and turbulence-closure models. MacCormack's two-step, unsplit, time-marching algorithm is used to solve the unsteady Navier-Stokes equations until steady-state solution is achieved. Cases discussed include (1) flat plate in supersonic free stream; (2) supersonic flow along an axial corner; (3) subsonic flow in an axial corner at M infinity = 0.95; and (4) supersonic flow in an axial corner at M infinity 1.5.

Vatsa, V. N.↗

Solution of three-dimensional time-dependent viscous flows

A procedure for solving three-dimensional, time-dependent turbulent flows is presented. The consistently split Linearized Block Implicit (LBI) scheme is used in conjunction with the QR Operator scheme to solve an approximate form of the Navier-Stokes equations in generalized nonorthogonal coordinates employing physical velocity components. Results of computations for both second order finite differences and the fourth order generalized Operator Compact Implicit (OCI) schemes are presented.

Weinberg, B. C.↗

The control network of Iapetus

A control network of the Saturnian satellite Iapetus has been established photogrammetrically from pictures taken by the two Voyager spacecraft. Coordinates of 62 control points have been computed and listed; pixel measurements of these points were made on 14 Voyager 1 and 66 Voyager 2 pictures. Some of these points are identified on the preliminary U.S. Geological Survey map of Iapetus and many are identified by name. The Voyager 1 and Voyager 2 pictures covered limited regions of the satellite's surface and contained no overlapping areas. The longitude system on Iapetus is defined by the crater Almeric; the 276 deg meridian passes through the center of this crater. The obliquity of Iapetus has been measured as 0.4 deg + or - 1.6 deg. The mean radius of Iapetus has been determined at 718 + or - 8 km.

Davies, M. E.↗

Measuring: The quantifying art

The photograph of a comet or asteroid against background stars contains information in analog form regarding the relative positions of comet and star images. Measuring is the procedure for covnerting this information to digital form for computation of the object's accurate equatorial coordinates at the time of observation. The use of one- and two-coordinate measuring engines to measure the positions of images which may be elongated, comatic, or both, as well as round images is discussed along with the training of novice measurers to cope with such images. The grading of plates for potential accuracy of positions, preparation of plates for measuring, the accuracy of measurement, effects of telescope focal length and reduction technique upon the accuracy of positions, and the procedures and checks which may prevent erroneous positions also are discussed.

Gibson, J.↗

A unified approach to computer analysis and modeling of spacecraft environmental interactions

A new, coordinated, unified approach to the development of spacecraft plasma interaction models is proposed. The objective is to eliminate the unnecessary duplicative work in order to allow researchers to concentrate on the scientific aspects. By streamlining the developmental process, the interchange between theories and experimentalists is enhanced, and the transfer of technology to the spacecraft engineering community is faster. This approach is called the UNIfied Spacecraft Interaction Model (UNISIM). UNISIM is a coordinated system of software, hardware, and specifications. It is a tool for modeling and analyzing spacecraft interactions. It will be used to design experiments, to interpret results of experiments, and to aid in future spacecraft design. It breaks a Spacecraft Ineraction analysis into several modules. Each module will perform an analysis for some physical process, using phenomenology and algorithms which are well documented and have been subject to review. This system and its characteristics are discussed.

Katz, I.↗

National plans for aircraft icing and improved aircraft icing forecasts and associated warning services

Recently, the United States has increased its activities related to aircraft icing in numerous fields: ice phobics, revised characterization of icing conditions, instrument development/evaluation, de-ice/anti-ice devices, simulated supercooled clouds, computer simulation and flight tests. The Federal Coordinator for Meteorology is involved in two efforts, one a National Plan on Aircraft Icing and the other a plan for Improved Aircraft Icing Forecasts and Associated Warning Services. These two plans will provide an approved structure for future U.S. activities related to aircraft icing. The recommended activities will significantly improve the position of government agencies to perform mandated activities and to enable U.S. manufacturers to be competitive in the world market.

Pass, Ralph P.↗

Recursive linearization of multibody dynamics equations of motion

The equations of motion of a multibody system are nonlinear in nature, and thus pose a difficult problem in linear control design. One approach is to have a first-order approximation through the numerical perturbations at a given configuration, and to design a control law based on the linearized model. Here, a linearized model is generated analytically by following the footsteps of the recursive derivation of the equations of motion. The equations of motion are first written in a Newton-Euler form, which is systematic and easy to construct; then, they are transformed into a relative coordinate representation, which is more efficient in computation. A new computational method for linearization is obtained by applying a series of first-order analytical approximations to the recursive kinematic relationships. The method has proved to be computationally more efficient because of its recursive nature. It has also turned out to be more accurate because of the fact that analytical perturbation circumvents numerical differentiation and other associated numerical operations that may accumulate computational error, thus requiring only analytical operations of matrices and vectors. The power of the proposed linearization algorithm is demonstrated, in comparison to a numerical perturbation method, with a two-link manipulator and a seven degrees of freedom robotic manipulator. Its application to control design is also demonstrated.

Lin, Tsung-Chieh↗

A method of smooth bivariate interpolation for data given on a generalized curvilinear grid

A method of locally bicubic interpolation is presented for data given at the nodes of a two-dimensional generalized curvilinear grid. The physical domain is transformed to a computational domain in which the grid is uniform and rectangular by a generalized curvilinear coordinate transformation. The metrics of the transformation are obtained by finite differences in the computational domain. Metric derivatives are determined by repeated application of the chain rule for partial differentiation. Given the metrics and the metric derivatives, the partial derivatives required to determine a locally bicubic interpolant can be estimated at each data point using finite differences in the computational domain. A bilinear transformation is used to analytically transform the individual quadrilateral cells in the physical domain into unit squares, thus allowing the use of simple formulas for bicubic interpolation.

Zingg, David W.↗

A 3D-CFD code for accurate prediction of fluid flows and fluid forces in seals

Current and future turbomachinery requires advanced seal configurations to control leakage, inhibit mixing of incompatible fluids and to control the rotodynamic response. In recognition of a deficiency in the existing predictive methodology for seals, a seven year effort was established in 1990 by NASA's Office of Aeronautics Exploration and Technology, under the Earth-to-Orbit Propulsion program, to develop validated Computational Fluid Dynamics (CFD) concepts, codes and analyses for seals. The effort will provide NASA and the U.S. Aerospace Industry with advanced CFD scientific codes and industrial codes for analyzing and designing turbomachinery seals. An advanced 3D CFD cylindrical seal code has been developed, incorporating state-of-the-art computational methodology for flow analysis in straight, tapered and stepped seals. Relevant computational features of the code include: stationary/rotating coordinates, cylindrical and general Body Fitted Coordinates (BFC) systems, high order differencing schemes, colocated variable arrangement, advanced turbulence models, incompressible/compressible flows, and moving grids. This paper presents the current status of code development, code demonstration for predicting rotordynamic coefficients, numerical parametric study of entrance loss coefficients for generic annular seals, and plans for code extensions to labyrinth, damping, and other seal configurations.

Athavale, M. M.↗

Twin Tail/Delta Wing Configuration Buffet Due to Unsteady Vortex Breakdown Flow

The buffet response of the twin-tail configuration of the F/A-18 aircraft; a multidisciplinary problem, is investigated using three sets of equations on a multi-block grid structure. The first set is the unsteady, compressible, full Navier-Stokes equations. The second set is the coupled aeroelastic equations for bending and torsional twin-tail responses. The third set is the grid-displacement equations which are used to update the grid coordinates due to the tail deflections. The computational model consists of a 76 deg-swept back, sharp edged delta wing of aspect ratio of one and a swept-back F/A-18 twin-tails. The configuration is pitched at 32 deg angle of attack and the freestream Mach number and Reynolds number are 0.2 and 0.75 x 10(exp 6) respectively. The problem is solved for the initial flow conditions with the twin tail kept rigid. Next, the aeroelastic equations of the tails are turned on along with the grid-displacement equations to solve for the uncoupled bending and torsional tails response due to the unsteady loads produced by the vortex breakdown flow of the vortex cores of the delta wing. Two lateral locations of the twin tail are investigated. These locations are called the midspan and inboard locations.

Kandil, Osama A.↗