Search NASA⌕ Search

SEARCH · Search NASA

Results for “Grid 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 343 records · Page 19

Operational Evolution of FTS3: A DevOps Driven Approach to Elastic Operations

The File Transfer Service (FTS3) is a distributed data movement service developed at CERN and widely used to transfer data across the Worldwide LHC Computing Grid (WLCG). At Fermilab, FTS3 supports data transfers for multiple experiments, including Intensity Frontier experiments such as DUNE, enabling reliable data movement between WebDAV endpoints in Europe and the Americas.​ At CHEP 2021, we reported on the initial containerized deployment of FTS3 on OKD, the community Kubernetes distribution of Red Hat OpenShift. In this work, we present the subsequent evolution of this deployment, focusing on new operational capabilities introduced to improve scalability, robustness, and long-term maintainability.​ We describe the adoption of more secure and reproducible container build workflows, the integration of DevOps-driven operational practices, and enhancements in monitoring and automation. A key new result is the introduction of horizontal scaling and elastic resource management, allowing FTS3 components to dynamically adapt to workload variations while maintaining service reliability. We also discuss improvements in fault tolerance and operational procedures derived from production experience.​ Finally, we summarize lessons learned from operating FTS3 as a Kubernetes-native service and outline how these developments have improved the resilience and efficiency of data movement operations at Fermilab.

Munoz Flores, Victor Leopoldo [Fermilab]↗

High y + Shear-Stress Turbulence Implementation for High Flux Isotope Reactor Narrow Channel Flows

The research objective of this work was to improve the engineering predictions of the turbulence characteristics of flows in curved narrow channels. Such channel flows are commonly encountered in nuclear research and test reactors, with one of them being the high-flux isotope reactor (HFIR). Research reactors bear high heat fluxes, and the proper computing of turbulence is paramount for safe and reliable reactor operation. The study builds on the results of a previous direct numerical simulation of turbulence to inform a well-known Reynolds-averaged Navier–Stokes shear-stress turbulence model and improves its accuracy in simulating parallel channel flows. A new formulation of the loss term in the dissipation conservation equation is suggested. Combined with high wall distance computational grids, the new implementation provides a fast-running flow solution, suitable for engineering purposes. Model generalization for parallel channel flows, in a broader range of frictional Reynolds numbers, is suggested by introducing a new form of the model constants.

CFD↗

Extended Conformal Mappings for Supersonic Aircraft Calculations

A gasdynamical analysis is presented for a three-dimensional, supersonic, inviscid, steady, shockless flow past an arbitrary airframe, using computational grids. The analysis, which includes special treatments for body points and bow-shock points, relies on the equations of motion written in terms of logarithm of pressure entropy and two angles reprising the velocity vector. A FORTRAN code was employed (sample geometry is shown). The flow analysis is considered reliable for aircraft cross-sections that are elliptical in shape.

Moretti, G.↗

Influence of nonconservative differencing on transonic streamline shapes

A computer program recently developed by South and Brandt (1976) which contained the Murman (1973) conservative finite-difference scheme is easily modified to use the Garabedian and Korn (1971) nonconservative finite difference scheme. This program solves the transonic small disturbance equation for only symmetric flow, but incorporates several iterative solution techniques. Results are presented for the case where the equally spaced computational grid extended to infinity in both the streamwise and normal directions. Streamline shapes are obtained along several grid lines by a streamwise integration of the normal component of the perturbation velocity. Comparison cases are run for a 10% thick parabolic arc airfoil at zero incidence for freestream Mach numbers of 0, 0.70, 0.84, and 0.95. It is shown that the use of a nonconservative finite-difference scheme in transonic flow calculations destroys the global mass balance when shocks are present. This lack of mass balance may prove to be more crucial in the case of an unconfined external flow.

Newman, P. A.↗

Numerical calculation of the transonic flow past a swept wing

A numerical method is presented for analyzing the transonic potential flow past a lifting, swept wing. A finite difference approximation to the full potential equation is solved in a coordinate system which is nearly conformally mapped from the physical space in planes parallel to the symmetry plane, and reduces the wing surface to a portion of one boundary of the computational grid. A coordinate invariant, rotated difference scheme is used, and the difference equations are solved by relaxation. The method is capable of treating wings of arbitrary planform and dihedral, although approximations in treating the tips and vortex sheet make its accuracy suspect for wings of small aspect ratio. Comparisons of calculated results with experimental data are shown for examples of both conventional and supercritical transport wings. Agreement is good for both types, but it was found necessary to account for the displacement effect of the boundary layer for the supercritical wing, presumably because of its greater sensitivity to changes in effective geometry.

Jameson, A.↗

Numerical solution of three-dimensional free turbulent shear flows

