Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Boundary Conditions”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,171 records · Page 65

Does fast magnetic reconnection exist?

The main features of the Priest-Forbes (1986) and Priest-Lee (1990) models of magnetic reconnection in astrophysical plasmas are discussed, and the Priest-Lee model is generalized to include inflow pressure gradients and thus different regimes of reconnection. It is shown that different scaling results can be obtained depending on the boundary conditions. These results are compared to the ones observed in the numerical experiments of Biskamp (1986) and Lee and Fu (1986). It is concluded that numerical experiments with suitably designed boundary conditions are likely to exhibit fast reconnection, and that such reconnection is a common process in astrophysical and space plasmas.

Priest, E. R.↗

Aeroheating Mapping to Thermal Model for Autonomous Aerobraking Capability

Thermal modeling has been performed to evaluate the potential for autonomous aerobraking of a spacecraft in the atmosphere of a planet. As part of this modeling, the aeroheating flux during aerobraking must be applied to the spacecraft solar arrays to evaluate their thermal response. On the Mars Reconnaissance Orbiter (MRO) mission, this was done via two separate thermal models and an extensive suite of mapping scripts. That method has been revised, and the thermal analysis of an aerobraking pass can now be accomplished via a single thermal model, using a new capability in the Thermal Desktop software. This capability, Boundary Condition Mapper, has the ability to input heating flux files that vary with time, position on the solar array, and with the skin temperature. A recently added feature to the Boundary Condition Mapper is that this module can also utilize files that describe the variation of aeroheating over the surface with atmospheric density (rather than time); this is the format of the MRO aeroheating files. This capability has allowed a huge streamlining of the MRO thermal process, simplifying the procedure for importing new aeroheating files and trajectory information. The new process, as well as the quantified time savings, is described.

Amundsen, Ruth M.↗

Numerical Models of Human Circulatory System under Altered Gravity: Brain Circulation

A computational fluid dynamics (CFD) approach is presented to model the blood flow through the human circulatory system under altered gravity conditions. Models required for CFD simulation relevant to major hemodynamic issues are introduced such as non-Newtonian flow models governed by red blood cells, a model for arterial wall motion due to fluid-wall interactions, a vascular bed model for outflow boundary conditions, and a model for auto-regulation mechanism. The three-dimensional unsteady incompressible Navier-Stokes equations coupled with these models are solved iteratively using the pseudocompressibility method and dual time stepping. Moving wall boundary conditions from the first-order fluid-wall interaction model are used to study the influence of arterial wall distensibility on flow patterns and wall shear stresses during the heart pulse. A vascular bed modeling utilizing the analogy with electric circuits is coupled with an auto-regulation algorithm for multiple outflow boundaries. For the treatment of complex geometry, a chimera overset grid technique is adopted to obtain connectivity between arterial branches. For code validation, computed results are compared with experimental data for steady and unsteady non-Newtonian flows. Good agreement is obtained for both cases. In sin-type Gravity Benchmark Problems, gravity source terms are added to the Navier-Stokes equations to study the effect of gravitational variation on the human circulatory system. This computational approach is then applied to localized blood flows through a realistic carotid bifurcation and two Circle of Willis models, one using an idealized geometry and the other model using an anatomical data set. A three- dimensional anatomical Circle of Willis configuration is reconstructed from human-specific magnetic resonance images using an image segmentation method. The blood flow through these Circle of Willis models is simulated to provide means for studying gravitational effects on the brain circulation under auto-regulation.

Kim, Chang Sung↗

SubTSBR to tackle high noise and outliers for data-driven discovery of differential equations

Data-driven discovery of differential equations has been an emerging research topic. We propose a novel algorithm subsampling-based threshold sparse Bayesian regression (SubTSBR) to tackle high noise and outliers. The subsampling technique is used for improving the accuracy of the Bayesian learning algorithm. It has two parameters: subsampling size and the number of subsamples. When the subsampling size increases with fixed total sample size, the accuracy of our algorithm goes up and then down. When the number of subsamples increases, the accuracy of our algorithm keeps going up. We demonstrate how to use our algorithm step by step and compare our algorithm with threshold sparse Bayesian regression (TSBR) for the discovery of differential equations. We show that our algorithm produces better results. We also discuss the merits of discovering differential equations from data and demonstrate how to discover models with random initial and boundary condition as well as models with bifurcations. The numerical examples are: (1) predator-prey model with noise, (2) shallow water equations with outliers, (3) heat diffusion with random initial and boundary condition, and (4) fish-harvesting problem with bifurcations.

