Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical computation”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS1979

HF-induced parametric decay in presence of field-aligned irregularities

A nonlinear theory for the saturation of the HF-induced parametric decay instability is formulated in terms of the normal modes of a field-aligned plasma wave duct. Coupling coefficients between modes are derived, and a saturation spectrum is computed numerically for a typical irregularity. The computed spectrum shows that unstable waves may propagate at larger angles to the magnetic field than would be the case in a regular medium. The relevance of this work to the experiment of Muldrew and Showen (1977), concerning the height of the HF-induced plasma line at Arecibo, is discussed. A possible explanation is also found for the 'decay line' feature in comparison with the 'broad bump' in the observed backscatter spectrum.

Rypdal, K.↗
NASA NTRS1985

Energy Stability of Thermocapillary Convection in Models of the Float Zone Process

The energy-stability of thermocapillary convection in models of the float-zone, crystal-growing process was studied. Stability limits, as functions of pertinent parameters, that will identify conditions which will not allow the existence of an undesirable oscillatory flow instability were determined. Such instabilities may occur in the space processing of semiconductor materials. The determination of the stability limits will involve two sets of numerical computations: (1) solution of the nonlinear governing equations together with the appropriate boundary conditions to determine the basic state (in general, velocity, pressure and temperature fields and the displacement of free surfaces and interfaces); and (2) solution of a nonlinear Euler-Lagrange systems for the energy-stability limit. Both computations, while difficult, should be within the scope of available computer capability and available concepts in numerical analysis. Finite-element methods are attractive candidates for the numerical work.

Neitzel, G. P.↗
NASA NTRS1988

Shape optimization of 3D continuum structures via force approximation techniques

The existing need to develop methods whereby the shape design efficiency can be improved through the use of high quality approximation methods is addressed. An efficient approximation method for stress constraints in 3D shape design problems is proposed based on expanding the nodal forces in Taylor series with respect to shape variations. The significance of this new method is shown through elementary beam theory calculations and via numerical computations using 3D solid finite elements. Numerical examples including the classical cantilever beam structure and realistic automotive parts like the engine connecting rod are designed for optimum shape using the proposed method. The numerical results obtained from these methods are compared with other published results, to assess the efficiency and the convergence rate of the proposed method.

Vanderplaats, Garret N.↗
NASA NTRS1980

Three-layer interactive method for computing supersonic laminar separated flows

An interactive model for numerical computation of complicated two-dimensional flowfields including regions of reversed flow is proposed. The present approach is one of dividing the flowfield into three regions, in each of which a simplified mathematical model is applied: (1) outer, supersonic flow for which the full potential equation (hyperbolic) is used; (2) viscous, laminar layer in which the compressible boundary-layer model (parabolic) is used; and (3) recirculating flow modeled by the incompressible Navier-Stokes equations (elliptic). For matching of the numerical solutions in the three layers, two interaction models are developed: one for pressure interaction, the other for interaction between the shear layer and the recirculating flow. The uniform solution for the whole flowfield is then obtained by iteration of the local solutions under the constraints imposed by matching. The three-layer interactive model is used for solution of the flowfield past an asymmetric cavity. The method is shown to be capable of dealing with backflow without encountering problems at separation, characteristic to the boundary-layer approach.

Brandeis, J.↗
NASA NTRS1992

Study of Fluid Experiment System (FES)/CAST/Holographic Ground System (HGS)

The use of holographic and schlieren optical techniques for studying the concentration gradients in solidification processes has been used by several investigators over the years. The HGS facility at MSFC has been primary resource in researching this capability. Consequently, scientific personnel have been able to utilize these techniques in both ground based research and in space experiments. An important event in the scientific utilization of the HGS facilities was the TGS Crystal Growth and the casting and solidification technology (CAST) experiments that were flown on the International Microgravity Laboratory (IML) mission in March of this year. The preparation and processing of these space observations are the primary experiments reported in this work. This project provides some ground-based studies to optimize on the holographic techniques used to acquire information about the crystal growth processes flown on IML. Since the ground-based studies will be compared with the space-based experimental results, it is necessary to conduct sufficient ground based studies to best determine how the experiment worked in space. The current capabilities in computer based systems for image processing and numerical computation have certainly assisted in those efforts. As anticipated, this study has certainly shown that these advanced computing capabilities are helpful in the data analysis of such experiments.

Workman, Gary L.↗
NASA NTRS1992

Study of FES/CAST/HGS

