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.

SourceResearch connectionFrom yearThrough year

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

At least 1,027 records · Page 57

All sources NASA NTRS DOE OSTI

NASA NTRS1976

Optimal reusable-tug and expendable-kickstage trajectories for high-energy planetary missions including correction for nodal precession

Equations are derived by using branched trajectory optimization techniques and the maximum principle to maximize the payload capability of a reusable tug/expendable kickstage vehicle configuration for planetary missions. The two stages and the payload are launched into a low earth orbit by a single space shuttle. The analysis includes correction for precession of the orbit. This correction is done by the tug. The tug propels the payload and the kickstage to an energy beyond earth escape and returns within a specified time to the precessed orbit. After separating from the tug, the kickstage accelerates the payload to the required injection conditions. Planetary injection conditions are specified by the mission energy and a fixed declination and right ascension of the outgoing asymptote. The multipoint boundary value problem resulting from the analysis is solved by a Newton-Raphson iteration technique. Partial derivatives of the boundary conditions are obtained by perturbing the initial conditions one at a time, integrating the trajectory and adjoint equations, and observing the changes in boundary conditions. Maximum payload capability is derived for two typical mission energies. In addition, the variations of several mission and stage parameters are also examined.

Borsody, J.↗

DOE OSTI2025

Investigating the relationship between simulation parameters and flow variables in simulating atmospheric gravity waves for wind energy applications

Wind farms, particularly offshore clusters, are becoming larger than ever before. Besides influencing the surface wind flow and the inflow for downstream wind farms, large wind farms can trigger atmospheric gravity waves in the inversion layer and the free atmosphere aloft. Wind-farm-induced gravity waves can cause adverse pressure gradients upstream of the wind farm, which contribute to the global blockage effect, and can induce favorable pressure gradients above and downstream of the wind farm that enhance wake recovery. Numerical modeling is a powerful means of studying these wind-farm-induced atmospheric gravity waves, but it comes with the challenge of handling spurious reflections of these waves from domain boundaries. Typically, approaches which employ radiation boundary conditions and forcing zones are used to avoid these reflections. However, the simulation setup of these approaches relies heavily on ad hoc processes. For instance, the widely used Rayleigh damping method requires ad hoc tuning to produce a setup that may only produce satisfactory results for a particular case. To provide more systematic guidance on setting up realistic simulations of atmospheric gravity waves, we conduct a large-eddy simulation (LES) study of flow over a 2D hill and through a wind farm canopy that explores the optimum domain size and damping layer setup depending on the fundamental parameters which determine the flow characteristics. In this work, we only consider linearly stratified conditions (i.e., no inversion layer), thereby focusing on internal gravity waves in the free atmosphere and their reflections from the domain boundaries. This type of flow is governed by a single Froude number, which dictates most of the internal wave properties, such as wavelength, amplitude, and direction. This, in turn, will dictate the optimum domain size and Rayleigh damping layer setup. We find the effective horizontal and vertical wavelengths (the representative wavelengths of the entire wave spectrum) to be the appropriate length scales to size the domain and damping layer thickness, and the optimal Rayleigh damping coefficient scales with the Brunt–Väisälä frequency. Considering Froude numbers seen in wind farm applications, we propose recommendations to limit the reflections to less than 10 % of the total upward-propagating wave energy. Typically, damping is done at the top boundary, but given the non-periodic lateral boundary conditions of practical wind farm simulation domains, we find that damping the inflow–outflow boundaries is of equal importance to damping the top boundary. The Brunt–Väisälä frequency-normalized damping coefficient should be between 1 and 10. The damping layer thickness should be at least one effective vertical wavelength; damping layers exceeding 1.5 times the vertical wavelength are found to be unnecessary. The domain length and height should accommodate at least one effective horizontal and vertical wavelength, respectively. Moreover, Rayleigh damping does not damp the waves completely, and the non-damped energy might accumulate over the simulation time.

17 WIND ENERGY↗

NASA NTRS2015

Validation of a Node-Centered Wall Function Model for the Unstructured Flow Code FUN3D