Numerical results are presented for three-dimensional compressible turbulent jet and wake flows. An alternating direction implicit numerical procedure is used to solve the finite-difference form of the parabolic elliptic Navier-Stokes equations. A coordinate transformation maps the boundaries at infinity into a finite computational domain in order to properly specify infinity boundary conditions as well as contain the downstream growth of the viscous flow field in a fixed computational grid. Turbulence closure is achieved through an algebraic mixing length eddy viscosity model. Numerical results for supersonic flow are presented for an axisymmetric jet, an elliptical jet, an elliptical wake, and two interacting rectangular jets. Experimental data were not available for comparison with the numerical results. However, the results compare well with empirical results for free shear flows.

Oh, Y. H.↗

Meteorite impact ejecta - Dependence of mass and energy lost on planetary escape velocity

The amounts of material and energy which escape a planet in a meteorite impact event is calculated as functions of impact and escape velocities. Results are obtained from the computed flow induced by the impact of iron and gabbroic anorthosite spheres onto a half-space of anorthosite at impact velocities of 5 to 45 km/sec. The impact-induced flows were determined by a numerical method using the mass, momentum, and energy conservation relations in finite-difference approximation, within an Eulerian computational grid. The impact velocities at which ejecta losses equal meteorite mass gains are found to be approximately 20, 35, and 45 km/sec for anorthosite objects and approximately 25, 35, and 40 km/sec for iron objects striking anorthosite surfaces for the gravity fields of the moon, Mercury and Mars.

Okeefe, J. D.↗

Numerical method for predicting flow characteristics and performance of nonaxisymmetric nozzles, theory

The theoretical foundation and formulation of a numerical method for predicting the viscous flowfield in and about isolated three dimensional nozzles of geometrically complex configuration are presented. High Reynolds number turbulent flows are of primary interest for any combination of subsonic, transonic, and supersonic flow conditions inside or outside the nozzle. An alternating-direction implicit (ADI) numerical technique is employed to integrate the unsteady Navier-Stokes equations until an asymptotic steady-state solution is reached. Boundary conditions are computed with an implicit technique compatible with the ADI technique employed at interior points of the flow region. The equations are formulated and solved in a boundary-conforming curvilinear coordinate system. The curvilinear coordinate system and computational grid is generated numerically as the solution to an elliptic boundary value problem. A method is developed that automatically adjusts the elliptic system so that the interior grid spacing is controlled directly by the a priori selection of the grid spacing on the boundaries of the flow region.

Thomas, P. D.↗

Utility of a finite element solution algorithm for initial-value problems

The Galerkin criterion within a finite element Weighted Residuals formulation is employed to establish an implicit solution algorithm for an initial-value partial differential equation. Numerical solutions of a transient parabolic and a hyperbolic equation, obtained using linear, quadratic and two cubic finite element basis functions, are employed to quantize accuracy and confirm and refine theoretical convergence rate estimates. The linear basis algorithm for the hyperbolic equation displays excellent accuracy on a coarse computational grid and a high-order convergence rate with discretization refinement. Good accuracy and a strong convergence rate in surface flux are determined for a nonhomogeneous Neumann boundary constraint applied to a parabolic equation. The results amply demonstrate the impact of the nondiagonal finite element initial-value matrix structure on solution accuracy and/or convergence rate.

Baker, A. J.↗

Large eddy simulation of turbulent channel flow: ILLIAC 5 calculation

The capabilities of large eddy simulation in the prediction and analyses of wall-bounded turbulent shear flows are demonstrated. The dynamical equations for large scale field motions are derived. The computational grid network is described and its relation to the observed physical length scales in the flow are discussed. Some aspects of the mechanics and structure of the flow are examined both in the vicinity of the wall and in regions away from the wall. An attempt is made to correlate numerical results with laboratory observations. Other significant observations and conclusions are presented.

Kim, J.↗

Analysis and calculation of macrosegregation in a casting ingot. MPS solidification model. Volume 1: Formulation and analysis

The physical and numerical formulation of a model for the horizontal solidification of a binary alloy is described. It can be applied in an ingot. The major purpose of the model is to calculate macrosegregation in a casting ingot which results from flow of interdendritic liquid during solidification. The flow, driven by solidification contractions and by gravity acting on density gradients in the interdendritic liquid, was modeled as flow through a porous medium. The symbols used are defined. The physical formulation of the problem leading to a set of equations which can be used to obtain: (1) the pressure field; (2) the velocity field: (3) mass flow and (4) solute flow in the solid plus liquid zone during solidification is presented. With these established, the model calculates macrosegregation after solidification is complete. The numerical techniques used to obtain solution on a computational grid are presented. Results, evaluation of the results, and recommendations for future development of the model are given. The macrosegregation and flow field predictions for tin-lead, aluminum-copper, and tin-bismuth alloys are included as well as comparisons of some of the predictions with published predictions or with empirical data.

