Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “embedded boundary”

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 199 records · Page 11

Transient, Small-Scale Field-Aligned Currents in the Plasma Sheet Boundary Layer During Storm Time Substorms

We report on field-aligned current observations by the four Magnetospheric Multiscale (MMS) spacecraft near the plasma sheet boundary layer (PSBL) during two major substorms on 23 June 2015. Small-scale field-aligned currents were found embedded in fluctuating PSBL flux tubes near the Separatrix region. We resolve, for the first time, short-lived earthward (downward) intense field-aligned current sheets with thicknesses of a few tens of kilometers, which are well below the ion scale, on flux tubes moving equatorward earth ward during outward plasma sheet expansion. They coincide with upward field-aligned electron beams with energies of a few hundred eV. These electrons are most likely due to acceleration associated with a reconnection jet or high-energy ion beam-produced disturbances. The observations highlight coupling of multiscale processes in PSBL as a consequence of magnetotail reconnection.

Nakamura, R.↗

Integrated Data Modeling and Simulation on the Joint Polar Satellite System Program

The Joint Polar Satellite System is a modern, large-scale, complex, multi-mission aerospace program, and presents a variety of design, testing and operational challenges due to: (1) System Scope: multi-mission coordination, role, responsibility and accountability challenges stemming from porous/ill-defined system and organizational boundaries (including foreign policy interactions) (2) Degree of Concurrency: design, implementation, integration, verification and operation occurring simultaneously, at multiple scales in the system hierarchy (3) Multi-Decadal Lifecycle: technical obsolesce, reliability and sustainment concerns, including those related to organizational and industrial base. Additionally, these systems tend to become embedded in the broader societal infrastructure, resulting in new system stakeholders with perhaps different preferences (4) Barriers to Effective Communications: process and cultural issues that emerge due to geographic dispersion and as one spans boundaries including gov./contractor, NASA/Other USG, and international relationships.

Roberts, Christopher J.↗

Global Sensitivity Analysis of Simulated Remote Sensing Polarimetric Observations Over Snow

This study presents a detailed theoretical assessment of the information content of passive polarimetric observations over snow scenes, using a global sensitivity analysis (GSA) method. Conventional sensitivity studies focus on varying a single parameter while keeping all other parameters fixed. In contrast, the GSA correctly addresses the covariance of state parameters across their entire parameter space, hence favoring a more correct interpretation of inversion algorithms and the optimal design of their state vectors. The forward simulations exploit a vector radiative transfer model to obtain the Stokes vector emerging at the top of the atmosphere for different solar zenith angles, when the bottom boundary consists of a vertically resolved snowpack of non-spherical grains. The presence of light-absorbing particulates (LAPs), either embedded in the snow or aloft in the atmosphere above in the form of aerosols, is also considered. The results are presented for a set of wavelengths spanning the visible (VIS), near-infrared (NIR), and shortwave infrared (SWIR) region of the spectrum. The GSA correctly captures the expected, high sensitivity of the reflectance to LAPs in the VIS–NIR and to grain size at different depths in the snowpack in the NIR–SWIR. With adequate viewing geometries, mono-angle measurements of total reflectance in the VIS–SWIR (akin to those of the Moderate Resolution Imaging Spectroradiometer, MODIS) resolve grain size in the top layer of the snowpack sufficiently well. The addition of multi-angle polarimetric observations in the VIS–NIR provides information on grain shape and microscale roughness. The simultaneous sensitivity in the VIS–NIR to both aerosols and snow-embedded impurities can be disentangled by extending the spectral range to the SWIR, which contains information on aerosol optical depth while remaining essentially unaffected when the same particulates are mixed with the snow. Multi-angle polarimetric observations can therefore (i) effectively partition LAPs between the atmosphere and the surface, which represents a notorious challenge for snow remote sensing based on measurements of total reflectance only and (ii) lead to better estimates of grain shape and roughness and, in turn, the asymmetry parameter, which is critical for the determination of albedo. The retrieval uncertainties are minimized when the degree of linear polarization is used in place of the polarized reflectance. The Sobol indices, which are the main metric for the GSA, were used to select the state parameters in retrievals performed on data simulated for multiple instrument configurations. Improvements in retrieval quality with the addition of measurements of polarization, multi-angle views, and different spectral channels reflect the information content, identified by the Sobol indices, relative to each configuration. The results encourage the development of new remote sensing algorithms that fully leverage multi-angle and polarimetric capabilities of modern remote sensors. They can also aid flight planning activities, since the optimal exploitation of the information content of multi-angle measurements depends on the viewing geometry. The better characterization of surface and atmospheric parameters in snow-covered regions advances research opportunities for scientists of the cryosphere and ultimately benefits albedo estimates in climate models.