97 MATHEMATICS AND COMPUTING↗

A Numerical Analysis of the Transient Response of an Ablation System Including Effects of Thermal Nonequilibrium, Mass Transfer and Chemical Kinetics

The differential equations governing the transient response of a one-dimensional ablative thermal protection system undergoing stagnation ablation are derived. These equations are for thermal nonequilibrium effects between the pyrolysis gases and the char layer and kinetically controlled chemical reactions and mass transfer between the pyrolysis gases and the char layer. The boundary conditions are written for the particular case of stagnation heating with surface removal by oxidation or sublimation and pyrolysis of the uncharred layer occurring in a plane. The governing equations and boundary conditions are solved numerically using the modified implicit method (Crank-Nicolson method). Numerical results are compared with exact solutions for a number of simplified cases. The comparison is favorable in each instance.

Clark, R. K.↗

The propagation of plane waves and higher order acoustic modes in multisectioned ducts

A theoretical and experimental study of acoustic propagation in an anechoically terminated multisectioned duct was performed. A unique source array consisting of two concentric rings of sources, providing phase and amplitude control in the radial as well as circumferential direction, was developed to generate plane waves and both spinning and non-spinning higher order modes. Measurement of attenuation and radial mode shapes were taken with various finite length liners inserted between the hard wall sections of a duct with an anechoic termination. A search technique was developed to find the complex eigenvalues for a liner under the assumption of a locally reacting boundary condition. The experimental results were compared with a theoretical analysis which includes the modal transmission and reflection at the interface of each duct section, and this indicated that the local reaction boundary condition is valid for three liner configurations.

Wyerman, B. R.↗

Finite strain calculations of continental deformation. I - Method and general results for convergent zones. II - Comparison with the India-Asia collision zone

The present investigation has the objective to perform numerical experiments on a rheologically simple continuum model for the continental lithosphere. It is attempted to obtain a better understanding of the dynamics of continental deformation. Calculations are presented of crustal thickness distributions, stress, strain, strain rate fields, latitudinal displacements, and finite rotations, taking into account as basis a model for continental collision which treats the litoshphere as a thin viscous layer subject to indenting boundary conditions. The results of this paper support the conclusions of England and McKenzie (1982) regarding the role of gravity in governing the deformation of a thin viscous layer subject to indenting boundary conditions. The results of the experiments are compared with observations of topography, stress and strain rate fields, and palaeomagnetic latitudinal displacements in Asia.

Houseman, G.↗

Large deflection random response of cross-ply laminated plates with elastically restrained edges and initial imperfections

The large deflection random response of symmetrically laminated cross-ply rectangular thin plates subjected to random excitation is studied. The out-of-plane boundary conditions are such that all the edges are rigidly supported against translation, but elastically restrained against rotation. The plate is also assumed to have a small initial imperfection. The assumed membrane boundary conditions are such that all the edges are free from normal and tangential forces in the plane of the plate. Mean-square deflections and mean-square strains are determined for a three-layered cross-ply laminate.

Prasad, C. B.↗

A new technique for simulating composite material

This project dealt with the development on new methodologies and algorithms for the multi-spectrum electromagnetic characterization of large scale nonmetallic airborne vehicles and structures. A robust, low memory, and accurate methodology was developed which is particularly suited for modern machine architectures. This is a hybrid finite element method that combines two well known numerical solution approaches. That of the finite element method for modeling volumes and the boundary integral method which yields exact boundary conditions for terminating the finite element mesh. In addition, a variety of high frequency results were generated (such as diffraction coefficients for impedance surfaces and material layers) and a class of boundary conditions were developed which hold promise for more efficient simulations. During the course of this project, nearly 25 detailed research reports were generated along with an equal number of journal papers. The reports, papers, and journal articles are listed in the appendices along with their abstracts.

