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 973 records · Page 54

Computation of three-dimensional compressible boundary layers to fourth-order accuracy on wings and fuselages

A solution method, fourth-order accurate in the body-normal direction and second-order accurate in the stream surface directions, to solve the compressible 3-D boundary layer equations is presented. The transformation used, the discretization details, and the solution procedure are described. Ten validation cases of varying complexity are presented and results of calculation given. The results range from subsonic flow to supersonic flow and involve 2-D or 3-D geometries. Applications to laminar flow past wing and fuselage-type bodies are discussed. An interface procedure is used to solve the surface Euler equations with the inviscid flow pressure field as the input to assure accurate boundary conditions at the boundary layer edge. Complete details of the computer program used and information necessary to run each of the test cases are given in the Appendix.

Iyer, Venkit↗

A spectrally accurate boundary-layer code for infinite swept wings

This report documents the development, validation, and application of a spectrally accurate boundary-layer code, WINGBL2, which has been designed specifically for use in stability analyses of swept-wing configurations. Currently, we consider only the quasi-three-dimensional case of an infinitely long wing of constant cross section. The effects of streamwise curvature, streamwise pressure gradient, and wall suction and/or blowing are taken into account in the governing equations and boundary conditions. The boundary-layer equations are formulated both for the attachment-line flow and for the evolving boundary layer. The boundary-layer equations are solved by marching in the direction perpendicular to the leading edge, for which high-order (up to fifth) backward differencing techniques are used. In the wall-normal direction, a spectral collocation method, based upon Chebyshev polynomial approximations, is exploited. The accuracy, efficiency, and user-friendliness of WINGBL2 make it well suited for applications to linear stability theory, parabolized stability equation methodology, direct numerical simulation, and large-eddy simulation. The method is validated against existing schemes for three test cases, including incompressible swept Hiemenz flow and Mach 2.4 flow over an airfoil swept at 70 deg to the free stream.

Pruett, C. David↗

Ignition and combustion of bulk metals under elevated, normal and reduced gravity conditions

This research effort is aimed at providing further insight into this multi-variable dependent phenomena by looking at the effects of gravity on the ignition and combustion behavior of metals. Since spacecraft are subjected to higher-than-1g gravity loads during launch and reentry and to zero-gravity environments while in orbit, the study of ignition and combustion of bulk metals at different gravitational potentials is of great practical concern. From the scientific standpoint, studies conducted under microgravity conditions provide simplified boundary conditions since buoyancy is removed, and make possible the identification of fundamental ignition mechanisms. The effect of microgravity on the combustion of bulk metals has been investigated by Steinberg, et al. on a drop tower simulator. However, no detailed quantitative work has been done on ignition phenomena of bulk metals at lower or higher-than-normal gravitational fields or on the combustion characteristics of metals at elevated gravity. The primary objective of this investigation is the development of an experimental system capable of providing fundamental physical and chemical information on the ignition of bulk metals under different gravity levels. The metals used in the study, iron (Fe), titanium (Ti), zirconium (Zr), magnesium (Mg), zinc (Zn), and copper (Cu) were selected because of their importance as elements of structural metals and their simple chemical composition (pure metals instead of multi-component alloys to avoid complication in morphology and spectroscopic studies). These samples were also chosen to study the two different combustion modes experienced by metals: heterogeneous or surface oxidation, and homogeneous or gas-phase reaction. The experimental approach provides surface temperature profiles, spectroscopic measurements, surface morphology, x-ray spectrometry of metals specimens and their combustion products, and high-speed cinematography of the heating, ignition and combustion stages of the metal specimen. This paper summarizes the results obtained to date from experiments conducted under normal and high-gravity conditions.

Abbud-Madrid, Angel↗

The Effect of Gravity on the Combustion of Bulk Metals

