Search NASA⌕ Search

SEARCH · Search NASA

Results for “numerics”

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
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Numerical modeling for an electric-field hyperthermia applicator

Hyperthermia, in conjunction with radiation and chemotherapy for treatment of cancers, is an area of current concern. Experiments have shown that hyperthermia can increase the potency of many chemotherapy drugs and the effectiveness of radiation for treating cancer. A combination of whole body or regional hyperthermia with chemotherapy or radiation should improve treatment results. Conventional methods for inducing whole body hyperthermia, such as exposing a patient in a radiant cabinet or under a hot water blanket, conduct heat very slowly from the skin to the body core. Thus a more efficient system, such as the three-plate electric-field hyperthermia applicator (EHA), is developed. This three-plate EHA has one top plate over and two lower plates beneath the patient. It is driven at 27.12 MHz with 500 Watts through a matching circuit. Using this applicator, a 50 kg pig was successfully heated to 42 C within 45 minutes. However, phantom and animal studies have indicated non-uniform heating near the side of the body. In addition, changes in the size and distance between the electrode plates can affect the heating (or electromagnetic field) pattern. Therefore, numerical models using the method of moments (MOM) or the finite difference time domain (FDTD) technique are developed to optimize the heating pattern of this EHA before it is used for human trials. The accuracy of the numerical modeling has been achieved by the good agreement between the MOM and FDTD results for the three-plate EHA without a biological body. The versatile FDTD technique is then applied to optimize the EHA design with a human body. Both the numerical and measured data in phantom blocks will be presented. The results of this study will be used to design an optimized system for whole body or regional hyperthermia.

Wu, Te-Kao↗
NASA NTRS1993

The use of transition region characteristics to improve the numerical simulation of heat transfer in bypass transitional flows

A method is presented for improving the numerical prediction of bypass transition heat transfer on a flat plate in a high-disturbance environment with zero or favorable pressure gradient. The method utilizes low Reynolds number k-epsilon turbulence models in combination with the characteristic parameters of the transition region. The parameters representing the characteristics of the transition region used are the intermittency, transition length and turbulent spot properties. An analysis is made of the transition length in terms of turbulent spot variables. The nondimensional spot formation rate, required for the prediction of the transition length, is shown by the analysis to be a function of the spot spreading angle, the dimensionless spot velocity ratio and the dimensionless spot area ratio. The intermittency form of the k-epsilon equations were derived from conditionally averaged equations which have been shown to be an improvement over global-time-averaged equations for the numerical calculation of the transition region. The numerical predictions are in general good agreement with the experimental data and indicate the potential use of the method in accelerating flows. Turbulence models of the k-epsilon type are known to underpredict the transition length. The present work demonstrates how incorporating transition region characteristics improves the ability of two-equation turbulence models to simulate bypass transition for flat plates with potential application to turbine vanes and blades.

Simon, Frederick F.↗
NASA NTRS1994

A new flux-conserving numerical scheme for the steady, incompressible Navier-Stokes equations

This paper is concerned with the continued development of a new numerical method, the space-time solution element (STS) method, for solving conservation laws. The present work focuses on the two-dimensional, steady, incompressible Navier-Stokes equations. Using first an integral approach, and then a differential approach, the discrete flux conservation equations presented in a recent paper are rederived. Here a simpler method for determining the flux expressions at cell interfaces is given; a systematic and rigorous derivation of the conditions used to simulate the differential form of the governing conservation law(s) is provided; necessary and sufficient conditions for a discrete approximation to satisfy a conservation law in E2 are derived; and an estimate of the local truncation error is given. A specific scheme is then constructed for the solution of the thin airfoil boundary layer problem. Numerical results are presented which demonstrate the ability of the scheme to accurately resolve the developing boundary layer and wake regions using grids which are much coarser than those employed by other numerical methods. It is shown that ten cells in the cross-stream direction are sufficient to accurately resolve the developing airfoil boundary layer.

Scott, James R.↗
NASA NTRS1994

Numerical solution of the two-dimensional time-dependent incompressible Euler equations

A numerical method is presented for solving the artificial compressibility form of the 2D time-dependent incompressible Euler equations. The approach is based on using an approximate Riemann solver for the cell face numerical flux of a finite volume discretization. Characteristic variable boundary conditions are developed and presented for all boundaries and in-flow out-flow situations. The system of algebraic equations is solved using the discretized Newton-relaxation (DNR) implicit method. Numerical results are presented for both steady and unsteady flow.

