Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical techniques”

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 1,261 records · Page 70
All sourcesNASA NTRSDOE OSTI
NASA NTRS1996

Artificial Boundary Conditions Based on the Difference Potentials Method

While numerically solving a problem initially formulated on an unbounded domain, one typically truncates this domain, which necessitates setting the artificial boundary conditions (ABC's) at the newly formed external boundary. The issue of setting the ABC's appears to be most significant in many areas of scientific computing, for example, in problems originating from acoustics, electrodynamics, solid mechanics, and fluid dynamics. In particular, in computational fluid dynamics (where external problems present a wide class of practically important formulations) the proper treatment of external boundaries may have a profound impact on the overall quality and performance of numerical algorithms. Most of the currently used techniques for setting the ABC's can basically be classified into two groups. The methods from the first group (global ABC's) usually provide high accuracy and robustness of the numerical procedure but often appear to be fairly cumbersome and (computationally) expensive. The methods from the second group (local ABC's) are, as a rule, algorithmically simple, numerically cheap, and geometrically universal; however, they usually lack accuracy of computations. In this paper we first present a survey and provide a comparative assessment of different existing methods for constructing the ABC's. Then, we describe a relatively new ABC's technique of ours and review the corresponding results. This new technique, in our opinion, is currently one of the most promising in the field. It enables one to construct such ABC's that combine the advantages relevant to the two aforementioned classes of existing methods. Our approach is based on application of the difference potentials method attributable to V. S. Ryaben'kii. This approach allows us to obtain highly accurate ABC's in the form of certain (nonlocal) boundary operator equations. The operators involved are analogous to the pseudodifferential boundary projections first introduced by A. P. Calderon and then also studied by R. T. Seeley. The apparatus of the boundary pseudodifferential equations, which has formerly been used mostly in the qualitative theory of integral equations and PDE'S, is now effectively employed for developing numerical methods in the different fields of scientific computing.

Tsynkov, Semyon V.↗
NASA NTRS1987

A finite volume Euler calculation of the aerodynamics of transonic airfoil-vortex interaction

Unsteady inviscid transonic airfoil-vortex interaction is numerically analyzed by solving the two-dimensional unsteady Euler equations in integral form using a finite volume scheme. The solution procedure is based on an explicit Runge-Kutta time-stepping scheme wherein the spatial terms are central-differenced and a combination of second- and fourth-differences in the flow variables is used to form the numerical dissipation terms to stabilize the scheme. A velocity decomposition technique is applied to alleviate the problem of vortex diffusion by the numerical dissipation terms and to treat the interaction of a Rankine vortex with an airfoil accurately. Results obtained are compared with available numerical data.

Damodaran, Murali↗
NASA NTRS1990

Microwave holographic metrology for antenna diagnosis

Advances in antenna diagnostic methodologies have been very significant in recent years. In particular, microwave holographic diagnostic techniques have been applied very successfully in improving the performance of reflector and array antennas. These techniques use the knowledge of the measured amplitude and phase of the antenna radiated fields and then take advantage of the existing Fourier transform relationships between the radiated fields and the effective aperture or current distribution to eventually determine the reflector surface or array excitation coefficients anomalies. In this paper an overview of the recent developments in applying microwave holography is presented. The theoretical, numerical and measurement aspects of this technique is detailed by providing representative results.

Rahmat-Samii, Y.↗
NASA NTRS1992

Homogeneous quantum electrodynamic turbulence

The electromagnetic field equations and Dirac equations for oppositely charged wave functions are numerically time-integrated using a spatial Fourier method. The numerical approach used, a spectral transform technique, is based on a continuum representation of physical space. The coupled classical field equations contain a dimensionless parameter which sets the strength of the nonlinear interaction (as the parameter increases, interaction volume decreases). For a parameter value of unity, highly nonlinear behavior in the time-evolution of an individual wave function, analogous to ideal fluid turbulence, is observed. In the truncated Fourier representation which is numerically implemented here, the quantum turbulence is homogeneous but anisotropic and manifests itself in the nonlinear evolution of equilibrium modal spatial spectra for the probability density of each particle and also for the electromagnetic energy density. The results show that nonlinearly interacting fermionic wave functions quickly approach a multi-mode, dynamic equilibrium state, and that this state can be determined by numerical means.

Shebalin, John V.↗
NASA NTRS2023

Phase-Space Energization of Ions in Oblique Shocks

Examining energization of kinetic plasmas in phase space is a growing topic of interest, owing to the wealth of data in phase space compared to traditional bulk energization diagnostics. Via the field-particle correlation (FPC) technique and using multiple means of numerically integrating the plasma kinetic equation, we have studied the energization of ions in phase space within oblique collisionless shocks. The perspective afforded to us with this analysis in phase space allows us to characterize distinct populations of energized ions. In particular, we focus on ions that reflect multiple times off the shock front through shock-drift acceleration, and how to distinguish these different reflected populations in phase space using the FPC technique. We further extend our analysis to simulations of three-dimensional shocks undergoing more complicated dynamics, such as shock ripple, to demonstrate the ability to recover the phase-space signatures of this energization process in a more general system. This work thus extends previous applications of the FPC technique to more realistic collisionless shock environments, providing stronger evidence of the technique's utility for simulation, laboratory, and spacecraft analysis.

Plasma astrophysics↗
NASA NTRS1986

Optimal control of quaternion propagation errors in spacecraft navigation

Optimal control techniques are used to drive the numerical error (truncation, roundoff, commutation) in computing the quaternion vector to zero. The normalization of the quaternion is carried out by appropriate choice of a performance index, which can be optimized. The error equations are derived from Friedland's (1978) theoretical development, and a matrix Riccati equation results for the computation of the gain matrix. Simulation results show that a high precision of the order of 10 to the -12th can be obtained using this technique in meeting the q(T)q=1 constraint. The performance of the estimator in the presence of the feedback control that maintains the normalization, is studied.

Vathsal, S.↗
NASA NTRS1987

Advances and trends in computational structural mechanics

The development status and applicational range of techniques in computational structural mechanics (CSM) are evaluated with a view to advances in computational models for material behavior, discrete-element technology, quality assessment, the control of numerical simulations of structural response, hybrid analysis techniques, techniques for large-scale optimization, and the impact of new computing systems on CSM. Primary pacers of CSM development encompass prediction and analysis of novel materials for structural components, computational strategies for large-scale structural calculations, and the assessment of response prediction reliability together with its adaptive improvement.

Noor, Ahmed K.↗
NASA NTRS1990

Upwind techniques for flows with multiple translational temperatures

The numerical simulation of high-speed flows out of chemical and thermodynamic equilibrium is considered. A general model for the species internal energy is utilized, which allows for multiple translational temperatures. A study of the characteristic wave speeds for the case of two translational temperatures, one associated with the average thermal speed of the heavy particles and the other associated with the translational temperature of free electrons, yields some very interesting results, with acoustic wave speeds of the symmetric form u + or - a no longer appearing. An asymptotic analysis in terms of small electron mass fractions yields approximate wave speed which recover the symmetric form, with the addition of a new wave speed associated with the electron energy equation.

Cinnella, Pasquale↗
NASA NTRS1993

An evaluation of excitation techniques for time domain based flutter data processing

Numerically simulated, wind tunnel and flight test experimental data were used to evaluate various types of system excitation sources for subcritical flutter testing. Emphasis was placed upon the determination of modal frequencies and damping from relatively short time series records. Parameter identification based upon digital time series models was used to predict modal characteristics from subcritical test data. Various types of excitation were considered. These were: 1) impulsive loading and the resulting transient free vibration response, 2) random forced response, in which the 'mechanical' forcing function could be directly measured, 3) sine sweep forced response, in which the 'mechanical' forcing function could be directly measured, and 4) forced response from unsteady aerodynamic excitation or turbulence, in which the excitation could only be indirectly inferred from other measurements.