The microgravity materials processing program has been instrumental in providing the crystal growth community with an experimental environment to better understand the phenomena associated with the growing of crystals. In many applications one may pursue the growth of large single crystals which cannot be grown on earth due to convective driven flows. A microgravity environment is characterized by neither convection of buoyancy. Consequently superior crystals are able to be grown in space. On the other hand, since neither convection nor buoyancy dominates the fluid flow in a microgravity environment, then lesser dominating phenomena can affect crystal growth, such as surface driven flows or diffusion limited solidification. In the case of experiments that are to be flown in space using the Fluid Experiments System (FES), diffusion limited growth should be the dominating phenomenon. The use of holographic and Schlieren optical techniques for studying the concentration gradients in solidification processes has been used by several investigators over the years. The Holographic Ground System (HGS) facility at MSFC has been a primary resource in researching this capability. Consequently scientific personnel have been able to utilize these techniques in both ground based research and in space experiments. An important event in the scientific utilization of the HGS facilities was the TGS (triglycine sulfate) Crystal Growth and the Casting and Solidification Technology (CAST) experiments that were flown on the International Microgravity Lab (IML) mission in March of this year. The preparation and processing of these space observations are the primary experiments reported in this work. This project provides some ground-based studies to optimize on the holographic techniques used to acquire information about the crystal growth processes flown on IML. Since the ground-based studies will be compared with the space-based experimental results, it is necessary to conduct sufficient ground based studies to best determine how the experiment in space worked. The current capabilities in computer based systems for image processing and numerical computation have certainly assisted in those efforts. As anticipated, this study has certainly shown that these advanced computing capabilities are helpful in the data analysis of such experiments.

Workman, Gary L.↗
NASA NTRS2013

Influence of Computational Drop Representation in LES of a Droplet-Laden Mixing Layer

Multiphase turbulent flows are encountered in many practical applications including turbine engines or natural phenomena involving particle dispersion. Numerical computations of multiphase turbulent flows are important because they provide a cheaper alternative to performing experiments during an engine design process or because they can provide predictions of pollutant dispersion, etc. Two-phase flows contain millions and sometimes billions of particles. For flows with volumetrically dilute particle loading, the most accurate method of numerically simulating the flow is based on direct numerical simulation (DNS) of the governing equations in which all scales of the flow including the small scales that are responsible for the overwhelming amount of dissipation are resolved. DNS, however, requires high computational cost and cannot be used in engineering design applications where iterations among several design conditions are necessary. Because of high computational cost, numerical simulations of such flows cannot track all these drops. The objective of this work is to quantify the influence of the number of computational drops and grid spacing on the accuracy of predicted flow statistics, and to possibly identify the minimum number, or, if not possible, the optimal number of computational drops that provide minimal error in flow prediction. For this purpose, several Large Eddy Simulation (LES) of a mixing layer with evaporating drops have been performed by using coarse, medium, and fine grid spacings and computational drops, rather than physical drops. To define computational drops, an integer NR is introduced that represents the ratio of the number of existing physical drops to the desired number of computational drops; for example, if NR=8, this means that a computational drop represents 8 physical drops in the flow field. The desired number of computational drops is determined by the available computational resources; the larger NR is, the less computationally intensive is the simulation. A set of first order and second order flow statistics, and of drop statistics are extracted from LES predictions and are compared to results obtained by filtering a DNS database. First order statistics such as Favre averaged stream-wise velocity, Favre averaged vapor mass fraction, and the drop stream-wise velocity, are predicted accurately independent of the number of computational drops and grid spacing. Second order flow statistics depend both on the number of computational drops and on grid spacing. The scalar variance and turbulent vapor flux are predicted accurately by the fine mesh LES only when NR is less than 32, and by the coarse mesh LES reasonably accurately for all NR values. This is attributed to the fact that when the grid spacing is coarsened, the number of drops in a computational cell must not be significantly lower than that in the DNS.

Bellan, Josette↗
NASA NTRS1972

Fluid dynamics computer programs for NERVA turbopump

During the design of the NERVA turbopump, numerous computer programs were developed for the analyses of fluid dynamic problems within the machine. Program descriptions, example cases, users instructions, and listings for the majority of these programs are presented.

Brunner, J. J.↗
NASA NTRS2011

Porting a Hall MHD Code to a Graphic Processing Unit

We present our experience porting a Hall MHD code to a Graphics Processing Unit (GPU). The code is a 2nd order accurate MUSCL-Hancock scheme which makes use of an HLL Riemann solver to compute numerical fluxes and second-order finite differences to compute the Hall contribution to the electric field. The divergence of the magnetic field is controlled with Dedner?s hyperbolic divergence cleaning method. Preliminary benchmark tests indicate a speedup (relative to a single Nehalem core) of 58x for a double precision calculation. We discuss scaling issues which arise when distributing work across multiple GPUs in a CPU-GPU cluster.

Dorelli, John C.↗
NASA NTRS2019

Dynamic Smagorinsky Modeled Large-Eddy Simulations of Turbulence Using Tetrahedral Meshes

