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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
NASA NTRS2010

Three-Dimensional Unsteady Simulation of a Modern High Pressure Turbine Stage Using Phase Lag Periodicity: Analysis of Flow and Heat Transfer

Unsteady three-dimensional RANS simulations have been performed on a highly loaded transonic turbine stage and results are compared to steady calculations as well as experiment. A low Reynolds number k- turbulence model is employed to provide closure for the RANS system. A phase-lag boundary condition is used in the periodic direction. This allows the unsteady simulation to be performed by using only one blade from each of the two rows. The objective of this paper is to study the effect of unsteadiness on rotor heat transfer and to glean any insight into unsteady flow physics. The role of the stator wake passing on the pressure distribution at the leading edge is also studied. The simulated heat transfer and pressure results agreed favorably with experiment. The time-averaged heat transfer predicted by the unsteady simulation is higher than the heat transfer predicted by the steady simulation everywhere except at the leading edge. The shock structure formed due to stator-rotor interaction was analyzed. Heat transfer and pressure at the hub and casing were also studied. Thermal segregation was observed that leads to the heat transfer patterns predicted by steady and unsteady simulations to be different.

Shyam, Vikram↗
DOE OSTI2024

Simulating Global Terrestrial Carbon and Nitrogen Biogeochemical Cycles With Implicit and Explicit Representations of Soil Microbial Activity

Abstract Nutrient limitation is widespread in terrestrial ecosystems. Accordingly, representations of nitrogen (N) limitation in land models typically dampen rates of terrestrial carbon (C) accrual, compared with C‐only simulations. These previous findings, however, rely on soil biogeochemical models that implicitly represent microbial activity and physiology. Here we present results from a biogeochemical model testbed that allows us to investigate how an explicit versus implicit representation of soil microbial activity, as represented in the MIcrobial‐MIneral Carbon Stabilization (MIMICS) and Carnegie‐Ames‐Stanford Approach (CASA) soil biogeochemical models, respectively, influence plant productivity, and terrestrial C and N fluxes at initialization and over the historical period. When forced with common boundary conditions, larger soil C pools simulated by the MIMICS model reflect longer inferred soil organic matter (SOM) turnover times than those simulated by CASA. At steady state, terrestrial ecosystems experience greater N limitation when using the MIMICS‐CN model, which also increases the inferred SOM turnover time. Over the historical period, however, warming‐induced acceleration of SOM decomposition over high latitude ecosystems increases rates of N mineralization in MIMICS‐CN. This reduces N limitation and results in faster rates of vegetation C accrual. Moreover, as SOM stoichiometry is an emergent property of MIMICS‐CN, we highlight opportunities to deepen understanding of sources of persistent SOM and explore its potential sensitivity to environmental change. Our findings underscore the need to improve understanding and representation of plant and microbial resource allocation and competition in land models that represent coupled biogeochemical cycles under global change scenarios.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2002

Initialization of Formation Flying Using Primer Vector Theory

In this paper, we extend primer vector analysis to formation flying. Optimization of the classical rendezvous or free-time transfer problem between two orbits using primer vector theory has been extensively studied for one spacecraft. However, an increasing number of missions are now considering flying a set of spacecraft in close formation. Missions such as the Magnetospheric MultiScale (MMS) and Leonardo-BRDF (Bidirectional Reflectance Distribution Function) need to determine strategies to transfer each spacecraft from the common launch orbit to their respective operational orbit. In addition, all the spacecraft must synchronize their states so that they achieve the same desired formation geometry over each orbit. This periodicity requirement imposes constraints on the boundary conditions that can be used for the primer vector algorithm. In this work we explore the impact of the periodicity requirement in optimizing each spacecraft transfer trajectory using primer vector theory. We first present our adaptation of primer vector theory to formation flying. Using this method, we then compute the AV budget for each spacecraft subject to different formation endpoint constraints.

Mailhe, Laurie↗
NASA NTRS1995

A mesoscale modeling study of the atmospheric circulation of high southern latitudes