Volakis, John L.↗

Improved Hydrology over Peatlands in a Global Land Modeling System

Peatlands of the Northern Hemisphere represent an important carbon pool that mainly accumulated since the last ice age under permanently wet conditions in specific geological and climatic settings. The carbon balance of peatlands is closely coupled to water table dynamics. Consequently, the future carbon balance over peatlands is strongly dependent on how hydrology in peatlands will react to changing boundary conditions, e.g. due to climate change or regional water level drawdown of connected aquifers or streams. Global land surface modeling over organic-rich regions can provide valuable global-scale insights on where and how peatlands are in transition due to changing boundary conditions. However, the current global land surface models are not able to reproduce typical hydrological dynamics in peatlands well. We implemented specific structural and parametric changes to account for key hydrological characteristics of peatlands into NASA's GEOS-5 Catchment Land Surface Model (CLSM, Koster et al. 2000). The main modifications pertain to the modeling of partial inundation, and the definition of peatland-specific runoff and evapotranspiration schemes. We ran a set of simulations on a high performance cluster using different CLSM configurations and validated the results with a newly compiled global in-situ dataset of water table depths in peatlands. The results demonstrate that an update of soil hydraulic properties for peat soils alone does not improve the performance of CLSM over peatlands. However, structural model changes for peatlands are able to improve the skill metrics for water table depth. The validation results for the water table depth indicate a reduction of the bias from 2.5 to 0.2 m, and an improvement of the temporal correlation coefficient from 0.5 to 0.65, and from 0.4 to 0.55 for the anomalies. Our validation data set includes both bogs (rain-fed) and fens (ground and/or surface water influence) and reveals that the metrics improved less for fens. In addition, a comparison of evapotranspiration and soil moisture estimates over peatlands will be presented, albeit only with limited ground-based validation data. We will discuss strengths and weaknesses of the new model by focusing on time series of specific validation sites.

Bechtold, M.↗

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↗

A Poisson–Nernst–Planck Model of Ion Transport and Interface Segregation in Metal–Insulator–Semiconductor Structures and Solar Cells

A numerical model that describes the transport of mobile ionic species in metal–insulator–semiconductor (MIS) and photovoltaic (PV) devices subject to temperature and voltage stress is presented. The finite element method (FEM) is used to solve the Nernst–Planck equation while imposing Poisson's equation self-consistently as a restriction for the electrostatic potential. This allows the contribution of the ionic species to the potential to be taken into account. Using a variational formulation eases the implementation of diverse boundary conditions, including the incorporation of segregation kinetics at the device interfaces. Segregation across the dielectric–semiconductor interface is relevant to modeling the electronic device degradation in systems where contamination reaches the semiconductor. The model in closed systems with no-flux boundary conditions is validated first. In the limiting case of low contamination levels with respect to the gate bias, the FEM solution matches analytically derived approximations. Then, the implementation is broadened to include an open boundary at the dielectric–semiconductor interface to account for leakage of ions. The predicted time dependence of the flatband voltage in Na-contaminated MIS test structures agrees well with measurements. This model successfully captures the role of long-range ion transport at concentrations of relevance to electronic and PV device instability and neuromorphic computing.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Chiral symmetry and scale invariance breaking in spin chains

The effects of the Dzyaloshinsky-Moriya interaction on finite size one-dimensional (1D) magnetic chains are investigated as a function of their length. The magnetic configuration the system adopts for varying boundary conditions are explored analytically, which leads to the appearance of chiral configurations that play a crucial role. The coercive and exchange bias fields show an unexpected chain length dependence, caused by the boundary conditions and by chiral symmetry breaking, which in turn leads to the breakdown of scale-invariance. Our treatment yields results in agreement with experimental evidence and ongoing research on phthalocyanine iron chains bonded to hydrogen.

36 MATERIALS SCIENCE↗

Electric field effects on galactic cosmic rays at the heliosphere boundary