remote sensing↗

An Embedded Statistical Method for Coupling Molecular Dynamics and Finite Element Analyses

The coupling of molecular dynamics (MD) simulations with finite element methods (FEM) yields computationally efficient models that link fundamental material processes at the atomistic level with continuum field responses at higher length scales. The theoretical challenge involves developing a seamless connection along an interface between two inherently different simulation frameworks. Various specialized methods have been developed to solve particular classes of problems. Many of these methods link the kinematics of individual MD atoms with FEM nodes at their common interface, necessarily requiring that the finite element mesh be refined to atomic resolution. Some of these coupling approaches also require simulations to be carried out at 0 K and restrict modeling to two-dimensional material domains due to difficulties in simulating full three-dimensional material processes. In the present work, a new approach to MD-FEM coupling is developed based on a restatement of the standard boundary value problem used to define a coupled domain. The method replaces a direct linkage of individual MD atoms and finite element (FE) nodes with a statistical averaging of atomistic displacements in local atomic volumes associated with each FE node in an interface region. The FEM and MD computational systems are effectively independent and communicate only through an iterative update of their boundary conditions. With the use of statistical averages of the atomistic quantities to couple the two computational schemes, the developed approach is referred to as an embedded statistical coupling method (ESCM). ESCM provides an enhanced coupling methodology that is inherently applicable to three-dimensional domains, avoids discretization of the continuum model to atomic scale resolution, and permits finite temperature states to be applied.

Saether, E.↗

NBODY - a multipurpose trajectory optimization computer program

Documentation of the NBODY trajectory optimization program is presented in the form of a mathematical development plus a user's manual. Optimal multistage-launch ascent trajectories may be determined by variational thrust steering during the upper phase. Optimal low-thrust interplanetary spacecraft trajectories may also be calculated with solar power or constant power, all-propulsion or embedded coast arcs, fixed or optimal thrust angles, and a variety of terminal end conditions. A hybrid iteration scheme solves the boundary-value problem, while either transversality conditions or a univariate search scheme optimize vehicle or trajectory parameters.

Strack, W. C.↗

The free-wake prediction of rotor hover performance using a vortex embedding method

A method is developed to predict the rotor hover performance. This method solves the compressible mass conservation equation much like current full potential codes and can therefore predict the transonic flows on a rotor. However, the newly developed approach also allows for the free convection of shed vorticity and permits the computation of the entire hover wake system. The method uses a vortex embedding scheme in potential flow and has been implemented in a computer code, HELIX -I. To predict power we implement a simple boundary layer and two different induced-drag integration schemes. The induced-drag is obtained from surface pressure integration and an energy flux integral. Comparisons between computations and experiment show good agreement for the prediction of power polars, surface pressure distribution, and tip vortex geometry.

Ramachandran, K.↗

Embedded function methods for compressible high speed turbulent flow

This is the final report on the work performed on the grant 'Embedded Function Methods for Compressible High Speed Turbulent Flow' carried out at Lehigh University during the contract period from September, 1987, to October of 1991. Work has continued at Lehigh on this project on an unfunded basis to the present. The original proposed work had two separate thrusts which were associated with developing embedded function methods in order to obviate the need to expend computational resources on turbulent wall layers in Navier Stokes and boundary-layer calculations. Previous work on the incompressible problem had indicated that this could be done successfully for two-dimensional and three-dimensional incompressible flows. The central objective here was to extend the basic approach to the high speed compressible problem.

Walker, J. D. A.↗

Estimation of 3-D Cloud Effects on TOMS Satellite Retrieval of Surface UV Irradiance

