SEARCH · Engineering Papers

Results for “Approximate boundary”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,099 records · Page 61

All sources NASA NTRS DOE OSTI

NASA NTRS1981

Development of a linear unsteady aerodynamic analysis for finite-deflection subsonic cascades

A linear unsteady potential flow analysis, which accounts for the effects of blade geometry and steady turning, is being developed to predict the aerodynamic response to blade vibrations in the fan, compressor or turbine stages of modern jet engines. In previous work numerical solutions were restricted to cascades of sharp-edged blades aligned with the mean flow. Under the present effort the solution procedure has been extended to treat blades with rounded or blunt edges. As part of this effort an analytical model-problem study has been conducted to clarify the behavior of first-order perturbation solutions in the vicinity of airfoil edges. Further, a numerical approximation using concepts from singular perturbation theory has been developed to resolve the unsteady boundary value problem for cascades of blunt leading-edged blades. Numerical results for NACA 0012 cascades, including detailed results in the vicinity of a blade leading edge, are presented and evaluated.

Verdon, J. M.↗

NASA NTRS1981

Shape determination and control for large space structures

An integral operator approach is used to derive solutions to static shape determination and control problems associated with large space structures. Problem assumptions include a linear self-adjoint system model, observations and control forces at discrete points, and performance criteria for the comparison of estimates or control forms. Results are illustrated by simulations in the one dimensional case with a flexible beam model, and in the multidimensional case with a finite model of a large space antenna. Modal expansions for terms in the solution algorithms are presented, using modes from the static or associated dynamic mode. These expansions provide approximated solutions in the event that a used form analytical solution to the system boundary value problem is not available.

Weeks, C. J.↗

NASA NTRS1983

An experimental study of transonic flow about a supercritical airfoil

A series of experiments was conducted on flow fields about two airfoil models whose sections are slight modifications of the original Whitcomb supercritical airfoil section. Data obtained include surface static-pressure distributions, far-wake surveys, oil-flow photographs, pitot-pressure surveys in the viscous regions, and holographic interferograms. These data were obtained for different combinations of lift coefficient and free-stream Mach number, which included both subcritical cases and flows with upper-surface shock waves. The availability of both pitot-pressure data and density data from interferograms allowed determination of flow-field properties in the vicinity of the trailing edge and in the wake without recourse to any assumptions about the local static pressure. The data show that significant static-pressure gradients normal to viscous layers exist in this region, and that they persist to approximately 10% chord downstream of the trailing edge. Comparisons are made between measured boundary-layer properties and results from boundary-layer computations that employed measured static-pressure distributions, as well as comparisons between data and results of airfoil flow-field computations.

Spaid, F. W.↗

NASA NTRS1981

Static shape determination and control for a large space antenna

An integral operator approach is used to derive solutions to static shape determination and control problems associated with large space structures. Problem assumptions include a linear self-adjoint system model, observations and control forces at discrete points, and quadratic performance criteria for the comparison of estimates or control forces. Results are illustrated by simulations with a finite element model of a large space antenna. Modal expansions for terms in the solution algorithms are presented, using modes from the static or associated dynamic model. These expansions provide approximate solutions in the event that a closed form analytical solution to the system boundary value problem is not available.

Weeks, C.↗

NASA NTRS1988

Predicting Temperatures In Ball Bearings

Computer simulations speed design studies. Analyses performed in two or three dimensions. Sizes and shapes of components approximated by zones or nodes connected by gridlines. From geometric information about grids and boundary conditions, properties of bearing and lubricant materials, and information supplied by users, thermal-analysis programs generate mathematical models for thermal transport. Thermal analysis of high-speed rolling contact bearings matured so much that computerized numerical simulations replace expensive time consuming full scale experiments.

Wagner, William R.↗

NASA NTRS1991

Lithospheric extension near Lake Mead, Nevada - A model for ductile flow in the lower crust

