Search NASA⌕ Search

SEARCH · Search NASA

Results for “boundary value problems”

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,369 records · Page 76

Elucidating the Onset of Plasticity in Sliding Contacts Using Differential Computational Orientation Tomography

Abstract Depending on the mechanical and thermal energy introduced to a dry sliding interface, the near-surface regions of the mated bodies may undergo plastic deformation. In this work, we use large-scale molecular dynamics simulations to generate “differential computational orientation tomographs” (dCOT) and thus highlight changes to the microstructure near tribological FCC alloy surfaces, allowing us to detect subtle differences in lattice orientation and small distances in grain boundary migration. The analysis approach compares computationally generated orientation tomographs with their undeformed counterparts via a simple image analysis filter. We use our visualization method to discuss the acting microstructural mechanisms in a load- and time-resolved fashion, focusing on sliding conditions that lead to twinning, partial lattice rotation, and grain boundary-dominated processes. Extracting and laterally averaging the color saturation value of the generated tomographs allows us to produce quantitative time- and depth-resolved maps that give a good overview of the progress and severity of near-surface deformation. Corresponding maps of the lateral standard deviation in the color saturation show evidence of homogenization processes occurring in the tribologically loaded microstructure, frequently leading to the formation of a well-defined separation between deformed and undeformed regions. When integrated into a computational materials engineering framework, our approach could help optimize material design for tribological and other deformation problems. Graphic Abstract .

36 MATERIALS SCIENCE↗

Numerical matching of the sheath and presheath solutions for a spherical probe in radial-motion theory

The theory of positive-ion collection by a probe immersed in a low-pressure plasma was reviewed and extended by Allen, Boyd, and Reynolds [Proc. Phys. Soc. 70, 297 (1957)]. For a given value of the ion current, the boundary values of the matched “nonneutral” or “sheath” solution V{sup ~}{sub nn}{sup (m)}(r; r{sub m}) were obtained from the “quasineutral” or “presheath” solution V{sup ~}{sub qn}(r) by choosing the small potential and electric-field values corresponding to some large “matching radius” r{sub m}. Here, a straightforward but efficient numerical method is presented for systematically determining an optimal value of the matching radius at which the presheath and sheath solutions are joined to yield the “matched” potential profile. Some suitable initial matching radius r{sub m1} is chosen and the related potential and electric-field values of the quasineutral solution are calculated. Using these as boundary conditions, Poisson's equation is integrated to yield the matched nonneutral solution including the corresponding potential at the probe surface. This procedure is repeated for increasing values r{sub m2}, r{sub m3},…. until the resulting potential at the probe surface becomes practically constant. The corresponding value of r{sub m} is taken as the “optimal” matching radius r{sub mo} at which the presheath and sheath solutions are ultimately joined to yield the “optimal” matched potential profile in the entire plasma-probe transition region. It is also shown that the Bohm criterion is inapplicable in the present problem.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

On the assumption of vanishing temperature fluctuations at the wall for heat transfer modeling

Boundary conditions for fluctuating wall temperature are required for near-wall heat transfer modeling. However, their correct specifications for arbitrary thermal boundary conditions are not clear. The conventional approach is to assume zero fluctuating wall temperature or zero gradient for the temperature variance at the wall. These are idealized specifications and the latter condition could lead to an ill posed problem for fully-developed pipe and channel flows. In this paper, the validity and extent of the zero fluctuating wall temperature condition for heat transfer calculations is examined. The approach taken is to assume a Taylor expansion in the wall normal coordinate for the fluctuating temperature that is general enough to account for both zero and non-zero value at the wall. Turbulent conductivity is calculated from the temperature variance and its dissipation rate. Heat transfer calculations assuming both zero and non-zero fluctuating wall temperature reveal that the zero fluctuating wall temperature assumption is in general valid. The effects of non-zero fluctuating wall temperature are limited only to a very small region near the wall.

Sommer, T. P.↗

Probabilistic Structural Analysis Program