Whitfield, David L.↗
NASA NTRS1993

Numerical investigation of multi-element airfoils

The flow over multi-element airfoils with flat-plate lift-enhancing tabs was numerically investigated. Tabs ranging in height from 0.25 percent to 1.25 percent of the reference airfoil chord were studied near the trailing edge of the main-element. This two-dimensional numerical simulation employed an incompressible Navier-Stokes solver on a structured, embedded grid topology. New grid refinements were used to improve the accuracy of the solution near the overlapping grid boundaries. The effects of various tabs were studied at a constant Reynolds number on a two-element airfoil with a slotted flap. Both computed and measured results indicated that a tab in the main-element cove improved the maximum lift and lift-to-drag ratio relative to the baseline airfoil without a tab. Computed streamlines revealed that the additional turning caused by the tab may reduce the amount of separated flow on the flap. A three-element airfoil was also studied over a range of Reynolds numbers. For the optimized flap rigging, the computed and measured Reynolds number effects were similar. When the flap was moved from the optimum position, numerical results indicated that a tab may help to reoptimize the airfoil to within 1 percent of the optimum flap case.

Cummings, Russell M.↗
NASA NTRS1992

Lumley's PODT definition of large eddies and a trio of numerical procedures

Lumley's 1967 Moscow paper provided, for the first time, a completely rational definition of the physically-useful term 'large eddy', popular for a half-century. The numerical procedures based upon his results are: (1) PODT (Proper Orthogonal Decomposition Theorem), which extracts the Large Eddy structure of stochastic processes from physical or computer simulation two-point covariances, and 2) LEIM (Large-Eddy Interaction Model), a predictive scheme for the dynamical large eddies based upon higher order turbulence modeling. Earlier Lumley's work (1964) forms the basis for the final member of the triad of numerical procedures: this predicts the global neutral modes of turbulence which have surprising agreement with both structural eigenmodes and those obtained from the dynamical equations. The ultimate goal of improved engineering design tools for turbulence may be near at hand, partly due to the power and storage of 'supermicrocomputer' workstations finally becoming adequate for the demanding numerics of these procedures.

Payne, Fred R.↗
NASA NTRS1994

Numerical simulation of the flow about the F-18 HARV at high angle of attack

This report summarizes research done over the past two years as part of NASA Grant NCC 2-729. This research has been aimed at validating numerical methods for computing the flow about the complete F-18 HARV at alpha = 30 deg and alpha = 45 deg. At 30 deg angle of attack, the flow about the F-18 is dominated by the formation, and subsequent breakdown, of strong vortices over the wing leading-edge extensions (LEX). As the angle of attack is increased to alpha = 45 deg, the fuselage forebody of the F-18 contains significant laminar and transitional regions which are not present at alpha = 30 deg. Further, the flow over the LEX at alpha = 45 deg is dominated by an unsteady shedding in time, rather than strong coherent vortices. This complex physics, combined with the complex geometry of a full aircraft configuration, provides a challenge for current computational fluid dynamics (CFD) techniques. The following sections present the numerical method and grid generation scheme that was used, a review of prior research done to numerically model the F-18 HARV, and a discussion of the current research. The current research is broken into two main topics: the effect of engine-inlet mass-flow rate on the F-18 vortex breakdown position, and the results using a refined F-18 computational model to compute the flow at alpha = 30 deg and alpha = 45 deg.

Murman, Scott M.↗
NASA NTRS1994

Numerical simulation of the flow about an F-18 aircraft in the high-alpha regime

The current research is aimed at developing and extending numerical methods to accurately predict the high Reynolds number flow about the NASA F-18 HARV at large angles of attack. The resulting codes are validated by comparison of the numerical results with in-flight aerodynamic measurements and flow visualization obtained on the HARV. Further, computations have been used to provide an analysis and numerical optimization of a pneumatic slot blowing concept, and a mechanical strake concept, for use as potential forebody flow control devices in improving high-alpha maneuverability.

Scott M Murman↗
NASA NTRS1994

Multiresolution representation and numerical algorithms: A brief review

In this paper we review recent developments in techniques to represent data in terms of its local scale components. These techniques enable us to obtain data compression by eliminating scale-coefficients which are sufficiently small. This capability for data compression can be used to reduce the cost of many numerical solution algorithms by either applying it to the numerical solution operator in order to get an approximate sparse representation, or by applying it to the numerical solution itself in order to reduce the number of quantities that need to be computed.

