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 541 records · Page 30

Melting and its relationship to impact crater morphology

Shock-melting features occur on planets at scales that range from micrometers to megameters. It is the objective of this study to determine the extent of thickness, volume geometry of the melt, and relationship with crater morphology. The variation in impact crater morphology on planets is influenced by a broad range of parameters: e.g., planetary density, thermal state, strength, impact velocity, gravitational acceleration. We modeled the normal impact of spherical projectiles on a semi-infinite planet over a broad range of conditions using numerical techniques.

Okeefe, John D.↗

Investigation of upwind, multigrid, multiblock numerical schemes for three dimensional flows. Volume 1: Runge-Kutta methods for a thin layer Navier-Stokes solver

A state-of-the-art computer code has been developed that incorporates a modified Runge-Kutta time integration scheme, upwind numerical techniques, multigrid acceleration, and multi-block capabilities (RUMM). A three-dimensional thin-layer formulation of the Navier-Stokes equations is employed. For turbulent flow cases, the Baldwin-Lomax algebraic turbulence model is used. Two different upwind techniques are available: van Leer's flux-vector splitting and Roe's flux-difference splitting. Full approximation multi-grid plus implicit residual and corrector smoothing were implemented to enhance the rate of convergence. Multi-block capabilities were developed to provide geometric flexibility. This feature allows the developed computer code to accommodate any grid topology or grid configuration with multiple topologies. The results shown in this dissertation were chosen to validate the computer code and display its geometric flexibility, which is provided by the multi-block structure.

Cannizzaro, Frank E.↗

Human impacts on soil carbon dynamics of deep-rooted Amazonian forests and effect of land use change on the carbon cycle in Amazon soils

The main objective of these NASA-funded projects is to improve our understanding of land-use impacts on soil carbon dynamics in the Amazon Basin. Soil contains approximately one half of tropical forest carbon stocks, yet the fate of this carbon following forest impoverishment is poorly studied. Our mechanistics approach draws on numerous techniques for measuring soil carbon outputs, inputs, and turnover time in the soils of adjacent forest and pasture ecosystems at our research site in Paragominas, state of Para, Brazil. We are scaling up from this site-specific work by analyzing Basin-wide patterns in rooting depth and rainfall seasonality, the two factors that we believe should explain much of the variation in tropical soil carbons dynamics. In this report, we summarize ongoing measurements at our Paragominas study site, progress in employing new field data to understand soil C dynamics, and some surprising results from our regional, scale-up work.

Nepstad, Daniel↗

Automatic computation of wing-fuselage intersection lines and fillet inserts with fixed-area constraint

Procedures for automatic computation of wing-fuselage juncture geometry are described. These procedures begin with a geometry in wave-drag format. First, an intersection line is computed by extrapolating the wing to the fuselage. Then two types of filleting procedures are described, both of which utilize a combination of analytical and numerical techniques appropriate for automatic calculation. An analytical technique for estimating the added volume due to the fillet is derived, and an iterative procedure for revising the fuselage to compensate for this additional volume is given. Sample results are included in graphical form.

Barger, Raymond L.↗

A new observational and numerical study of tidal interactions in M81-M82-NGC3077 system

A nearby system of interacting galaxies M81-M82-NGC3077 triplet (D = 3.3 Mpc; Freeman & Madore 1988) has been studied using multi-wavelength observations and numerical simulations to obtain a comprehensive understanding on the dynamics and the consequences of tidal interactions in a group environment. The VLA 12-field Mosaic H I observations of 2 x 1.5 deg. region have revealed a vast array of H I filaments which suggests that the severity and extent of tidal disruptions far exceed the previous estimates. A tidal remnant of the former H I disk of M82 extending up to 30 kpc (in projection) is identified for the first time, and the pervasive effects of the tidal disruption are traced into the inner disk by optical and CO observations, including a kinematic trace of a large scale bar potential (Yun, Ho, & Lo 1992). The H I disk of M81 is traced out to 40 kpc in radius, and a large scale (l approx. 20 kpc) velocity anomaly ('High Velocity Trough'), which may be a remnant of a gaseous collision, is found within the disk of M81. The large H I bridge between M81 and NGC 3077 (van der Hulst 1979) is also found to extend approx. greater than 50 kpc further, bending around NGC 3077, toward M82. The total H I detected in this experiment, 5.6 x 10(exp 9) solar mass, represents the majority of the single-dish flux (Appleton, Davies, & Stephenson 1981) and suggests that the bulk of H I found in the region belongs to the three galaxies and the tidal filaments. The impact and details of the tidal interactions have been further examined through the use of numerical techniques. The 'restricted 3-body' approach was used to simulate the observed distribution of tidal H I streamers connecting the three galaxies, and the success of the simulation is further strengthened by the accurate predictions on the gas kinematics.

Yun, M. S.↗

Boundary conditions for direct computation of aerodynamic sound generation

