Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Boundary Conditions”

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,135 records · Page 63
All sourcesNASA NTRSDOE OSTI
NASA NTRS1998

Transition in Turbine Flows

We have further developed our capabilities to analyze transition in turbine boundary layers from first principles by integrating the nonlinear parabolized stability equations (PSE) with improved initial and boundary conditions. With modified iteration schemes, we are able to proceed deeper into the transition region where skin friction coefficient and heat transfer coefficient significantly increase. Initial and boundary conditions at elevated turbulence levels can be derived by receptivity analysis. Test runs for ERCOFTAC test case T3A at 2.4\% turbulence level provide results in good agreement with the experimental data. The sharper minimum of the skin coefficient also shown by DNS results is likely due to the missing intermittency. The method has been applied to various experimentally studied turbine blades (UTRC, VKI, Zierke, Langston, Hippensteele, and others). The PSE results, though physically reasonable, do not agree as well as expected with the experimental findings. We have, therefore, performed an extensive search for the reasons of the seemingly systematic deviations. A first source of uncertainty has been found in the often insufficient documentation of the experiments (e.g. on blockage by end-wall boundary layers). However, variation of the relevant parameters does not lead to more satisfactory agreement. A second reason has been found in the "standard procedure" which considers a 2D flow at midspan and uses a panel code and subsequent boundary-layer code to obtain the laminar basic flow for the transition analysis. Comparison with the pressure distribution obtained with a 3D design code (RVC3D) shows significant three-dimensionality of the flow (e.g. in the UTRC experiments). The spanwise variation has been neglected in our original PSE code. To overcome this problem, we have developed the PSE/3D for fully 3D boundary layers to account for streamwise and spanwise variations. Since the design code does not provide the boundary-layer flow with sufficient resolution, we have generated the Euler solution and employed a 3D boundary-layer code to obtain the viscous basic flow. Although only the linear stability level of PSE/3D has been implemented so far, the discrepancies with the experiments change but do not disappear. We still find deviations between the computed and experimental variations of C(sub f), and St along the blade for laminar flow. The main reason can be seen by comparing the solution of the boundary-layer code with the viscous results of the design code. The conventional boundary-layer solution exhibits an asymptotic behavior appropriate in external aerodynamics but does not match the steep gradients of the inviscid flow through the passage and consequently provides biased results for C(sub f), and St. An attempt is currently being made to correct this deficiency. Before attempting to perform the transition analysis for the viscous flow provided by the design code, we have analyzed the implementation and "best possible" results. Code and results exhibit flaws that may negatively affect the design and are intolerable for transition analysis. Therefore, we have decided to develop a new code to obtain a reliable basis for stability and transition studies. We expect to report improved results by the time of the meeting.

Herbert, Thorwald↗
NASA NTRS1985

Application of Runge Kutta time marching scheme for the computation of transonic flows in turbomachines

Numerical solutions of the unsteady Euler equations are obtained using the classical fourth order Runge Kutta time marching scheme. This method is fully explicit and is applied to the governing equations in the finite volume, conservation law form. In order to determine the efficiency of this scheme for solving turbomachinery flows, steady blade-to-blade solutions are obtained for compressor and turbine cascades under subsonic and transonic flow conditions. Computed results are compared with other numerical methods and wind tunnel measurements. The present study also focuses on other important numerical aspects influencing the performance of the algorithm and the solution accuracy such as grid types, boundary conditions, and artificial viscosity. For this purpose, H, O, and C type computational grids as well as characteristic and extrapolation type boundary conditions are included in the solution procedure.

Subramanian, S. V.↗
NASA NTRS1981

The modification stage of basin formation - Conditions of ring formation

The impact cratering process may be divided into three stages, including a short high pressure phase of initial contact between the projectile and target, a longer cratering flow phase during which material is both ejected and displaced to form the transient crater, and a still longer modification phase during which the transient crater is modified into the final observed form. The present investigation is concerned with the nature of the modification stage as constrained by the initial boundary condition which comprises the transient crater and its surroundings, and the final boundary condition which comprises the structure of observed craters and basins. A brief description is provided of an interpretation of basin ring structure based on geophysical, geological, and dimensional analysis of craters and basins on the moon, Mercury, and Mars, on structural observations of terrestrial impact and explosion craters, and on an empirical crater restoration model. A scenario of basin formation is constructed, and relations of this scenario to observed morphology are discussed.