In this paper, the implementation of two wall function models in the Reynolds averaged Navier-Stokes (RANS) computational uid dynamics (CFD) code FUN3D is described. FUN3D is a node centered method for solving the three-dimensional Navier-Stokes equations on unstructured computational grids. The first wall function model, based on the work of Knopp et al., is used in conjunction with the one-equation turbulence model of Spalart-Allmaras. The second wall function model, also based on the work of Knopp, is used in conjunction with the two-equation k-! turbulence model of Menter. The wall function models compute the wall momentum and energy flux, which are used to weakly enforce the wall velocity and pressure flux boundary conditions in the mean flow momentum and energy equations. These wall conditions are implemented in an implicit form where the contribution of the wall function model to the Jacobian are also included. The boundary conditions of the turbulence transport equations are enforced explicitly (strongly) on all solid boundaries. The use of the wall function models is demonstrated on four test cases: a at plate boundary layer, a subsonic di user, a 2D airfoil, and a 3D semi-span wing. Where possible, different near-wall viscous spacing tactics are examined. Iterative residual convergence was obtained in most cases. Solution results are compared with theoretical and experimental data for several variations of grid spacing. In general, very good comparisons with data were achieved.

Carlson, Jan-Renee↗

NASA NTRS1995

Kinetic structure of rotational discontinuities: Implications for the magnetopause

Magnetic field rotations in the high ion beta magnetosheath that are part of the magnetopause structure are expected to have only a small normal component. We have studied the properties of rotational discontinuities (RDs) under these conditions, viewed as the limit of weak intermediate shocks (ISs), by performing hybrid simulations with a reflecting wall boundary condition (piston method). With this dynamic formation, the sense and size of rotation are not arbitrarily predetermined, but rather evolve from the given upstream (magnetosheath) and downstream (magnetospheric) boundary conditions, similar to what takes place at the magnetopause. This work focuses on several aspects: the observed minimum shear of RDs, their width, their internal signature, and their relation to ISs in isotropic plasmas. Our simulation results are in agreement with the minimum shear observations, that is, the RDs choose the sense of rotation that corresponds to the minimum angle between the upstream and downstream field vector. The RDs are stable, with a unique scale size. Typical gradient scale half widths are one to four ion inertial lengths with a total width up to ten times of that, in agreement with magnetopause observations. We develop a generalized fluid theory of RDs and discuss the characteristic internal signatures of the rotational layer, comparing the kinetic simulation results to predictions from the generalized fluid theory. The results show that ion inertia, anisotropic pressure, finite Larmor radius effects, nonzero ion heat flux, and reflected ions all contribute to the signatures of RDs on kinetic scales. The RDs may have upstream or downstream wave trains, which become weak for high ion beta and small normal components of the magnetic field. We explain the presence and direction of wave trains in terms of the kinetic properties of the Alfven/ion-cyclotron mode. Away from the RD limit there is a smooth transition to weak intermediate shocks, which have small jumps close to expected Rankine-Hugoniot values. Apart from that, there are few kinetic plasma signatures that distinguish RDs from their neighboring ISs. However, noncoplanar ISs evolve in time into thin RDs. Using the properties of RDs and ISs, we make specific suggestions how these discontinuities can be distinguished observationally in the case of an isotropic plasma.

Krauss-Varban, D.↗

NASA NTRS2005

Navier-Stokes Simulation of Airconditioning Facility of a Large Modem Computer Room

NASA recently assembled one of the world's fastest operational supercomputers to meet the agency's new high performance computing needs. This large-scale system, named Columbia, consists of 20 interconnected SGI Altix 512-processor systems, for a total of 10,240 Intel Itanium-2 processors. High-fidelity CFD simulations were performed for the NASA Advanced Supercomputing (NAS) computer room at Ames Research Center. The purpose of the simulations was to assess the adequacy of the existing air handling and conditioning system and make recommendations for changes in the design of the system if needed. The simulations were performed with NASA's OVERFLOW-2 CFD code which utilizes overset structured grids. A new set of boundary conditions were developed and added to the flow solver for modeling the roomls air-conditioning and proper cooling of the equipment. Boundary condition parameters for the flow solver are based on cooler CFM (flow rate) ratings and some reasonable assumptions of flow and heat transfer data for the floor and central processing units (CPU) . The geometry modeling from blue prints and grid generation were handled by the NASA Ames software package Chimera Grid Tools (CGT). This geometric model was developed as a CGT-scripted template, which can be easily modified to accommodate any changes in shape and size of the room, locations and dimensions of the CPU racks, disk racks, coolers, power distribution units, and mass-storage system. The compute nodes are grouped in pairs of racks with an aisle in the middle. High-speed connection cables connect the racks with overhead cable trays. The cool air from the cooling units is pumped into the computer room from a sub-floor through perforated floor tiles. The CPU cooling fans draw cool air from the floor tiles, which run along the outside length of each rack, and eject warm air into the center isle between the racks. This warm air is eventually drawn into the cooling units located near the walls of the room. One major concern is that the hot air ejected to the middle isle might recirculate back into the cool rack side and cause thermal short-cycling. The simulations analyzed and addressed the following important elements of the computer room: 1) High-temperature build-up in certain regions of the room; 2) Areas of low air circulation in the room; 3) Potential short-cycling of the computer rack cooling system; 4) Effectiveness of the perforated cooling floor tiles; 5) Effect of changes in various aspects of the cooling units. Detailed flow visualization is performed to show temperature distribution, air-flow streamlines and velocities in the computer room.

