SEARCH · Engineering Papers

Results for “Approximate 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.

SourceResearch connectionFrom yearThrough year

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

At least 1,045 records · Page 58

All sources NASA NTRS DOE OSTI

NASA NTRS1981

A comparative study of nonreflecting far-field boundary condition procedures for unsteady transonic flow computation

Various nonreflecting far-field boundary condition procedures are compared by implementing them in the computer code LTRAN2. This code solves the implicit finite-difference representation of the small-disturbance equations for transonic flows about airfoils. The first- and second-approximate nonreflecting conditions, as proposed by Engquist and Majda, are compared with the condition derived from the full-characteristic equation. The far-field boundary conditions and the description of the algorithm for implementing these conditions in LTRAN2 are discussed. Various cases are computed and compared with results from the older, more conventional procedures. One concludes that the full-characteristic equation produces the most effective results, thus allowing the far-field boundary to be located closer to the airfoil; this decreases the computer time required to obtain the solution because fewer mesh points are required.

Kwak, D.↗

NASA NTRS1990

Global configuration of a magnetic cloud

A magnetic cloud associated with a 2N flare on January 1, 1978 was observed by IMP-8, Helios A, Helios B, and Voyager 2. The variation of the magnetic field observed at each spacecraft is represented to good approximation by Lundquist's solution for a cylindrically symmetric force-free magnetic field with constant alpha. A least-squares fit of Lundquist's solution to the data from each spacecraft gives the local orientation of the axis of the magnetic cloud. The times of the estimated boundaries of the magnetic cloud at each spacecraft, together with the speeds of the boundaries and the spacecraft position, give the positions of the boundaries at a given time. From these results the magnetic cloud is determined to resemble a flux rope whose minor radius is approximately 0.15 AU at 1 AU, and whose radius of curvature at 1 AU is approximately 1/3 AU.

Burlaga, L. F.↗

NASA NTRS2018

Radiative Transfer in a Discrete Random Medium Adjacent to a Half-Space with a Rough Interface

For a macroscopically plane-parallel discrete random medium, the boundary conditions for the specific coherency dyadic at a rough interface are derived. The derivation is based on a modification of the Twersky approximation for a scattering system consisting of a group of particles and the rough surface, and reduces to the solution of the scattering problem for a rough surface illuminated by a plane electromagnetic wave propagating in a discrete random medium with non-scattering boundaries. In a matrix-form setting, the boundary conditions for the specific coherency dyadic imply the boundary conditions for specific intensity column vectors which in turn, yield the expressions for the reflection and transmission matrices. The derived expressions are shown to be identical to those obtained by applying a phenomenological approach based on a facet model to the solution of the scattering problem for a rough surface illuminated by a plane electromagnetic wave.

Adrian Doicu↗

NASA NTRS1988

Boundary identification for 2-D parabolic problems arising in thermal testing of materials

Problems on the identification of two-dimensional spatial domains arising in the detection and characterization of structural flaws in materials are considered. For a thermal diffusion system with external boundary input, observations of the temperature on the surface are used in an output least square approach. Parameter estimation techniques based on the method of mappings are discussed, and approximation schemes are developed based on a finite-element Galerkin approach. Theoretical convergence results for computational techniques are given, and the results are applied to the identification of two kinds of boundary shapes.

Banks, H. T.↗

NASA NTRS2004

Adjoint Formulation for an Embedded-Boundary Cartesian Method

Many problems in aerodynamic design can be characterized by smooth and convex objective functions. This motivates the use of gradient-based algorithms, particularly for problems with a large number of design variables, to efficiently determine optimal shapes and configurations that maximize aerodynamic performance. Accurate and efficient computation of the gradient, however, remains a challenging task. In optimization problems where the number of design variables dominates the number of objectives and flow- dependent constraints, the cost of gradient computations can be significantly reduced by the use of the adjoint method. The problem of aerodynamic optimization using the adjoint method has been analyzed and validated for both structured and unstructured grids. The method has been applied to design problems governed by the potential, Euler, and Navier-Stokes equations and can be subdivided into the continuous and discrete formulations. Giles and Pierce provide a detailed review of both approaches. Most implementations rely on grid-perturbation or mapping procedures during the gradient computation that explicitly couple changes in the surface shape to the volume grid. The solution of the adjoint equation is usually accomplished using the same scheme that solves the governing flow equations. Examples of such code reuse include multistage Runge-Kutta schemes coupled with multigrid, approximate-factorization, line-implicit Gauss-Seidel, and also preconditioned GMRES. The development of the adjoint method for aerodynamic optimization problems on Cartesian grids has been limited. In contrast to implementations on structured and unstructured grids, Cartesian grid methods decouple the surface discretization from the volume grid. This feature makes Cartesian methods well suited for the automated analysis of complex geometry problems, and consequently a promising approach to aerodynamic optimization. Melvin e t al. developed an adjoint formulation for the TRANAIR code, which is based on the full-potential equation with viscous corrections. More recently, Dadone and Grossman presented an adjoint formulation for the Euler equations. In both approaches, a boundary condition is introduced to approximate the effects of the evolving surface shape that results in accurate gradient computation.

