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 289 records · Page 16

Vector potential and surface magnetic field in magnetoelectric antiferromagnetic materials

A general formula for the average vector potential of bulk periodic systems is proposed and shown to set the boundary conditions at magnetic interfaces. Here, for antiferromagnetic materials, the study reveals a unique relation between the macroscopic potential and the orientation-dependent magnetic quadrupole as a result of the different crystalline and magnetic symmetries. In particular, at surfaces and interfaces of a truncated bulk without inversion and time-reversal symmetries, the average vector potential exhibits a discontinuity which results in an interfacial magnetic field. In general, however, due to the surface and interface electronic and atomic relaxations, additional magnetization may result. For the experimentally observed magnetoelectric antiferromagnets, in particular, our symmetry analysis suggests that the relaxation effects could well be a system response to the presence of such a potential discontinuity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Object-oriented analysis as a foundation for building climate storylines of compounding short-term drought and crop heat stress

Introduction: Crops are vulnerable to precipitation and heat extremes during late spring through summer. Methods: We analyzed for a north-central U.S. region short-term drought and agricultural heat stress during April-May-June-July. We used the 4-km Parameter Elevation Regression on Independent Slopes Model (PRISM) for observations, aggregated to a 25-km grid, and two 25-km Regional Climate Model version 4 (RegCM4) simulns used either GFDL- or MPI-GCM boundary conditions. We chose 1981-2000 as our contemporary time period, and 2041- 2060 as our scenario time period, which used the Representative Concentration Pathway 8.5 emissions scenario. We used object-oriented analysis to identify events of interest in observations and simulations by identifying objects in a space-time domain that meet specified criteria, such as exceeding a heat-stress temperature threshold. The event diagnosis allowed analysis of compound events, occurring when temperature and drought objects overlap. Results: Identified objects yielded events that can undermine agricultural productivity and which are thus relevant to decision makers, making them building blocks for possible climate storylines. The observations and simulations showed similar spatial distributions of event frequencies across the analysis region. However, the simulations attained this distribution by having fewer events that tend to cover larger areas compared to observed events, suggesting that the effective resolution of the simulations was coarser than their 25-km grids. Short-term drought frequency increased and heat-stress frequency decreased in transitioning to the scenario climate. When compounding occurred heat-stress events generally preceded the short-term drought events. The overlapping, compound events tended to be more extreme compared to non-overlapping events of either type. Discussion: The information yielded projected changes in these agriculturally motivated events. One prominent conditional behavior emerging from the work was that a heat-stress event should be a warning to watch for potential drought, as both could compound each other to more intense levels.

54 ENVIRONMENTAL SCIENCES↗

Break-up characteristics of Chena River Basin

The author has identified the following significant results. Utilization of the color additive viewer was made to delineate snow boundary conditions from imagery taken during the fall snow accumulation period. By using filters that falsely colored the vegetation red, the snow cover appeared white.

Carlson, R. F.↗

Transition to turbulence in plane channel flow

A numerical simulation of the final stages of transition to turbulence in plane channel flow is reported. Three dimensional, incompressible Navier-Stokes equations are numerically integrated to obtain the time-evolution of two and three dimensional finite amplitude disturbances. Computations are performed on the CYBER-203 vector processor for a 32x51x32 grid. Results are presented for no-slip boundary conditions at the solid walls as well as for periodic suction-blowing to simulate active control of transition by mass transfer. Solutions indicate that the method is capable of simulating the complex character of vorticity dynamics during the various stages of transition and final breakdown. In particular, evidence points to the formation of a lambda-shape vortex and the subsequent system of horseshoe vortices inclined to the main flow direction as the main elements of transition. Calculations involving suction-blowing indicate that interference with a wave of suitable phase and amplitude reduces the disturbance growth rates.

Biringen, S.↗

Transition to turbulence in plane channel flows