Maples, A. L.↗

Shockless design and analysis of transonic blade shapes

A fast computer program was developed to eliminate the shocks by slightly altering portions of the contour of a given airfoil in the cascade. The program can be used in two basic modes: (1) An analysis for steady, transonic, potential flow through a given planar cascade of airfoils and (2) a design for converting a given cascade into a shockless transonic cascade. The design mode can automatically be followed by the analysis mode, which confirms that the flow field is shock free. The program generates its own multilevel boundary conforming computational grids and solves a full potential equation in a fully conservative form. The shockless design is performed by implementing Sobieczky's fictitious-gas elliptic continuation concept.

Dulikravich, D. S.↗

A potential-flow/boundary-layer method for calculating subsonic and transonic airfoil flow with trailing-edge separation

The development of a potential-flow/boundary-layer method for calculating subsonic and transonic turbulent flow past airfoils with trailing-edge separation is reported. A moment-of-momentum integral boundary-layer method is used which employs the law-of-the-wall/law-of-the-wake velocity profile and a two-layer eddy-viscosity model and ignores the laminar sublayer. All integrals across the boundary layer are obtained in closed form. Separation is assumed to occur when the shearing-stress velocity vanishes. A closed-form solution is derived for separated-flow regions where the shearing stress is negligible. In the potential-flow method, the exact form of the airfoil boundary condition is used, but it is applied at the chord line rather than the airfoil surface. This allows the accurate computation of flow about airfoils at large angles of attack but permits the use of body-oriented Cartesian computational grids. The governing equation for the perturbation velocity potential contains several terms in addition to the classical small-disturbance terms.

Barnwell, R. W.↗

Finite element analysis of transonic flows in highly staggered cascades

Finite element method is employed in the solution of full-potential equation for analyzing steady, transonic flows through a series of highly staggered cascades. Convergence characteristics and efficiency of the developed numerical procedure are discussed. Accuracy of the mathematical model is investigated in terms of the employed computational grid. In particular, the importance of accurate representation of the trailing edge together with the shock on the obtained flow results are studied. Numerical results include choked flows, for a cascade of NACA 0012 airfoils with pitch/chord length = 1, and stagger angle of 45 deg.

Akay, H. U.↗

Transonic flow computations in cascades using finite element method

Analysis of transonic flow through a cascade of airfoils is investigated using the finite element method. Development of a computational grid suitable for complex flow structures and different types of boundary conditions is presented. An efficient pseudo-time integration scheme is developed for the solution of equations. Modeling of the shock and the convergence characteristics of the developed scheme are discussed. Numerical results include a 45 deg staggered cascade of NACA 0012 airfoils with inlet flow Mach number of 0.8 and angles of attack 1, 0, and 1 deg.

Akay, H. U.↗

Computation of three-dimensional inviscid flow over hypersonic missile configurations using the GIM code

A three-dimensional computational technique was used to obtain flowfield solutions to the Euler equations over selected hypersonic missile configurations. The General Interpolants Method (GIM) computer code was used with interpolation functions in an algebraic approach to generate a discrete computational grid for each configuration. The spatial marching version of the GIM code, which treats the parabolized Navier-Stokes (PNS) equations or the Euler equations with a shock capturing, 'MacCormack-like' scheme, was used to advance the solution hyperbolically over each configuration. The inviscid flowfield solutions over the two three-dimensional missile configurations, calculated using the GIM hyperbolic scheme, are presented here. The flow field over a wing/body configuration at zero degree angle of attack is presented. Flow over the fuselage of a tactical missile, termed the TAME 10, at both zero degree and 7.5 degree angles of attack is presented. In addition, an inviscid, two-dimensional analysis of an inlet configuration designed to mount on the TAME 10 is included. Contour maps of velocity and pressure are included for each configuration. Comparison of calculation and data show good agreement.

Xiques, K. E.↗

Approximate boundary condition procedure for the two-dimensional numerical solution of vortex wakes

Research on efficient computational methods for general vorticity fields has been conducted in connection with a need for basic research on vortex-dominated flows. The present investigation is concerned with the evolution of vortex wakes behind aircraft wings. An efficient procedure is presented for the calculation of the boundary values used in the numerical solution of the unsteady, incompressible, two-dimensional Navier-Stokes equations for an unbounded flow field. The extent of the computational grid can be reduced compared to methods utilizing standard boundary conditions, without loss of accuracy. The efficiencies realized make it feasible to calculate the vortex wake development for realistic wing configurations, including the merging of multiple vortices, for Reynolds numbers of about 10,000 based on wing chord.

Weston, R. P.↗