Nemec, Marian↗

NASA NTRS2012

The First in situ Observation of Kelvin-Helmholtz Waves at High-Latitude Magnetopause during Strongly Dawnward Interplanetary Magnetic Field Conditions

We report the first in situ observation of high-latitude magnetopause (near the northern duskward cusp) Kelvin-Helmholtz waves (KHW) by Cluster on January 12, 2003, under strongly dawnward interplanetary magnetic field (IMF) conditions. The fluctuations unstable to Kelvin-Helmholtz instability (KHI) are found to propagate mostly tailward, i.e., along the direction almost 90 deg. to both the magnetosheath and geomagnetic fields, which lowers the threshold of the KHI. The magnetic configuration across the boundary layer near the northern duskward cusp region during dawnward IMF is similar to that in the low-latitude boundary layer under northward IMF, in that (1) both magnetosheath and magnetospheric fields across the local boundary layer constitute the lowest magnetic shear and (2) the tailward propagation of the KHW is perpendicular to both fields. Approximately 3-hour-long periods of the KHW during dawnward IMF are followed by the rapid expansion of the dayside magnetosphere associated with the passage of an IMF discontinuity that characterizes an abrupt change in IMF cone angle, Phi = acos (B(sub x) / absolute value of Beta), from approx. 90 to approx. 10. Cluster, which was on its outbound trajectory, continued observing the boundary waves at the northern evening-side magnetopause during sunward IMF conditions following the passage of the IMF discontinuity. By comparing the signatures of boundary fluctuations before and after the IMF discontinuity, we report that the frequencies of the most unstable KH modes increased after the discontinuity passed. This result demonstrates that differences in IMF orientations (especially in f) are associated with the properties of KHW at the high-latitude magnetopause due to variations in thickness of the boundary layer, and/or width of the KH-unstable band on the surface of the dayside magnetopause.

solar wind↗

NASA NTRS2006

Adjoint Sensitivity Computations for an Embedded-Boundary Cartesian Mesh Method and CAD Geometry

Cartesian-mesh methods are perhaps the most promising approach for addressing the issues of flow solution automation for aerodynamic design problems. In these methods, the discretization of the wetted surface is decoupled from that of the volume mesh. This not only enables fast and robust mesh generation for geometry of arbitrary complexity, but also facilitates access to geometry modeling and manipulation using parametric Computer-Aided Design (CAD) tools. Our goal is to combine the automation capabilities of Cartesian methods with an eficient computation of design sensitivities. We address this issue using the adjoint method, where the computational cost of the design sensitivities, or objective function gradients, is esseutially indepeudent of the number of design variables. In previous work, we presented an accurate and efficient algorithm for the solution of the adjoint Euler equations discretized on Cartesian meshes with embedded, cut-cell boundaries. Novel aspects of the algorithm included the computation of surface shape sensitivities for triangulations based on parametric-CAD models and the linearization of the coupling between the surface triangulation and the cut-cells. The objective of the present work is to extend our adjoint formulation to problems involving general shape changes. Central to this development is the computation of volume-mesh sensitivities to obtain a reliable approximation of the objective finction gradient. Motivated by the success of mesh-perturbation schemes commonly used in body-fitted unstructured formulations, we propose an approach based on a local linearization of a mesh-perturbation scheme similar to the spring analogy. This approach circumvents most of the difficulties that arise due to non-smooth changes in the cut-cell layer as the boundary shape evolves and provides a consistent approximation tot he exact gradient of the discretized abjective function. A detailed gradient accurace study is presented to verify our approach. Thereafter, we focus on a shape optimization problem for an Apollo-like reentry capsule. The optimization seeks to enhance the lift-to-drag ratio of the capsule by modifyjing the shape of its heat-shield in conjunction with a center-of-gravity (c.g.) offset. This multipoint and multi-objective optimization problem is used to demonstrate the overall effectiveness of the Cartesian adjoint method for addressing the issues of complex aerodynamic design. This abstract presents only a brief outline of the numerical method and results; full details will be given in the final paper.