Pinkelman, J. K.↗
NASA NTRS2000

Magnetic Susceptibility Effects and Lorentz Damping in Diamagnetic Fluids

A great number of crystals (semi-conductor and protein) grown in space are plagued by convective motions which contribute to structural flaws. The character of these instabilities is not well understood but is associated with density variations in the presence of residual gravity and g-jitter. Both static and dynamic (rotating or travelling wave) magnetic fields can be used to reduce the effects of convection in materials processing. In semi-conductor melts, due to their relatively high electrical conductivity, the induced Lorentz force can be effectively used to curtail convective effects. In melts/solutions with reduced electrical conductivity, such as aqueous solutions used in solution crystal growth, protein crystal growth and/or model fluid experiments for simulating melt growth, however, the variation of the magnetic susceptibility with temperature and/or concentration can be utilized to better damp fluid convection than the Lorentz force method. This paper presents a comprehensive, comparative numerical study of the relative damping effects using static magnetic fields and gradients in a simple geometry subjected to a thermal gradient. The governing equations are formulated in general terms and then simplified for the numerical calculations. Operational regimes, based on the best damping technique for different melts/solutions are identified based on fluid properties. Comparisons are provided between the numerical results and available results from experiments in surveyed literature.

Ramachandran, Narayanan↗
NASA NTRS1982