Results obtained from a numerical simulation of the final stages of transition to turbulence in plane channel flow are described. Three dimensional, incompressible Navier-Stokes equations are numerically integrated to obtain the time evolution of two and three dimensional finite amplitude disturbances. Computations are performed on CYBER-203 vector processor for a 32x51x32 grid. Results are presented for no-slip boundary conditions at the solid walls as well as for periodic suction blowing to simulate active control of transition by mass transfer. Solutions indicate that the method is capable of simulating the complex character of vorticity dynamics during the various stages of transition and final breakdown. In particular, evidence points to the formation of a lambda-shape vortex and the subsequent system of horseshoe vortices inclined to the main flow direction as the main elements of transition. Calculations involving periodic suction-blowing indicate that interference with a wave of suitable phase and amplitude reduces the disturbance growth rates.

Biringen, S.↗

An explicit example of Hopf bifurcation in fluid mechanics

It is observed that a complete and explicit example of Hopf bifurcation appears not to be known in fluid mechanics. Such an example is presented for the rotating Benard problem with free boundary conditions on the upper and lower faces, and horizontally periodic solutions. Normal modes are found for the linearization, and the Veronis computation of the wave numbers is modified to take into account the imposed horizontal periodicity. An invariant subspace of the phase space is found in which the hypotheses of the Joseph-Sattinger theorem are verified, thus demonstrating the Hopf bifurcation. The criticality calculations are carried through to demonstrate rigorously, that the bifurcation is subcritical for certain cases, and to demonstrate numerically that it is subcritical for all the cases in the paper.

Kloeden, P.↗

Pliocene Model Intercomparison Project (PlioMIP): Experimental Design and Boundary Conditions (Experiment 2)

The Palaeoclimate Modelling Intercomparison Project has expanded to include a model intercomparison for the mid-Pliocene warm period (3.29 to 2.97 million yr ago). This project is referred to as PlioMIP (the Pliocene Model Intercomparison Project). Two experiments have been agreed upon and together compose the initial phase of PlioMIP. The first (Experiment 1) is being performed with atmosphere only climate models. The second (Experiment 2) utilizes fully coupled ocean-atmosphere climate models. Following on from the publication of the experimental design and boundary conditions for Experiment 1 in Geoscientific Model Development, this paper provides the necessary description of differences and/or additions to the experimental design for Experiment 2.

Haywood, A. M.↗

Spurious fields in time domain computations of scattering problems

In this paper two-dimensional electromagnetic scattering problems with a time-periodic incident field are considered. The scatterer is a perfect conductor, and an artificial boundary condition is used. The large time behavior of solutions, depending on (divergence-free) initial conditions, is characterized. It turns out that in addition to the expected time-periodic solution the limiting solution may also contain a spurious stationary field. The source of the stationary field is explained and equations describing it are obtained. Several avoidance strategies are discussed, and numerical comparisons of these techniques are given.

Kangro, Urve↗

A four-dimensional primitive equation model for coupled coastal-deep ocean studies

A prototype four dimensional continental shelf/deep ocean model is described. In its present form, the model incorporates the effects of finite amplitude topography, advective nonlinearities, and variable stratification and rotation. The model can be forced either directly by imposed atmospheric windstress and surface pressure distributions, and energetic mean currents imposed by the exterior oceanic circulation; or indirectly by initial distributions of shoreward propagation mesoscale waves and eddies. To avoid concerns over the appropriate specification of 'open' boundary conditions on the cross-shelf and seaward model boundaries, a periodic channel geometry (oriented along-coast) is used. The model employs a traditional finite difference expansion in the cross-shelf direction, and a Fourier (periodic) representation in the long-shelf coordinate.

Haidvogel, D. B.↗

Computation of Transonic Nozzle Sound Transmission and Rotor Problems by the Dispersion-Relation-Preserving Scheme

The transonic nozzle transmission problem and the open rotor noise radiation problem are solved computationally. Both are multiple length scales problems. For efficient and accurate numerical simulation, the multiple-size-mesh multiple-time-step Dispersion-Relation-Preserving scheme is used to calculate the time periodic solution. To ensure an accurate solution, high quality numerical boundary conditions are also needed. For the nozzle problem, a set of nonhomogeneous, outflow boundary conditions are required. The nonhomogeneous boundary conditions not only generate the incoming sound waves but also, at the same time, allow the reflected acoustic waves and entropy waves, if present, to exit the computation domain without reflection. For the open rotor problem, there is an apparent singularity at the axis of rotation. An analytic extension approach is developed to provide a high quality axis boundary treatment.