Nemec, Marian↗

NASA NTRS2000

2D Quantum Mechanical Study of Nanoscale MOSFETs

With the onset of quantum confinement in the inversion layer in nanoscale MOSFETs, behavior of the resonant level inevitably determines all device characteristics. While most classical device simulators take quantization into account in some simplified manner, the important details of electrostatics are missing. Our work addresses this shortcoming and provides: (a) a framework to quantitatively explore device physics issues such as the source-drain and gate leakage currents, DIBL, and threshold voltage shift due to quantization, and b) a means of benchmarking quantum corrections to semiclassical models (such as density-gradient and quantum-corrected MEDICI). We have developed physical approximations and computer code capable of realistically simulating 2-D nanoscale transistors, using the non-equilibrium Green's function (NEGF) method. This is the most accurate full quantum model yet applied to 2-D device simulation. Open boundary conditions and oxide tunneling are treated on an equal footing. Electrons in the ellipsoids of the conduction band are treated within the anisotropic effective mass approximation. We present the results of our simulations of MIT 25, 50 and 90 nm "well-tempered" MOSFETs and compare them to those of classical and quantum corrected models. The important feature of quantum model is smaller slope of Id-Vg curve and consequently higher threshold voltage. Surprisingly, the self-consistent potential profile shows lower injection barrier in the channel in quantum case. These results are qualitatively consistent with ID Schroedinger-Poisson calculations. The effect of gate length on gate-oxide leakage and subthreshold current has been studied. The shorter gate length device has an order of magnitude smaller current at zero gate bias than the longer gate length device without a significant trade-off in on-current. This should be a device design consideration.

Svizhenko, Alexei↗

NASA NTRS2001

2-D Modeling of Nanoscale MOSFETs: Non-Equilibrium Green's Function Approach

We have developed physical approximations and computer code capable of realistically simulating 2-D nanoscale transistors, using the non-equilibrium Green's function (NEGF) method. This is the most accurate full quantum model yet applied to 2-D device simulation. Open boundary conditions and oxide tunneling are treated on an equal footing. Electrons in the ellipsoids of the conduction band are treated within the anisotropic effective mass approximation. Electron-electron interaction is treated within Hartree approximation by solving NEGF and Poisson equations self-consistently. For the calculations presented here, parallelization is performed by distributing the solution of NEGF equations to various processors, energy wise. We present simulation of the "benchmark" MIT 25nm and 90nm MOSFETs and compare our results to those from the drift-diffusion simulator and the quantum-corrected results available. In the 25nm MOSFET, the channel length is less than ten times the electron wavelength, and the electron scattering time is comparable to its transit time. Our main results are: (1) Simulated drain subthreshold current characteristics are shown, where the potential profiles are calculated self-consistently by the corresponding simulation methods. The current predicted by our quantum simulation has smaller subthreshold slope of the Vg dependence which results in higher threshold voltage. (2) When gate oxide thickness is less than 2 nm, gate oxide leakage is a primary factor which determines off-current of a MOSFET (3) Using our 2-D NEGF simulator, we found several ways to drastically decrease oxide leakage current without compromising drive current. (4) Quantum mechanically calculated electron density is much smaller than the background doping density in the poly silicon gate region near oxide interface. This creates an additional effective gate voltage. Different ways to. include this effect approximately will be discussed.

Svizhenko, Alexei↗

NASA NTRS1985

Vortex boundary-layer interactions

The interaction of a turbulent boundary layer on a flat plate with a strong artificially generated longitudinal vortex which may or may not actually enter the boundary layer is studied. The vortices are generated by a delta wing suspended ahead of the test plate, so that the configuration is approximately that of a close coupled carnard with zero main-wing sweep and an invisible body. All necessary configuration and parametric checks are completed, and data acquisition and analysis on the first configuration chosen for detailed study, in which the vortex starts to merge with the boundary layer a short distance downstream of the leading edge of the test plate, are nearly complete.

Bradshaw, P.↗