A numerical scheme suitable for the computation of both the near field acoustic sources and the far field sound produced by turbulent free shear flows utilizing the Navier-Stokes equations is presented. To produce stable numerical schemes in the presence of shear, damping terms must be added to the boundary conditions. The numerical technique and boundary conditions are found to give stable results for computations of spatially evolving mixing layers.

Colonius, Tim↗

Three-dimensional flow phenomena in a transonic, high-through-flow, axial-flow compressor stage

A detailed aerodynamic study of a transonic, high-through-flow, single stage compressor is presented. The compressor stage was comprised of a low-aspect-ratio rotor combined alternately with two different stator designs. Both experimental and numerical studies are conducted to understand the details of the complex flow field present in this stage. Aerodynamic measurements using high-frequency, Kulite pressure transducers and conventional probes are compared with results from a three-dimensional viscous flow analysis. A steady multiple blade row approach is used in the numerical technique to examine the detailed flow structure inside the rotor and the stator passages. The comparisons indicate that many flow field features are correctly captured by viscous flow analysis, and therefore unmeasured phenomena can be studied with some level of confidence.

Copenhaver, William W.↗

A note on the effects of viscosity on the stability of a trailing-line vortex

The linear stability of the Batchelor (1964) vortex is examined with emphasis on new viscous modes recently found numerically by Khorrami (1991). Unlike the previously reported inviscid modes of instability, these modes are destabilized by viscosity and exhibit small growth rates at large Reynolds numbers. The analysis presented here uses a combination of asymptotic and numerical techniques. The results confirm the existence of the additional modes of instability due to viscosity.

Duck, Peter W.↗

Source term evaluation for combustion modeling

A modification is developed for application to the source terms used in combustion modeling. The modification accounts for the error of the finite difference scheme in regions where chain-branching chemical reactions produce exponential growth of species densities. The modification is first applied to a one-dimensional scalar model problem. It is then generalized to multiple chemical species, and used in quasi-one-dimensional computations of shock-induced combustion in a channel. Grid refinement studies demonstrate the improved accuracy of the method using this modification. The algorithm is applied in two spatial dimensions and used in simulations of steady and unsteady shock-induced combustion. Comparisons with ballistic range experiments give confidence in the numerical technique and the 9-species hydrogen-air chemistry model.

Sussman, Myles A.↗

Fuel optimal reorientation of axisymmetric spin-stabilized satellites

The fuel optimal spin axis reorientation problem for large and small maneuver times is solved using a new numerical technique. The approach, which employs the adaptive grid bisection search, is shown to be useful in determining the fuel optimal impulsive control strategy, particularly when rapid maneuvers are desired. For large enough maneuver times, the method produces the two-impulse reorientation maneuver. For even larger maneuver times, the existence of multiple optimal pulsing schemes is demonstrated.

Redmond, Jim↗

Computation of unsteady shock-induced combustion using logarithmic species conservation equations

Numerical simulations are used to investigate periodic combustion instabilities observed in ballistic-range experiments of blunt bodies flying at supersonic speeds through hydrogen-air mixtures. The computations are validated by comparing experimental shadowgraphs with shadowgraphs created from the computed flowfields and by comparing the experimentally measured instability frequencies with computed frequencies. The numerical simulations use a logarithmic transformation of the species conservation equations as a way to reduce the grid requirements for computing shock-induced combustion. The transformation is applied to the Euler equations coupled to a detailed hydrogen-air chemical reaction mechanism with 13 species and 33 reactions. The resulting differential equations are solved using a finite volume formulation and a two-step predictor-corrector scheme to advance the solution in time. Results are presented and compared for both a flux-vector splitting scheme and an upwind TVD scheme. The computations add insight to the physical processes observed in the experiments and the numerical methods needed to simulate them. The usefulness of the ballistic-range experiments for the validation of numerical techniques and chemical kinetic models is also demonstrated.

Wilson, Gregory J.↗

Two-dimensional hybrid continuum/particle approach for rarefied flows

A hybrid numerical technique previously developed for one-dimensional rarefied gas flows is generalized to two dimensions. The method is based on the fact that the flowfield that develops near a body in a rarefied gas typically contains local regions of continuum, transitional and free-molecular flow. By utilizing the solution technique most appropriate for each region and coupling the techniques in an interface region where both are applicable, more computationally efficient solutions can be obtained, or equivalently, more complex flowfields can be analyzed. The present method combines finite difference solution of the Navier-Stokes equations in the continuum regions, with Direct Simulation Monte Carlo in the more rarefied regions. The two schemes are coupled interactively via a general conservative flux boundary condition. The method is tested by application to the model problem of pressure-driven rarefied flow through a slit. Results show the hybrid scheme offers a speedup of a factor of nearly two for the nominal conditions considered, due to the decrease in the size of the Direct Simulation domain.

Wadsworth, Dean C.↗

A numerical model for predicting crack path and modes of damage in unidirectional metal matrix composites