Croft, S. K.↗
DOE OSTI2026

An implicit-explicit time splitting strategy for the far SOL plasma fluid model with DG-FEM discretization

We consider a far scrape-off layer (SOL) plasma fluid model of ions that is governed by a Braginskiitype model: a one-dimensional, nonlinear system of advection-diffusion equations coupled with a diffusion equation for neutral particles. Our motivation for studying this system arises from the coupling between the edge plasma and radio-frequency (RF) heating, where solving a far SOL plasma fluid model provides critical insights into edge plasma dynamics. Numerical simulations of plasma fluid models require advanced computational techniques to achieve both efficiency and accuracy, especially when resolving the boundary layer in magnetically confined plasmas. In this work, we propose an implicit-explicit time operator splitting strategy that allows for an efficient solution algorithm, where the diffusive terms are treated semi-implicitly requiring only a linear solve, while the advection part is handled explicitly using a strong-stability-preserving Runge-Kutta (SSP-RK3) scheme. This leads to a fully decoupled system in which the diffusion and advection sub-problems can be solved separately, simplifying the overall solution procedure and allowing for efficient parallelization, which is particularly relevant for exploring the impact of RF heating on the SOL plasma. The main challenge of the discretization is due to the strong coupling between diffusion and advection, particularly through the boundary conditions. This makes implementation of such a scheme in an accurate and stable manner nontrivial. We discuss in detail how to split the equations and manage boundary conditions to maintain stability and well-posedness for each subsystem. We also describe a spatial discretization approach, based on the discontinuous Galerkin finite element method (DG-FEM) and present numerical results for a one-dimensional system.

Burkovska, Olena [ORNL] (ORCID:0000000163101130)↗
NASA NTRS1979

The evolution of soft X-ray-emitting flare loops

We have constructed a numerical model for a cooling flare loop in which the complete set of single-fluid equations in a one-dimensional geometry (i.e., parallel to the magnetic field) is solved. Both evaporative and static boundary conditions for the chromosphere-corona interface have been developed. This model is used to investigate the effects of initial temperature and density, loop geometry, and boundary conditions on the form of the plasma evolution and the soft X-ray emission. The results are then compared with Skylab S-056 observations of the 1973 August 9 flare. For this comparison, and under the present assumptions, we conclude that even highly compact flares must have a multiloop structure similar to large flares, and that both radiative and conductive cooling are necessary to explain the observations. The data appear to be consistent with the predicted emission from a combination of evaporative cooling loops.

Antiochos, S. K.↗
NASA NTRS1982

Lubrication of rigid ellipsida solids

The influence of geometry on the isothermal hydrodynamic film separating two rigid solids was investigated. The minimum film thickness is derived for fully flooded conjunctions by using the Reynolds boundary conditions. It was found that the minimum film thickness had the same speed, viscosity, and load dependence as Kapitza' classical solution. However, the incorporation of Reynolds boundary conditions resulted in an additional geometry effect. Solutions using the parabolic film approximation are compared by using the exact expression for the film in the analysis. Contour plots are known that indicate in detail the pressure developed between the solids.

Hamrock, B. J.↗
NASA NTRS1986

Support of long-wavelength loads on Venus and implications for internal structure

The development of a thin-shell formulation to compute stress in a planetary lithosphere is discussed. The formulation uses spherical harmonic representations of the gravity field and topography as boundary conditions to produce lateral density anomalies at two depths. The stress fields for displacement, crustal thickness anomaly, and mantle density anomaly, which are required for the boundary conditions, are calculated. The computed stresses are compared to the tectonic features of Venus. The data reveal that the highland areas are formed in a state of isostasy or uplift and are supported by dynamic processes in the mantle. Lithospheric models predict low-density mantle beneath Alta, Beta, Eisila, Thetis, and Thetus Regiones, and thick crust accompanied by higher-density mantle under Ishtar Terra and Ovda Regio.

Banerdt, W. B.↗
NASA NTRS1990

Navier-Stokes computations about complex configurations including a complete F-16 aircraft

Transonic Navier-Stokes (TNS) code solutions gathered from the literature for three-dimensional geometries, including two different wings and a complete F-16A aircraft, are presently discussed. The TNS codes use a zonal grid approach whose number of zones vary from four to 54. The Euler equations are solved in zones away from no-slip surfaces, and the thin-layer TNS equations are solved in all zones immediately adjacent to no-slip surfaces. In the case of 'corner' zones possessing no-slip boundary conditions on two different surfaces, a thin-layer formulation along both directions is employed. Employing these features, a zonal construction with the requisite set of boundary conditions can be devised for almost any application.