NASA NTRS1990

A computational efficient modelling of laminar separation bubbles

In predicting the aerodynamic characteristics of airfoils operating at low Reynolds numbers, it is often important to account for the effects of laminar (transitional) separation bubbles. Previous approaches to the modelling of this viscous phenomenon range from fast but sometimes unreliable empirical correlations for the length of the bubble and the associated increase in momentum thickness, to more accurate but significantly slower displacement-thickness iteration methods employing inverse boundary-layer formulations in the separated regions. Since the penalty in computational time associated with the more general methods is unacceptable for airfoil design applications, use of an accurate yet computationally efficient model is highly desirable. To this end, a semi-empirical bubble model was developed and incorporated into the Eppler and Somers airfoil design and analysis program. The generality and the efficiency was achieved by successfully approximating the local viscous/inviscid interaction, the transition location, and the turbulent reattachment process within the framework of an integral boundary-layer method. Comparisons of the predicted aerodynamic characteristics with experimental measurements for several airfoils show excellent and consistent agreement for Reynolds numbers from 2,000,000 down to 100,000.

Dini, Paolo↗

NASA NTRS1994

Ultraviolet illuminated molecular cloud boundaries: Extended (C II) 158 micrometer emission toward L1630

The construction of a large scale map of the 158 micrometer C+ line toward the L1630/Orion B molecular cloud, covering an approximately 35' by 45' area which includes the NGC 2024 H II region, zeta-Ori, the reflection nebula NGC 2023, and the Horsehead nebula, is reported. Emission in the C II line is very widespread. The line was detected at levels in excess of a few 0.0001 erg/sq cm/s/sr over almost the entire mapped region. Extended emission associated with the NGC 2024 H II region and its envelope accounts for more than half of the C II flux. Over this approximately 1.5 by 2.5 pc region, the amount of gas phase carbon in the form of C+ is comparable to the amount of carbon in CO. This result, together with the C II distribution implies that C II emission arises on the surface of clumps throughout the cloud rather than in a single layer at the H II region boundary. Away from the H II region, most of the C II emission comes from the western edge of the L1630 cloud and probably results from excitation by external OB stars. The overall extent of the C II emission is comparable to that of millimeter molecular lines but the distributions are different in detail. The difference in C II and molecular line distributions, in particular, the larger extent of the C II emission west of NGC 2024 implies large variations in the ratio of the C II and J = 1 towards 0 intensities. Models of photon dominated regions can explain the relation between C II and CO intensities only if the cloud edges and cloud interior are considered separately. A method for using C II and radio continuum emission to characterize the relationship between OB stars and photon dominated regions is proposed.

Jaffe, D. T.↗

NASA NTRS1976

Interplanetary boundary layers at 1 AU

The structure and nature of discontinuities in the interplanetary magnetic field at 1 AU in the period March 18, 1971 to April 9, 1971, is determined by using high-resolution magnetic field measurements from Explorer 34. The discontinuities that were selected for this analysis occurred under a variety of interplanetary conditions at an average rate of 0.5/hr. This set does not include all discontinuities that were present, but the sample is large and it is probably representative. Both tangential and rotational discontinuities were identified, the ratio of TD's to RD's being approximately 3 to 1. Tangential discontinuities were observed every day, even among Alfvenic fluctuations. The structure of most of the boundary layers was simple and ordered, i.e., the magnetic field usually changed smoothly and monotonically from one side of the boundary layer to the other.

Burlaga, L. F.↗

NASA NTRS1981

Dynamic behavior of an unsteady trubulent boundary layer

Experiments on an unsteady turbulent boundary layer are reported in which the upstream portion of the flow is steady (in the mean) and in the downstream region, the boundary layer sees a linearly decreasing free stream velocity. This velocity gradient oscillates in time, at frequencies ranging from zero to approximately the bursting frequency. For the small amplitude, the mean velocity and mean turbulence intensity profiles are unaffected by the oscillations. The amplitude of the periodic velocity component, although as much as 70% greater than that in the free stream for very low frequencies, becomes equal to that in the free stream at higher frequencies. At high frequencies, both the boundary layer thickness and the Reynolds stress distribution across the boundary layer become frozen. The behavior at higher amplitude is quite similar. At sufficiently high frequencies, the boundary layer thickness remains frozen at the mean value over the oscillation cycle, even though flow reverses near the wall during a part of the cycle.

Parikh, P. G.↗

