Search NASA⌕ Search

Engineering topics

Molvik, G. A.

Publications and source records attributed to Molvik, G. A..

12 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

Global aspects of the motion of the heliospheric termination shock: A gasdynamic model

The heliospheric termination shock is expected to move in response to variation in upstream solar wind conditions. Using numerical techniques, we extend an earlier strictly one-dimensional analytic gas dynamic model of shock motion to two dimensions, to investigate the qualitative features of global behavior of shock motion, and the consequences of latitudinal variation in dynamic pressure. The boundary conditions of the calculation are given by the solar wind parameters as a function of latitude and time on an inner spherical boundary, and a constant pressure (roughly simulating the effect of the local interstellar medium) on an outer boundary. Density variations, specified at the inner boundary as a function of time, are convected into the termination shock. Immediately after the interaction, the shock moves with speeds given by the earlier analytic model. However, as the termination shock propagates outward (or inward), it begins to slow down. After about 2 to 10 years, depending on details of boundary conditions, the signal from the shock interaction has reached the outer boundary and propagates inward to the position of the termination shock, strongly affecting the behavior of the shock. Assuming no further disturbances in the solar wind, the termination shock will reach its new equilibrium after some tens of years. In reality, large-scale variations in solar wind dynamic pressure occur on time scales short in comparison with the eleven year solar cycle, so that one expects that the termination shock is never in an equilibrium position, but rather oscillates inward and outward; this oscillation will vary with heliographic latitude. The effects of a variety of types of solar wind disturbances are investigated and summarized.

Ratkiewicz, R.↗
NASA NTRS1995

Heliospheric Termination Shock Motion Due to Fluctuations in the Solar Wind Upstream Conditions: Spherically Symmetric Model

Large-scale fluctuations in the solar wind plasma upstream of the heliospheric termination shock (TS) will cause inward and outward motions of the shock. Using numerical techniques, we extend an earlier strictly one-dimensional (planar) analytic gas dynamic model to spherical symmetry to investigate the features of global behavior of shock motion. Our starting point is to establish a steady numerical solution of the gasdynamic equations describing the interaction between the solar wind and the interstellar medium. We then introduce disturbances of the solar wind dynamic pressure at an inner boundary, and follow the subsequent evolution of the system, especially the motion of the termination shock. Our model solves spherically symmetric gasdynamic equations as an initial-boundary value problem. The equations in conservative form are solved using a fully implicit Total Variation Diminishing (TVD) upwind scheme with Roe-type Riemann solver. Boundary conditions are given by the solar wind parameters on an inner spherical boundary, where they are allowed to vary with time for unsteady calculations, and by a constant pressure (roughly simulating the effect of the local interstellar medium) on an outer boundary. We find that immediately after the interaction, the shock moves with speeds given by the earlier analogous analytic models. However, as the termination shock propagates it begins to slow down, seeking a new equilibrium position. In addition, the disturbance transmitted through the TS, either a shock or rarefaction wave, will encounter the heliopause boundary and be reflected back. The reflected signal will encounter the TS, causing it to oscillate. The phenomenon may be repeated for a number of reflections, resulting in a "ringing" of the outer heliosphere.

Ratkiewicz, R.↗
NASA NTRS1995

Computer Simulation of Motion of Heliospheric Termination Shock

A time-dependent computer model is used to study a non-equilibrium structure of the heliosphere, resulting from the interaction between the solar wind and the interstellar plasma. An interaction of the heliospheric termination shock with various disturbances such as density or velocity jumps or interplanetary shocks in upstream solar wind as well as pressure jumps in local interstellar medium is investigated. Consequences of latitudinal variation in dynamic pressure are also discussed. The initial boundary conditions of unsteady calculations are given by the solar wind parameters as functions of time and/or latitude on an inner boundary, and/or by the local interstellar pressure as a function of time (roughly simulating the effect of the local interstellar medium) on an outer boundary.

Ratkiewicz, R. E.↗
NASA NTRS1994

Shock Excursion Due to Fluctuations in the Solar Wind Upstream Conditions

