Search NASA⌕ Search

SEARCH · Search NASA

Results for “explicit solution”

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 883 records · Page 49

Explicit large time-step schemes for the shallow water equations

Modifications to explicit finite difference schemes for solving the shallow water equations for meteorological applications by increasing the time step for the fast gravity waves are analyzed. Terms associated with the gravity waves in the shallow water equations are treated on a coarser grid than those associated with the slow Rossby waves, which contain much more of the available energy and must be treated with higher accuracy, enabling a several-fold increase in time step without degrading the accuracy of the solution. The method is presented in Cartesian and spherical coordinates for a rotating earth, using generalized leapfrog, frozen coefficient, and Fourier filtering finite difference schemes. Computational results verify the numerical stability of the approach.

Turkel, E.↗

A multistage time-stepping scheme for the Navier-Stokes equations

A class of explicit multistage time-stepping schemes is used to construct an algorithm for solving the compressible Navier-Stokes equations. Flexibility in treating arbitrary geometries is obtained with a finite-volume formulation. Numerical efficiency is achieved by employing techniques for accelerating convergence to steady state. Computer processing is enhanced through vectorization of the algorithm. The scheme is evaluated by solving laminar and turbulent flows over a flat plate and an NACA 0012 airfoil. Numerical results are compared with theoretical solutions or other numerical solutions and/or experimental data.

Swanson, R. C.↗

A multistage time-stepping scheme for the Navier-Stokes equations

A class of explicit multistage time-stepping schemes is used to construct an algorithm for solving the compressible Navier-Stokes equations. Flexibility in treating arbitrary geometries is obtained with a finite-volume formulation. Numerical efficiency is achieved by employing techniques for accelerating convergence to steady state. Computer processing is enhanced through vectorization of the algorithm. The scheme is evaluated by solving laminar and turbulent flows over a flat plate and an NACA 0012 airfoil. Numerical results are compared with theoretical solutions or other numerical solutions and/or experimental data.

Swanson, R. C.↗

A study of the effects of numerical dissipation on the calculation of supersonic separated flows

An extensive investigation of the effect of numerical dissipation on the calculation of supersonic, separated flow over a rearward-facing step is carried out. The complete two-dimensional Navier-Stokes equations are solved by means of MacCormack's standard explicit, unsplit, time-dependent, finite difference method. A fourth-order numerical dissipation term is added explicitly. The magnitude of this term is progressively varied, and its consequences on the flowfield calculations are identified and studied. For a cold-wall, heat transfer case, numerical dissipation had a major effect on the results, particularly in the separated region. However, rather dramatically for an adiabatic wall case, numerical dissipation had virtually no effect on the results. The role of grid size on both the influence of numerical dissipation, and on the overall accuracy of the separated flow solutions is discussed.

Kuruvila, G.↗

High-speed compressible flow and other advection-dominated problems of fluid dynamics

Finite element methods are described for modeling high speed compressible flows with strong advection, problems important to aerodynamics. The situations are characterized by high pressure and temperature gradients, transients and the appearance of discontinuities, factors which require mesh refinement during computations. Techniques are developed for temporal and spatial discretization of a model problem. Several observations are made regarding the explicit and implicit features of the calculations, the use of the Lax-Wendroff scheme to produce a mass-matrix for obtaining accurate results for transients, methods of performing stability analyses, and simplification techniques. Examples are provided of solving the nonlinear shallow-water equations and describing compressible flows, particularly transonic flows. Domain splitting is defined for improving the calculations at each time step and in different parts of the flow regime while simultaneously advancing the calculations towards a solution.

Zienkiewicz, O. C.↗

Application of Krylov exponential propagation to fluid dynamics equations

This paper presents an application of matrix exponentiation via Krylov subspace projection, to the solution of fluid dynamics problems. The main idea is to approximate the operation exp(A)v by means of a projection-like process onto a Krylov subspace. This results in a computation of an exponential matrix vector product similar to the one above but of a much smaller size. Time integration schemes can then be devised to exploit this basic computational kernel. The motivation of this approach is to provide time-integration schemes that are essentially of an explicit nature but which have good stability properties.

Saad, Y.↗

Parallel Methods on Large-Scale Structural Analysis and Physics Applications; Symposium, Hampton, VA, Feb. 5, 6, 1991, Selected Papers