In recent years, metal combustion studies at the University of Colorado have focused on the effects of gravity (g) on the ignition and burning behavior of bulk metals. The impetus behind this effort is the understanding of the ignition conditions and flammability properties of structural metals found in oxygen (O2) systems for space applications. Since spacecraft are subjected to higher-than-1g loads during launch and reentry and to a zero-gravity environment while in orbit, the study of ignition and combustion of bulk metals at different gravitational accelerations is of great practical concern. From the scientific standpoint, studies conducted under low gravity conditions provide simplified boundary conditions, since buoyancy is removed, and make possible the identification of fundamental ignition and combustion mechanisms. This investigation is intended to provide experimental verification of the influence of natural convection on the burning behavior of metals. In addition, the study offers the first findings of the influence of gravity on ignition of bulk metals and on the combustion mechanism and structure of metal-oxygen, vapor-phase diffusion flames in a buoyancy-free environment. Titanium (Ti) and magnesium (Mg) metals were chosen because of their importance as elements of structural materials and their simple chemical composition-pure metals instead of multicomponent alloys to simplify chemical and spectroscopic analyses. In addition, these elements present the two different combustion modes observed in metals: heterogeneous or surface burning (for Ti) and homogeneous or gas-phase reaction (for Mg). Finally, Mg, Ti, and their oxides exhibit a wide range of thermophysical and chemical properties. Metal surface temperature profiles, critical and ignition temperatures, propagation rates, burning times, and spectroscopic measurements are obtained under normal and reduced gravity. Visual evidence of all phenomena is provided by high-speed photography.

Branch, Melvyn C.↗

A Sectored-One-Dimensional Model for Simulating Combustion Instabilities in Premix Combustors

A one-dimensional, CFD based combustor simulation has been developed that exhibits self-excited, thermoacoustic oscillations in premixed combustor geometries that typically have large, abrupt changes in cross sectional area. The combustor geometry is approximated by dividing it into a finite number of one-dimensional sectors. Within each sector, the equations of motion are integrated numerically, along with a species transport and a reaction equation. Across the sectors, mass and energy are conserved, and momentum loss is prescribed using appropriately compatible boundary conditions that account for the area change. The resulting simulation and associated boundary conditions essentially represent a one-dimensional, multi-block technique. Details of the simulation code are presented herein. Results are then shown comparing experimentally observed and simulated operation of a particular combustor rig that exhibited different instabilities at different operating points. It will be shown that the simulation closely matched the rig data in oscillation amplitudes, frequencies, and operating points at which the instabilities occurred. Finally, advantages and limitations of the simulation technique are discussed.

Paxson, Daniel E.↗

A Semianalytical Ion Current Model for Radio Frequency Driven Collisionless Sheaths

We propose a semianalytical ion dynamics model for a collisionless radio frequency biased sheath. The model uses bulk plasma conditions and electrode boundary condition to predict ion impact energy distribution and electrical properties of the sheath. The proposed model accounts for ion inertia and ion current modulation at bias frequencies that are of the same order of magnitude as the ion plasma frequency. A relaxation equation for ion current oscillations is derived which is coupled with a damped potential equation in order to model ion inertia effects. We find that inclusion of ion current modulation in the sheath model shows marked improvements in the predictions of sheath electrical properties and ion energy distribution function.

Bose, Deepak↗

Melt Flow Control in the Directional Solidification of Binary Alloys

Our main project objectives are to develop computational techniques based on inverse problem theory that can be used to design directional solidification processes that lead to desired temperature gradient and growth conditions at the freezing front at various levels of gravity. It is known that control of these conditions plays a significant role in the selection of the form and scale of the obtained solidification microstructures. Emphasis is given on the control of the effects of various melt flow mechanisms on the local to the solidification front conditions. The thermal boundary conditions (furnace design) as well as the magnitude and direction of an externally applied magnetic field are the main design variables. We will highlight computational design models for sharp front solidification models and briefly discuss work in progress toward the development of design techniques for multi-phase volume-averaging based solidification models.

Zabaras, Nicholas↗

Enhanced Mixing in a Rectangular Duct