Harten, Amiram↗
NASA NTRS1994

Numerical aperture limits on efficient ball lens coupling of laser diodes to single-mode fibers with defocus to balance spherical aberration

The potential capabilities and limitations of single ball lenses for coupling laser diode radiation to single-mode optical fibers have been analyzed; parameters important to optical communications were specifically considered. These parameters included coupling efficiency, effective numerical apertures, lens radius, lens refractive index, wavelength, magnification in imaging the laser diode on the fiber, and defocus to counterbalance spherical aberration of the lens. Limiting numerical apertures in object and image space were determined under the constraint that the lens perform to the Rayleigh criterion of 0.25-wavelength (Strehl ratio = 0.80). The spherical aberration-defocus balance to provide an optical path difference of 0.25 wavelength units was shown to define a constant coupling efficiency (i.e., 0.56). The relative numerical aperture capabilities of the ball lens were determined for a set of wavelengths and associated fiber-core diameters of particular interest for single-mode fiber-optic communication. The results support general continuing efforts in the optical fiber communications industry to improve coupling links within such systems with emphasis on manufacturing simplicity, system packaging flexibility, relaxation of assembly alignment tolerances, cost reduction of opto-electronic components and long term reliability and stability.

Wilson, R. Gale↗
NASA NTRS1995

The space-time solution element method: A new numerical approach for the Navier-Stokes equations

This paper is one of a series of papers describing the development of a new numerical method for the Navier-Stokes equations. Unlike conventional numerical methods, the current method concentrates on the discrete simulation of both the integral and differential forms of the Navier-Stokes equations. Conservation of mass, momentum, and energy in space-time is explicitly provided for through a rigorous enforcement of both the integral and differential forms of the governing conservation laws. Using local polynomial expansions to represent the discrete primitive variables on each cell, fluxes at cell interfaces are evaluated and balanced using exact functional expressions. No interpolation or flux limiters are required. Because of the generality of the current method, it applies equally to the steady and unsteady Navier-Stokes equations. In this paper, we generalize and extend the authors' 2-D, steady state implicit scheme. A general closure methodology is presented so that all terms up through a given order in the local expansions may be retained. The scheme is also extended to nonorthogonal Cartesian grids. Numerous flow fields are computed and results are compared with known solutions. The high accuracy of the scheme is demonstrated through its ability to accurately resolve developing boundary layers on coarse grids. Finally, we discuss applications of the current method to the unsteady Navier-Stokes equations.

Scott, James R.↗
NASA NTRS1995

Multiple burn fuel-optimal orbit transfers: Numerical trajectory computation and neighboring optimal feedback guidance

This report describes current work in the numerical computation of multiple burn, fuel-optimal orbit transfers and presents an analysis of the second variation for extremal multiple burn orbital transfers as well as a discussion of a guidance scheme which may be implemented for such transfers. The discussion of numerical computation focuses on the use of multivariate interpolation to aid the computation in the numerical optimization. The second variation analysis includes the development of the conditions for the examination of both fixed and free final time transfers. Evaluations for fixed final time are presented for extremal one, two, and three burn solutions of the first variation. The free final time problem is considered for an extremal two burn solution. In addition, corresponding changes of the second variation formulation over thrust arcs and coast arcs are included. The guidance scheme discussed is an implicit scheme which implements a neighboring optimal feedback guidance strategy to calculate both thrust direction and thrust on-off times.

Chuang, C.-H.↗
NASA NTRS1995

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↗
NASA NTRS1995

Numerical simulation of the flow about the F-18 HARV at high angle of attack

This research has been aimed at validating numerical methods for computing the flow about the complete F-18 HARV at alpha = 30 deg and alpha = 45 deg. At 30 deg angle of attack, the flow about the F-18 is dominated by the formation, and subsequent breakdown, of strong vortices over the wing leading-edge extensions (LEX). As the angle of attack is increased to alpha = 45 deg, the fuselage forebody of the F-18 contains significant laminar and transitional regions which are not present at alpha = 30 deg. Further, the flow over the LEX at alpha = 45 deg is dominated by an unsteady shedding in time, rather than strong coherent vortices. This complex physics, combined with the complex geometry of a full-aircraft configuration, provides a challenge for current computational fluid dynamics (CFD) techniques. The following sections present the numerical method and grid generation scheme that was used, a review of prior research done to numerically model the F-18 HARV, and a discussion of the current research. The current research is broken into three main topics; the effect of engine-inlet mass-flow rate on the F-18 vortex breakdown position, the results using a refined F-18 computational model to compute the flow at alpha = 30 deg and alpha = 45 deg, and research done using the simplified geometry of an ogive-cylinder configuration to investigate the physics of unsteady shear-layer shedding. The last section briefly summarizes the discussion.