The meteorology of high southern latitudes during winter is simulated using a cloud-free version of The Pennsylvania State University-National Center for Atmospheric Research Mesoscale Model version 4 (MM4) with a 100-km horizontal resolution. Comparisons between idealized simulations of Antarctic with MM4 and with the mesoscale model of Parish and Waight reveal that both models produce similarly realistic velocity fields in the boundary layer. The latter model tends to produce slightly faster drainage winds over East Antarctica. The intensity of the katabatic winds produced by MM4 is sensitive to parameterizations of boundary layer fluxes. Two simulations performed with MM4 using analyses from the European Center for Medium-Range Weather Forecasts (ECMWF) for June 1988 as initial and boundary conditions. A simulation of the period from 000 UTC 2 June to 0000 UTC 8 June produces realistic synoptic phenomena including ridge development over East Antarctica, frontogenesis over the Amundsen Sea, and a katabatic surge over the Ross Ice Shelf. The simulated time-averaged fields for June 1988, particularly that of a 500-hPa height, are in good agreement with time-averaged fields analyzed by the ECMWF. The results of the simulations provide detailed features of the Antarctic winter boundary layer along the steeply sloping terrain. Highest boundary layer wind speeds averaged over the month-long simulation are approximately 20 m/s. The lack of latent heating in the simulations apparently results in some bias in the results. In particular, the cloud-free version of MM4 underpredicts the intensity of lows in the sea level pressure field.

Hines, K. M.↗
DOE OSTI2025

Elasto-viscoplastic fast Fourier transform modeling framework for assessing microstructural effects on stress intensity factors characterizing fracture toughness

A large-strain elasto-viscoplastic fast Fourier transform (LS-EVPFFT) model with non-periodic (NP) velocity-based boundary conditions is adapted to simulate the sensitivity of stress intensity factors on microstructure for 304L stainless steel. The material was characterized via electron backscattered diffraction (EBSD) serial-sectioning to obtain a measured 3-D microstructural cell to perform simulations. The NP-LS-EVPFFT model, including the simulation setup and boundary conditions, was verified using a crystal plasticity finite element (CPFE) model. To this end, the generation of meshes of notched specimens was developed, which involved creating Python scripts for mesh “cutting” in Abaqus, and Sculpt scripts in Cubit for meshing of the measured microstructural cell processed with DREAM.3D. The complexity of the mesh preparation highlighted the advantages of the FFT-based model, which circumvents the mesh generation process. Given the efficiency of the FFT-based model, statistical distribution of stress intensity factors in function of crystal orientation at the crack tip, grain structure, and crystallographic texture surrounding the crack tip were predicted. Further, the distributions reveal about 10% variation of stress intensity factors with microstructure with the most significant sensitivity found to be the crystal orientation at the crack tip. The methodology developed in this work is discussed as a practical simulation tool for predicting the sensitivity of stress intensity factors on microstructural variability in metallic materials.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Assessing boundary condition and parametric uncertainty in numerical-weather-prediction-modeled, long-term offshore wind speed through machine learning and analog ensemble

To accurately plan and manage wind power plants, not only does the time-varying wind resource at the site of interest need to be assessed but also the uncertainty connected to this estimate. Numerical weather prediction (NWP) models at the mesoscale represent a valuable way to characterize the wind resource offshore, given the challenges connected with measuring hub-height wind speed. The boundary condition and parametric uncertainty associated with modeled wind speed is often estimated by running a model ensemble. However, creating an NWP ensemble of long-term wind resource data over a large region represents a computational challenge. Here, we propose two approaches to temporally extrapolate wind speed boundary condition and parametric uncertainty using a more convenient setup in which a mesoscale ensemble is run over a short-term period (1 year), and only a single model covers the desired long-term period (20 year). We quantify hub-height wind speed boundary condition and parametric uncertainty from the short-term model ensemble as its normalized across-ensemble standard deviation. Then, we develop and apply a gradient-boosting model and an analog ensemble approach to temporally extrapolate such uncertainty to the full 20-year period, for which only a single model run is available. As a test case, we consider offshore wind resource characterization in the California Outer Continental Shelf. Both of the proposed approaches provide accurate estimates of the long-term wind speed boundary condition and parametric uncertainty across the region (R 2 >0.75), with the gradient-boosting model slightly outperforming the analog ensemble in terms of bias and centered root-mean-square error. At the three offshore wind energy lease areas in the region, we find a long-term median hourly uncertainty between 10 % and 14 % of the mean hub-height wind speed values. Finally, we assess the physical variability in the uncertainty estimates. In general, we find that the wind speed uncertainty increases closer to land. Also, neutral conditions have smaller uncertainty than the stable and unstable cases, and the modeled wind speed in winter has less boundary condition and parametric sensitivity than summer.