Holst, Terry L.↗
DOE OSTI2024

Topanga: A kinetic ion plasma code for large-scale ionospheric simulations on magnetohydrodynamic timescales

Topanga is a kinetic ion code developed for simulating large-scale plasma phenomena in the Earth's ionosphere on magnetohydrodynamic timescales. It is a domain-decomposed parallel code that runs on high-performance computing platforms. Features of Topanga include spherical geometry for simplified boundary conditions and computational efficiency; a hybrid plasma model with inertia-less fluid electrons, kinetic ions, and an electric field specified via an Ohm's law; a Maxwell-FDTD (finite difference time domain) plasma model which retains the displacement current in Maxwell's equations and models electron currents in the ionosphere with a tensor conductivity; sponge-layer boundary conditions for absorption of electromagnetic and plasma waves incident on the domain boundaries; and a novel mixed-implicit algorithm for evolving the EM fields inside the Maxwell-FDTD region that is stable over many orders of magnitude in the electron–ion collision frequency. We verify the numerical methods used in Topanga on a pair of test problems. The first test involves modeling a three-dimensional collisionless shock using the hybrid set of equations. The second test involves modeling a spherical TEM mode in vacuum using the Maxwell-FDTD set of equations. Finally, we demonstrate how using the combined set of hybrid and Maxwell-FDTD equations to model the Starfish Prime high-altitude nuclear test recovers a “missing” EM signal on the ground that is not present when using only the hybrid set of equations. The magnitude of this signal in the simulation containing the Maxwell-FDTD region agrees well with the E3a portion of the magnetohydrodynamic electromagnetic pulse from Starfish Prime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2000

Asymptotic Solutions for Optical Properties of Large Particles with Strong Absorption

For scattering calculations involving nonspherical particles such as ice crystals, we show that the transverse wave condition is not applicable to the refracted electromagnetic wave in the context of geometric optics when absorption is involved. Either the TM wave condition (i.e., where the magnetic field of the refracted wave is transverse with respect to the wave direction) or the TE wave condition (i.e., where the electric field is transverse with respect to the propagating direction of the wave) may be assumed for the refracted wave in an absorbing medium to locally satisfy the electromagnetic boundary condition in the ray tracing calculation. The wave mode assumed for the refracted wave affects both the reflection and refraction coefficients. As a result, a nonunique solution for these coefficients is derived from the electromagnetic boundary condition. In this study we have identified the appropriate solution for the Fresnel reflection/refraction coefficients in light scattering calculation based on the ray tracing technique. We present the 3 x 2 refraction or transmission matrix that completely accounts for the inhomogeneity of the refracted wave in an absorbing medium. Using the Fresnel coefficients for an absorbing medium, we derive an asymptotic solution in an analytical format for the scattering properties of a general polyhedral particle. Numerical results are presented for hexagonal plates and columns with both preferred and random orientations. The asymptotic theory can produce reasonable accuracy in the phase function calculations in the infrared window region (wavelengths near 10 micron) if the particle size (in diameter) is on the order of 40 micron or larger. However, since strong absorption is assumed in the computation of the single-scattering albedo in the asymptotic theory, the single scattering albedo does not change with variation of the particle size. As a result, the asymptotic theory can lead to substantial errors in the computation of single-scattering albedo for small and moderate particle sizes. However, from comparison of the asymptotic results with the FDTD solution, it is expected that a convergence between the FDTD results and the asymptotic theory results can be reached when the particle size approaches 200 micron. We show that the phase function at side-scattering and backscattering angles is insensitive to particle shape if the random orientation condition is assumed. However, if preferred orientations are assumed for particles, the phase function has a strong dependence on scattering azimuthal angle. The single-scattering albedo also shows very strong dependence on the inclination angle of incident radiation with respect to the rotating axis for the preferred particle orientations.

Yang, Ping↗
DOE OSTI2021

(U) A Theoretical Study of Asay Foil Trajectories (V.02)

