Search NASA⌕ Search

SEARCH · Search NASA

Results for “Periodic 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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

Validation of the SINDA/FLUINT code using several analytical solutions

The Systems Improved Numerical Differencing Analyzer and Fluid Integrator (SINDA/FLUINT) code has often been used to determine the transient and steady-state response of various thermal and fluid flow networks. While this code is an often used design and analysis tool, the validation of this program has been limited to a few simple studies. For the current study, the SINDA/FLUINT code was compared to four different analytical solutions. The thermal analyzer portion of the code (conduction and radiative heat transfer, SINDA portion) was first compared to two separate solutions. The first comparison examined a semi-infinite slab with a periodic surface temperature boundary condition. Next, a small, uniform temperature object (lumped capacitance) was allowed to radiate to a fixed temperature sink. The fluid portion of the code (FLUINT) was also compared to two different analytical solutions. The first study examined a tank filling process by an ideal gas in which there is both control volume work and heat transfer. The final comparison considered the flow in a pipe joining two infinite reservoirs of pressure. The results of all these studies showed that for the situations examined here, the SINDA/FLUINT code was able to match the results of the analytical solutions.

Keller, John R.↗
NASA NTRS2004

BASS Code Development

The BASS computational aeroacoustic code solves the fully nonlinear Euler equations in the time domain in two-dimensions. The acoustic response of the stator is determined simultaneously for the first three harmonics of the convected vortical gust of the rotor. The spatial mode generation, propagation and decay characteristics are predicted by assuming the acoustic field away from the stator can be represented as a uniform flow with small harmonic perturbations superimposed. The computed field is then decomposed using a joint temporal-spatial transform to determine the wave amplitudes as a function of rotor harmonic and spatial mode order. This report details the following technical aspects of the computations and analysis. 1) the BASS computational technique; 2) the application of periodic time shifted boundary conditions; 3) the linear theory aspects unique to rotor-stator interactions; and 4) the joint spatial-temporal transform. The computational results presented herein are twofold. In each case, the acoustic response of the stator is determined simultaneously for the first three harmonics of the convected vortical gust of the rotor. The fan under consideration here like modern fans is cut-off at +, and propagating acoustic waves are only expected at 2BPF and 3BPF. In the first case, the computations showed excellent agreement with linear theory predictions. The frequency and spatial mode order of acoustic field was computed and found consistent with linear theory. Further, the propagation of the generated modes was also correctly predicted. The upstream going waves propagated from the domain without reflection from the in ow boundary. However, reflections from the out ow boundary were noticed. The amplitude of the reflected wave was approximately 5% of the incident wave. The second set of computations were used to determine the influence of steady loading on the generated noise. Toward this end, the acoustic response was determined with three steady loading conditions: design, low-loading, high-loading. The overall trend showed significant (approximately 10 dB) increases in the generated noise for the highly loaded stator.

Sawyer, Scott↗
NASA NTRS1983

Wave scattering and guidance by dielectric waveguides with periodic surfaces

The extended-boundary-condition method is used to solve the wave scattering and guidance by a stratified medium with a periodic surface. First, the transition matrix for two sets of Flocquet waves incident from above and below a periodic surface separating the two media is derived. The results are applied to the diffraction of an incident plane wave by a multilayered medium having a periodic surface. The theoretical results can easily be applied to the guidance problem for which the complex guiding constant of the periodic structure is obtained. Numerical results for both the scattering and the guidance problems are obtained. The theory is compared with experiment for metallic gratings with surface-plasmon excitations. The calculated guiding constants of a thin film with a periodic sinusoidal surface are also compared with other numerical methods.

Chuang, S. L.↗
NASA NTRS1993

A new conformal absorbing boundary condition for finite element meshes and parallelization of FEMATS