17 WIND ENERGY↗
DOE OSTI2025

Characterizing Near-Surface Moisture Increase during the Clear-Sky Afternoon-to-Evening Transition Using a Single-Column Model

The afternoon-to-evening transition is the period when the atmospheric boundary layer transitions from convective to stable conditions. One noticeable feature during this transition period is the rapid increase in water vapor concentration near the surface. However, the mechanism behind the increase in water vapor remains poorly understood. This study investigated the processes contributing to the water vapor increase and the impacts of the land cover and horizontal advection on the water vapor increase using a single-column model for the clear-sky condition. Numerical experiments were conducted on three cases at the Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) site. Evapotranspiration was found to be the main moisture source to the water vapor increase, and that change of the near-surface vertical temperature gradient from negative to positive is the trigger for this increase during the afternoon-to-evening transition. This is because the near-surface turbulence divergence term decreases due to the change in the buoyancy profile caused by vertical temperature gradient change. The impact of horizontal advection on water vapor varies, and it can either lead to an increase or a decrease in water vapor, depending on the spatial horizontal water vapor gradient. This study also found that land cover can influence the timing of the water vapor increase since different land covers may have different Bowen ratios. Water vapor increases earlier under conditions with smaller Bowen ratios compared to those with larger values.

Atmosphere-land interaction↗
DOE OSTI2021

Injection data analysis using material balance time for CO 2 storage capacity estimation in deep closed saline aquifers

Estimating the ultimate storage capacity of deep saline aquifers is important to address the formation potential to store the envisioned large volumes of CO 2 . Injection data (i.e. injection rate, bottomhole pressure, and cumulative injected volume of CO 2 ) are routinely recorded during storage operations. These data contain valuable information on the subsurface (e.g. the reservoir pore volume and the formation storage capacity) that can be extracted. In this paper, we present a two-step graphical technique to infer the pore volume and the ultimate storage capacity of closed saline aquifers by analyzing the available injection data. First, the pore volume is inferred through adapting the concept of the material balance time. Material balance time is an approximate superposition time function developed to interpret production data from oil and gas wells operating at variable pressure/rate conditions during the boundary-dominated flow period. Using material balance techniques, the ultimate storage capacity is then estimated through linear extrapolation of the average pressure trend to the maximum allowable pressure the formation can withstand. The average pressure is not available in practice, but is can be obtained from the injection data. Two approaches are presented in this study to calculate the average pressure; namely the rigorous and the approximate approaches. Unlike the rigorous approach, the approximate approach does not require a prior knowledge of some reservoir properties (e.g. relative permeability, absolute permeability, formation porosity and thickness) to calculate the average pressure. To investigate its potential and reliability in analyzing CO 2 injection data, the proposed technique is applied to four synthetic cases representing different well operating conditions. Results indicate that the approximate approach consistently overestimates the actual (simulated) storage capacity as compared to the rigorous approach. The agreement - between the inferred and the simulated reservoir pore volume, and between the analytical and numerical estimates of storage capacity - validates the potential application of the technique to CO 2 storage in closed saline aquifers. The technique is further substantiated through application to a field data set utilized from a commercial-scale geological storage (CGS) project. Finally, field data interpretation shows that the proposed technique can be utilized to identify the degree of hydraulic continuity and reservoir compartmentalization within a target formation by interpreting the corresponding pressure and rate responses.

02 PETROLEUM↗
DOE OSTI2020

Dislocation-graphene interactions in Cu/graphene composites and the effect of boundary conditions: a molecular dynamics study

Here, in this paper, systematic molecular dynamics simulations were employed to study dislocation-graphene interactions in Cu/graphene systems using two carefully selected configurations: (i) edge dislocations in the pileup interacting with graphene during nanoindentation and (ii) complex dislocations interacting with graphene in compression of layered nanopillars. The intrinsic and extrinsic size effects were investigated with respect to varying the Cu lamella thickness and pillar sizes, revealing an anomalous extrinsic size effect: the smaller, the weaker. To understand which boundary conditions are physically consistent, both free and periodic graphene were considered. It was seen that edge dislocations in the pileup continuously transmitted across the free graphene through Cu/graphene interfacial sliding and graphene reorientation, whereas transmission was very difficult for periodic graphene as its in-plane deformation was required. In compression of Cu/graphene nanopillars, the strengthening effect of free graphene was found to be less obvious than that of the periodic case, which could be attributed to free graphene edges acting as dislocation sources after the compressed Cu overflowed the graphene sheet from all directions. We therefore conclude that the strengthening effect of periodic graphene inclusions appears to be overestimated, such that free graphene should be more appropriate.