To improve surface UV irradiance retrieval from the Total Ozone Mapping Spectrometer (TOMS) we simulate errors of the TOMS cloud correction algorithm for summertime broken cloud conditions. Cloud scenes (50 km by 50 km) are modeled by a normal random (Gaussian) field with a fixed lower boundary and conservative scattering. The model relates stochastic field characteristics with the cloud amount, mean cloud diameter and aspect ratio. Clouds are embedded into Rayleigh atmosphere with standard ozone profile. Radiative transfer calculations of the radiance at the top of the atmosphere and irradiance at the surface were performed using 3-D Monte Carlo (MC) code. The results are averaged over the satellite field of view on the surface (50 km by 50 km) and compared with TOMS predicted surface irradiance for the same scene reflectance. The TOMS algorithm assumes horizontally homogeneous Cl-type cloud between 3 km and 5.5 km. The effective optical depth is determined by fitting observed (MC) radiance at 380 nm. Having the same radiance at the satellite the homogeneous and broken cloud models predict different average irradiances at the surface. This is due to the differences in Bidirectional Reflection Distribution Function (BRDF) for homogeneous and broken cloud scenes with the same hemispherical albedo. For typical TOMS observational geometry at mid-latitudes the simulated single pixels errors may be as large as +/- 20%. Qualitatively these errors are due to the dominance of the non-horizontal cloud surfaces, which are not accounted for in the homogeneous cloud model. However, due to high variability of the real cloud shapes and types it is unclear how these single pixel errors would affect TOMS time-integrated UV exposure over extended periods (weeks to months) for different regions.

Krotkov, Nickolay A.↗

Turboelectric Distributed Propulsion Engine Cycle Analysis for Hybrid-Wing-Body Aircraft

Meeting NASA's N+3 goals requires a fundamental shift in approach to aircraft and engine design. Material and design improvements allow higher pressure and higher temperature core engines which improve the thermal efficiency. Propulsive efficiency, the other half of the overall efficiency equation, however, is largely determined by the fan pressure ratio (FPR). Lower FPR increases propulsive efficiency, but also dramatically reduces fan shaft speed through the combination of larger diameter fans and reduced fan tip speed limits. The result is that below an FPR of 1.5 the maximum fan shaft speed makes direct drive turbines problematic. However, it is the low pressure ratio fans that allow the improvement in propulsive efficiency which, along with improvements in thermal efficiency in the core, contributes strongly to meeting the N+3 goals for fuel burn reduction. The lower fan exhaust velocities resulting from lower FPRs are also key to meeting the aircraft noise goals. Adding a gear box to the standard turbofan engine allows acceptable turbine speeds to be maintained. However, development of a 50,000+ hp gearbox required by fans in a large twin engine transport aircraft presents an extreme technical challenge, therefore another approach is needed. This paper presents a propulsion system which transmits power from the turbine to the fan electrically rather than mechanically. Recent and anticipated advances in high temperature superconducting generators, motors, and power lines offer the possibility that such devices can be used to transmit turbine power in aircraft without an excessive weight penalty. Moving to such a power transmission system does more than provide better matching between fan and turbine shaft speeds. The relative ease with which electrical power can be distributed throughout the aircraft opens up numerous other possibilities for new aircraft and propulsion configurations and modes of operation. This paper discusses a number of these new possibilities. The Boeing N2 hybrid-wing-body (HWB) is used as a baseline aircraft for this study. The two pylon mounted conventional turbofans are replaced by two wing-tip mounted turboshaft engines, each driving a superconducting generator. Both generators feed a common electrical bus which distributes power to an array of superconducting motor-driven fans in a continuous nacelle centered along the trailing edge of the upper surface of the wing-body. A key finding was that traditional inlet performance methodology has to be modified when most of the air entering the inlet is boundary layer air. A very thorough and detailed propulsion/airframe integration (PAI) analysis is required at the very beginning of the design process since embedded engine inlet performance must be based on conditions at the inlet lip rather than freestream conditions. Examination of a range of fan pressure ratios yielded a minimum Thrust-specific-fuel-consumption (TSFC) at the aerodynamic design point of the vehicle (31,000 ft /Mach 0.8) between 1.3 and 1.35 FPR. We deduced that this was due to the higher pressure losses prior to the fan inlet as well as higher losses in the 2-D inlets and nozzles. This FPR is likely to be higher than the FPR that yields a minimum TSFC in a pylon mounted engine. 1