NASA/NESSUS 6.2c is a general-purpose, probabilistic analysis program that computes probability of failure and probabilistic sensitivity measures of engineered systems. Because NASA/NESSUS uses highly computationally efficient and accurate analysis techniques, probabilistic solutions can be obtained even for extremely large and complex models. Once the probabilistic response is quantified, the results can be used to support risk-informed decisions regarding reliability for safety-critical and one-of-a-kind systems, as well as for maintaining a level of quality while reducing manufacturing costs for larger-quantity products. NASA/NESSUS has been successfully applied to a diverse range of problems in aerospace, gas turbine engines, biomechanics, pipelines, defense, weaponry, and infrastructure. This program combines state-of-the-art probabilistic algorithms with general-purpose structural analysis and lifting methods to compute the probabilistic response and reliability of engineered structures. Uncertainties in load, material properties, geometry, boundary conditions, and initial conditions can be simulated. The structural analysis methods include non-linear finite-element methods, heat-transfer analysis, polymer/ceramic matrix composite analysis, monolithic (conventional metallic) materials life-prediction methodologies, boundary element methods, and user-written subroutines. Several probabilistic algorithms are available such as the advanced mean value method and the adaptive importance sampling method. NASA/NESSUS 6.2c is structured in a modular format with 15 elements.

Pai, Shantaram S.↗

Effect of radiation and magnetohydrodynamic free convection boundary layer flow on a solid sphere with Newtonian heating in a micropolar fluid

In this paper, the effect of radiation on magnetohydrodynamic free convection boundary layer flow on a solid sphere with Newtonian heating in a micropolar fluid, in which the heat transfer from the surface is proportional to the local surface temperature, is considered. The transformed boundary layer equations in the form of nonlinear partial differential equations are solved numerically using an implicit finite difference scheme known as the Keller-box method. Numerical solutions are obtained for the local wall temperature and the local skin friction coefficient, as well as the velocity, angular velocity and temperature profiles. The features of the flow and heat transfer characteristics for various values of the Prandtl number Pr, micropolar parameter K, magnetic parameter M, radiation parameter N{sub R}, the conjugate parameter γ and the coordinate running along the surface of the sphere, x are analyzed and discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Systematic over-estimation of lattice thermal conductivity in materials with electrically-resistive grain boundaries

Reducing the thermal conductivity κ of a material via nano-structuring to create small grain sizes is one of the most common strategies to improve thermoelectric materials. In such polycrystalline materials heat carrying phonons are scattered at the grain boundaries, which directly improves the thermoelectric quality factor and ultimately the figure-of-merit zT. In some cases, however, such as in Mg 3 Sb 2 , SnSe, and Mg 2 Si an opposite trend is found where higher lattice thermal conductivity reported in small grain polycrystalline material than in large grain or single crystal materials. This unphysical result indicates a problem with the conventional use of the Wiedemann–Franz law. Here, we trace this problematic finding to the electrical resistance at the grain boundaries, which leads to an overestimation of the phonon or lattice contribution of the thermal conductivity κ L . In materials with significant grain boundary electrical resistance, the estimated electronic contribution to the thermal conductivity LσT is low because the measured electrical conductivity σ is low. However within the grain electrons may still be transporting more heat than the total conductivity suggests, leading to an overestimation of κ L if the conventional κ L = κ – LσT is used with the measured values of κ and σ. Here, the overestimation of κL in small-grain samples is shown to be pervasive across a broad range of thermoelectric materials, including Mg 3 Sb 2 , Mg 2 Si, PbTe, PbSe, SnSe, (Hf,Zr)CoSb, CoSb 3 , and Bi 2 Te 3 alloys, and a correction is necessary to properly understand and predict their charge and heat transport.

36 MATERIALS SCIENCE↗

Benchmark study of the two-dimensional Hubbard model with auxiliary-field quantum Monte Carlo method

Ground-state properties of the Hubbard model on a two-dimensional square lattice are studied by the auxiliary-field quantum Monte Carlo method. Accurate results for energy, double occupancy, effective hopping, magnetization, and momentum distribution are calculated for interaction strengths of U/t from 2 to 8, for a range of densities including half-filling and n = 0.3,0.5,0.6, 0.75, and 0.875. At half-filling, the results are numerically exact. Away from half-filling, the constrained path Monte Carlo method is employed to control the sign problem. Our results are obtained with several advances in the computational algorithm, which are described in this paper in detail. We discuss the advantages of generalized Hartree-Fock trial wave functions and its connection to pairing wave functions, as well as the interplay with different forms of Hubbard-Stratonovich decompositions. We study the use of different twist angle sets when applying the twist averaged boundary conditions. We propose the use of quasirandom sequences, which improves the convergence to the thermodynamic limit over pseudorandom and other sequences. With it and a careful finite size scaling analysis, we are able to obtain accurate values of ground-state properties in the thermodynamic limit. Detailed results for finite-sized systems up to 16 × 16 are also provided for benchmark purposes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Effect of velocity gradients on measurements of turbulent shear stress