Some of the progress toward the development and parallelization of an improved version of the finite element code FEMATS is described. This is a finite element code for computing the scattering by arbitrarily shaped three dimensional surfaces composite scatterers. The following tasks were worked on during the report period: (1) new absorbing boundary conditions (ABC's) for truncating the finite element mesh; (2) mixed mesh termination schemes; (3) hierarchical elements and multigridding; (4) parallelization; and (5) various modeling enhancements (antenna feeds, anisotropy, and higher order GIBC).

Chatterjee, A.↗
DOE OSTI2020

Axial magnetic flux amplification in Hall-magnetohydrodynamic simulations of externally magnetized z-pinches

It is shown that externally magnetized gas-puff z-pinches under compression can exhibit robust axial flux amplification under suitable axial boundary conditions. This effect relies upon the Hall term in the generalized Ohm's law to generate azimuthal currents in the presence of a driving axial electric field. Under dynamic compression, the total current tends to flow in a mostly force-free boundary layer, separating the predominately azimuthal field outside the layer and the predominately axial field inside the layer. The effect only appears to occur if the axial boundaries allow for outflow or absorption of electromagnetic energy. The effect is mitigated by imposing either periodic or conducting axial boundary conditions. A semi-analytic equilibrium analysis agrees with steady-state solutions of the time-dependent electron-magnetohydrodynamic equations and provides an estimate of the scaling of the boundary layer as well as suggesting a scenario for the formation of the boundary layer. When operative, the effect can significantly impede plasma compression due to the increase in axial flux that diffuses into the pre-compressed plasma or through the presence of conductors that inhibit movement of the generated axial flux. Several facilities have noted unusual implosion behavior in z-pinch experiments with applied axial fields that does not appear to be explainable within the standard magnetohydrodynamic model. In conclusion, it is suggested that these experiments can be explained by the axial flux amplification and concentration effect

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2005

Intermittent Behavior of the Separated Boundary Layer along the Suction Surface of a Low Pressure Turbine Blade under Periodic Unsteady Flow Conditions

The paper experimentally and theoretically studies the effects of periodic unsteady wake flow and aerodynamic characteristics on boundary layer development, separation and re-attachment along the suction surface of a low pressure turbine blade. The experiments were carried out at Reynolds number of 110,000 (based on suction surface length and exit velocity). For one steady and two different unsteady inlet flow conditions with the corresponding passing frequencies, intermittency behaviors were experimentally and theoretically investigated. The current investigation attempts to extend the intermittency unsteady boundary layer transition model developed in previously to the LPT cases, where separation occurs on the suction surface at a low Reynolds number. The results of the unsteady boundary layer measurements and the intermittency analysis were presented in the ensemble-averaged and contour plot forms. The analysis of the boundary layer experimental data with the flow separation, confirms the universal character of the relative intermittency function which is described by a Gausssian function.

Oeztuerk, B↗
NASA NTRS2009

Numerical Study of Outlet Boundary Conditions for Unsteady Turbulent Internal Flows Using the NCC

This paper presents the results of studies on the outlet boundary conditions for turbulent internal flow simulations. Several outlet boundary conditions have been investigated by applying the National Combustion Code (NCC) to the configuration of a LM6000 single injector flame tube. First of all, very large eddy simulations (VLES) have been performed using the partially resolved numerical simulation (PRNS) approach, in which both the nonlinear and linear dynamic subscale models were employed. Secondly, unsteady Reynolds averaged Navier- Stokes (URANS) simulations have also been performed for the same configuration to investigate the effects of different outlet boundary conditions in the context of URANS. Thirdly, the possible role of the initial condition is inspected by using three different initial flow fields for both the PRNS/VLES simulation and the URANS simulation. The same grid is used for all the simulations and the number of mesh element is about 0.5 million. The main purpose of this study is to examine the long-time behavior of the solution as determined by the imposed outlet boundary conditions. For a particular simulation to be considered as successful under the given initial and boundary conditions, the solution must be sustainable in a physically meaningful manner over a sufficiently long period of time. The commonly used outlet boundary condition for steady Reynolds averaged Navier-Stokes (RANS) simulation is a fixed pressure at the outlet with all the other dependent variables being extrapolated from the interior. The results of the present study suggest that this is also workable for the URANS simulation of the LM6000 injector flame tube. However, it does not work for the PRNS/VLES simulation due to the unphysical reflections of the pressure disturbances at the outlet boundary. This undesirable situation can be practically alleviated by applying a simple unsteady convection equation for the pressure disturbances at the outlet boundary. The numerical results presented in this paper suggest that this unsteady convection of pressure disturbances at the outlet works very well for all the unsteady simulations (both PRNS/VLES and URANS) of the LM6000 single injector flame tube.

Liu, Nan-Suey↗
DOE OSTI2023

Technical descriptions of the experimental dynamical downscaling simulations over North America by the CAM–MPAS variable-resolution model

Abstract. Comprehensive assessment of climate datasets is important for communicating model projections and associated uncertainties to stakeholders. Uncertainties can arise not only from assumptions and biases within the model but also from external factors such as computational constraint and data processing. To understand sources of uncertainties in global variable-resolution (VR) dynamical downscaling, we produced a regional climate dataset using the Model for Prediction Across Scales (MPAS; dynamical core version 4.0) coupled to the Community Atmosphere Model (CAM; version 5.4), which we refer to as CAM–MPAS hereafter. This document provides technical details of the model configuration, simulations, computational requirements, post-processing, and data archive of the experimental CAM–MPAS downscaling data. The CAM–MPAS model is configured with VR meshes featuring higher resolutions over North America as well as quasi-uniform-resolution meshes across the globe. The dataset includes multiple uniform- (240 and 120 km) and variable-resolution (50–200, 25–100, and 12–46 km) simulations for both the present-day (1990–2010) and future (2080–2100) periods, closely following the protocol of the North American Coordinated Regional Climate Downscaling Experiment. A deviation from the protocol is the pseudo-warming experiment for the future period, using the ocean boundary conditions produced by adding the sea surface temperature and sea-ice changes from the low-resolution version of the Max Planck Institute Earth System Model (MPI-ESM-LR) in the Coupled Model Intercomparison Project Phase 5 to the present-day ocean state from a reanalysis product. Some unique aspects of global VR models are evaluated to provide background knowledge to data users and to explore good practices for modelers who use VR models for regional downscaling. In the coarse-resolution domain, strong resolution sensitivity of the hydrological cycles exists over the tropics but does not appear to affect the midlatitude circulations in the Northern Hemisphere, including the downscaling target of North America. The pseudo-warming experiment leads to similar responses of large-scale circulations to the imposed radiative and boundary forcings in the CAM–MPAS and MPI-ESM-LR models, but their climatological states in the historical period differ over various regions, including North America. Such differences are carried to the future period, suggesting the importance of the base state climatology. Within the refined domain, precipitation statistics improve with higher resolutions, and such statistical inference is verified to be negligibly influenced by horizontal remapping during post-processing. Limited (≈50 % slower) throughput of the current code is found on a recent many-core/wide-vector high-performance computing system, which limits the lengths of the 12–46 km simulations and indirectly affects sampling uncertainty. Our experience shows that global and technical aspects of the VR downscaling framework require further investigations to reduce uncertainties for regional climate projection.

58 GEOSCIENCES↗
NASA NTRS1975

Nonlinear vibration of cylindrical shells

The large amplitude vibration of a thin-walled cylindrical shell is analyzed. A perturbation method is used to solve the steady-state forced vibration problem. The simply-supported boundary conditions and the circumferential periodicity condition are satisfied. The resulting solution indicates that in addition to the fundamental mode, the response contains asymmetric modes as well as axisymmetric modes with the frequency twice that of the fundamental mode. Vibrations involving a single driven mode response are investigated. The results indicate that the nonlinearity is either softening or hardening depending on the mode. The vibrations involving both a driven mode and a companion mode are also investigated. An experimental investigation is conducted, and the results are in qualitative agreement with the theory. 'Nonstationary' response is detected at some frequencies for large amplitude response in which the amplitude drifts from one value to another. Various nonlinear phenomena are observed and comparisons with the theoretical results are made.

Chen, J. C.↗
NASA NTRS1978

Numerical experiments on the climatic sensitivity of an atmospheric hydrologic cycle

It is shown for an intermediate numerical model that fractional cloudiness and relative humidity decrease with increasing temperature. The fractional cloudiness decreases at a rate about 1 per deg K. This occurs in spite of an increase in the evaporation, water transport, condensation, precipitation and cloud water content with increasing temperature. These results are quite similar to those found from models with more highly parameterized clouds, notably the NCAR model. The fractional cloudiness in this model is measured by the fractional coverage of total cloud water and the fractional coverage of positive condensation, in addition to the relative humidity. It is also shown that some of the characteristics of a temperate climate can be simulated in an intermediate numerical model with periodic, antisymmetric and symmetric boundary conditions on an f plane. Intermediate models of this sort may therefore be useful to investigate general questions about the earth's hydrologic cycle on climatic space and time scales

Roads, J. O.↗
NASA NTRS1985

Direct numerical simulation of equilibrium turbulent boundary layers

This paper describes the simulation of turbulent boundary layers by direct numerical solution of the three-dimensional, time-dependent Navier-Stokes equations, using a spectral method. The flow is incompressible, with Re sub delta = 1000. The equations are written in the self-similar coordinate system and periodic streamwise and spanwise boundary conditions are imposed. A family of nine 'equilibrium' boundary layers, from the strongly accelerated 'sink' flow to Stratford's separating boundary layer is treated. Good general agreement with experiments is observed. The effects of pressure gradients on the structures and statistics, both in the wall and wake regions are discussed.

Spalart, P. R.↗
NASA NTRS1992

The program FANS-3D (finite analytic numerical simulation 3-dimensional) and its applications

In this study, the program named FANS-3D (Finite Analytic Numerical Simulation-3 Dimensional) is presented. FANS-3D was designed to solve problems of incompressible fluid flow and combined modes of heat transfer. It solves problems with conduction and convection modes of heat transfer in laminar flow, with provisions for radiation and turbulent flows. It can solve singular or conjugate modes of heat transfer. It also solves problems in natural convection, using the Boussinesq approximation. FANS-3D was designed to solve heat transfer problems inside one, two and three dimensional geometries that can be represented by orthogonal planes in a Cartesian coordinate system. It can solve internal and external flows using appropriate boundary conditions such as symmetric, periodic and user specified.

Bravo, Ramiro H.↗
NASA NTRS1993

Simulated rarefied aerodynamics of the Magellan spacecraft during aerobraking

Aerodynamic loads upon the Magellan spacecraft during aerobraking through the atmosphere of Venus are computed at off-design attitudes with a direct simulation Monte Carlo (DSMC) particle method. Simulated rarefied flows at nominal altitudes near 140 km and an entry speed of 8.6 km/s were compared to simulated and analytic free molecular results. Aerodynamic moments, forces, and heating for rarefied entry at all attitudes were 7-10 percent below free molecular results. All moments acted to restore the vehicle to its nominal zero-pitch, zero-yaw attitude. Suggested canting of the solar panels is an innovative configuration to assess gas-surface interaction during aerobraking. The resulting roll torques about the central body-axis as predicted in rarefied flow simulations were nearly twice that predicted for free molecular flow, although differences became less distinct for thermal accommodation coefficients well below unity. Roll torques increased dramatically with reduced accommodation coefficients employed in the simulation. In the DSMC code, periodic free-molecule boundary conditions and a coarse computational grid and body resolution served to minimize the simulation size and cost while retaining solution validity.

Haas, Brian L.↗
NASA NTRS2006

Using global aerosol models and satellite data for air quality studies: Challenges and data needs

Aerosol particles, also known as PM2.5 (particle diameter less than 2.5 pm) and PM10 (particle diameter less than 10 pm), are one of the key atmospheric components that determines air quality. Yet, air quality forecasts for PM are still in their infancy and remain a challenging task. It is difficult to simply relate PM levels to local meteorological conditions, and large uncertainties exist in regional air quality model emission inventories and initial and boundary conditions. Especially challenging are periods when a significant amount of aerosol comes from outside the regional modeling domain through long-range transport. In the past few years, NASA has launched several satellites with global aerosol measurement capabilities, providing large-scale chemical weather pictures. NASA has also supported development of global models which simulate atmospheric transport and transformation processes of important atmospheric gas and aerosol species. I will present the current modeling and satellite capabilities for PM2.5 studies, the possibilities and challenges in using satellite data for PM2.5 forecasts, and the needs of future remote sensing data for improving air quality monitoring and modeling.

Chin, Mian↗
NASA NTRS2015

Investigation of Micro-Scale Architectural Effects on Damage of Composites

This paper presents a three-dimensional, energy based, anisotropic, stiffness reduction, progressive damage model for composite materials and composite material constituents. The model has been implemented as a user-defined constitutive model within the Abaqus finite element software package and applied to simulate the nonlinear behavior of a damaging epoxy matrix within a unidirectional composite material. Three different composite microstructures were considered as finite element repeating unit cells, with appropriate periodicity conditions applied at the boundaries. Results representing predicted transverse tensile, longitudinal shear, and transverse shear stress-strain curves are presented, along with plots of the local fields indicating the damage progression within the microstructure. It is demonstrated that the damage model functions appropriately at the matrix scale, enabling localization of the damage to simulate failure of the composite material. The influence of the repeating unit cell geometry and the effect of the directionality of the applied loading are investigated and discussed.

damage↗
DOE OSTI2024

Parametric Instability of Alfvén Waves and Wave Packets in Periodic and Open Systems

The parametric decay instability of Alfvén waves has been widely studied, but few investigations have examined wave packets of finite size and the effect of different boundary conditions on the growth rate. In this paper, we perform a linear analysis of circular and arc-polarized wave trains and wave packets in periodic and open boundary systems in a low- β plasma. We find that both types of wave are 3–5 times more stable in open boundary conditions compared to periodic. Additionally, once the wave packet width ℓ becomes smaller than the system size L , the growth rate decreases nearly with a power law γ ∝ ℓ / L . This study demonstrates that the stability of a pump wave cannot be separated from the laboratory settings, and that the growth rate of daughter waves depends on the conditions downstream and upstream of the pump wave and on the fraction of volume it fills. Our results can explain simulations and experiments of localized Alfvén waves. They also suggest that Alfvénic fluctuations in the solar wind, including sharp impulses known as switchbacks, can be more stable than traditional theory suggests depending on wind conditions.

Alfven waves↗
NASA NTRS1990

Spurious frequencies as a result of numerical boundary treatments

The stability theory for finite difference Initial Boundary-Value approximations to systems of hyperbolic partial differential equations states that the exclusion of eigenvalues and generalized eigenvalues is a sufficient condition for stability. The theory, however, does not discuss the nature of numerical approximations in the presence of such eigenvalues. In fact, as was shown previously, for the problem of vortex shedding by a 2-D cylinder in subsonic flow, stating boundary conditions in terms of the primitive (non-characteristic) variables may lead to such eigenvalues, causing perturbations that decay slowly in space and remain periodic time. Characteristic formulation of the boundary conditions avoided this problem. A more systematic study of the behavior of the (linearized) one-dimensional gas dynamic equations under various sets of oscillation-inducing legal boundary conditions is reported.

Abarbanel, Saul↗
DOE OSTI2022

Discrete dislocation dynamics for crystal RVEs. Part 1: Periodic network kinematics

A novel implementation of the dislocation flux boundary condition in discrete dislocation dynamics is presented. The continuity of the individual dislocation loops in a periodic representative crystal volume (RVE) is enforced across the boundary of the RVE with the help of a dual topological description for representing dislocation line kinematics in two equivalent spaces representing the deforming crystal, the RVE and the unbounded crystal spaces. The former describes the motion of the dislocations in the simulated crystal RVE whereas the latter represents the motion of dislocations in an infinite space containing all replicas of the RVE. A mapping between the two spaces forms the basis of the implementation of flux boundary condition. The implementation details are discussed in the context of statistical homogeneity of bulk crystals undergoing macroscopically homogeneous plastic deformation. In this case, the boundary nodes associated with dislocation segments bear no relevance in the motion of the dislocations. Finally, some test cases are presented and discussed to establish the proposed approach.

42 ENGINEERING↗