An experimental investigation of the mixing of non-reacting opposed rows of jets injected normal to a confined rectangular crossflow has been conducted. Planar Mie-scattering was used to measure the time-average concentration distribution of the jet fluid in planes perpendicular to the duct axis. Particular emphasis was placed on the study of closely spaced orifice configurations applicable to the mixing zone of an RQL combustor. Baseline studies were performed of mixing under "ideal" conditions, i.e., plenum fed jets injecting into a crossflow uniform in velocity and turbulence intensity. In addition, more practical ("non-ideal") issues encountered during hardware design were also studied. As in other studies, mixing effectiveness, determined using a spatial unmixedness parameter based on the variance of mean jet concentration distributions, was found to be optimum when the spacing-to-duct-height ratio was inversely proportional to the square root of the jet-to-mainstream momentum-flux ratio. This relationship is suitable for design under ideal flow conditions. Inlet flow boundary conditions of the jet and approach flow (mainstream) were found to strongly influence mixing performance, but no attempt was made to determine optimum performance under non-ideal conditions. The tests performed do offer some guidance as to expected mixing behavior for several common variables likely to be imposed by hardware constraints. Additionally, in this study it was found that for rows of orifices with opposite centerlines inline, mixing was similar for blockages up to 89 percent (previous crossflow mixing studies concerned with dilution zone configurations, blockages were typically less than 50 percent). Lower levels of unmixedness were obtained as a function of downstream location when axial injection length was minimized. Mixing may be enhanced if orifice centerlines of opposed rows are staggered, but blockage must be =50 percent in this configuration. Round hole and "square" orifice shapes had similar performance. Other variations in orifice shape did not substantially augment overall mixing performance. Furthermore an isothermal mixing data set was generated and used by CFDRC as input to a NO(x) inference code.

Liscinsky, D. S.↗

A Computational Model for Rotor-Fuselage Interactional Aerodynamics

A novel unsteady rotor-fuselage interactional aerodynamics model has been developed. This model loosely couples a Generalized Dynamic Wake Theory (GDWT) to a thin-layer Navier-Stokes solution procedure. This coupling is achieved using an unsteady pressure jump boundary condition in the Navier-Stokes model. The new unsteady pressure jump boundary condition models each rotor blade as a moving pressure jump which travels around the rotor azimuth and is applied between two adjacent planes in a cylindrical, non-rotating grid. Comparisons are made between measured and predicted time-averaged and time-accurate rotor inflow ratios. Additional comparisons are made between measured and predicted unsteady surface pressures on the top centerline and sides of the fuselage.

D. Douglas Boyd, Jr.↗

The Land Surface Data Toolkit (LDT v7.2) - A Data Fusion Environment for Land Data Assimilation Systems

The effective applications of land surface models (LSMs) and hydrologic models pose a varied set of data input and processing needs, ranging from ensuring consistency checks to more derived data processing and analytics. This article describes the development of the Land surface Data Toolkit (LDT), which is an integrated framework designed specifically for processing input data to execute LSMs and hydrological models. LDT not only serves as a preprocessor to the NASA Land Information System (LIS), which is an integrated framework designed for multi-model LSM simulations and data assimilation (DA) integrations, but also as a land-surface-based observation and DA input processor. It offers a variety of user options and inputs to processing datasets for use within LIS and stand-alone models. The LDT design facilitates the use of common data formats and conventions. LDT is also capable of processing LSM initial conditions and meteorological boundary conditions and ensuring data quality for inputs to LSMs and DA routines. The machine learning layer in LDT facilitates the use of modern data science algorithms for developing data-driven predictive models. Through the use of an object-oriented framework design, LDT provides extensible features for the continued development of support for different types of observational datasets and data analytics algorithms to aid land surface modeling and data assimilation.

droughts and floods↗

Test and Analysis Correlation of Standard and Hybrid Standard/Thin-ply Composite Notched Test Specimens