Felder, James L.↗

Flow and heat transfer measurements along a cooled supersonic diffuser

An experimental investigtion was conducted to ascertain the mean flowfield, including shock wave structure, separated flow regions, turbulent boundary-layer growth, static pressure variations, wall heat transfer, and shear stresses in a second-throat, axisymmetric, supersonic diffuser with wall cooling. The diffuser inlet Mach number of the heated air flow was 3.76, the stagnation pressure was 6.8 atm, the ratio of wall to total gas temperature was 0.44, and the diffuser discharged to the atmosphere. The complex flowfield involved deceleration and acceleration regions, supersonic and embedded subsonic regions, and strong viscous regions with relatively large radial and axial variations. The heat transfer and wall static pressure distributions were remarkably similar, and heat transfer rates were high locally at oblique shock/turbulent boundary-layer interactions, in the pseudoshock region, and in the separation region in the diffuser outlet section.

Back, L. H.↗

The effect of closed boundary conditions on a stationary dynamo

One of two boundary conditions generally assumed in solutions of the dynamo equation is related to the disappearance of the azimuthal field at the boundary. Parker (1984) points out that for the realization of this condition the field must escape freely through the surface. Escape requires that the field be detached from the gas in which it is embedded. In the case of the sun, this can be accomplished only through reconnection in the tenuous gas above the visible surface. Parker concludes that the observed magnetic activity on the solar surface permits at most three percent of the emerging flux to escape. He arrives at the conclusion that, instead of B(phi) = 0, the partial derivative of B(phi) to r is equal to zero. The present investigation is concerned with the effect of changing the boundary condition according to Parker's conclusion. Implications for the solar convection zone are discussed.

Choudhuri, A. R.↗

Analysis of twisting stiffness for a multifiber composite layer

The twisting stiffness of a rectangular cross section consisting of a single row of solid circular cross-section fibers embedded in a matrix is analyzed. The problem is formulated as a Dirichlet torsion problem of a multielement region and solved by the boundary-point least-squares method. Numerical results for a single-fiber square cross section compare favorably with previous relaxation-method results. New numerical results for three and five-fiber composites suggest that the torsional rigidity of a multifiber composite can be approximated from the torsional rigidities of single and three-fiber models.

Bert, C. W.↗

Relaminarization of turbulent flow on a flat plate by localized surface heating

Relaminarization of a turbulent boundary layer in air on a flat plate is demonstrated experimentally using localized surface heating. Localized heating is achieved by electrically heating a wire embedded in a thermally insulated substrate (Space Shuttle Tile) on the surface. The stability of the flow downstream of the applied control point increases with decreasing stream temperature in the flow direction. The mean and perturbation velocity profiles without control show that the flow is turbulent. With control, these profiles sequentially change from intermittently turbulent to a fully laminar state. In the relaminarization stage, the turbulent energy is dissipated by molecular transport due to viscous and conductivity mechanisms. The new profile adjusts to a lower Reynolds number based on the momentum thickness than that of the previous turbulent state.

Maestrello, Lucio↗

Toward Understanding Wake Vortices and Atmospheric Turbulence Interactions Using Large-Eddy Simulation

The vortices produced by an aircraft in flight are a complex phenomena created from a 'sheet of vorticity' leaving the trailing edge of the aircraft surfaces. This sheet tends to roll-up into two counter-rotating vortices. After a few spans downstream of the aircraft, the roll-up process is complete and the vortex pair may be characterized in a simple manner for modeling purposes. Our research will focus on what happens to these post roll-up vortices in the vicinity of an airport terminal. As the aircraft wake vortices descend, they are transported by the air mass which they are embedded and are decayed by both internal and external processes. In the vicinity of the airport, these external influences are usually due to planetary boundary layer (PBL) turbulence. Using large-eddy simulation (LES), one may simulate a variety of PBL conditions. In the LES method, turbulence is generated in the PBL as a response to surface heat flux, horizontal pressure gradient, wind shear, and/or stratification, and may produce convective or unstably stratified, neutral, or stably stratified PBL's. Each of these PBL types can occur during a typical diurnal cycle of the PBL. Thus it is important to be able to characterize these conditions with the LES method. Once this turbulent environment has been generated, a vortex pair will be introduced and the interactions are observed. The objective is to be able to quantify the PBL turbulence vortex interaction and be able to draw some conclusions of vortex behavior from the various scale interactions. This research is ongoing, and we will focus on what has been accomplished to date and the future direction of this research. We will discuss the model being used, show results that validate its use in the PBL, and present a nested-grid method proposed to analyze the entire PBL and vortex pair simultaneously.