Finite element analysis of the effect of a non-planar solid-liquid interface on the lateral solute segregation during unidirectional solidification

The effect of solid-liquid interface shape on lateral solute segregation during steady-state unidirectional solidification of a binary mixture is calculated under the assumption of no convection in the liquid. A finite element technique is employed to compute the concentration field in the liquid and the lateral segregation in the solid with a curved boundary between the liquid and solid phases. The computational model is constructed assuming knowledge of the solid-liquid interface shape; no attempt is made to relate this shape to the thermal field. The influence of interface curvature on the lateral compositional variation is investigated over a range of system parameters including diffusivity, growth speed, distribution coefficient, and geometric factors of the system. In the limiting case of a slightly nonplanar interface, numerical results from the finite element technique are in good agreement with the analytical solutions of Coriell and Sekerka obtained by using linear theory. For the general case of highly non-planar interface shapes, the linear theory fails and the concentration field in the liquid as well as the lateral solute segregation in the solid can be calculated by using the finite element method.

Carlson, F. M.↗
NASA NTRS2016

Development of a Perfectly Matched Layer Technique for a Discontinuous-Galerkin Spectral-Element Method

The numerical simulation of many aerodynamic non-periodic flows of practical interest involves discretized computational domains that often must be artificially truncated. Appropriate boundary conditions are required at these truncated domain boundaries, and ideally, these boundary conditions should be perfectly "absorbing" or "nonreflecting" so that they do not contaminate the flow field in the interior of the domain. The proper specification of these boundaries is critical to the stability, accuracy, convergence, and quality of the numerical solution, and has been the topic of considerable research. The need for accurate boundary specification has been underscored in recent years with efforts to apply higher-fidelity methods (DNS, LES) in conjunction with high-order low-dissipation numerical schemes to realistic flow configurations. One of the most popular choices for specifying these boundaries is the characteristics-based boundary condition where the linearized flow field at the boundaries are decomposed into characteristic waves using either one-dimensional Riemann or other multi-dimensional Riemann approximations. The values of incoming characteristics are then suitably modified. The incoming characteristics are specified at the in flow boundaries, and at the out flow boundaries the variation of the incoming characteristic is zeroed out to ensure no reflection. This, however, makes the problem ill-posed requiring the use of an ad-hoc parameter to allow small reflections that make the solution stable. Generally speaking, such boundary conditions work reasonably well when the characteristic flow direction is normal to the boundary, but reflects spurious energy otherwise. An alternative to the characteristic-based boundary condition is to add additional "buffer" regions to the main computational domain near the artificial boundaries, and solve a different set of equations in the buffer region in order to minimize acoustic reflections. One approach that has been used involves modeling the pressure fluctuations as acoustic waves propagating in the far-field relative to a single noise-source inside the buffer region. This approach treats vorticity-induced pressure fluctuations the same as acoustic waves. Another popular approach, often referred to as the "sponge layer," attempts to dampen the flow perturbations by introducing artificial dissipation in the buffer region. Although the artificial dissipation removes all perturbations inside the sponge layer, incoming waves are still reflected from the interface boundary between the computational domain and the sponge layer. The effect of these refkections can be somewhat mitigated by appropriately selecting the artificial dissipation strength and the extent of the sponge layer. One of the most promising variants on the buffer region approach is the Perfectly Matched Layer (PML) technique. The PML technique mitigates spurious reflections from boundaries and interfaces by dampening the perturbation modes inside the buffer region such that their eigenfunctions remain unchanged. The technique was first developed by Berenger for application to problems involving electromagnetic wave propagation. It was later extended to the linearized Euler, Euler and Navier-Stokes equations by Hu and his coauthors. The PML technique ensures the no-reflection property for all waves, irrespective of incidence angle, wavelength, and propagation direction. Although the technique requires the solution of a set of auxiliary equations, the computational overhead is easily justified since it allows smaller domain sizes and can provide better accuracy, stability, and convergence of the numerical solution. In this paper, the PML technique is developed in the context of a high-order spectral-element Discontinuous Galerkin (DG) method. The technique is compared to other approaches to treating the in flow and out flow boundary, such as those based on using characteristic boundary conditions and sponge layers. The superiority of the current PML technique over other approaches is demonstrated for a range of test cases, viz., acoustic pulse propagation, convective vortex, shear layer flow, and low-pressure turbine cascade flow. The paper is structured as follows. We first derive the PML equations from the non{linear Euler equations. A short description of the higher-order DG method used is then described. Preliminary results for the four test cases considered are then presented and discussed. Details regarding current work that will be included in the final paper are also provided.

Spectral↗
NASA NTRS1990

Visualization of fluid dynamics at NASA Ames

Some of the hardware and software tools and techniques in use at NASA's numerical aerodynamic simulation facility for the analysis of computational fluid dynamics are described. The visualization process can be illustrated by video tapes and stereo pictures. Although these visualization tools have dramatically improved the ability to conduct research in fluid dynamics, a comparison of the current environment for analysis with an 'ideal' environment illustrates that there are still major improvements that should be made. The most time-consuming task in future analyses of the increasingly complex computer simulations will be the extraction and clear display of the key features. In addition, the interface between the workstation and the scientist should be improved significantly. Current research on techniques for creating these improvements is described.

Watson, V.↗
NASA NTRS2000

Modeling of Ultrasonically Generated Liquid-Liquid Dispersions During Controlled Directional Solidification

There are innumerable two-component systems in which two very different liquid phases co-exist in equilibrium over a range of temperature and composition, e.g., oil and water, salt fluxes and solders, aluminum and lead. Often it is of practical concern to fabricate a solid component consisting of a uniform dispersion of one phase in the other. Unfortunately, uniform microstructural development during solidification of two immiscible liquids is hampered by inherent, often large, density differences between the phases that lead to severe segregation. Uniformity is also compromised by preferential wetting and coalescence phenomena. It is, however, well known that ultrasonic energy can initiate and maintain a fine liquid-liquid dispersion. The work presented here extends that observation by application of ultrasonic energy to promote uniform phase incorporation during controlled directional solidification. To this end experiments with the transparent organic, immiscible, succinonitrile-glycerol system were conducted and the numerous processing parameters associated with this technique were evaluated in view of optimizing dispersion uniformity. In view of the initial experimental results a model that predicts the dispersed liquid droplet size as a function of material properties, sample geometry, and applied energy has been developed. In the mathematical model we consider the ultrasonic field in an experimental ampoule of length L and diameter D induced by a probe having a vibration frequency of f=2OKhz (circular frequency omega = 2 pi f). The amplitude is adjustable from A=65 to 13Omicrons. The probe tip diameter is d, the liquid has a density of p, in which the speed of sound and surface tension are, respectively, c and sigma. The mathematical model and numerical investigation for the experiments [1] is done using the following assumptions: (i) The droplet size is small in comparison to the sound wave length; (ii) The forces between droplets are neglected (relative concentration is small); (iii) The droplet is stable if the kinetic energy, E(sub K), of the liquid motion due to ultrasonic field influence is less then the binding energy, E(sub S), due to the surface tension (it is easy to show that the surface energy of two droplets resulting from one is larger by about a factor of two.); (iv) The stability limit is characterized by E(sub S) to approx. E(sub K).

Grugel, R. N.↗
NASA NTRS1991

Evolution and dynamics of shear-layer structures in near-wall turbulence

Near-wall flow structures in turbulent shear flows are analyzed, with particular emphasis on the study of their space-time evolution and connection to turbulence production. The results are obtained from investigation of a database generated from direct numerical simulation of turbulent channel flow at a Reynolds number of 180 based on half-channel width and friction velocity. New light is shed on problems associated with conditional sampling techniques, together with methods to improve these techniques, for use both in physical and numerical experiments. The results clearly indicate that earlier conceptual models of the processes associated with near-wall turbulence production, based on flow visualization and probe measurements need to be modified. For instance, the development of asymmetry in the spanwise direction seems to be an important element in the evolution of near-wall structures in general, and for shear layers in particular. The inhibition of spanwise motion of the near-wall streaky pattern may be the primary reason for the ability of small longitudinal riblets to reduce turbulent skin friction below the value for a flat surface.

Johansson, Arne V.↗
NASA NTRS1992

Implementation of Vigneron's streamwise pressure gradient approximation in the PNS equations

For single sweep parabolized Navier-Stokes solvers, the streamwise pressure gradient must be modified in the subsonic region to eliminate numerical instabilities. The accuracy of this modification on the solution of the parabolized Navier-Stokes equations with Vigneron's technique is shown to depend on how the numerical approximation of the pressure gradient is formed. A simple test case of supersonic laminar flow over a flat plate is computed with two different numerical methods for solving the PNS equations. Significant errors in the temperature profile and skin friction coefficient are demonstrated using a fully conservative differencing treatment of Vigneron's splitting for the pressure gradient typically used in parabolized Navier-Stokes solvers. The physical reason for this error is discussed. An alternate formulation is demonstrated which minimizes these errors.

Morrison, Joseph H.↗
NASA NTRS1983

Refinement and application of acoustic impulse technique to study nozzle transmission characteristics

An improved acoustic impulse technique was developed and was used to study the transmission characteristics of duct/nozzle systems. To accomplish the above objective, various problems associated with the existing spark-discharge impulse technique were first studied. These included (1) the nonlinear behavior of high intensity pulses, (2) the contamination of the signal with flow noise, (3) low signal-to-noise ratio at high exhaust velocities, and (4) the inability to control or shape the signal generated by the source, specially when multiple spark points were used as the source. The first step to resolve these problems was the replacement of the spark-discharge source with electroacoustic driver(s). These included (1) synthesizing on acoustic impulse with acoustic driver(s) to control and shape the output signal, (2) time domain signal averaging to remove flow noise from the contaminated signal, (3) signal editing to remove unwanted portions of the time history, (4) spectral averaging, and (5) numerical smoothing. The acoustic power measurement technique was improved by taking multiple induct measurements and by a modal decomposition process to account for the contribution of higher order modes in the power computation. The improved acoustic impulse technique was then validated by comparing the results derived by an impedance tube method. The mechanism of acoustic power loss, that occurs when sound is transmitted through nozzle terminations, was investigated. Finally, the refined impulse technique was applied to obtain more accurate results for the acoustic transmission characteristics of a conical nozzle and a multi-lobe multi-tube supressor nozzle.

Salikuddin, M.↗
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.↗