The boundary condition for galactic cosmic rays at the outer edge of the modulation region was examined. With few exceptions authors have in the past uncritically used a time-independent, spherically-symmetric condition. It was found that, depending on the mechanism by which the heliospheric magnetic field connects to the general galactic magnetic field, the induced electric field may have important effects on the boundary condition. Sample calculations are presented which illustrate the potentially large, nonspherically-symmetric effects which may be expected.

Jokipii, J. R.↗

Numerical simulation of the atmospheric response ito the time-varying El Nino SST anomalies during May 1982 through October 1983

An attempt is made to simulate the atmospheric circulation anomalies corresponding to the observed SST anomalies in the Pacific Ocean for the 18-month period of May 1982 through October 1983. A GCM is first integrated for 25 months with monthly climatological boundary conditions of SST, soil moisture, sea, ice, and albedo. Starting from day 165 of this 'control' integration, which corresponds to May 1, the 18-month integration is carried out with the same boundary conditions except that the observed monthly SST anomalies for May 1982-October 1983 are added to the climatological values in the Pacific from 40 S to 60 N. The evolution of the model-simulated circulation and rainfall anomalies are compared to actual observations for the same period, and remarkable agreement is found.

Fennessy, M. J.↗

Ground test program for a full-size solar dynamic heat receiver

Test hardware, facilities, and procedures were developed to conduct ground testing of a full size, solar dynamic heat receiver in a partially simulated, low Earth orbit environment. The heat receiver was designed to supply 102 kW of thermal energy to a helium and xenon gas mixture continuously over a 94 minute orbit, including up to 36 minutes of eclipse. The purpose of the test program was to quantify the receiver thermodynamic performance, its operating temperatures, and thermal response to changes in environmental and power module interface boundary conditions. The heat receiver was tested in a vacuum chamber using liquid nitrogen cold shrouds and an aperture cold plate. Special test equipment were designed to provide the required ranges in interface boundary conditions that typify those expected or required for operation as part of the solar dynamic power module on the Space Station Freedom. The support hardware includes an infrared quartz lamp heater with 30 independently controllable zones and a closed Brayton cycle engine simulator to circulate and condition the helium xenon gas mixture. The test article, test support hardware, facilities, and instrumentation developed to conduct the ground test program are all described.

Sedgwick, L. M.↗

Ground test program for a full-size solar dynamic heat receiver

Test hardware, facilities, and procedures were developed to conduct ground testing of a full-size, solar dynamic heat receiver in a partially simulated, low earth orbit environment. The heat receiver was designed to supply 102 kW of thermal energy to a helium and xenon gas mixture continuously over a 94 minute orbit, including up to 36 minutes of eclipse. The purpose of the test program was to quantify the receiver thermodynamic performance, its operating temperatures, and thermal response to changes in environmental and power module interface boundary conditions. The heat receiver was tested in a vacuum chamber using liquid nitrogen cold shrouds and an aperture cold plate. Special test equipment was designed to provide the required ranges in interface boundary conditions that typify those expected or required for operation as part of the solar dynamic power module on the Space Station Freedom. The support hardware includes an infrared quartz lamp heater with 30 independently controllable zones and a closed-Brayton cycle engine simulator to circulate and condition the helium-xenon gas mixture. The test article, test support hardware, facilities, and instrumentation developed to conduct the ground test program are all described.

Sedgwick, L. M.↗

Is the earth's magnetotail balloon unstable?

In the past, the onset of magnetospheric substorms has been attributed to the plasma tearing mode instabilities. This paper investigates the ideal MHD ballooning instability of the near-and middle-tail magnetosphere region, as a first step toward determining whether it could trigger the tearing mode, by using the energy principle to investigate whether standard 2D tail models with the 'hard' ionospheric boundary condition are unstable to ballooning instability. Numerical results are presented for compressible ballooning modes that are symmetric about the center of the current sheet. It is shown that, for such a hard boundary condition, no reasonable magnetotail configuration exists that would be unstable to compressible symmetric ballooning but stable against interchange.

Lee, D.-Y.↗