Experiments were carried out to evaluate the effects of normal velocity gradients on hot wire measurements in a subsonic boundary layer of the same size as the flow investigated by Johnson and Rose (1975). Both hot wire and film anemometers were used to measure the turbulent properties of the boundary layer. A special X-wire probe with one wire vertical and the other at an angle of about 40 deg to the flow was used to demonstrate the gradient effects. The results indicate that major errors are encountered when mean and turbulent velocity gradients exist along the length of hot wire sensors, the problem being more pronounced at high speeds. Although the split film sensor results show a significant improvement over the X-wire sensor, further reduction in the space resolution of sensors by roughly an order of magnitude would appear to be necessary to reduce the error to acceptable values near the wall.

Sandborn, V. A.↗

H2O-H2SO4 system in Venus' clouds and OCS, CO, and H2SO4 profiles in Venus' troposphere

A coupled problem of diffusion and condensation is solved for the H2SO4-H2O system in Venus' cloud layer. The position of the lower cloud boundary and profiles of the H2O and H2SO4 vapor mixing ratios and of the H2O/H2SO4 ratio of sulfuric acid aerosol and its flux are calculated as functions of the column photochemical production rate of sulfuric acid, Phi(sub H2SO4). Variations of the lower cloud boundary are considered. Our basic model, which is constrained to yield f(sub H2O) (30 km) = 30 ppm (Pollack et al. 1993), predicts the position of the lower cloud boundary at 48.4 km coinciding with the mean Pioneer Venus value, the peak H2SO4 mixing ratio of 5.4 ppm, and the H2SO4 production rate Phi(sub H2SO4) = 2.2 x 10(exp 12)/sq cm/s. The sulfur to sulfuric acid mass flux ratio in the clouds is 1:27 in this model, and the mass loading ratio may be larger than this value if sulfur particles are smaller than those of sulfuric acid. The model suggests that the extinction coefficient of sulfuric acid particles with radius 3.7 micrometers (mode 3) is equal to 0.3/km in the middle cloud layer. The downward flux of CO is equal to 1.7 x 10(exp 12)/sq cm/s in this model. Our second model, which is constrained to yield f(sub H2O) = 10 ppm at the lower cloud boundary, close to the value measured by the Magellan radiooccultations, predicts the position of this boundary to be at 46.5 km, which agrees with the Magellan data; f(sub H2O) (30 km) = 90 ppm, close to the data of Moroz et al. (1983) at this altitude; Phi(sub H2SO4) = 6.4 x 10(exp 12)/sq cm/s; and Phi(sub co) = 4.2 x 10(exp 12)/sq cm/s. The S/H2SO4 flux mass ratio is 1:18, and the extinction coefficient of the mode 3 sulfuric acid particles is equal to 0.9 km in the middle cloud layer. A strong gradient of the H2SO4 vapor mixing ratio near the bottom of the cloud layer drives a large upward flux of H2SO4, which condenses and forms the excessive downward flux of liquid sulfuric acid, which is larger by a factor of 4-7 than the flux in the middle cloud layer. This is the mechanism of formation of the lower cloud layer. Variations of the lower cloud layer are discussed. Our modeling of the OCS and CO profiles in the lower atmosphere measured by Pollack et al. (1993) provides a reasonable explanation of these data and shows that the rate coefficient of the reaction SO3 + CO yields CO2 + SO2 is equal to 10(exp -11) exp(-(13,100 +/- 1000)/T) cc/s. The main channel of the reaction between SO3 and OCS is CO2 + (SO)2, and its rate coefficient is equal to 10(exp -11) exp(-(8900 +/- 500)/T) cc/s. In the conditions of Venus' lower atmosphere, (SO)2 is removed by the reaction (SO)2 + OCS yields CO + S2 + SO2. The model predicts an OCS mixing ratio of 28 ppm near the surface.

Krasnopolsky, V. A.↗

Bidirectional reflectance of oceanic waters: A comparison of modeled and measured upward radiance fields