NASA NTRS1994

ROSAT observations of cataclysmic variables: A search for the boundary layer emission

Many cataclysmic variables (CVs) and low-mass X-ray binaries show He II line emission at wavelength 1640 and wavelength 4686. As recombination lines of an ion requiring 55 eV for ionization, these lines cannot be attributed to gravitational heating in the disk but could be due to reprocessed 55-280 eV X-rays which may be produced in a boundary layer at the disk/star surface. We undertook a search for this soft component using the Positon Sensitve Proportional Counter (PSPC) on ROSAT. The photon fluxes we observe are not sufficient to produce the He II wavelength 1640 fluxes observed by IUE. If we assume that the X-ray luminosity is due to emission from a boundary layer, then the observed quiescent luminosity is one-quarter of the total accretion luminosity as determined by optical and ultraviolet observations rather than the one-half that is predicted by some theories. An additional very soft component (e.g. an approximately equal 10-20 eV blackbody) could provide sufficient photons in the 0.055-0.100 keV band to produce the observed He II fluxes and contribute the additional luminosity from the boundary layer while remaining hidden to the PSPC and IUE. Alternatively, the high He II line fluxes could be evidence for high-temperature collisionally ionized material. We also find evidence for line emission at approximately equal 0.9 keV consistent with L-shell emission from highly ionized iron.

Vrtilek, S. D.↗

NASA NTRS1997

Analytical Studies of Boundary Layer Generated Aircraft Interior Noise

An analysis is made of the "interior noise" produced by high, subsonic turbulent flow over a thin elastic plate partitioned into "panels" by straight edges transverse to the mean flow direction. This configuration models a section of an aircraft fuselage that may be regarded as locally flat. The analytical problem can be solved in closed form to represent the acoustic radiation in terms of prescribed turbulent boundary layer pressure fluctuations. Two cases are considered: (i) the production of sound at an isolated panel edge (i.e., in the approximation in which the correlation between sound and vibrations generated at neighboring edges is neglected), and (ii) the sound generated by a periodic arrangement of identical panels. The latter problem is amenable to exact analytical treatment provided the panel edge conditions are the same for all panels. Detailed predictions of the interior noise depend on a knowledge of the turbulent boundary layer wall pressure spectrum, and are given here in terms of an empirical spectrum proposed by Laganelli and Wolfe. It is expected that these analytical representations of the sound generated by simplified models of fluid-structure interactions can used to validate more general numerical schemes.

Howe, M. S.↗

NASA NTRS1981

Dynamic behavior of an unsteady turbulent boundary layer

Experiments on an unsteady turbulent boundary layer are reported in which the upstream portion of the flow is steady (in the mean) and in the downstream region, the boundary layer sees a linearly decreasing free stream velocity. This velocity gradient oscillates in time, at frequencies ranging from zero to approximately the bursting frequency. For the small amplitude, the mean velocity and mean turbulence intensity profiles are unaffected by the oscillations. The amplitude of the periodic velocity component, although as much as 70 percent greater than that in the free stream for very low frequencies, becomes equal to that in the free stream at higher frequencies. At high frequencies, both the boundary layer thickness and the Reynolds stress distribution across the boundary layer become frozen. The behavior at higher amplitude is quite similar. At sufficiently high frequencies, the boundary layer thickness remains frozen at the mean value over the oscillation cycle, even though flow reverses near the wall during a part of the cycle. Previously announced in STAR as N81-28392

Parikh, P. G.↗

NASA NTRS1977

Introductory Remarks

Advanced transonic configurations, such as the supercritical wing, are inherently more sensitive to Reynolds number than earlier configurations because the pressure recovery gradients imposed on the boundary layer are generally steeper. The results of two-dimensional supercritical airfoil investigations and theoretical calculations show this effect. A technique for approximately simulating full-scale Reynolds number characteristics at present wind-tunnel Reynolds numbers for near-cruise conditions is described. The transition strip, which in the past has been located near the leading edge of the wing, is rearward so that the relative displacement thickness of the boundary layer at the trailing edge of the wing is the same as might be expected on a full-scale configuration with the transition near the leading edge. Two-dimensional wind-tunnel results indicate that the technique provides a very good simulation of airfoil characteristics at full-scale Reynolds number. The need for sorting the effects of Reynolds number and aeroelasticity, which can be done in the NTF, is also addressed.

Whitcomb, R.↗