Source record↗

NASA NTRS1991

Thermal and solutal convection with conduction effects inside a rectangular enclosure

We numerically investigate the effects of various boundary conditions on the flow field characteristics of the physical vapor transport process. We use a prescribed temperature profile as boundary condition on the enclosure walls, and we consider parametric variations applicable to ground-based and space microgravity conditions. For ground-based applications, density gradients in the fluid phase generate buoyancy-driven convection which in turn disrupts the uniformity of the mass flux at the interface depending on the orientation. Heat conduction in the crystal can affect the fluid flow near the interface of the crystal. When considering isothermal source and sink at the interfaces, we observe a diffusive mode and three modes (i.e., thermal, solutal, and thermo-solutal). The convective modes show opposing flow field trends between thermal and solutal convection; theoretically, these trends can be used to achieve a uniform mass flux near the crystal. However, under the physical conditions chosen, the mathematical condition necessary for uniform mass flux cannot be satisfied because of thermodynamic restrictions. When a longitudinal thermal gradient is prescribed on the boundary of the crystal, a non-uniform interface temperature results, which induces a symmetrical fluid flow near the interface for the vertical case. For space microgravity applications, we show that the flow field is dominated by the Stefan wind and a uniform mass flux results at the interface.

Mennetrier, Christophe↗

NASA NTRS1974

Numerical and experimental studies of the natural convection within a horizontal cylinder

Numerical solutions are obtained for the quasi-compressible Navier-Stokes equations governing the time-dependent natural convection within a horizontal cylinder. The early flow development and wall heat transfer are obtained after a uniformly cold wall is imposed as a boundary condition on the cylinder. Results are also obtained for a time-varying cold wall as a boundary condition with windward explicit differencing used for the numerical solutions. The viscous truncation error associated with this scheme is controlled so that first-order accuracy is maintained in time and space. Experiments within a small-scale instrumented horizontal cylinder revealed the time development of the temperature distribution across the boundary layer as well as the decay of wall heat transfer with time. Agreement between temperature distributions obtained experimentally and numerically was generally good. The time decay of the dimensionless ratio of the Nusselt number to the one-fourth power of the Grashof number is found both numerically and experimentally, and good agreement is obtained between these two results over most of the cylinder wall.

Stewart, R. B.↗

NASA NTRS1991

Optimal guidance of aero-assisted transfer vehicles based on matched asymptotic expansions

This paper addresses a number of issues related to the matched asymptotic expansions analysis of skip trajectories, or any class of problems that give rise to inner layers that are not associated directly with satisfying boundary conditions. The procedure for matching inner and outer solutions, and using the composite solution to satisfy boundary conditions is rigorously followed in developing a complete algebraic solution to the problem of inclination change with minimum energy loss. Repeated solution of these algebraic equations along the trajectory, treating each current state as an initial state, constitute a feedback guidance algorithm.

Calise, A. J.↗

NASA NTRS1994

The field of a screened magnetic dipole

The purpose of this note is to quantitatively study the asymptotic behavior of the dipole magnetic field in the tail region of a paraboloidal or cylindrical model of the magnetosphere, assuming the complete screening of the internal field by magnetopause currents. This screening assumption is equivalent to imposing the boundary condition that the normal component of the magnetic field is zero at the magnetopause. With this boundary condition, the screened dipole field falls off exponentially with distance down the tail, in sharp constrast to the bare dipole field. Analytic expressions for a cylindrical and paraboloidal magnetopause are given.

Greene, J. M.↗

NASA NTRS2003