Hybrid composite laminates, a combination of standard-and thin-ply layers, are structurally efficient and can offer opportunities for weight savings over composite laminates composed of only standard-thickness plies. Three hybrid laminates and three standard laminates were tested and analyzed to provide a baseline comparison for not only weight savings and strength benefits of the hybrid stacking sequences, but also to provide fidelity for progressive failure modeling techniques using the two different laminate designs. The standard and hybrid composite specimens were designed to have the same thickness and in-plane stiffness properties, as closely as possible, and were loaded with axial tension and axial compression. Small notched specimens were tested under compression and tension loading, and larger intermediate specimens were tested under tension loading, only. The finite element modeling represented the as-designed panel dimensions as well as the real-life testing conditions through various boundary conditions. The progressive failure modeling and analysis technique described, exhibited very good correlation with the hybrid test specimens, and therefore is appropriate for use in future investigations into more appealing stacking sequences. The progressive failure modeling can also be used to determine the proper usage of thin-ply materials in order to determine the most ideal strength benefits over their standard-ply counterparts.

Alana Zahn↗

A compact solution to computational acoustics

This paper demonstrates that the linearized, dimensional Euler equations for acoustic computation can be accurately solved as a set of decoupled first-order wave equations, and that if ordered properly, this system of simple waves has unambiguous, easily implemented boundary conditions, allowing waves of same group speeds to pass through numerical boundaries or comply with wall conditions. Thus, the task of designing a complex multi-dimensional scheme with approximate far-field boundary conditions reduces to the design of higher order schemes for the one-dimensional simple wave equation. A compact finite-difference scheme and a characteristically exact but numerically n(th) order accurate boundary condition are introduced for solving the first order wave equation. Spanning a three-point two-level stencil, this low-dispersion implicit scheme has a third order spatial accuracy when used on nonuniform meshes, fourth order accurate on uniform meshes, and a temporal accuracy of second order due to the choice of trapezoidal integration for algorithmic simplicity. The robustness and accuracy of the scheme are demonstrated through a series of numerical experiments and comparisons with published results. When tested on the one-dimensional wave equation on a uniform grid, this scheme allows a Gaussian wave packet to pass through any finite domain with low numerical dispersion characteristic of a spatially fourth-order scheme and reflections at numerical boundaries maintained below truncation error. On highly stretched and irregular grids, only mild dispersions are found in the solution while solutions by other methods fail or are severely distorted. Yet, this scheme is no more sophisticated to solve or implement than the Crank-Nicolson scheme. This scheme has been tested on four categories of the ICASE/LaRC benchmark problems, which include propagation of acoustic and convective waves in Cartesian and cylindrical domains, reflection of acoustic wave at stationary/moving boundaries, and sound generation by gust-blade interaction.

Fung, K.-Y.↗

Morning boundary layer conditions for shallow to deep convective cloud evolution during the dry season in the central Amazon

Observations of the boundary layer (BL) processes are analyzed statistically for dry seasons of 2 years and in detail, as case studies, for 4 shallow convective days (ShCu) and 4 shallow-to-deep convective days (ShDeep) using a suite of ground-based measurements from the Observation and Modeling of the Green Ocean Amazon (GoAmazon 2014/5) Experiment. The BL stages in ShDeep days, from the nighttime to the cloudy mixing layer stage, are then described in comparison with ShCu days. Atmospheric thermodynamics and dynamics, environmental profiles, and surface turbulent fluxes were employed to compare these two distinct situations for each stage of the BL evolution. Particular attention is given to the morning transition stage, in which the BL changes from stable to unstable conditions in the early morning hours. Results show that the decrease in time duration of the morning transition on ShDeep days is associated with high humidity and well-established vertical wind shear patterns. Higher humidity since nighttime not only contributes to lowering the cloud base during the rapid growth of the BL but also contributes to the balance between radiative cooling and turbulent mixing during nighttime, resulting in higher sensible heat flux in the early morning. The sensible heat flux promotes rapid growth of the well-mixed layer, thus favoring the deeper BL starting from around 08:00 LST (UTC-4 h). Under these conditions, the time duration of morning transition is used to promote convection, having an important effect on the convective BL strength and leading to the formation of shallow cumulus clouds and their subsequent evolution into deep convective clouds. Statistical analysis was used to validate the conceptual model obtained from the case studies. Despite the case-to-case variability, the statistical analyses of the processes in the BL show that the described processes are very representative of cloud evolution during the dry season.