We consider the trajectory of a generic Asay foil ejecta momentum diagnostic sensor, for scenarios where ejecta are produced at a planar surface and fly ballistically through a perfect vacuum to the sensor. To do so, we build upon a previously established mathematical formalism derived for the analytic study of stationary sensors (i.e., piezopins). First, we derive the momentum conservation equation for the problem, in a form amenable to accelerating sensors, in terms of a generic ejecta source areal mass function (“source model”). This defines an integro-differential equation (IDE) for the foil trajectory. When ejecta production is instantaneous - as is generally assumed in momentum diagnostic data analyses - the IDE leads to an implicit and easily calculable closed-form solution for the foil trajectory in a perfect system, as long as the ejecta particle velocity distribution is twice-integrable. General properties of the instant-production trajectory solution indicate the existence of a boundary condition the particle velocity distribution must satisfy in order for the analytically predicted foil trajectories to be compatible with certain features commonly observed in foil data. This boundary condition is identical to one derived previously from a consideration of piezopin data. Armed with the analytic solution for instant production, we also consider various techniques used to extract time-dependent cumulative ejecta masses from foil trajectories, and derive an expression for the error imposed by using an approximated equation of motion. This analytic trajectory solution furthermore makes it possible to examine the common practice of presenting inferred cumulative ejecta masses as a function of a normalized implied velocity; we derive conditions under which this methodology is and is not meaningful. We also propose a strategy for extending the instant-production trajectory solution to time-dependent source functions.

42 ENGINEERING↗
DOE OSTI2021

Final report of activities for the LDRD-express project #223796 titled: “Fluid models of charged species transport: numerical methods with mathematically guaranteed properties”, PI: Ignacio Tomas, Co-PI: John Shadid

This report summarizes the findings and outcomes of the LDRD-express project with title “Fluid models of charged species transport: numerical methods with mathematically guaranteed properties”. The primary motivation of this project was the computational/mathematical exploration of the ideas advanced aiming to improve the state-of-the-art on numerical methods for the one-fluid Euler-Poisson models and gain some understanding on the Euler-Maxwell model. Euler-Poisson and Euler-Maxwell, by themselves are not the most technically relevant PDE plasma-models. However, both of them are elementary building blocks of PDE-models used in actual technical applications and include most (if not all) of their mathematical difficulties. Outside the classical ideal MHD models, rigorous mathematical and numerical understanding of one-fluid models is still a quite undeveloped research area, and the treatment/understanding of boundary conditions is minimal (borderline non-existent) at this point in time. This report focuses primarily on bulk-behaviour of Euler-Poisson’s model, touching boundary conditions only tangentially.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1975

Seasonal changes in the atmospheric heat balance simulated by the GISS general circulation model

Tests of the ability of numerical general circulation models to simulate the atmosphere have focussed so far on simulations of the January climatology. These models generally present boundary conditions such as sea surface temperature, but this does not prevent testing their ability to simulate seasonal changes in atmospheric processes that accompany presented seasonal changes in boundary conditions. Experiments to simulate changes in the zonally averaged heat balance are discussed since many simplified models of climatic processes are based solely on this balance.

Stone, P. H.↗
NASA NTRS1978

Effect of geometry on hydrodynamic film thickness

The influence of geometry on the isothermal hydrodynamic film separating two rigid solids was investigated. Pressure-viscosity effects were not considered. The minimum film thickness is derived for fully flooded conjunctions by using the Reynolds boundary conditions. It was found that the minimum film thickness had the same speed, viscosity, and load dependence as Kapitza's classical solution. However, the incorporation of Reynolds boundary conditions resulted in an additional geometry effect. Solutions using the parabolic film approximation are compared with those using the exact expression for the film in the analysis. Contour plots are shown that indicate in detail the pressure developed between the solids.

Brewe, D. E.↗
NASA NTRS1992

Kolmogorov Behavior of Near-Wall Turbulence and Its Application in Turbulence Modeling

The near-wall behavior of turbulence is re-examined in a way different from that proposed by Hanjalic and Launder and followers. It is shown that at a certain distance from the wall, all energetic large eddies will reduce to Kolmogorov eddies (the smallest eddies in turbulence). All the important wall parameters, such as friction velocity, viscous length scale, and mean strain rate at the wall, are characterized by Kolmogorov microscales. According to this Kolmogorov behavior of near-wall turbulence, the turbulence quantities, such as turbulent kinetic energy, dissipation rate, etc. at the location where the large eddies become Kolmogorov eddies, can be estimated by using both direct numerical simulation (DNS) data and asymptotic analysis of near-wall turbulence. This information will provide useful boundary conditions for the turbulent transport equations. As an example, the concept is incorporated in the standard k-epsilon model which is then applied to channel and boundary flows. Using appropriate boundary conditions (based on Kolmogorov behavior of near-wall turbulence), there is no need for any wall-modification to the k-epsilon equations (including model constants). Results compare very well with the DNS and experimental data.