Small variations in gravity anomalies and topographic elevation observed in areas that have undergone highly variable amounts of upper crustal thinning can be satisfactorily explained by ductile flow of lower crustal material under the proper conditions. The boundary between the unextended Colorado Plateau and a strongly extended domain in the Basin and Range Province in the Lake Mead (Nevada) region is examined. Finite element modeling of Newtonian flow and power law creep shows that flow over the length scale of the eastern Basin and Range (500) km or more) corresponding to upper crustal extension by a factor of 1.4-3 over 10 million years requires effective viscosities less than 10 to the 18th - 10 to the 20th Pa s for ductile channels 10-25 km thick. Modeling suggests that these effective viscosities may be sustained by lower crustal material deforming at laboratory-derived power law creep rates. The longer-scale flow may require elevated crustal temperatures (more than 700 C), depending on the composition and material properties assumed. Under the boundary conditions assumed in this study the linear viscous flow models yield a satisfactory approximation to deformation by power law creep. This work suggests that flow in the lower crust may be a viable mechanism for producing small variations in total crustal thickness between strongly extended and less extended regions, and thereby explaining the relative uniformity in gravity and topography between such regions.

Kruse, Sarah↗

NASA NTRS1992

The gabbro-eclogite phase transition and the elevation of mountain belts on Venus

Among the four mountain belts surrounding Lakshmi Planum, Maxwell Montes is the highest and stands up to 11 km above the mean planetary radius and 7 km above Lakshmi Planum. The bulk composition and radioactive heat production of the crust on Venus, where measured, are similar to those of terrestrial tholeiitic basalt. Because the thickness of the low-density crust may be limited by the gabbro-garnet granulite-eclogite phase transitions, the 7-11 km maximum elevation of Maxwell Montes is difficult to understand except in the unlikely situation that the crust contains a large volume of magma. A possible explanation is that the base of the crust is not in phase equilibrium. It has been suggested that under completely dry conditions, the gabbro-eclogite phase transition takes place by solid-state diffusion and may require a geologically significant time to run to completion. Solid-state diffusion is a strongly temperature-dependent process. In this paper we solve the thermal evolution of the mountain belt to attempt to constrain the depth of the gabbro-eclogite transition and thus to assess this hypothesis quantitatively. The one-dimensional heat equation is solved numerically by a finite difference approximation. The deformation of the horizontally shortening crustal and mantle portions of the thermal boundary layer is assumed to occur by pure shear, and therefore the vertical velocity is given by the product of the horizontal strain rate and depth.

Namiki, Noriyuki↗

NASA NTRS1993

On the Daubechies-based wavelet differentiation matrix

The differentiation matrix for a Daubechies-based wavelet basis is constructed and superconvergence is proven. That is, it will be proven that under the assumption of periodic boundary conditions that the differentiation matrix is accurate of order 2M, even though the approximation subspace can represent exactly only polynomials up to degree M-1, where M is the number of vanishing moments of the associated wavelet. It is illustrated that Daubechies-based wavelet methods are equivalent to finite difference methods with grid refinement in regions of the domain where small-scale structure is present.

Jameson, Leland↗

NASA NTRS1994

Parallel inhomogeneity and the Alfven resonance. 1: Open field lines

In light of a recent demonstration of the general nonexistence of a singularity at the Alfven resonance in cold, ideal, linearized magnetohydrodynamics, we examine the effect of a small density gradient parallel to uniform, open ambient magnetic field lines. To lowest order, energy deposition is quantitatively unaffected but occurs continuously over a thickened layer. This effect is illustrated in a numerical analysis of a plasma sheet boundary layer model with perfectly absorbing boundary conditions. Consequences of the results are discussed, both for the open field line approximation and for the ensuing closed field line analysis.

Hansen, P. J.↗

NASA NTRS1995

A similarity hypothesis for the two-point correlation tensor in a temporally evolving plane wake

The analysis demonstrated that the governing equations for the two-point velocity correlation tensor in the temporally evolving wake admit similarity solutions, which include the similarity solutions for the single-point moment as a special case. The resulting equations for the similarity solutions include two constants, beta and Re(sub sigma), that are ratios of three characteristic time scales of processes in the flow: a viscous time scale, a time scale characteristic of the spread rate of the flow, and a characteristic time scale of the mean strain rate. The values of these ratios depend on the initial conditions of the flow and are most likely measures of the coherent structures in the initial conditions. The occurrences of these constants in the governing equations for the similarity solutions indicates that these solutions, in general, will only be the same for two flows if these two constants are equal (and hence the coherent structures in the flows are related). The comparisons between the predictions of the similarity hypothesis and the data presented here and elsewhere indicate that the similarity solutions for the two-point correlation tensors provide a good approximation of the measures of those motions that are not significantly affected by the boundary conditions caused by the finite extent of real flows. Thus, the two-point similarity hypothesis provides a useful tool for both numerical and physical experimentalist that can be used to examine how the finite extent of real flows affect the evolution of the different scales of motion in the flow.