A finite element-based numerical technique has been developed to simulate damage growth in unidirectional composites. This technique incorporates elastic-plastic analysis, micromechanics analysis, failure criteria, and a node splitting and node force relaxation algorithm to create crack surfaces. Any combination of fiber and matrix properties can be used. One of the salient features of this technique is that damage growth can be simulated without pre-specifying a crack path. In addition, multiple damage mechanisms in the forms of matrix cracking, fiber breakage, fiber-matrix debonding and plastic deformation are capable of occurring simultaneously. The prevailing failure mechanism and the damage (crack) growth direction are dictated by the instantaneous near-tip stress and strain fields. Once the failure mechanism and crack direction are determined, the crack is advanced via the node splitting and node force relaxation algorithm. Simulations of the damage growth process in center-slit boron/aluminum and silicon carbide/titanium unidirectional specimens were performed. The simulation results agreed quite well with the experimental observations.

Bakuckas, J. G.↗

A study of the influence of the data acquisition system sampling rate on the accuracy of measured acceleration loads for transport aircraft

A research effort was initiated at National Aeronautics and Space Administration (NASA) Langley Research Center (LaRC), to describe the relationship between the sampling rate and the accuracy of acceleration loads obtained from the data acquisition system of a transport aircraft. An accelerometer was sampled and digitized at a rate of 100 samples per second onboard a NASA Boeing 737 (B-737) flight research aircraft. Numerical techniques were used to reconstruct 2.5 hours of flight data into its original input waveform and then re-sample the waveform into rates of 4, 8, 16, and 32 samples per second. Peak-between-means counting technique and power spectral analysis were used to evaluate each sampling rate using the 32 samples per second data as the comparison. This paper presents the results from these methods and includes in appendix A, the peak-between-means counting results used in a general fatigue analysis for each of the sampling rates.

Whitehead, Julia H.↗

Flux-split algorithms for hypersonic flows

This study reviews numerical techniques for the simulation of fluid flows spanning the range of reactive regimes, from local chemical equilibrium to full thermochemical non-equilibrium. In particular, characteristic-based algorithms are considered, of the flux-vector and the flux-difference type, which are becoming ever more popular due to their accurate rendition of physically complicated, shock-wave dominated flows. Consideration is given to the problems associated with modeling the thermo-chemical behavior of reactive mixtures of thermally perfect gases, and a few topics of current research are outlined.

Cinnella, P.↗

Airborne gravimetry, altimetry, and GPS navigation errors

Proper interpretation of airborne gravimetry and altimetry requires good knowledge of aircraft trajectory. Recent advances in precise navigation with differential GPS have made it possible to measure gravity from the air with accuracies of a few milligals, and to obtain altimeter profiles of terrain or sea surface correct to one decimeter. These developments are opening otherwise inaccessible regions to detailed geophysical mapping. Navigation with GPS presents some problems that grow worse with increasing distance from a fixed receiver: the effect of errors in tropospheric refraction correction, GPS ephemerides, and the coordinates of the fixed receivers. Ionospheric refraction and orbit error complicate ambiguity resolution. Optimal navigation should treat all error sources as unknowns, together with the instantaneous vehicle position. To do so, fast and reliable numerical techniques are needed: efficient and stable Kalman filter-smoother algorithms, together with data compression and, sometimes, the use of simplified dynamics.

Colombo, Oscar L.↗

Investigation of high-speed free shear flows using improved pressure-strain correlated Reynolds stress turbulence model

A high-speed shear layer is studied using compressibility corrected Reynolds stress turbulence model which employs newly developed model for pressure-strain correlation. MacCormack explicit prediction-corrector method is used for solving the governing equations and the turbulence transport equations. The stiffness arising due to source terms in the turbulence equations is handled by a semi-implicit numerical technique. Results obtained using the new model show a sharper reduction in growth rate with increasing convective Mach number. Some improvements were also noted in the prediction of the normalized streamwise stress and Reynolds shear stress. The computed results are in good agreement with the experimental data.

Tiwari, S. N.↗

Thermocapillary flow with evaporation and condensation and its effect on liquid retention in low-G fluid acquisition devices

The steady motion, thermal and free surface behavior of a volatile, wetting liquid in microgravity are studied using scaling and numerical techniques. The objective is to determine whether the thermocapillary and two-phase convection arising from thermodynamic nonequilibrium along the porous surfaces of spacecraft liquid acquisition devices could cause the retention failures observed with liquid hydrogen and heated vapor pressurant. Why these devices seem immune to retention loss when pressurized with heated helium or heated directly through the porous structure was also examined. Results show that highly wetting fluids exhibit large negative and positive dynamic pressure gradients towards the meniscus interline when superheated and subcooled, respectively. With superheating, the pressure variation and recoil force arising from liquid/vapor phase change exert the same influence on surface morphology and promote retention. With subcooling, however, the pressure distribution produces a suction that degrades mechanical equilibrium of the surface. This result indicates that thermocapillary-induced deformation arising from subcooling and condensation is the likely cause for retention loss. In addition, increasing the level of nonequilibrium by reducing accommodation coefficient suppresses deformation and explains why this failure mode does not occur in instances of direct screen heating or pressurization with a heated inert gas.

Schmidt, George R.↗