Murman, Scott M.↗
NASA NTRS1995

A numerical study of fundamental shock noise mechanisms

The results of this thesis demonstrate that direct numerical simulation can predict sound generation in unsteady aerodynamic flows containing shock waves. Shock waves can be significant sources of sound in high speed jet flows, on helicopter blades, and in supersonic combustion inlets. Direct computation of sound permits the prediction of noise levels in the preliminary design stage and can be used as a tool to focus experimental studies, thereby reducing cost and increasing the probability of a successfully quiet product in less time. This thesis reveals and investigates two mechanisms fundamental to sound generation by shocked flows: shock motion and shock deformation. Shock motion is modeled by the interaction of a sound wave with a shock. During the interaction, the shock wave begins to move and the sound pressure is amplified as the wave passes through the shock. The numerical approach presented in this thesis is validated by the comparison of results obtained in a quasi-one dimensional simulation with linear theory. Analysis of the perturbation energy demonstrated for the first time that acoustic energy is generated by the interaction. Shock deformation is investigated by the numerical simulation of a ring vortex interacting with a shock. This interaction models the passage of turbulent structures through the shock wave. The simulation demonstrates that both acoustic waves and contact surfaces are generated downstream during the interaction. Analysis demonstrates that the acoustic wave spreads cylindrically, that the sound intensity is highly directional, and that the sound pressure level increases significantly with increasing shock strength. The effect of shock strength on sound pressure level is consistent with experimental observations of shock noise, indicating that the interaction of a ring vortex with a shock wave correctly models a dominant mechanism of shock noise generation.

Meadows, Kristine R.↗
NASA NTRS1994

Electron transport and energy degradation in the ionosphere: Evaluation of the numerical solution, comparison with laboratory experiments and auroral observations

Auroral electron transport calculations are a critical part of auroral models. We evaluate a numerical solution to the transport and energy degradation problem. The numerical solution is verified by reproducing simplified problems to which analytic solutions exist, internal self-consistency tests, comparison with laboratory experiments of electron beams penetrating a collision chamber, and by comparison with auroral observations, particularly the emission ratio of the N2 second positive to N2(+) first negative emissions. Our numerical solutions agree with range measurements in collision chambers. The calculated N(2)2P to N2(+)1N emission ratio is independent of the spectral characteristics of the incident electrons, and agrees with the value observed in aurora. Using different sets of energy loss cross sections and different functions to describe the energy distribution of secondary electrons that emerge from ionization collisions, we discuss the uncertainties of the solutions to the electron transport equation resulting from the uncertainties of these input parameters.

Lummerzheim, D.↗
NASA NTRS1992

A numerical study of three-dimensional surface tension driven convection with fre surface deformation

The steady three-dimensional thermocapillary motion with a deformable free surface is studied numerically in both normal and zero gravity environments. Flow configurations consist of a square cavity heated from the side. In the analysis, the free surface is allowed to deform and the grid distribution is adapted to the surface deformation. The divergence-free condition is satisfied by using a dual time-stepping approach in the numerical scheme. Convective flux derivatives are evaluated using a third-order accurate upwind-biased flux-split differencing technique. The numerical solutions at the midplane of the square cavity are compared with the results from two-dimensional calculations. In addition, numerial results for cases under zero and normal gravity conditions are compared. Significantly different flow structures and surface deformation have been observed. The comparison of calculated results will be compared with experimental data in the updated version of this paper.

Hsieh, Kwang-Chung↗
NASA NTRS1995

Numerical Stability In Hyperbolic Boundary-Value Problems

Technical memorandum discusses stability of numerical solutions involving semidiscrete approximations to hyperbolic partial differential equations in initial-and-boundary-value problems. Topic of practical significance because hyperbolic partial differential equations arise in mathematical modeling of waves and blasts. Solutions often needed over restricted regions of unbounded spaces. Outer boundaries artificial, introduced only to limit domains of numerical solutions. Conditions at such artificial boundaries cause numerical instabilities that degrade computed solutions.

Warming, Robert F.↗