Tam, Christopher K. W.↗

The use of NASA GEOS Global Analysis in MM5/WRF Initialization: Current Studies and Future Applications

An effort has been made at NASA/GSFC to use the Goddard Earth Observing system (GEOS) global analysis in generating the initial and boundary conditions for MM5/WRF simulation. This linkage between GEOS global analysis and MM5/WRF models has made possible for a few useful applications. As one of the sample studies, a series of MM5 simulations were conducted to test the sensitivity of initial and boundary conditions to MM5 simulated precipitation over the eastern; USA. Global analyses horn different operational centers (e.g., NCEP, ECMWF, I U ASA/GSFCj were used to provide first guess field and boundary conditions for MM5. Numerical simulations were performed for one- week period over the eastern coast areas of USA. the distribution and quantities of MM5 simulated precipitation were compared. Results will be presented in the workshop. In addition,other applications from recent and future studies will also be addressed.

Pu, Zhao-Xia↗

Helfrich-Hurault elastic instabilities driven by geometrical frustration

Here, the Helfrich-Hurault (HH) elastic instability is a well-known mechanism behind patterns that form as a result of strain upon liquid crystal systems with periodic ground states. In the HH model, layered structures undulate and buckle in response to local, geometric incompatibilities in order to maintain the preferred layer spacing. Classic HH systems include cholesteric liquid crystals under electromagnetic field distortions and smectic liquid crystals under mechanical strains, where both materials are confined between rigid substrates. However, richer phenomena are observed when undulation instabilities occur in the presence of deformable interfaces and variable boundary conditions. Understanding how the HH instability is affected by deformable surfaces is imperative for applying the instability to a broader range of materials. In this review, the HH mechanism is reexamined and special focus is given to how the boundary conditions influence the response of lamellar systems to geometrical frustration. Lamellar liquid crystals confined within a spherical shell geometry are used as the model system. Made possible by the relatively recent advances in microfluidics within the past 15 years, liquid crystal shells are composed entirely of fluid interfaces and have boundary conditions that can be dynamically controlled at will. Past and recent work that exemplifies how topological constraints, molecular anchoring conditions, and boundary curvature can trigger the HH mechanism in liquid crystals with periodic ground states is examined. The review ends by identifying similar phenomena across a wide variety of materials, both biological and synthetic. The fact that the HH mechanism is a generic and often overlooked response of periodic materials to geometrical frustration is highlighted.

42 ENGINEERING↗

(abstract) Role of Long Equatorial Wave Reflection in the Low-Frequency Variability Observed During the TOPEX/POSEIDON Period: Data Analysis and Coupled Model Study

The analysis of TOPEX/POSEIDON sea level and ERS-1 zonal wind stress data gives evidence that first mode meridonal long Rossby waves fully reflected into Kelvin waves at the equatorial western Pacific boundary at different periods including the January-June 1994 period. The evolution of the conditions (zonal wind stress, sea surface temperature and sea level anomalies) in the central Pacific suggests that this reflection may have played a role in the reversal of cold to warm anomalies observed in the central Pacific in June-July 1994. To investigate the actual role of reflected Kelvin waves observed by TOPEX/POSEIDON, a simple ocean-atmosphere coupled model is designed. Starting from these initial conditions, several simulations are run for testing the role western bundary reflection evidenced in TOPEX/POSEIDON data. Coupled mechanisms are suggested to explain the low-frequency variability observed during this period.

Rossby waves Kelvin waves zonal wind stress TOPEX/↗

Large Eddy Simulation in a Channel with Exit Boundary Conditions

The influence of the exit boundary conditions (vanishing first derivative of the velocity components and constant pressure) on the large eddy simulation of the fully developed turbulent channel flow has been investigated for equidistant and stretched grids at the channel exit. Results show that the chosen exit boundary conditions introduce some small disturbance which is mostly damped by the grid stretching. The difference between the fully developed turbulent channel flow obtained with LES with periodicity condition and the inlet and exit and the LES with fully developed flow at the inlet and the exit boundary condition is less than 10% for equidistant grids and less than 5% for the case grid stretching. The chosen boundary condition is of interest because it may be used in complex flows with backflow at exit.

Cziesla, T.↗

Numerical simulations of plasma double layers

The results of analytical studies of quasi-static electric fields along geomagnetic field lines are discussed. The calculations were targeted at the structure, generation mechanisms and stability parameters. The field consists of two oppositely charged layers, either weakly or strongly charged, with an electric field between. Existence conditions are defined for the double layer field and balancing requirements are explored. Details of the simulation techniques, i.e., particle in cell and Vlasov simulations, for studying the double layer are outlined, noting that both periodic and quasi-periodic simulations are used. Solutions to Poisson's equation for fixed and floating point boundary conditions are generated. Finally, attention is also given to oblique and two-dimensional magnetic double layers.

Goertz, C. K.↗

Time-dependent transonic flow solutions for axial turbomachinery

Three-dimensional unsteady transonic flow through an axial turbomachine stage is described in terms of a pair of two-dimensional formulations pertaining to orthogonal surfaces, namely, a blade-to-blade surface and a hub-to-casing surface. The resulting systems of nonlinear, inviscid, compressible equations of motion are solved by an explicit finite-difference technique. The blade-to-blade program includes the periodic interaction between rotor and stator blade rows. Treatment of the boundary conditions and of the blade slipstream motion by a characteristic type procedure is discussed in detail. Harmonic analysis of the acoustic far field produced by the blade row interaction, including an arbitrary initial transient, is outlined. Results from the blade-to-blade program are compared with experimental measurements of the rotating pressure field at the tip of a high-speed fan. The hub-to-casing program determines circumferentially averaged flow properties on a meridional plane. Blade row interactions are neglected in this formulation, but the force distributions over the entire blade surface for both the rotor and stator are obtained. Results from the hub-to-casing program are compared with a relaxation method solution for a subsonic rotor. Results are also presented for a quiet fan stage which includes transonic flow in both the rotor and stator and a normal shock in the stator.

Erdos, J.↗

PROTEUS two-dimensional Navier-Stokes computer code, version 1.0. Volume 2: User's guide

A new computer code was developed to solve the two-dimensional or axisymmetric, Reynolds averaged, unsteady compressible Navier-Stokes equations in strong conservation law form. The thin-layer or Euler equations may also be solved. Turbulence is modeled using an algebraic eddy viscosity model. The objective was to develop a code for aerospace applications that is easy to use and easy to modify. Code readability, modularity, and documentation were emphasized. The equations are written in nonorthogonal body-fitted coordinates, and solved by marching in time using a fully-coupled alternating direction-implicit procedure with generalized first- or second-order time differencing. All terms are linearized using second-order Taylor series. The boundary conditions are treated implicitly, and may be steady, unsteady, or spatially periodic. Simple Cartesian or polar grids may be generated internally by the program. More complex geometries require an externally generated computational coordinate system. The documentation is divided into three volumes. Volume 2 is the User's Guide, and describes the program's general features, the input and output, the procedure for setting up initial conditions, the computer resource requirements, the diagnostic messages that may be generated, the job control language used to run the program, and several test cases.

Towne, Charles E.↗

PROTEUS two-dimensional Navier-Stokes computer code, version 1.0. Volume 3: Programmer's reference

A new computer code was developed to solve the 2-D or axisymmetric, Reynolds-averaged, unsteady compressible Navier-Stokes equations in strong conservation law form. The thin-layer or Euler equations may also be solved. Turbulence is modeled using an algebraic eddy viscosity model. The objective was to develop a code for aerospace applications that is easy to use and easy to modify. Code readability, modularity, and documentation were emphasized. The equations are written in nonorthogonal body-fitted coordinates, and solved by marching in time using a fully-coupled alternating-direction-implicit procedure with generalized first- or second-order time differencing. All terms are linearized using second-order Taylor series. The boundary conditions are treated implicitly, and may be steady, unsteady, or spatially periodic. Simple Cartesian or polar grids may be generated internally by the program. More complex geometries require an externally generated computational coordinate system. The documentation is divided into three volumes. Volume 3 is the Programmer's Reference, and describes the program structure, the FORTRAN variables stored in common blocks, and the details of each subprogram.

Towne, Charles E.↗