Ewing, D. W.↗

NASA NTRS1996

The Sensitivity of Large-Eddy Simulation to Local and Nonlocal Drag Coefficients at the Lower Boundary

It was found that the homogeneity of the surface drag coefficient plays an important role in the large scale structure of turbulence in large-eddy simulation of the convective atmospheric boundary layer. Particularly when a ground surface temperature was specified, large horizontal anisotropies occurred when the drag coefficient depended upon local velocities and heat fluxes. This was due to the formation of streamwise roll structures in the boundary layer. In reality, these structures have been found to form when shear is approximately balanced by buoyancy. The present cases, however, were highly convective. The formation was caused by particularly low values of the drag coefficient at the entrance to thermal plume structures.

Schowalter, D. G.↗

NASA NTRS1997

Lateral Shock of the R Aquarii Jet

The R Aqr jet was observed with the VLA B-configuration at two epochs separated by approximately 13.2 yr. Comparison of the resulting 6 cm continuum images show that the radio jet has undergone a lateral counterclockwise rotation of approximately 6 deg-12 deg on the plane of the sky. The model of jet parcels on independent trajectories is difficult to reconcile with these observations and leads us to consider a path-oriented jet (i.e., younger parcels follow the same path as older parcels). Comparison of the most recent radio image with a nearly contemporaneous HST/FOC ultraviolet image at approximately 2330 Angstroms suggests that the ultraviolet emission lies along the leading side of the rotating radio jet. In conjunction with a proper motion analysis of the jet material that yields empirical space-velocity and resulting acceleration-magnitude relationships as a function of distance from the central source, we evaluate the observational results in terms of a schematic model in which the jet emission consists of plane-parallel isothermal shocks along the leading edge of rotation. In such a radiating shock, the ultraviolet-emitting region is consistent with the adiabatic region in the form of a high-temperature, low-density sheath that surrounds the cooled postshock radio-emitting region. Within the context of the schematic model, we obtain the temperatures, densities, and pressures within the preshock, adiabatic, and postshock regions as a function of distance from the central source; the physical parameters so derived compare favorably to previously published estimates. We obtain a total jet mass of 3.1 x 10(exp -5) solar mass and an age of approximately 115 yr. We evaluate the model in the context of its density-boundary condition, its applicability to an episodic or quasi-continuous jet, and angular momentum considerations.

Hollis, J. M.↗

NASA NTRS1958

Stability of Cylindrical and Conical Shells of Circular Cross Section, with Simultaneous Action of Axial Compression and External Normal Pressure

We consider in this report the determination of the upper limit of critical loads in the case of simultaneous action of a compressive force, uniformly distributed over plane cross sections, and of isotropic external normal pressure on cylindrical or conical shells of circular cross section. As a starting point we use the differential equations for neutral equilibrium of conical shells which have been used for the solution of the problem of stability of conical shells under torsion and under axial compression; upon solution of the problem it is possible to satisfy all boundary conditions, in contrast to the report where no attention is paid to the fulfillment of the boundary conditions, and to the report where only part of the boundary conditions are satisfied by solution of the problem according to Galerkin's method. Approximate formulas are used for the determination of the critical external normal pressure with simultaneous action of longituninal compression. Let us note that the formulas suggested in reference 5 are not well founded and may lead, in a number of cases, to a substantial mistake in the magnitude of the critical load.

Mushtari, K. M.↗

DOE OSTI2022

MFIX-Exa: Performance prediction of multiphase energy conversion devices