Recent advances in parallel methods and algorithms integrated into large-scale codes are presented. Consideration is given to problem decomposition (substructuring), efficient matrix solution algorithms for shared memory architectures, dynamic and transient analysis algorithms for shared memory architectures, and algorithms for distributed and massively parallel architectures. Particular attention is given to partitioning of unstructured problems for parallel processing, parallel-vector computation for linear-structural analysis and nonlinear unconstraint optimization problems, a parallel-vector equation solver for unsymmetric matrices on supercomputers, parallel nonlinear finite element dynamic response, multigrid algorithms for solving structural mechanics problems on supercomputers, structural analysis on massively parallel computers, explicit finite element methods with contact-impact on SIMD computers, and the impact of mapping and sparsity on parallelized finite element method modules.

Storaasli, Olaf O.↗

A well-posed numerical method to track isolated conformal map singularities in Hele-Shaw flow

We present a new numerical method for calculating an evolving 2D Hele-Shaw interface when surface tension effects are neglected. In the case where the flow is directed from the less viscous fluid into the more viscous fluid, the motion of the interface is ill-posed; small deviations in the initial condition will produce significant changes in the ensuing motion. This situation is disastrous for numerical computation, as small round-off errors can quickly lead to large inaccuracies in the computed solution. Our method of computation is most easily formulated using a conformal map from the fluid domain into a unit disk. The method relies on analytically continuing the initial data and equations of motion into the region exterior to the disk, where the evolution problem becomes well-posed. The equations are then numerically solved in the extended domain. The presence of singularities in the conformal map outside of the disk introduces specific structures along the fluid interface. Our method can explicitly track the location of isolated pole and branch point singularities, allowing us to draw connections between the development of interfacial patterns and the motion of singularities as they approach the unit disk. In particular, we are able to relate physical features such as finger shape, side-branch formation, and competition between fingers to the nature and location of the singularities. The usefulness of this method in studying the formation of topological singularities (self-intersections of the interface) is also pointed out.

Baker, Gregory↗

Numerical solution of two-dimensional turbulent blunt body flows with an impinging shock

An implicit finite-difference method has been developed to compute two-dimensional, turbulent, blunt body flows with an impinging shock wave. The full time-averaged Navier-Stokes equations are solved with algebraic eddy viscosity and turbulent Prandtl number models employed for shear stress and heat flux. The irregular-shaped bow shock is treated as a discontinuity across which the Rankine-Hugoniot equations are applied. A Type III turbulent shock interference flow field has been computed and the numerical results compare favorably with existing experimental data. In addition, comparisons are made between the present implicit code and a previous explicit code.

Tannehill, J. C.↗

Application of Krylov exponential propagation to fluid dynamics equations

An application of matrix exponentiation via Krylov subspace projection to the solution of fluid dynamics problems is presented. The main idea is to approximate the operation exp(A)v by means of a projection-like process onto a krylov subspace. This results in a computation of an exponential matrix vector product similar to the one above but of a much smaller size. Time integration schemes can then be devised to exploit this basic computational kernel. The motivation of this approach is to provide time-integration schemes that are essentially of an explicit nature but which have good stability properties.

Saad, Youcef↗

Preliminary Work for Examining the Scalability of Reinforcement Learning

Researchers began studying automated agents that learn to perform multiple-step tasks early in the history of artificial intelligence (Samuel, 1963; Samuel, 1967; Waterman, 1970; Fikes, Hart & Nilsonn, 1972). Multiple-step tasks are tasks that can only be solved via a sequence of decisions, such as control problems, robotics problems, classic problem-solving, and game-playing. The objective of agents attempting to learn such tasks is to use the resources they have available in order to become more proficient at the tasks. In particular, each agent attempts to develop a good policy, a mapping from states to actions, that allows it to select actions that optimize a measure of its performance on the task; for example, reducing the number of steps necessary to complete the task successfully. Our study focuses on reinforcement learning, a set of learning techniques where the learner performs trial-and-error experiments in the task and adapts its policy based on the outcome of those experiments. Much of the work in reinforcement learning has focused on a particular, simple representation, where every problem state is represented explicitly in a table, and associated with each state are the actions that can be chosen in that state. A major advantage of this table lookup representation is that one can prove that certain reinforcement learning techniques will develop an optimal policy for the current task. The drawback is that the representation limits the application of reinforcement learning to multiple-step tasks with relatively small state-spaces. There has been a little theoretical work that proves that convergence to optimal solutions can be obtained when using generalization structures, but the structures are quite simple. The theory says little about complex structures, such as multi-layer, feedforward artificial neural networks (Rumelhart & McClelland, 1986), but empirical results indicate that the use of reinforcement learning with such structures is promising. These empirical results make no theoretical claims, nor compare the policies produced to optimal policies. A goal of our work is to be able to make the comparison between an optimal policy and one stored in an artificial neural network. A difficulty of performing such a study is finding a multiple-step task that is small enough that one can find an optimal policy using table lookup, yet large enough that, for practical purposes, an artificial neural network is really required. We have identified a limited form of the game OTHELLO as satisfying these requirements. The work we report here is in the very preliminary stages of research, but this paper provides background for the problem being studied and a description of our initial approach to examining the problem. In the remainder of this paper, we first describe reinforcement learning in more detail. Next, we present the game OTHELLO. Finally we argue that a restricted form of the game meets the requirements of our study, and describe our preliminary approach to finding an optimal solution to the problem.