Large-scale fluctuations in the solar wind upstream of the termination shock will cause inward and outward motions of the shock. In earlier work, Barnes analyzed such motion by calculating of the response of a planar gasdynamic shock to upstream disturbances. We now generalize this analysis to the case of a spherically symmetric shock. Our procedure is first to solve numerically the set of gasdynamic equations describing the interaction between the solar wind and the interstellar medium to establish a dynamic equilibrium. The next step is to impose upstream fluctuations of the solar wind dynamical pressure on this equilibrium state at an inner boundary, and then to follow the subsequent shock motion.

Ratkiewicz, Romana E.↗
NASA NTRS1992

Two Algorithms For Hypersonic Computations

Two new algorithms developed containing several desirable features for computation of viscous, hypersonic flows about three-dimensional geometries. Employ upwind-differencing methods, having unsurpassed abiltiy to capture flow-field discontinuities without need for any user-specified smoothing terms. Particularly important in hypersonic computations where many complicated wave structures exist. Both algorithms based on finite-volume formulation to ensure schemes, including boundary conditions, fully conservative.

Molvik, G. A.↗
NASA NTRS1992

Conservative Grid-Interface Algorithm For Computing Flows

Best features of structured- and unstructured-grid methods combined. Gaps and overlaps between zonal grids eliminated by grid-interface algorithm, which generates single interfacial grid and corrects fluxes of flow quantities accordingly. Incorporated into two three-dimensional Navier-Stokes finite-volume codes and tested in computations of incompressible and compressible flows about simple bodies. Good numerical results obtained. General enough to be incorporated into other finite-volume codes without restrictions on complexities of shapes of bodies and zonal interfaces.

Klopfer, G. H.↗
NASA NTRS1991

Conservative multizonal interface algorithm for the 3-D Navier-Stokes equations

A conservative zonal interface algorithm using features of both structured and unstructured mesh CFD technology is presented. The flow solver within each of the zones is based on structured mesh CFD technology. The interface algorithm was implemented into two three-dimensional Navier-Stokes finite volume codes and was found to yield good results.

Klopfer, G. H.↗
NASA NTRS1989

Experimental and computational results for 5 degree blunt cones with shock generators at high velocity

Experiments and computations have been performed under laminar conditions in air on 5-deg blunt cones at velocities of 5 km/s and 6 km/s and at Reynolds numbers of 100,000 and 1 million. The computations were performed using ideal-, equilibrium- and nonequilibrium-chemistry models for air. At the conditions of the tests, the aerodynamic coefficients are sensitive to the real-gas effects present, and both experimental and computational aerodynamic coefficients show real-gas and nonlinear effects. The nonequilibrium computations show that a large amount of oxygen is dissociated in the blunt nose region of the flow and much of the oxygen remains dissociated over the entire length of the body, providing an insight into the source of the observed effects in the aerodynamic coefficients. The experimental and computational shock-shapes are in good agreement.

Strawa, A. W.↗
NASA NTRS1985

Numerical solution of strong radiation gasdynamic interactions in a hydrogen-seedant mixture

Theoretical and experimental results are provided for the absorption of laser energy by a gas flow. A two-dimensional numerical model described the temperature field, fluid dynamics, variable absorptivity and gas dynamics. The model was tested using a hydrogen-seedant mixture injected with 10.6 micron laser radiation. The laser mix was selected to allow reduced temperature gradients and energy source terms. It was found possible, by varying parameters such as the laser power and flash pressure, to precisely position the peak absorptance region in higher speed flows. More refined models are being tested to achieve accurate absorptance predictions in flows below Mach 0.0001.

Merkle, C. L.↗
NASA NTRS1984

A two-dimensional analysis of laser heat addition in converging nozzles

The two-dimensional equations of motion describing the interaction between a laser beam and a flowing gas are considered. An implicit numerical scheme is used to solve these equations for unchoked flow through a converging-diverging nozzle. Separate grids are used for the fluid dynamics and the radiation equations. The effects of beam focusing and cross-beam intensity profiles are included. The calculations are based upon real gas properties for all quantities except the gas absorptivity, which is taken as a constant. The solutions contain the expected hot central core region with cool gas near the walls. This results in steep temperature gradients in both the streamwise and cross-stream directions. The absorption zone acts as a blockage in the nozzle causing a nonuniform velocity profile at the inlet and an overall decrease in mass flow. The absorption region also forces the streamlines to move away from the axis of symmetry, although this effect is not strong.

Molvik, G. A.↗