MFIX-Exa targets the high-fidelity CFD-DEM model in which particles are unresolved by the fluid grid (typically using a fluid mesh approximately twice the particle diameter) but discrete particle dynamics including collisions are resolved with a simple linear-spring dashpot model. The fluid field is solved with a modern low Mach number formation of a cell-centered, nodal-pressure approximate projection method using a Godunov scheme . Physics capabilities include: complex geometries via embedded boundaries (EBs), open and closed system ideal gas equations of state, species transport and heterogeneous chemi- cal reactions. MFIX-Exa is built on the AMReX software framework, the ECP Block-Structured Adaptive Mesh Refinement (AMR) Co-Design Center, which allows the code to be portable and performant.Further integration into the ECP ecosystem includes HYPRE linear solvers and Ascent for in situ visualization. Current project efforts focus on scaling up realistic simulations to the KPP of the challenge problem: NETs 50kW chemical looping reactor, shakedown and scaling on JSLE and OLCF TDS systems and reaching for any remaining performance improvements.

Musser, Jordan↗

DOE OSTI2023

Physics-informed machine learning of redox flow battery based on a two-dimensional unit cell model

In this paper, we present a physics-informed neural network (PINN) approach for predicting the performance of an all-vanadium redox flow battery, with its physics constraints enforced by a two-dimensional (2D) mathematical model. The 2D model, which includes 6 governing equations and 24 boundary conditions, provides a detailed representation of the electrochemical reactions, mass transport and hydrodynamics occurring inside the redox flow battery. To solve the 2D model with the PINN approach, a composite neural network is employed to approximate species concentration and potentials; the input and output are normalized according to prior knowledge of the battery system; the governing equations and boundary conditions are first scaled to an order of magnitude around 1, and then further balanced with a self-weighting method. Our numerical results show that the PINN is able to predict cell voltage correctly, but the prediction of potentials shows a constant-like shift. To fix the shift, the PINN is enhanced by further constrains derived from the current collector boundary. Finally, we show that the enhanced PINN can be even further improved if a small number of labeled data is available.

25 ENERGY STORAGE↗

NASA NTRS1976

Rocket engine thrust chamber wall temperature distribution calculation and analysis

An analytical computational concept is presented which predicts the temperature profiles along a regeneratively cooled thrust chamber wall on the hot gas side and on the coolant side, and also the coolant bulk temperature profile. The computational model is based upon a coupling of the boundary layer heat transfer process with the heat transfer process through the chamber wall and the coolant flow heat absorption. The calculation is started with approximate temperature distributions for the hot gas side wall and the coolant flow. The iteration process of the computer program is terminated when the total heat transfer rates from the hot gas boundary layer to the wall and from the wall to the coolant are equal. The computer program for the integration of regenerative cooling process to a thrust chamber is kept general such that this program can be used with any boundary layer analysis computer program for temperature profile and heat transfer studies. A sample application of this concept is shown by using a boundary layer analysis program for the RL10 rocket engine thrust chamber.

Saha, H.↗

NASA NTRS1987

A model of the wall boundary layer for ducted propellers

The objective of the present study is to include a representation of a wall boundary layer in an existing finite element model of the propeller in the wind tunnel environment. The major consideration is that the new formulation should introduce only modest alterations in the numerical model and should still be capable of producing economical predictions of the radiated acoustic field. This is accomplished by using a stepped approximation in which the velocity profile is piecewise constant in layers. In the limit of infinitesimally thin layers, the velocity profile of the stepped approximation coincides with that of the continuous profile. The approach described here could also be useful in modeling the boundary layer in other duct applications, particularly in the computation of the radiated acoustic field for sources contained in a duct.

Eversman, Walter↗

NASA NTRS1975

A numerical method for solving the Navier-Stokes equations with application to shock-boundary layer interactions

A numerical method for solving the compressible form of the unsteady Navier-Stokes equations is described. This method was originally presented in 1970 and has since been modified during the development of computer programs at Ames for implementing models that account for the effects of turbulence in shock-induced separated flows. Although this paper does not describe the turbulence models themselves, a complete description of the basic numerical method is given with emphasis on the choice of a computational mesh for high Reynolds number flows, finite-difference approximations for mixed partial derivatives, extension of the Courant-Friedrichs-Lewy stability condition for viscous flows, mesh boundary conditions, and numerical smoothing for strong shock-wave calculations.

Maccormack, R. W.↗