DeCroix, David↗

The Stress-strain Behavior of Polymer-Nanotube Composites from Molecular Dynamics Simulations

Stress-strain curves of polymer-carbon nanotube composites are derived from molecular dynamics simulations of a single-walled carbon nanotube embedded in polyethylene. A comparison is made between the response to mechanical loading of a composite with a long, continuous nanotube (replicated via periodic boundary conditions) and the response of a composite with a short, discontinuous nanotube. Both composites are mechanically loaded in the direction of and transverse to the NT axis. The long-nanotube composite shows an increase in the stiffness relative to the polymer and behaves anisotropically under the different loading conditions. The short-nanotube composite shows no enhancement relative to the polymer, most probably because of its low aspect ratio. The stress-strain curves are compared with rule-of-mixtures predictions.

Frankland, S. J. V.↗

Laminarization of Turbulent Boundary Layer on Flexible and Rigid Surfaces

An investigation of the control of turbulent boundary layer flow over flexible and rigid surfaces downstream of a concave-convex geometry has been made. The concave-convex curvature induces centrifugal forces and a pressure gradient on the growth of the turbulent boundary layer. The favorable gradient is not sufficient to overcome the unfavorable; thus, the net effect is a destabilizing, of the flow into Gortler instabilities. This study shows that control of the turbulent boundary layer and structural loading can be successfully achieved by using localized surface heating because the subsequent cooling and geometrical shaping downstream over a favorable pressure gradient is effective in laminarization of the turbulence. Wires embedded in a thermally insulated substrate provide surface heating. The laminarized velocity profile adjusts to a lower Reynolds number, and the structure responds to a lower loading. In the laminarization, the turbulent energy is dissipated by molecular transport by both viscous and conductivity mechanisms. Laminarization reduces spanwise vorticity because of the longitudinal cooling gradient of the sublayer profile. The results demonstrate that the curvature-induced mean pressure gradient enhances the receptivity of the flow to localized surface heating, a potentially viable mechanism to laminarize turbulent boundary layer flow; thus, the flow reduces the response of the flexible structure and the resultant sound radiation.

Maestrello, Lucio↗

Scattering and distortion of the unsteady motion on transversely sheared mean flows

It is shown that the pressure and velocity fluctuations of the unsteady motion on a transversely sheared mean flow can be expressed entirely in terms of the derivatives of two potential functions. One of these is a convected quantity that can be specified as a boundary condition and is related to a transverse component of the upstream velocity field. The other can be determined by solving an inhomogeneous wave equation whose source term is also a convected quantity that can be specified as a boundary condition in any given problem. The general theory is used to study the interaction of an unsteady flow with a semi-infinite plate embedded in a shear layer. The acoustic field produced by this interaction is calculated in the limits of low and high frequency. The results are compared with experimental one-third octave sound pressure level radiation patterns. The agreement is found to be excellent, especially in the low frequency range, where the mean-flow and convective effects are shown to have a strong influence on the directivity of the sound.

Goldstein, M. E.↗

Finite element analysis and optimization of composite structures

Linearly elastic fiber reinforced composite discs and laminates in plane stress with variable local orientation and concentration of one or two fiber fields embedded in the matrix material, are considered. The thicknesses and the domain of the discs or laminates are assumed to be given, together with prescribed boundary conditions and in-plane loading along the edge. The problem under study consists in determining throughout the structural domain the optimum orientations and concentrations of the fiber fields in such a way as to maximize the integral stiffness of the composite disc or laminate under the seven loading. Minimization of the integral stiffness can also be performed. The optimization is performed subject to a prescribed bound on the total cost or weight of the composite that for given unit cost factors or specific weights determines the amounts of fiber and matrix materials in the structure. Examples are presented.

Thomsen, Jan↗