Shih, Tsan-Hsing↗
NASA NTRS1996

Comparative Study of Advanced Turbulence Models for Turbomachinery

A computational study has been undertaken to study the performance of advanced phenomenological turbulence models coded in a modular form to describe incompressible turbulent flow behavior in two dimensional/axisymmetric and three dimensional complex geometry. The models include a variety of two equation models (single and multi-scale k-epsilon models with different near wall treatments) and second moment algebraic and full Reynolds stress closure models. These models were systematically assessed to evaluate their performance in complex flows with rotation, curvature and separation. The models are coded as self contained modules that can be interfaced with a number of flow solvers. These modules are stand alone satellite programs that come with their own formulation, finite-volume discretization scheme, solver and boundary condition implementation. They will take as input (from any generic Navier-Stokes solver) the velocity field, grid (structured H-type grid) and computational domain specification (boundary conditions), and will deliver, depending on the model used, turbulent viscosity, or the components of the Reynolds stress tensor. There are separate 2D/axisymmetric and/or 3D decks for each module considered. The modules are tested using Rocketdyn's proprietary code REACT. The code utilizes an efficient solution procedure to solve Navier-Stokes equations in a non-orthogonal body-fitted coordinate system. The differential equations are discretized over a finite-volume grid using a non-staggered variable arrangement and an efficient solution procedure based on the SIMPLE algorithm for the velocity-pressure coupling is used. The modules developed have been interfaced and tested using finite-volume, pressure-correction CFD solvers which are widely used in the CFD community. Other solvers can also be used to test these modules since they are independently structured with their own discretization scheme and solver methodology. Many of these modules have been independently tested by Professor C.P. Chen and his group at the University of Alabama at Huntsville (UAH) by interfacing them with own flow solver (MAST).

Hadid, Ali H.↗
NASA NTRS1972

Numerical solutions of the Navier-Stokes equations for the supersonic laminar flow over a two-dimensional compression corner

Numerical solutions have been obtained for the supersonic, laminar flow over a two-dimensional compression corner. These solutions were obtained as steady-state solutions to the unsteady Navier-Stokes equations using the finite difference method of Brailovskaya, which has second-order accuracy in the spatial coordinates. Good agreement was obtained between the computed results and wall pressure distributions measured experimentally for Mach numbers of 4 and 6.06, and respective Reynolds numbers, based on free-stream conditions and the distance from the leading edge to the corner. In those calculations, as well as in others, sufficient resolution was obtained to show the streamline pattern in the separation bubble. Upstream boundary conditions to the compression corner flow were provided by numerically solving the unsteady Navier-Stokes equations for the flat plate flow field, beginning at the leading edge. The compression corner flow field was enclosed by a computational boundary with the unknown boundary conditions supplied by extrapolation from internally computed points.

Carter, J. E.↗
NASA NTRS1998

Inner Core Rotation from Geomagnetic Westward Drift and a Stationary Spherical Vortex in Earth's Core

The idea that geomagnetic westward drift indicates convective leveling of the planetary momentum gradient within Earth's core is pursued in search of a differentially rotating mean state, upon which various oscillations and secular effects might be superimposed. The desired state conforms to roughly spherical boundary conditions, minimizes dissipative interference with convective cooling in the bulk of the core, yet may aid core cooling by depositing heat in the uppermost core and lower mantle. The variational calculus of stationary dissipation applied to a spherical vortex within the core yields an interesting differential rotation profile, akin to spherical Couette flow bounded by thin Hartmann layers. Four boundary conditions are required. To concentrate shear induced dissipation near the core-mantle boundary, these are taken to be: (i) no-slip at the core-mantle interface; (ii) geomagnetically estimated bulk westward flow at the base of the core-mantle boundary layer; (iii) no-slip at the inner-outer core interface; and, to describe magnetic locking of the inner core to the deep outer core; (iv) hydrodynamically stress-free at the inner-outer core boundary. By boldly assuming the axial core angular momentum anomaly to be zero, the super-rotation of the inner core relative to the mantle is calculated to be at most 1.5 deg./yr.

Voorhies, Coerte V.↗