Clouse, Jeff↗

A comparison of upwind schemes for computation of three-dimensional hypersonic real-gas flows

The method of Suresh and Liou (1992) is extended, and the resulting explicit noniterative upwind finite-volume algorithm is applied to the integration of 3D parabolized Navier-Stokes equations to model 3D hypersonic real-gas flowfields. The solver is second-order accurate in the marching direction and employs flux-limiters to make the algorithm second-order accurate, with total variation diminishing in the cross-flow direction. The algorithm is used to compute hypersonic flow over a yawed cone and over the Ames All-Body Hypersonic Vehicle. The solutions obtained agree well with other computational results and with experimental data.

Gerbsch, R. A.↗

Relativistic Navigation: A Theoretical Foundation

We present a theoretical foundation for relativistic astronomical measurements in curved space-time. In particular, we discuss a new iterative approach for describing the dynamics of an isolated astronomical N-body system in metric theories of gravity. To do this, we generalize the Fock-Chandrasekhar method of the weak-field and slow-motion approximation (WFSMA) and develop a theory of relativistic reference frames (RF's) for a gravitationally bounded many-extended-body problem. In any proper RF constructed in the immediate vicinity of an arbitrary body, the N-body solutions of the gravitational field equations are formally presented as a sum of the Riemann-flat inertial space-time, the gravitational field generated by the body itself, the unperturbed solutions for each body in the system transformed to the coordinates of this proper RF, and the gravitational interaction term. We develop the basic concept of a general WFSMA theory of the celestial RF's applicable to a wide class of metric theories of gravity and an arbitrary model of matter distribution. We apply the proposed method to general relativity. Celestial bodies are described using a perfect fluid model; as such, they possess any number of internal mass and current multipole moments that explicitly characterize their internal structures. The obtained relativistic corrections to the geodetic equations of motion arise because of a coupling of the bodies' multiple moments to the surrounding gravitational field. The resulting relativistic transformations between the different RF's extend the Poincare group to the motion of deformable self-gravitating bodies. Within the present accuracy of astronomical measurements we discuss the properties of the Fermi-normal-like proper RF that is defined in the immediate vicinity of the extended compact bodies. We further generalize the proposed approximation method and include two Eddington parameters (gamma, Beta). This generalized approach was used to derive the relativistic equations of satellite motion in the vicinity of the extended bodies. Anticipating improvements in radio and laser tracking technologies over the next few decades, we apply this method to spacecraft orbit determination. We emphasize the number of feasible relativistic gravity tests that may be performed within the context of the parameterized WFSMA. Based on the planeto-centric equations of motion of a spacecraft around the planet, we suggested a new null test of the Strong Equivalence Principle (SEP). The experiment to measure the corresponding SEP violation effect could be performed with the future Mercury Orbiter mission. We discuss other relativistic effects, including the perihelion advance and the redshift and geodetic precession of the orbiter's orbital plane about Mercury, as well as the possible future implementation of the proposed formalism in software codes developed for solar-system orbit determination. All the important calculations are completely documented, and the references contain an extensive list of cited literature.

Turyshev, Slava G.↗

Micromechanics analysis of space simulated thermal deformations and stresses in continuous fiber reinforced composites

Space simulated thermally induced deformations and stresses in continuous fiber reinforced composites were investigated with a micromechanics analysis. The investigation focused on two primary areas. First, available explicit expressions for predicting the effective coefficients of thermal expansion (CTEs) for a composite were compared with each other, and with a finite element (FE) analysis, developed specifically for this study. Analytical comparisons were made for a wide range of fiber/matrix systems, and predicted values were compared with experimental data. The second area of investigation focused on the determination of thermally induced stress fields in the individual constituents. Stresses predicted from the FE analysis were compared to those predicted from a closed-form solution to the composite cylinder (CC) model, for two carbon fiber/epoxy composites. A global-local formulation, combining laminated plate theory and FE analysis, was used to determine the stresses in multidirectional laminates. Thermally induced damage initiation predictions were also made.

Bowles, David E.↗

Computations of Internal and External Axisymmetric Nozzle Aerodynamics at Transonic Speeds

Computational Fluid Dynamics (CFD) analyses of axisymmetric circular-arc boattail nozzles have been completed in support of NASA's Next Generation Launch Technology Program to investigate the effects of high-speed nozzle geometries on the nozzle internal flow and the surrounding boattail regions. These computations span the very difficult transonic flight regime, with shock-induced separations and strong adverse pressure gradients. External afterbody and internal nozzle pressure distributions computed with the Wind code are compared with experimental data. A range of turbulence models were examined in Wind, including an Explicit Algebraic Stress model (EASM). Computations on two nozzle geometries have been completed at freestream Mach numbers ranging from 0.6 to 0.9, driven by nozzle pressure ratios (NPR) ranging from 2.9 to 5. Results obtained on converging-only geometry indicate reasonable agreement to experimental data, with the EASM and Shear Stress Transport (SST) turbulence models providing the best agreement. Calculations completed on a converging-diverging geometry involving large-scale internal flow separation did not converge to a true steady-state solution when run with variable timestepping (steady-state). Calculations obtained using constant timestepping (time-accurate) indicate less variations in flow properties compared with steady-state solutions. This failure to converge to a steady-state solution was found to be the result of difficulties in using variable time-stepping with large-scale separations present in the flow. Nevertheless, time-averaged boattail surface pressure coefficient and internal nozzle pressures show fairly good agreement with experimental data. The SST turbulence model demonstrates the best over-all agreement with experimental data.

Dalbello, Teryn↗

The Effect of the Density Ratio on the Nonlinear Dynamics of the Unstable Fluid Interface

Here we report multiple harmonic theoretical solutions for a complete system of conservation laws, which describe the large-scale coherent dynamics in RTI and RMI for fluids with a finite density ratio in the general three-dimensional case. The analysis yields new properties of the bubble front dynamics. In either RTI or RMI, the obtained dependencies of the bubble velocity and curvature on the density ratio differ qualitatively and quantitatively from those suggested by the models of Sharp (1984), Oron et al. (2001), and Goncharov (2002). We show explicitly that these models violate the conservation laws. For the first time, our theory reveals an important qualitative distinction between the dynamics of the RT and RM bubbles.

Abarzhi, S. I.↗

Dispersion of ensembles of non-interacting particles

The dynamics of an ensemble of noninteracting particles dispersing from a common origin and moving in a common force field with an initial distribution of momenta is analyzed using an approach where the particles are considered as a continuum described by a phase-space distribution function. General solutions are obtained for both the distribution function and the associated spatial density function. The linear case of small departures from circular orbits in an axisymmetric gravitational field is treated along with the specific case of particle dispersion from an object in a circular orbit in the same type of field. Numerical results are presented for the latter case, and consideration is given to the inverse problem of determining the initial time and velocity distribution from knowledge of the ensemble structure at a later time. Explicit results are provided for the case of an ellipsoidal distribution of initial momenta, and a numerical procedure is indicated for treating more general cases.

Heard, W. B.↗

Statistical mechanics of dilute liquid mixtures of He3 in He4

The Woo-Tan-Massey quantum-mechanical theory of dilute mixtures of He-3 in liquid He-4 is extended to incorporate superfluid flow. Then statistical-mechanical principles are employed to find explicit formulas for thermodyanmic and hydrodynamic properties of mixtures, such as internal energy, free energy, entropy, pressure, chemical potential, and normal-fluid and superfluid densities. This paper is the first in a planned series that applies methods originally developed for pure liquid He-4 to dilute mixtures as a means of establishing a comprehensive theory that includes a theoretical bridge between a microscopic description of the mixtures and a hydrodynamic two-fluid model. Ultimately effects of applied electric and magnetic fields on dilute solutions are also planned to be treated.

Jackson, H. W.↗