54 ENVIRONMENTAL SCIENCES↗

Scale setting of SU⁡(𝑁) Yang–Mills theory, topology and large-𝑁 volume independence

We set the scale of SU⁡(𝑁) Yang-Mills theories for 𝑁 =3, 5, 8 and in the large-𝑁 limit via gradient flow, as a first step towards the computation of the large-𝑁 Λ-parameter using step scaling. We adopt twisted boundary conditions to achieve large-𝑁 volume reduction and the Parallel Tempering on Boundary Conditions algorithm to tame topological freezing. This setup allows accurate determinations of the gradient-flow scales down to lattice spacings as fine as ∼0.025 fm for all the explored values of 𝑁, a regime that has never been reached with ergodic algorithms. Moreover, we are able to precisely estimate the finite-size systematics related to topological freezing, and to show the suppression of finite-volume effects expected by virtue of large-𝑁 twisted volume reduction.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Can changes in deformation regimes be inferred from crystallographic preferred orientations in polar ice?

Creep due to ice flow is generally thought to be the main cause for the formation of crystallographic preferred orientations (CPOs) in polycrystalline anisotropic ice. However, linking the development of CPOs to the ice flow history requires a proper understanding of the ice aggregate's microstructural response to flow transitions. In this contribution the influence of ice deformation history on the CPO development is investigated by means of full-field numerical simulations at the microscale. We simulate the CPO evolution of polycrystalline ice under combinations of two consecutive deformation events up to high strain, using the code VPFFT (visco-plastic fast Fourier transform algorithm) within ELLE. A volume of ice is first deformed under coaxial boundary conditions, which results in a CPO. The sample is then subjected to different boundary conditions (coaxial or non-coaxial) in order to observe how the deformation regime switch impacts the CPO. The model results indicate that the second flow event tends to destroy the first, inherited fabric with a range of transitional fabrics. However, the transition is slow when crystallographic axes are critically oriented with respect to the second imposed regime. Therefore, interpretations of past deformation events from observed CPOs must be carried out with caution, particularly in areas with complex deformation histories.

54 ENVIRONMENTAL SCIENCES↗

Numerical calculations of flow fields

The solutions to the equations of motion for inviscid fluid flow around a pointed elliptic cone at incidence are presented. The numerical method used, MacCormack's second order preferential predictor-corrector finite difference approximation, is applied to the fluid flow equations derived in conservation-law form. The entropy boundary condition, hitherto unused for elliptic cone problems, is investigated and compared to reflection boundary condition solutions. The stagnation streamline movement of the inclined elliptic cone is noted and surface pressure coefficients are plotted. Also presented are solutions for an elliptic cone and a circular cone at zero incidence and a circular cone at a small angle of attack. Comparisons are made between these present solutions and previously published theory.

Anderson, D. M.↗

Structure of the conical flowfield about external axial corners

A theoretical and experimental investigation of the flowfield about symmetrical and asymmetrical external axial corners formed by the juncture of swept compressive wedges is presented. A study is made of the nature of inviscid cross flow stagnation points and their dependence on certain flow parameters is derived. The importance of imposing the physically relevant boundary condition for numerical simulations is discussed. For the asymmetrical wedge configuration, the proper boundary condition along the axial corner is shown to be a conical analog of Prandtl-Meyer flow. Further insight into the flow structure is revealed by the experimental results.

Salas, M. D.↗

The application of MINIQUASI to thermal program boundary and initial value problems

The feasibility of applying the solution techniques of Miniquasi to the set of equations which govern a thermoregulatory model is investigated. For solving nonlinear equations and/or boundary conditions, a Taylor Series expansion is required for linearization of both equations and boundary conditions. The solutions are iterative and in each iteration, a problem like the linear case is solved. It is shown that Miniquasi cannot be applied to the thermoregulatory model as originally planned.

Source record↗