36 MATERIALS SCIENCE↗
NASA NTRS1986

Coupled aerodynamic and acoustical predictions for turboprops

To predict the noise fields for proposed turboprop airplanes, an existing turboprop noise code by Farassat has been modified to accept blade pressure inputs from a three-dimensional aerodynamic code. A Euler-type code can handle the nonlinear transonic flow of these high-speed, highly swept blades. This turbofan code was modified to allow the calculation mesh to extend to about twice the blade radius and to apply circumferential periodicity rather than solid-wall boundary conditions on the blade in the region between the blade tip and the outer shroud. Outputs were added for input to the noise prediction program and for color contour plots of various flow variables. The Farassat input subroutines were modified to read files of blade coordinates and predicted surface pressures. Aerodynamic and acoustic results are shown for the SR-3 model blade. Comparison of the acoustic predicted results with measured data show good agreement.

Clark, Bruce J.↗
NASA NTRS1992

Influences of CO2 sublimation/condensation processes on the long-term evolution of the Martian atmosphere

The accumulation process of frozen CO2 within the Martian regolith is analyzed. The boundary conditions are taken to be the periodically variable temperature at the Mars surface and the constant heat flux at the bottom of the regolith layers, some hundred meters thick. The mean value of temperature as well as its amplitude are assumed to be dependent on latitude and could be variable with the geological time scale. The flux of gaseous CO2 within the regolith is controlled by the mechanism of gas transport through the porous medium. An appropriate initial geometry of distribution of pores is assumed. The porosity and flux change due to condensation/sublimation processes. The equations of heat and mass transfer are solved numerically for a given latitude. Next, the local (over the parallel of latitude belt) annual balance of CO2 is calculated. The procedure is repeated over the whole Martian surface and next the total annual net flux is calculated.

Kossacki, K.↗
DOE OSTI2024

Phase Pattern Formation in Grain Boundaries

Here, in this Letter, we derive conditions that predict the existence of two-phase periodic-pattern grain boundary structures that are stable against coarsening. While previous research has established that elastic effects can lead to phase pattern formation on crystal surfaces, the possibility of stable grain boundary structures composed of alternating grain boundary phases has not been previously analyzed. Our theory identifies the specific combination of grain boundary and materials properties that enable the emergence of patterned grain boundary states and shows that the dislocation content of grain boundary phase junctions, absent in surface phenomena, weakens the stability of the patterned structures. The predictions of the theory are tested using a model copper grain boundary that exhibits multiple phases and two-phase pattern formation. We discuss how, similarly to surfaces, elastic effects associated with grain boundary phase junctions have profound implications for how grain boundary phases transform.

36 MATERIALS SCIENCE↗
NASA NTRS2001

A Regional Model Study of Synoptic Features Over West Africa

Synoptic weather features over West Africa were studied in simulations by the regional simulation model (RM) at the NASA/Goddard Institute for Space Studies. These pioneering simulations represent the beginning of an effort to adapt regional models for weather and climate prediction over West Africa. The RM uses a cartesian grid with 50 km horizontal resolution and fifteen vertical levels. An ensemble of four simulations was forced with lateral boundary conditions from ECMWF global analyses for the period 8-22 August 1988. The simulated mid-tropospheric circulation includes the skillful development and movement of several African wave disturbances. Wavelet analysis of mid-tropospheric winds detected a dominant periodicity of about 4 days and a secondary periodicity of 5-8 days. Spatial distributions of RM precipitation and precipitation time series were validated against daily rain gauge measurements and ISCCP satellite infrared cloud imagery. The time-space distribution of simulated precipitation was made more realistic by combining the ECMWR initial conditions with a 24-hr spin-up of the moisture field and also by damping high frequency gravity waves by dynamic initialization. Model precipitation "forecasts" over the Central Sahel were correlated with observations for about three days, but reinitializing with observed data on day 5 resulted in a dramatic improvement in the precipitation validation over the remaining 9 days. Results imply that information via the lateral boundary conditions is not always sufficient to minimize departures between simulated and actual precipitation patterns for more than several days. In addition, there was some evidence that the new initialization may increase the simulations' sensitivity to the quality of lateral boundary conditions.