The bidirectional reflectance of oceanic waters is conveniently described in a normalized way by forming the ratio of the upwelling irradiance E(sub u) to any upwelling radiance L(sub u)(theta prime, phi). This ratio, Q (theta prime, theta(sub 0), (phi(sub 0) - phi), where theta prime, phi are the nadir and azimuth angles for the upward radiance and theta(sub 0), phi(sub 0) are the zenith and azimuth angles of the Sun, has been determined from measurements at sea and computed via Monte Carlo simulations using the inherent optical properties measured in the field and appropriate boundary conditions (clear sky, no wind, varying Sun angle). Experimental ad computed Q values are in excellent agreement. This successful comparison confirms the importance of the bidirectional character of ocean reflectance, already pointed out from a purely numerical approach without field validation, and corroborates the extended range of the Q variations. The later point is of importance when interpreting the marine signals detected by an ocean color satellite-borne sensor. The validation is extended by considering the historical data for the radiance distributions in Lake Pend Oreille determined at various depths. The closure issue in ocean optics is examined by solving the direct problem of radiative transfer and through a model-data comparison in terms of radiance field.

Morel, Andre↗

The interaction between giant gaseous protoplanets and the primitive solar nebula

The manner in which a giant gaseous protoplanet becomes embedded in the primitive solar nebula determines surface boundary conditions which must be used in studying the evolution of such objects. On the one hand, if the system resembles a contact binary system, then the envelope of the protoplanet should approach the entropy of the surrounding nebula. On the other hand angular momentum transfer by resonance and tidal effects between the nebula and the protoplanet may cause the nebula to exhibit a zone of avoidance near the protoplanet, thus inhibiting exchange of material. This problem has been studied with a computer program developed by D. N. C. Lin which simulates disk hydrodynamics by particle motions with dissipation. These studies suggest that for expected values of the protoplanet/protosun mass ratios, significant inhibition of mass exchange is likely, so that it is a reasonable next step to undertake protoplanet evolution studies with the imposition of minimum protoplanet surface temperatures.

Cameron, A. G. W.↗

Boundary-layer receptivity due to a wall suction and control of Tollmien-Schlichting waves

A numerical study of the generation of Tollmien-Schlichting (T-S) waves due to the interaction between a small free-stream disturbance and a small localized suction slot on an otherwise flat surface was carried out using finite difference methods. The nonlinear steady flow is of the viscous-inviscid interactive type while the unsteady disturbed flow is assumed to be governed by the Navier-Stokes equations linearized about this flow. Numerical solutions illustrate the growth or decay of T-S waves generated by the interaction between the free-stream disturbance and the suction slot, depending on the value of the scaled Strouhal number. An important result of this receptivity problem is the numerical determination of the amplitude of the T-S waves and the demonstration of the possible active control of the growth of T-S waves.

Bodonyi, R. J.↗

Boundary-layer receptivity due to a wall suction and control of Tollmien-Schlichting waves

A numerical study of the generation of Tollmien-Schlichting (T-S) waves due to the interaction between a small free-stream disturbance and a small localized suction slot on an otherwise flat surface was carried out using finite difference methods. The nonlinear steady flow is of the viscous-inviscid interactive type while the unsteady disturbed flow is assumed to be governed by the Navier-Stokes equations linearized about this flow. Numerical solutions illustrate the growth or decay of T-S waves generated by the interaction between the free-stream disturbance and the suction slot, depending on the value of the scaled Strouhal number. An important result of this receptivity problem is the numerical determination of the amplitude of the T-S waves and the demonstration of the possible active control of the growth of T-S waves.

Bodonyi, R. J.↗

An analysis of the vertical structure equation for arbitrary thermal profiles

The vertical structure equation is a singular Sturm-Liouville problem whose eigenfunctions describe the vertical dependence of the normal modes of the primitive equations linearized about a given thermal profile. The eigenvalues give the equivalent depths of the modes. The spectrum of the vertical structure equation and the appropriateness of various upper boundary conditions, both for arbitrary thermal profiles were studied. The results depend critically upon whether or not the thermal profile is such that the basic state atmosphere is bounded. In the case of a bounded atmosphere it is shown that the spectrum is always totally discrete, regardless of details of the thermal profile. For the barotropic equivalent depth, which corresponds to the lowest eigen value, upper and lower bounds which depend only on the surface temperature and the atmosphere height were obtained. All eigenfunctions are bounded, but always have unbounded first derivatives. It was proved that the commonly invoked upper boundary condition that vertical velocity must vanish as pressure tends to zero, as well as a number of alternative conditions, is well posed. It was concluded that the vertical structure equation always has a totally discrete spectrum under the assumptions implicit in the primitive equations.

Cohn, Stephen E.↗

An analysis of the vertical structure equation for arbitrary thermal profiles

The vertical structure equation is a singular Sturm-Liouville problem whose eigenfunctions describe the vertical dependence of the normal modes of the primitive equations linearized about a given thermal profile. The eigenvalues give the equivalent depths of the modes. The spectrum of the vertical structure equation and the appropriateness of various upper boundary conditions, both for arbitrary thermal profiles were studied. The results depend critically upon whether or not the thermal profile is such that the basic state atmosphere is bounded. In the case of a bounded atmosphere it is shown that the spectrum is always totally discrete, regardless of details of the thermal profile. For the barotropic equivalent depth, which corresponds to the lowest eigen value, upper and lower bounds which depend only on the surface temperature and the atmosphere height were obtained. All eigenfunctions are bounded, but always have unbounded first derivatives. It was proved that the commonly invoked upper boundary condition that vertical velocity must vanish as pressure tends to zero, as well as a number of alternative conditions, is well posed. It was concluded that the vertical structure equation always has a totally discrete spectrum under the assumptions implicit in the primitive equations.

Cohn, Stephen E.↗

Calculated blade response at high tip-speed ratios

This paper presents the initial results of a digital computer study of helicopter rotor-blade-motion stability for a broad range of forward speeds encompassing proposed conventional and compound helicopter designs. The analysis treats the general case of nonlinear, coupled flapping and lagging motions of hinged rigid blades. The results of the study to date indicate that blade-motion stability boundaries cannot be specified for a given blade design in terms of a fixed value of rotor tip-speed ratio. The stability boundaries can shift significantly within the desired operational speed range, with the stability limits depending upon the rotor loading, the mag­nitude of the disturbance encountered by the blade, and the blade position at the instant the external disturbance is encountered. The sample methods suggested for improving blade-motion stability include increasing the effec­tive hinge spring restraint by incorporating pitch-flap coupling or by using hingeless or teetering rotor systems. Reducing rotor loading also has a beneficial effect on the blade-motion stability. The results indi­cate that methods used to deal with the blade-motion stability problem can be expected to diminish the overall vibration levels of the helicopter rotor system across the entire design speed range.

TIP SPEED↗

A New Higher-Order Composite Theory for Analysis and Design of High Speed Tilt-Rotor Blades

A higher-order theory is developed to model composite box beams with arbitrary wall thicknesses. The theory, based on a refined displacement field, represents a three-dimensional model which approximates the elasticity solution. Therefore, the cross-sectional properties are not reduced to one-dimensional beam parameters. Both inplane and out-of-plane warping are automatically included in the formulation. The model accurately captures the transverse shear stresses through the thickness of each wall while satisfying all stress-free boundary conditions. Several numerical results are presented to validate the present theory. The developed theory is then used to model the load carrying member of a tilt-rotor blade which has thick-walled sections. The composite structural analysis is coupled with an aerodynamic analysis to compute the aeroelastic stability of the blade. Finally, a multidisciplinary optimization procedure is developed to improve the aerodynamic, structural and aeroelastic performance of the tilt-rotor aircraft. The Kreisselmeier-Steinhauser function is used to formulate the multiobjective function problem and a hybrid approximate analysis is used to reduce the computational effort. The optimum results are compared with the baseline values and show significant improvements in the overall performance of the tilt-rotor blade.

McCarthy, Thomas Robert↗

Finite analytic numerical solution of axisymmetric Navier-Stokes and energy equations

Convective heat transfer for steady state laminar flow in axisymmetric coordinates is considered. Numerical solutions for flow pattern and temperature distribution are obtained by the finite analytic numerical method applied to the Navier-Stokes equations expressed in terms of vorticity and stream function, and the energy equation. The finite analytic numerical method differs from other numerical methods in that it utilizes a local analytic solution in an element of the problem to construct the total numerical solution. Finite analytic solutions of vorticity, stream function, temperature and heat transfer coefficients for flow with Reynolds number of 5.0, 100.0, 1000.0, and 2000.0, and Prandtl number of 0.1, 1.0, and 10.0 with uniform grid sizes are reported for axisymmetric pipe with sudden expansion and contraction. The wall temperature is considered to be isothermal and differs from the inlet temperature. It is shown that the finite analytic solution is stable, converges rapidly, and simulates the convection of fluid flow accurately since the local analytic solution is capable of simulating automatically the influence of skewed convection through the element boundary on the interior nodal values thereby minimizing the false numerical diffusion.

Chen, C.-J.↗