Eddy-resolving numerical computations of turbulent flows are emerging as viable alternatives to Reynolds Averaged Navier-Stokes (RANS) calculations for flows with an intrinsically steady mean state due to the advances in large-scale parallel computing. In these computations, medium to large turbulent eddies are resolved by the numerics while the smaller or subgrid scales are either modeled or taken care of by the inherent numerical dissipation. To advance the state of the art of unstructured-mesh turbulence simulation capabilities, large eddy simulations (LES) using the dynamic Smagorinsky model (DSM) on tetrahedral meshes are carried out with the space-time conservation element, solution element (CESE) method. In contrast to what has been reported in the literature, the present implementation of dynamic models allows for active backscattering without any ad-hoc limiting of the eddy viscosity calculated from the subgrid-scale model. For the benchmark problems involving compressible isotropic turbulence decay as well as the shock/turbulent boundary layer interaction benchmark problems, no numerical instability associated with kinetic energy growth is observed and the volume percentage of the backscattering portion accounts for about 38-40% of the simulation domain. A slip-wall model in conjunction with the implemented DSM is used to simulate a relatively high Reynolds number Mach 2.85 turbulent boundary layer over a 30° ramp with several tetrahedral meshes and a wall-normal spacing of either Δγ& = 10 or Δγ& = 20. The computed mean wall pressure distribution, separation region size, mean velocity profiles, and Reynolds stress agree reasonably well with experimental data.

CESE method↗
NASA NTRS2019

Converging-Diverging Nozzle Thruster Code for Nuclear or Chemical Rocket Performance Computations

This report describes development of a variable area ratio C-D nozzle numerical computational code "THRST", applicable to the evaluation of nuclear or chemical rocket performance. The propellant for the nuclear rockets is para-hydrogen with gas properties supplied by a sub-routine, while properties for chemical thrusters need to be user input. Detailed documentation of input and output variables and sample output listings enable the report to function as a "User Manual".

Juhasz, Albert J.↗
NASA NTRS1973

Dynamic loads due to moving liquid.

The marker-and-cell method has been extended to simulate the flow and velocity fields of a moving liquid in an axisymmetric container. Dynamic loads induced by the moving liquid are calculated from the pressure distribution along the container wall as functions of time. Using viscosity and equivalent gravitational acceleration as parameters, several cases of the transient flow of a liquid in a container are studied. Subsequently, a three-dimensional marker-and-cell formulation is obtained and a numerical scheme to implement the formulation is developed. The resulting computer program is used to simulate the three-dimensional flow of a liquid in a rectangular container. The study shows that dynamic loads induced by a moving liquid can be computed numerically.

Feng, G. C.↗
NASA NTRS1989

Overview of optimal trajectories for flight in a windshear

Optimal flight trajectories for the B-727, B-737, and B-747 aircraft in the presence of wind shear are studied. The takeoff problem and the abort landing problem are considered with reference to flight in a vertical plane. In the former, optimal trajectories are computed by minimizing the peak deviation of the absolute path inclination from a reference value; in the latter, optimal trajectories are computed by minimizing the peak value of the altitude drop. Numerical computations show that, for both the problems under consideration, the optimal trajectories of the three aircraft show the same qualitative behavior. Hence, it appears that the near-optimal guidance schemes developed for the B-727 can be extended to the other two aircraft, albeit with some quantitative modification.

Miele, A.↗
NASA NTRS1988

A numerical method for computing unsteady 2-D boundary layer flows

A numerical method for computing unsteady two-dimensional boundary layers in incompressible laminar and turbulent flows is described and applied to a single airfoil changing its incidence angle in time. The solution procedure adopts a first order panel method with a simple wake model to solve for the inviscid part of the flow, and an implicit finite difference method for the viscous part of the flow. Both procedures integrate in time in a step-by-step fashion, in the course of which each step involves the solution of the elliptic Laplace equation and the solution of the parabolic boundary layer equations. The Reynolds shear stress term of the boundary layer equations is modeled by an algebraic eddy viscosity closure. The location of transition is predicted by an empirical data correlation originating from Michel. Since transition and turbulence modeling are key factors in the prediction of viscous flows, their accuracy will be of dominant influence to the overall results.

Krainer, Andreas↗
NASA NTRS1991

Computation of turbulent flow in a thin liquid layer of fluid involving a hydraulic jump

Numerically computed flow fields and free surface height distributions are presented for the flow of a thin layer of liquid adjacent to a solid horizontal surface that encounters a hydraulic jump. Two kinds of flow configurations are considered: two-dimensional plane flow and axisymmetric radial flow. The computations used a boundary-fitted moving grid method with a k-epsilon model for the closure of turbulence. The free surface height was determined by an optimization procedure which minimized the error in the pressure distribution on the free surface. It was also checked against an approximate procedure involving integration of the governing equations and use of the MacCormack predictor-corrector method. The computed film height also compared reasonably well with previous experiments. A region of recirculating flow was found to be present adjacent to the solid boundary near the location of the jump, which was caused by a rapid deceleration of the flow.

Rahman, M. M.↗
NASA NTRS1990

AI/OR computational model for integrating qualitative and quantitative design methods

A theoretical framework for integrating qualitative and numerical computational methods for optimally-directed design is described. The theory is presented as a computational model and features of implementations are summarized where appropriate. To demonstrate the versatility of the methodology we focus on four seemingly disparate aspects of the design process and their interaction: (1) conceptual design, (2) qualitative optimal design, (3) design innovation, and (4) numerical global optimization.

Agogino, Alice M.↗