Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “periodic 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.

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

At least 19 records

Irregular, long-period boundary oscillations beyond approximately 100 R(sub e): GEOTAIL plasma observations

Near the tail boundary beyond about 100 Re, GEOTAIL often measures irregular, long-period oscillations in plasma velocity and density. Flow speed and density oscillate between magnetosheath values and values an order of magnitude less. The oscillations can persist for days. A typical oscillation lasts 100 minutes, but the range is large. The oscillations are highly asymmetric in that the increasing phase of the oscillation is an order of magnitude faster than the decreasing phase. This asymmetry shows that they are a distinct class of oscillations, not previously explicitly reported, and that they are not mere consequences of tail flapping in a variable solar wind. The changes in flow direction through an oscillation imply that the oscillation results from a motion of the boundary toward and away from the spacecraft with an amplitude between 5 and 10 R(sub e). A consideration of options suggests that the most plausible cause of these oscillations is the 'breathing' of the magnetotail that attends the substorm cycle.

Siscoe, G. L.↗

Heat transfer in a periodic boundary layer near an axisymmetric stagnation point on a circular cylinder

An analysis is presented to investigate the time-mean characteristics of the laminar boundary layer near an axisymmetric stagnation point when the velocity of the oncoming flow relative to the body oscillates. Different solutions are obtained for the small and high values of the reduced frequency parameter. Numerical solutions for the velocity and temperature functions are presented and the wall values of the velocity gradients are tabulated.

Gorla, R. S. R.↗

GRID3C: Computer program for generation of C type multilevel, three dimensional and boundary conforming periodic grids

A fast computer program, GRID3C, was developed for accurately generating periodic, boundary conforming, three dimensional, consecutively refined computational grids applicable to realistic axial turbomachinery geometries. The method is based on using two functions to generate two dimensional grids on a number of coaxial axisymmetric surfaces positioned between the centerbody and the outer radial boundary. These boundary fitted grids are of the C type and are characterized by quasi-orthogonality and geometric periodicity. The built in nonorthogonal coordinate stretchings and shearings cause the grid clustering in the regions of interest. The stretching parameters are part of the input to GRID3C. In its present version GRID3C can generate and store a maximum of four consecutively refined three dimensional grids. The output grid coordinates can be calculated either in the Cartesian or in the cylindrical coordinate system.

Dulikravich, D. S.↗

Fast generation of three-dimensional computational boundary-conforming periodic grids of C-type

A fast computer program, GRID3C, was developed to generate multilevel three dimensional, C type, periodic, boundary conforming grids for the calculation of realistic turbomachinery and propeller flow fields. The technique is based on two analytic functions that conformally map a cascade of semi-infinite slits to a cascade of doubly infinite strips on different Riemann sheets. Up to four consecutively refined three dimensional grids are automatically generated and permanently stored on four different computer tapes. Grid nonorthogonality is introduced by a separate coordinate shearing and stretching performed in each of three coordinate directions. The grids are easily clustered closer to the blade surface, the trailing and leading edges and the hub or shroud regions by changing appropriate input parameters. Hub and duct (or outer free boundary) have different axisymmetric shapes. A vortex sheet of arbitrary thickness emanating smoothly from the blade trailing edge is generated automatically by GRID3C. Blade cross sectional shape, chord length, twist angle, sweep angle, and dihedral angle can vary in an arbitrary smooth fashion in the spanwise direction.

Dulikravich, D. S.↗

Four-Dimensional Oceanic and Atmosperic Data Assimilation with Tropical Rainfall Measuring Mission Data