Sheet Flows, Avalanches, and Dune Evolution on Earth and Mars

This investigation is a collaboration between researchers at Cornell University, the University of Florida, and the University of Rennes 1, France. Flow modeling at Cornell University focused on mechanisms for the suspension and transport of wind-blown sand that are important in both terrestrial and Martian environments. These mechanisms include the saltation (or jumping) of grains, collisions between grains, and the interaction of grains with the velocity fluctuations of the turbulent wind. Of particular interest are sheet flows; these are relatively thin, highly concentrated regions of grains flowing near the ground under the influence of a strong turbulent wind. In them, the grains are suspended by interparticle collisions. Sheet flows may be relatively rare events, but they have the capacity to move great amounts of sand. In order to describe sheet flows, a turbulent mixture theory was formulated for particles in a fluid in which fluctuations in the volume fiaction of the particles take place on the scale of the turbulent eddies. Ensemble averaged equations for particle and fluid mass, momentum, and energy and fluid rate of dissipation were expressed in terms of Farve (concentration) averaged velocities and the associated velocity fluctuations. Correlations that describe the turbulent suspension of particles and dissipation of turbulent energy of both phases due to fluid particle interactions were modeled and boundary conditions at the bed and at the upper surface of the collisional flow were formulated. The boundary conditions at the upper surface were tested in a numerical simulation developed at the University of Florida. Steady and unsteady solutions for steady and unsteady fully-developed flows were obtained over a range of wind speeds fiom the lowest for which collisional between particles occurred to at which turbulent suspension is found to dominate collisional suspension. Below the value of the wind speed at which collisions between particles were unimportant, numerical solutions were obtained for the velocity distribution function and the resulting fields of concentration, particle and gas mean velocity, and particle shear stress for the steady two-dimensional saltation of spherical sand particles driven by a turbulent wind over a bed characterized by a simple relationship (the splash function) between the properties of incoming particles and those of the rebounding particles and other particles ejected fiom the bed. At the University of Rennes 1, experiments devoted to the characterization of the splash function for beds consisting of either random or ordered arrays of spheres in two- dimensions were completed. These indicated the role played by the packing geometry in the rebound and ejection of grains. Preliminary experiments on response of a three- dimensional collision bed to a collision with a single particle were performed. Data was taken with a single camera focused on the plane of collision. Here, for example, the decrease of the effective coefficient of restitution of the bed with an increase of the angle of incidence of the incoming particle has been measured. Other experiments on avalanches at Rennes studied the properties of the flows of particles that are responsible for the motion of the leeward side of a dune. In these, the dependence of the initiation of avalanches on the packing and depth of the particles was measured. Particle migration was studied in inclined flows of a binary mixture of disks and the mechanisms of diffision and segregation were isolated and characterized. The influence of side wall on dense, rapid inclined flows was measured and shown to be the reason why the angle of the free surface in such flows can exceed the static angle of repose. Future research will be devoted to a better understanding the transition between saltating (collisionless) and collisional flows as the wind speed the increases. This will involve the understanding of the evolution of the splash function as clisions with the bed become more numerous, more frequent, and more violent.

Source record↗

NASA NTRS2006

Rotating Rake Mode Measurements Over Passive Treatment in a Ducted Fan

The NASA Glenn Research Center s Rotating Rake mode measurement system has been successful in measuring the modal content propagating in hardwall ducts. This paper proposes an extension of the Rotating Rake measurement and analysis technique to treated sections by developing basis functions based on wall impedance boundary conditions for flow conditions (i.e., constant duct area and Mach number) where the closed form analytical solution exists. Analytical equations developed to estimate mode power are incorporated. This method is verified by decomposing and analyzing radial pressure profiles generated numerically by the Eversman propagation code. Several modes, frequencies and impedances are evaluated. Data from a low-speed ducted fan with several different impedance conditions was acquired and reduced to determine the best fit to the data. Using the impedance boundary conditions result in better mode measurement solutions.

Sutliff, Daniel L.↗

NASA NTRSDate not supplied

Ares: A Coupling Methodology for Ablation Modeling

To enable modeling of complex and coupled ablation problems, a multi-physics framework is developed. A methodology for modeling shape change in coupled systems is presented. The approach taken to model gas-surface interactions and translate coupled surface phenomena to physically meaningful boundary conditions in the distinct solvers is discussed. Particular emphasis is placed on the nature of coupled boundary conditions pertaining to surface energy and mass balances as well as surface chemistry modeling. The developed methodology is used to simulate a shear test in arc-jet conditions in order to assess the validity of the coupled approach as well as the implementation of the relevant physical processes.