Druyan, Leonard M.↗
DOE OSTI2025

The effect of African Easterly Wave suppression by periodicity on Atlantic tropical cyclones

Research has shown that suppressing African Easterly Waves (AEWs) does not reduce basin-wide North Atlantic tropical cyclone (TC) frequency but can enhance TC environmental favorability. We investigated the AEW-TC relationship further by examining the effects of suppressing the two AEW periodicities individually on TC activity. Using regional model simulations, AEWs were prescribed or suppressed in the 2–6 d or 6–10 d ranges through the lateral boundary conditions. Seasonal TC frequency increased significantly when either AEW periodicity was suppressed, with a larger increase when the 2–6 d waves were suppressed. We also found that suppressing the 2–6 d waves increased mid-tropospheric moisture by up to 8%, as well as overall atmospheric instability, near the western coast of northern Africa. Furthermore, the convective disturbances that developed into TCs exhibited stronger rotation, increased ascending motion, and higher rainfall. Our results suggest that reduced 2–6 d AEW activity may trigger a more active TC season.

Danso, Derrick K↗
NASA NTRS2008

Improving Air Quality Forecasts with AURA Observations

Past studies have identified model initial and boundary conditions as sources of reducible errors in air-quality simulations. In particular, improving the initial condition improves the accuracy of short-term forecasts as it allows for the impact of local emissions to be realized by the model and improving boundary conditions improves long range transport through the model domain, especially in recirculating anticyclones. During the August 2006 period, we use AURA/OMI ozone measurements along with MODIS and CALIPSO aerosol observations to improve the initial and boundary conditions of ozone and Particulate Matter. Assessment of the model by comparison of the control run and satellite assimilation run to the IONS06 network of ozonesonde observations, which comprise the densest ozone sounding campaign ever conducted in North America, to AURA/TES ozone profile measurements, and to the EPA ground network of ozone and PM measurements will show significant improvement in the CMAQ calculations that use AURA initial and boundary conditions. Further analyses of lightning occurrences from ground and satellite observations and AURA/OMI NO2 column abundances will identify the lightning NOx signal evident in OMI measurements and suggest pathways for incorporating the lightning and NO2 data into the CMAQ simulations.

Newchurch, M. J.↗
NASA NTRS1988

Scattering of electromagnetic waves from a periodic surface with random roughness

Equations for the scattering of electromagnetic waves from a randomly perturbed periodic surface have been formulated using the extended boundary condition method and solved using the small perturbation method. Surface currents and scattered fields are solved for up to the second order. The results indicate that as the correlation length of the random roughness increases, the bistatic scattering patterns of the scattered fields show several beams associated with each Bragg diffraction direction of the periodic surface. The beam shape becomes broader with smaller correlation length. Results obtained using the Kirchhoff approximation are found to agree well with the present results for the hh and vv polarized backscattering coefficients for small angles of incidence.

Yueh, H. A.↗
DOE OSTI2020

Evaluation of Constitutive Relationship for Silicon Carbide Cladding Material

A curved representative volume element (RVE) was employed to analyze SiC/SiC tubular composites. By imposing periodicity on the tubular model through appropriate boundary conditions, the elastic constants of the composite was evaluated. It was found that as the radius of the tubular composite increases, the more the flat RVE model and curved RVE model disagreed.

Thandaga Nagaraju, Hemanth↗
NASA NTRS1980

Numerical calculation of steady inviscid full potential compressible flow about wind turbine blades

An exact nonlinear mathematical model that accounts for three-dimensional cascade effects about the inner portions of the rotor blades and compressibility effects about the tip regions of the blades was derived. An artificially time dependent version was iteratively solved by a finite volume technique involving an artificial viscosity and a three-level consecutive mesh refinement. The exact boundary conditions were applied by generating a boundary conforming periodic computation mesh.

Dulikravich, D. S.↗