An oceanic data assimilation system which allows to utilize the forthcoming Tropical Rainfall Measuring Mission (TRMM) data has been developed and applied to the Pacific Ocean to produce the velocity field. The assimilated data will be indispensable to examine the effects of rainfall and its variability on the structure and circulation of the tropical oceans and to assess the impact of global warming due to the increase of carbon dioxide on the ocean circulation system and the marine pollution caused by oil spill and ocean damping of radionuclide. The data will also provide the verification for the oceanic and ocean-atmosphere coupled General Circulation Models (GCM's). The system consists of oceanic GCM, analysis scheme and data. In the system the flow field has been determined to be physically consistent with the observed density field and the sea surface winds derived from the Special Sensor Microwave Imagery (SSM/I) data which drive the ocean current. The time integration has been performed for five years until the flow field near the surface attained the steady state starting from the rest ocean with observed temperature and salinity fields, and the SSM/I surface wind velocity. The resultant flow field showed high producibility of the system. Especially the flow near the ocean surface agreed well with available observed data. The system, for the first time, succeeded to produce the eastward subtropical current which has been discovered in the joint investigation on Kuroshio current (CSK) in the 1960s. To verify the quality of the flow field a trajectory analysis has been carried out and compared with the Algos buoy data. BRIEF DESCRIPTION OF THE DATA ASSIMILATION SYSTEM ## Oceanic GCM and analysis scheme--The basic equations are much the same as used for the GCM's, except for the Newtonian damping terms introduced into the prediction equations for the potential temperature and salinity to maintain these fields as observed. The C grid of 2'lat. by 2'long. in horizontal and the 11 vertical levels are applied to the entire Pacific Ocean. At the east and west ocean boundaries the periodic boundary conditions are applied creating fictitious ocean there. The SMAC Method is used to increase the accuracy of mass conservation. * Data--The JODC temperature and salinity data obtained from 1906 to 1988 are used in the system between Long.100'E. and 60'W. The surface wind data are derived from the SSM/I data by Dr-R. Atlas of NASA/GSFC. The data set contains every 6 hours data from July 1987 to June 1989 on the grid of 2'lat. by 2.5'long. The averaged for the whole period and then interpolated into the 2'lat. by 2'long. grid data are used to force the system. The sea bottom topography data was based on the General Bathymetric Chart of the Ocean (GEBCO) supplied by the Canadian Hydrographic Service under contract with the International Hydrographic Organization and International Oceanographic Commission of UNESCO.

Takano, Kenji↗

Simulation of turbulent shear flows at Stanford and NASA-Ames - What can we do and what have we learned?

The capabilities and limitations of large eddy simulation (LES) and full turbulence simulation (FTS) are outlined. It is pointed out that LES, although limited at the present time by the need for periodic boundary conditions, produces large-scale flow behavior in general agreement with experiments. What is more, FTS computations produce small-scale behavior that is consistent with available experiments. The importance of the development work being done on the National Aerodynamic Simulator is emphasized. Studies at present are limited to situations in which periodic boundary conditions can be applied on boundaries of the computational domain where the flow is turbulent.

Reynolds, W. C.↗

Molecular Modeling for Calculation of Mechanical Properties of Epoxies with Moisture Ingress

Atomistic models of epoxy structures were built in order to assess the effect of crosslink degree, moisture content and temperature on the calculated properties of a typical representative generic epoxy. Each atomistic model had approximately 7000 atoms and was contained within a periodic boundary condition cell with edge lengths of about 4 nm. Four atomistic models were built with a range of crosslink degree and moisture content. Each of these structures was simulated at three temperatures: 300 K, 350 K, and 400 K. Elastic constants were calculated for these structures by monitoring the stress tensor as a function of applied strain deformations to the periodic boundary conditions. The mechanical properties showed reasonably consistent behavior with respect to these parameters. The moduli decreased with decreasing crosslink degree with increasing temperature. The moduli generally decreased with increasing moisture content, although this effect was not as consistent as that seen for temperature and crosslink degree.

Clancy, Thomas C.↗

Direct numerical simulation of instabilities in parallel flow with spherical roughness elements

Results from a direct numerical simulation of laminar flow over a flat surface with spherical roughness elements using a spectral-element method are given. The numerical simulation approximates roughness as a cellular pattern of identical spheres protruding from a smooth wall. Periodic boundary conditions on the domain's horizontal faces simulate an infinite array of roughness elements extending in the streamwise and spanwise directions, which implies the parallel-flow assumption, and results in a closed domain. A body force, designed to yield the horizontal Blasius velocity in the absence of roughness, sustains the flow. Instabilities above a critical Reynolds number reveal negligible oscillations in the recirculation regions behind each sphere and in the free stream, high-amplitude oscillations in the layer directly above the spheres, and a mean profile with an inflection point near the sphere's crest. The inflection point yields an unstable layer above the roughness (where U''(y) is less than 0) and a stable region within the roughness (where U''(y) is greater than 0). Evidently, the instability begins when the low-momentum or wake region behind an element, being the region most affected by disturbances (purely numerical in this case), goes unstable and moves. In compressible flow with periodic boundaries, this motion sends disturbances to all regions of the domain. In the unstable layer just above the inflection point, the disturbances grow while being carried downstream with a propagation speed equal to the local mean velocity; they do not grow amid the low energy region near the roughness patch. The most amplified disturbance eventually arrives at the next roughness element downstream, perturbing its wake and inducing a global response at a frequency governed by the streamwise spacing between spheres and the mean velocity of the most amplified layer.

Deanna, R. G.↗

On the spline-based wavelet differentiation matrix

The differentiation matrix for a spline-based wavelet basis is constructed. Given an n-th order spline basis it is proved that the differentiation matrix is accurate of order 2n + 2 when periodic boundary conditions are assumed. This high accuracy, or superconvergence, is lost when the boundary conditions are no longer periodic. Furthermore, it is shown that spline-based bases generate a class of compact finite difference schemes.

Jameson, Leland↗

On the Generation of Turbulent Inflow Conditions for Boundary Layer Simulations

Turbulent flows that exhibit inhomogeneities in the streamwise direction pose a particular challenge to numerical simulation approaches due to the need to prescribe time-dependent turbulent inflow conditions. In most cases the flow downstream is more or less 'driven' by the conditions at the inlet, making it necessary to specify realistic turbulent fluctuations that are in equilibrium with the assumed mean flow. This requirement often dictates that the inflow data should satisfy the Navier-Stokes equations, which in turn implies that an independent simulation be used to generate them. Detailed simulations for the purpose of creating inflow conditions can be costly and thus certain levels of approximation are desirable. In this paper we shall focus on an approximate yet accurate method for generating inflow conditions for spatially-developing boundary layer simulations. The proposed method is essentially a simplification of the method of Spalart and Leonard (1985), who devised an ingenious transformation that allows for the calculation of spatially evolving boundary layers in conjunction with periodic boundary conditions applied in the streamwise direction. While this method is elegant and highly accurate, it is more complicated than is necessary for the purpose of generating inflow data. A few key approximations are used in this work to arrive at a 'modified Spalart method' that is very easy to implement and efficient to use. The new method is shown to yield results that compare well with the computations of Spalart (1988). When used as a means of generating inflow data, the modified Spalart method is shown to be superior to existing approaches.

Lund, T. S.↗

Model Predictive Optimal Control of a Time-Delay Distributed-Parameter Systems

This paper presents an optimal control method for a class of distributed-parameter systems governed by first order, quasilinear hyperbolic partial differential equations that arise in many physical systems. Such systems are characterized by time delays since information is transported from one state to another by wave propagation. A general closed-loop hyperbolic transport model is controlled by a boundary control embedded in a periodic boundary condition. The boundary control is subject to a nonlinear differential equation constraint that models actuator dynamics of the system. The hyperbolic equation is thus coupled with the ordinary differential equation via the boundary condition. Optimality of this coupled system is investigated using variational principles to seek an adjoint formulation of the optimal control problem. The results are then applied to implement a model predictive control design for a wind tunnel to eliminate a transport delay effect that causes a poor Mach number regulation.

Nguyen, Nhan↗

Simulation of the turbulent Rayleigh-Benard problem using a spectral/finite difference technique

The three-dimensional, incompressible Navier-Stokes and energy equations with the Bousinesq assumption have been directly simulated at a Rayleigh number of 3.8 x 10 to the 5th power and a Prandtl number of 0.76. In the vertical direction, wall boundaries were used and in the horizontal, periodic boundary conditions were used. A spectral/finite difference numerical method was used to simulate the flow. The flow at these conditions is turbulent and a sufficiently fine mesh was used to capture all relevant flow scales. The results of the simulation are compared to experimental data to justify the conclusion that the small scale motion is adequately resolved.

Eidson, T. M.↗

High Fidelity CFD Analysis and Validation of Rotorcraft Gearbox Aerodynamics Under Operational and Oil-Out Conditions

This document represents the evolving formal documentation of the NPHASE-PSU computer code. Version 3.15 is being delivered along with the software to NASA in 2013.Significant upgrades to the NPHASE-PSU have been made since the first delivery of draft documentation to DARPA and USNRC in 2006. These include a much lighter, faster and memory efficient face based front end, support for arbitrary polyhedra in front end, flow-solver and back-end, a generalized homogeneous multiphase capability, and several two-fluid modelling and algorithmic elements. Specific capability installed for the NASA Gearbox Windage Aerodynamics NRA are included in this version: Hybrid Immersed Overset Boundary Method (HOIBM) [Noack et. al (2009)] Periodic boundary conditions for multiple frames of reference, Fully generalized immersed boundary method, Fully generalized conjugate heat transfer, Droplet deposition, bouncing, splashing models, and, Film transport and breakup.

Gears↗

Filament Wound Composite Analysis Using the NASA Multiscale Analysis Tool (NASMAT) and Finite Element Analysis

Fiber reinforced composite materials, owing to their tailorable thermomechanical and functional properties, allow one to produce a structure that is stronger, stiffer, and lighter than its metal counterpart while performing the same function, yielding a more efficient structure. This not only allows for the improvement of current technologies like aircraft structures, but also enables new technologies like gaseous hydrogen storage for mobility applications, which are otherwise impractical when manufactured using traditional metals due to weight and space restrictions or material embrittlement. However, the use of composites imposes greater design and manufacturing challenges on an engineer, since they are heterogenous, having a distinct structure across multiple length scale, behave generally anisotropically at the structural level and require complex manufacturing and processing methods. Capturing this complex behavior requires detailed numerical simulations, including the modeling of microstructural features like undulations, voids, and fiber alignment. In this paper, multiple repeating unit cells (RUCs), representing filament wound composites, are developed (via a script provided in the Appendix) and analyzed. The refinement of these RUCs is varied, and the analyses are performed using both the Abaqus finite element software and the NASA Multiscale Analysis Tool (NASMAT). A study is undertaken to compare the predicted effective elastic properties of the wound RUC to a laminate representation of the wound RUC, which neglects the undulations. Additionally, two different sets of periodic boundary conditions (PBCs) have been examined. One approximates the real boundary conditions using a standard approach and the other represents the PBCs exactly through the use of an offset. Lastly, a comparison of the local elastic stress fields is made among the models and approaches. Since wound structures are often approximated as laminated structures, it is important to understand the degree to which this assumption is valid, namely by first comparing the elastic constants and local elastic fields. This will provide, on the one hand, information concerning the bulk mechanical behavior and, on the other hand, insights concerning local load distributions and likely damage initiation sites.

Composites↗

On the differentiation matrix for Daubechies-based wavelets on an interval

The differentiation matrix for a Daubechies-based wavlet basis defined on an interval will be constructed. It will be shown that the differentiation matrix based on the currently available boundary constructions does not maintain the superconvergence encountered under periodic boundary conditions.

Jameson, Leland↗

Direct Numerical Simulation of Supersonic Turbulent Boundary Layers

The talk will present some initial results from the direct numerical simulation (DNS) of compressible turbulent boundary layers. We solve numerically the compressible Navier-Stokes equations using a method based on Spalart's transformation for the incompressible turbulent boundary layer. This allows the spatially developing boundary layer to be transformed to a calculation with periodic boundary conditions in the streamwise and spanwise directions. The equations are solved using Fourier expansions in the horizontal directions and B-splines in the wall-normal direction. The first simulation is at Mach 2.5 with a momentum thickness Reynolds number based on wall viscosity of R(sub theta(sup 1)) = 825. We are examining the physics of the compressible boundary layer using turbulence statistics and budget equations. The turbulence statistics include: rms (root mean square) and mean profiles, energy spectra, and two-point correlations. It is found that there are large density gradients which require significantly more resolution than the incompressible case.

Guarini, Stephen↗

Partition function for a two dimensional plasma in the random phase approximation

The partition function for a two-dimensional plasma is evaluated within the random phase approximation. The periodic boundary conditions are fully taken into account by including the periodic image interactions. In the guiding-center limit, the negative temperature threshold energy is evaluated, and a value different from previous calculations results. When an identical random phase evaluated, and a value different from previous calculations results. When an identical random phase evaluation is applied to the finite gyroradius plasma, the Salzberg-Prager-May equation of state is recovered.

Seyler, C. E., Jr.↗

Partition function for a two-dimensional plasma in the random-phase approximation

The partition function for a two-dimensional plasma is evaluated within the random-phase approximation. The periodic boundary conditions are taken fully into account by including the periodic image interactions. In the guiding-center limit, the 'negative temperature' threshold energy is evaluated, and a value different from previous calculations results. When an identical random-phase evaluation is applied to the finite-gyroradius plasma, the Salzberg-Prager-May equation of state is recovered.

Seyler, C. E., Jr.↗