EDL↗

NASA NTRSDate not supplied

Ares: A Coupling Methodology for Ablation Modeling

EDL↗

NASA NTRS1990

Effects of shock on the stability of hypersonic boundary layers

A set of linearized shock boundary conditions is derived, which is then imposed at the shock to account for the interaction of the shock wave with the boundary/shock layer instability wave; these boundary conditions are used to study the effect of shock on hypersonic boundary layer stability under the assumption of quasi-parallel flow. The result show that the shock has little effect on the boundary layer instability (subsonic first and second mode disturbances) when the shock is located outside the boundary layer edge. When the shock is located near the boundary layer edge, it exerts a stabilizing influence on the first and second modes. The shock also induces unstable supersonic modes with oscillatory structure in the shock layer, but these modes grow slower than the subsonic modes.

Chang, Chau-Lyan↗

NASA NTRS1972

The phase discrepancy between temperature and density in the thermosphere

The observations of thermospheric densities by satellite drag analysis and the thermospheric temperatures by back scatter radar measurements show a phase difference. The dependence of this phase difference on various factors is briefly reviewed. It is shown that by choosing certain boundary conditions at 120 km the observations of both temperature and densities become consistent. The possible range of boundary conditions is analysed. Furthermore it is shown that the phase discrepancy is largely due to the higher harmonic components of both the temperature and the density variation in the thermosphere.

Blum, P. W.↗

NASA NTRS1977

Compressible laminar boundary layers with suction on swept and tapered wings

In this paper we present a numerical method for solving the compressible laminar boundary-layer equations with suction on swept and tapered wings. The method employs an efficient two-point finite-difference method to solve the governing equations, and a very convenient similarity transformation which removes the wall normal velocity as a boundary condition and places it into the governing equations as a parameter. In this way the awkward nonlinear boundary condition which couples all the variables is avoided. To test and demonstrate the method, we present a sample calculation for a typical laminar-flow-control (LFC) wing.

Kaups, K.↗

NASA NTRS1984

Approximate-factorization algorithms - Theory and applications in viscous-flow computations

A systematic development of implicit approximate-factorization algorithms in delta form for both unsteady and steady viscous flow is presented. The algorithms are cast in conservation-law form and simplified by using a thin-layer approximation to the governing equations. The implementation of implicit surface viscous boundary conditions is discussed in detail, and an example is presented illustrating the advantage of using the implicit boundary conditions. Three-dimensional results from the steady form of the algorithm are presented and compared with experimental data.

Chaussee, D. S.↗

NASA NTRS1993

Numerical modeling of the interstellar medium in galactic disks

We have been developing detailed hydrodynamic models of the global interstellar medium in the hope of understanding the mass and volume occupied by various phases, as well as their structure and kinematics. In our model, the gas is modeled by one fluid while representative Pop 1 stars are modeled by a second fluid. The two fluids are coupled in that the gas forms into stars at a rate given by a Schmidt law while stellar mass loss returns matter into the gas phase (on a time scale of 100 Myr). Also, the stars heat the gas through stellar winds and the gas cools through optically thin radiation. The time behavior of these two fluids is studied in two spatial dimensions with the Eulerian finite difference numerical hydrodynamic code Zen. The two spatial dimensions are along the plane of a disk (x, total length of 2 kpc) and perpendicular to the disk (z, total height of +/- 15 kpc) and a galactic gravitational field in the z direction, typical of that at the solar circle, is imposed upon the simulation; self-gravity and rotation are absent. For the boundary conditions, outflow is permitted at the top and bottom of the grid (z = +/- 15 kpc) while periodic boundary conditions are imposed upon left and right sides of the grid. As initial conditions, we assumed a gaseous distribution like that seen for the H1 by earlier researchers, although the results are insensitive to the initial conditions. We have run simulations in which the heating due to stars, parameterized as a stellar wind velocity, a, is varied from low (a = 150 km/s), to intermediate (a = 300 km/s), to high (a = 600 km/s). Since the intermediate case is roughly equivalent to the Galactic energy injection rate from supernovae, this summary will concentrate on results from this simulation.

Rosen, A.↗