Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Boundary value problem”

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,027 records · Page 57

On the closedness and geometry of tensor network state sets

Tensor network states (TNS) are a powerful approach for the study of strongly correlated quantum matter. The curse of dimensionality is addressed by parametrizing the many-body state in terms of a network of partially contracted tensors. These tensors form a substantially reduced set of effective degrees of freedom. In practical algorithms, functionals like energy expectation values or overlaps are optimized over certain sets of TNS. Concerning algorithmic stability, it is important whether the considered sets are closed because, otherwise, the algorithms may approach a boundary point that is outside the TNS set and tensor elements diverge. Here we discuss the closedness and geometries of TNS sets, and we propose regularizations for optimization problems on non-closed TNS sets. We show that sets of matrix product states (MPS) with open boundary conditions, tree tensor network states, and the multiscale entanglement renormalization ansatz are always closed, whereas sets of translation-invariant MPS with periodic boundary conditions (PBC), heterogeneous MPS with PBC, and projected entangled pair states are generally not closed. The latter is done using explicit examples like the W state, states that we call two-domain states, and fine-grained versions thereof.

97 MATHEMATICS AND COMPUTING↗

The initial-value problem for viscous channel flows

Plane viscous channel flows are perturbed and the ensuing initial-value problems are investigated in detail. Unlike traditional methods where traveling wave normal modes are assumed for solution, this works offers a means whereby completely arbitrary initial input can be specified without having to resort to eigenfunction expansions. The full temporal behavior, including both early time transients and the long time asymptotics, can be determined for any initial disturbance. Effects of three-dimensionality can be assessed. The bases for the analysis are: (a) linearization of the governing equations; (b) Fourier decomposition in the spanwise and streamwise directions of the flow; and (c) direct numerical integration of the resulting partial differential equations. All of the stability data that are known for such flows can be reproduced. Also, the optimal initial condition can be determined in a straight forward manner and such optimal conditions clearly reflect transient growth data that is easily determined by a rational choice of a basis for the initial conditions. Although there can be significant transient growth for subcritical values of the Reynolds number using this approach it does not appear possible that arbitrary initial conditions will lead to the exceptionally large transient amplitudes that have been determined by optimization of normal modes. The approach is general and can be applied to other classes of problems where only a finite discrete spectrum exists, such as the boundary layer for example.

Criminale, W. O.↗

The causes for geographical variations in OS187/OS186 at the Cretaceous-Tertiary boundary

Researchers at Yale has approached the problem of the osmium isotopic composition of marine deposits formed in contact with both oxidized and reduced bottom waters. The measured (187) Os/(186) Os ratios of modern bulk sediment can be explained using mixing equations involving continental detrital, volcaniclastic, cosmogenic and hydrogeneous components. These studies show that sediments deposited under reducing marine conditions contain a hydrogenous component which is enriched in Re and has a radiogenic (187) Os/(186) Os ratio. The presence of such a hydrogenous component in the marine fish clay at Stevns Klint can account for the elevation of its (187) Os/(186) Os ration above the expected meteoritic value. Mass balance considerations require the Re/Os ratio of the phase precipitated from the terminal Cretaceous sea at Stevns Klint to have been about one tenth the value observed in contemporary deposits in the Black Sea, assuming Re has not been lost (or Os gained) subsequent to precipitation. In continental sections, the elevation of the (187) Os/(186) Os ratio in boundary layers may be due to precipitation from continental waters of crustally-derived radiogenic osmium either contemporaneous with the meteoritic (or mantle) osmium deposition or later during diagenesis.

Turekian, K. K.↗

Bayesian inference in geomagnetism

The inverse problem in empirical geomagnetic modeling is investigated, with critical examination of recently published studies. Particular attention is given to the use of Bayesian inference (BI) to select the damping parameter lambda in the uniqueness portion of the inverse problem. The mathematical bases of BI and stochastic inversion are explored, with consideration of bound-softening problems and resolution in linear Gaussian BI. The problem of estimating the radial magnetic field B(r) at the earth core-mantle boundary from surface and satellite measurements is then analyzed in detail, with specific attention to the selection of lambda in the studies of Gubbins (1983) and Gubbins and Bloxham (1985). It is argued that the selection method is inappropriate and leads to lambda values much larger than those that would result if a reasonable bound on the heat flow at the CMB were assumed.

Backus, George E.↗

N-point correlation functions in the CfA and SSRS redshift distribution of galaxies

Using counts in cells, we estimate the volume-average N-point galaxy correlation functions for N = 2, 3, and 4, in redshift samples of the CfA and SSRS catalogs. Volume-limited samples of different sizes are used to study the uncertainties at different scales, the shot noise, and the problem with the boundaries. The hierarchical constants S3 and S4 agree well in all samples in CfA and SSRS, with average S3 = 194 +/- 0.07 and S4 = 4.56 +/- 0.53. We compare these results with estimates obtained from angular catalogs and recent analysis over IRAS samples. The amplitudes SJ seem larger in real space than in redshift space, although the values from the angular analysis correspond to smaller scales, where we might expect larger nonperturbative effects. It is also found that S3 and S4 are smaller for IRAS than for optical galaxies. This, together with the fact that IRAS galaxies have smaller amplitude for the above correlation functions, indicates that the density fluctuations of IRAS galaxies cannot be simply proportional to the density fluctuations of optical galaxies, i.e., biasing has to be nonlinear between them.

Gaztanaga, Enrique↗

An Airfoil Spanning an Open Jet

Proceeding from the fundamental problem on the mutual relation of a wing and free boundaries the distribution of the circulation is determined for an airfoil spanning an open jet of rectangular section at different aspect ratios, and then for an open jet of circular section. The solution is obtained by means of a Fourier series and computations have been performed for different values of the variables. The second part describes the experiments performed for the purpose of proving the theory. The results confirm the theory. In conclusion it defines the induced drag of a wing extending across an open jet and compares it with the drag of a monoplane having a span equal to the jet width at equal total lift.

Stuper, J↗

LSENS: A General Chemical Kinetics and Sensitivity Analysis Code for homogeneous gas-phase reactions. Part 1: Theory and numerical solution procedures

LSENS, the Lewis General Chemical Kinetics and Sensitivity Analysis Code, has been developed for solving complex, homogeneous, gas-phase chemical kinetics problems and contains sensitivity analysis for a variety of problems, including nonisothermal situations. This report is part 1 of a series of three reference publications that describe LENS, provide a detailed guide to its usage, and present many example problems. Part 1 derives the governing equations and describes the numerical solution procedures for the types of problems that can be solved. The accuracy and efficiency of LSENS are examined by means of various test problems, and comparisons with other methods and codes are presented. LSENS is a flexible, convenient, accurate, and efficient solver for chemical reaction problems such as static system; steady, one-dimensional, inviscid flow; reaction behind incident shock wave, including boundary layer correction; and perfectly stirred (highly backmixed) reactor. In addition, the chemical equilibrium state can be computed for the following assigned states: temperature and pressure, enthalpy and pressure, temperature and volume, and internal energy and volume. For static problems the code computes the sensitivity coefficients of the dependent variables and their temporal derivatives with respect to the initial values of the dependent variables and/or the three rate coefficient parameters of the chemical reactions.

Radhakrishnan, Krishnan↗

LSENS: A General Chemical Kinetics and Sensitivity Analysis Code for homogeneous gas-phase reactions. Part 3: Illustrative test problems

LSENS, the Lewis General Chemical Kinetics and Sensitivity Analysis Code, has been developed for solving complex, homogeneous, gas-phase chemical kinetics problems and contains sensitivity analysis for a variety of problems, including nonisothermal situations. This report is part 3 of a series of three reference publications that describe LSENS, provide a detailed guide to its usage, and present many example problems. Part 3 explains the kinetics and kinetics-plus-sensitivity analysis problems supplied with LSENS and presents sample results. These problems illustrate the various capabilities of, and reaction models that can be solved by, the code and may provide a convenient starting point for the user to construct the problem data file required to execute LSENS. LSENS is a flexible, convenient, accurate, and efficient solver for chemical reaction problems such as static system; steady, one-dimensional, inviscid flow; reaction behind incident shock wave, including boundary layer correction; and perfectly stirred (highly backmixed) reactor. In addition, the chemical equilibrium state can be computed for the following assigned states: temperature and pressure, enthalpy and pressure, temperature and volume, and internal energy and volume. For static problems the code computes the sensitivity coefficients of the dependent variables and their temporal derivatives with respect to the initial values of the dependent variables and/or the three rate coefficient parameters of the chemical reactions.

Bittker, David A.↗

LSENS, a general chemical kinetics and sensitivity analysis code for homogeneous gas-phase reactions. 2: Code description and usage

LSENS, the Lewis General Chemical Kinetics Analysis Code, has been developed for solving complex, homogeneous, gas-phase chemical kinetics problems and contains sensitivity analysis for a variety of problems, including nonisothermal situations. This report is part 2 of a series of three reference publications that describe LSENS, provide a detailed guide to its usage, and present many example problems. Part 2 describes the code, how to modify it, and its usage, including preparation of the problem data file required to execute LSENS. Code usage is illustrated by several example problems, which further explain preparation of the problem data file and show how to obtain desired accuracy in the computed results. LSENS is a flexible, convenient, accurate, and efficient solver for chemical reaction problems such as static system; steady, one-dimensional, inviscid flow; reaction behind incident shock wave, including boundary layer correction; and perfectly stirred (highly backmixed) reactor. In addition, the chemical equilibrium state can be computed for the following assigned states: temperature and pressure, enthalpy and pressure, temperature and volume, and internal energy and volume. For static problems the code computes the sensitivity coefficients of the dependent variables and their temporal derivatives with respect to the initial values of the dependent variables and/or the three rate coefficient parameters of the chemical reactions. Part 1 (NASA RP-1328) derives the governing equations describes the numerical solution procedures for the types of problems that can be solved by lSENS. Part 3 (NASA RP-1330) explains the kinetics and kinetics-plus-sensitivity-analysis problems supplied with LSENS and presents sample results.

Radhakrishnan, Krishnan↗

LSENS, A General Chemical Kinetics and Sensitivity Analysis Code for Homogeneous Gas-Phase Reactions: Code Description and Usage - Part 2

LSENS, the Lewis General Chemical Kinetics and Sensitivity Analysis Code, has been developed for solving complex, homogeneous, gas-phase chemical kinetics problems and contains sensitivity analysis for a variety of problems, including nonisothermal situations. This report is part II of a series of three reference publications that describe LSENS, provide a detailed guide to its usage, and present many example problems. Part II describes the code, how to modify it, and its usage, including preparation of the problem data file required to execute LSENS. Code usage is illustrated by several example problems, which further explain preparation of the problem data file and show how to obtain desired accuracy in the computed results. LSENS is a flexible, convenient, accurate, and efficient solver for chemical reaction problems such as static system; steady, one-dimensional, inviscid flow; reaction behind incident shock wave, including boundary layer correction; and perfectly stirred (highly backmixed) reactor. In addition, the chemical equilibrium state can be computed for the following assigned states: temperature and pressure, enthalpy and pressure, temperature and volume, and internal energy and volume. For static problems the code computes the sensitivity coefficients of the dependent variables and their temporal derivatives with respect to the initial values of the dependent variables and/or the three rate coefficient parameters of the chemical reactions. Part I (NASA RP-1328) derives the governing equations and describes the numerical solution procedures for the types of problems that can be solved by LSENS. Part III (NASA RP-1330) explains the kinetics and kinetics-plus-sensitivity-analysis problems supplied with LSENS and presents sample results.

Radhakrishnan, Krishnan↗

Mixing and Transition Control Studied

Considerable progress in understanding nonlinear phenomena in both unbounded and wallbounded shear flow transition has been made through the use of a combination of high- Reynolds-number asymptotic and numerical methods. The objective of this continuing work is to fully understand the nonlinear dynamics so that ultimately (1) an effective means of mixing and transition control can be developed and (2) the source terms in the aeroacoustic noise problem can be modeled more accurately. Two important aspects of the work are that (1) the disturbances evolve from strictly linear instability waves on weakly nonparallel mean flows so that the proper upstream conditions are applied in the nonlinear or wave-interaction streamwise region and (2) the asymptotic formulations lead to parabolic problems so that the question of proper out-flow boundary conditions--still a research issue for direct numerical simulations of convectively unstable shear flows--does not arise. Composite expansion techniques are used to obtain solutions that account for both mean-flow-evolution and nonlinear effects. A previously derived theory for the amplitude evolution of a two-dimensional instability wave in an incompressible mixing layer (which is in quantitative agreement with available experimental data for the first nonlinear saturation stage for a plane-jet shear layer, a circular-jet shear layer, and a mixing layer behind a splitter plate) have been extended to include a wave-interaction stage with a three-dimensional subharmonic. The ultimate wave interaction effects can either give rise to explosive growth or an equilibrium solution, both of which are intimately associated with the nonlinear self-interaction of the three dimensional component. The extended theory is being evaluated numerically. In contrast to the mixing-layer situation, earlier comparisons of theoretical predictions based on asymptotic methods and experiments in wall-bounded shear-flow transition have been somewhat lacking in one aspect or another. The current work strongly suggests that the main weakness is the underlying asymptotic representation of the linear "part" of the problem and not the explicit modeling of the nonlinear/wave-interaction effects. Consequently, the long-wave-length/high-Reynolds-number asymptotic limit for the Blasius boundary-layer stability problem was reexamined, and a new dispersion relationship for the instability waves that is uniformly valid for both the upper- and lower branch regions to the required order of approximation was obtained. A comparison with numerical results, obtained by solving the Orr-Sommerfeld stability problem, shows that the asymptotic formula provides surprisingly good results, even for values of the frequency parameter usually encountered in experimental investigations. This is particularly evident in the dynamically important upper-branch region, where much of the nonlinear interactions in transition experiments are believed to take place. The result is important in that it can be used to greatly improve the accuracy of weakly nonlinear critical-layer-based theories, and a consistent nonlinear theory is currently under evaluation.

Source record↗

Spatially quasi-periodic water waves of finite depth

We present a numerical study of spatially quasi-periodic gravity-capillary waves of finite depth in both the initial value problem and travelling wave settings. We adopt a quasi-periodic conformal mapping formulation of the Euler equations, where one-dimensional quasi-periodic functions are represented by periodic functions on a higher-dimensional torus. We compute the time evolution of free surface waves in the presence of a background flow and a quasi-periodic bottom boundary and observe the formation of quasi-periodic patterns on the free surface. Two types of quasi-periodic travelling waves are computed: small-amplitude waves bifurcating from the zero-amplitude solution and larger-amplitude waves bifurcating from finite-amplitude periodic travelling waves. We derive weakly nonlinear approximations of the first type and investigate the associated small-divisor problem. We find that waves of the second type exhibit striking nonlinear behaviour, e.g. the peaks and troughs are shifted non-periodically from the corresponding periodic waves due to the activation of quasi-periodic modes.

Science & Technology - Other Topics↗

A new solution to the accretion problem.

The Bondi-Hoyle theory of line-accretion is briefly reviewed, and it is shown that there are an infinity of steady-state solutions satisfying the requirements originally imposed as boundary conditions. Questions are raised as to the applicability of these solutions to physical reality. The existence of an entirely different type of steady-state solution that implies slow velocity in the accretion stream beyond the neutral point is demonstrated. Accurate numerical values of this solution are obtained for a number of cases. The braking action given by this solution, for any selected cut-off distance, is far stronger than in the Bondi-Hoyle solution. Consideration of the physical state in the accretion stream shows that for this slow solution the appropriate cutoff range would be only of the order of 100 AU rather than of interstellar-distance order.

Lyttleton, R. A.↗

Corrections for attached sidewall boundary-layer effects in 2-dimensional airfoil testing

The problems of sidewall boundary-layer effects in airfoil testing is treated by considering the changes in the flow area due to boundary-layer thinning under the influence of the airfoil flowfield. Using von Karman's momentum integral equation, it is shown that the sidewall boundary-layer thickness in the region of the airfoil can reduce to about half the undisturbed value under the conditions prevailing in testing of supercritical airfoils. A Mach number correction due to this increased width of the flow passage is proposed. Using the small disturbance approximation, the effect of the sidewall boundary-layers is shown to be equivalent to a change in the test Mach number and also in the airfoil thickness. Comparison of the results of this approach with other similarity rules and correlation of the experimental data demonstrate the applicability of the analysis presented from low speeds to transonic speeds.

Murthy, A. V.↗

Numerical simulations of turbulent thermal convection with differential rotation

The solar nebula, from which the planets in our solar system formed, featured a disk of gas and dust grains in rapid, differential rotation, and at some stage was likely to have been unstable to thermal convection. This situation is suspected by many to lead to significant turbulent Reynolds stress production and angular momentum transport in such systems, and estimates of transport rates have been attempted from unsubstantiated phenomenological models. In order to determine the circumstances and physical conditions under which our own planetary system formed and to explain recent observations of young stellar systems, it is necessary to develop realistic models of heat and angular momentum transport for such flows. Developing an understanding of complicated flows featuring thermal convection, rotation, and shear is also of wide interest in stellar astrophysics and in planetary and terrestrial atmospheric studies. The ultimate objective is to develop workable models based on the numerical simulations for constructing global solar nebula models; viz., relatively simple prescriptions for heat and angular momentum fluxes from given system parameters (e.g., ratios of rotation, shear, and convective lapse rates) are characterized, quantified, and developed. Toward this end, our program has been to attempt to understand the behavior of the direct numerical simulations of Boussinesq convection, which, despite the complexity of the results, is still an overly simplified approximation to the real system and should be more amenable to analysis. These results are also intended to be tested against turbulence models, especially those designed for atmospheric boundary layers, and may provide a basis for subgrid-scale models. In order to make the numerical simulations more realistic with regard to the solar nebula problem, a fully compressible code that will allow incorporation of large density stratifications and realistic thermodynamic and radiative properties is developed. In order to explore the properties of these flows at the very high values of Re found in natural systems and the very low values of Pr found in most astrophysical contexts, we will need to employ large-eddy simulations for which we want to determine the most appropriate subgrid-scale model to incorporate.

Cabot, William H.↗

A review and analysis of neural networks for classification of remotely sensed multispectral imagery

A literature survey and analysis of the use of neural networks for the classification of remotely sensed multispectral imagery is presented. As part of a brief mathematical review, the backpropagation algorithm, which is the most common method of training multi-layer networks, is discussed with an emphasis on its application to pattern recognition. The analysis is divided into five aspects of neural network classification: (1) input data preprocessing, structure, and encoding; (2) output encoding and extraction of classes; (3) network architecture, (4) training algorithms; and (5) comparisons to conventional classifiers. The advantages of the neural network method over traditional classifiers are its non-parametric nature, arbitrary decision boundary capabilities, easy adaptation to different types of data and input structures, fuzzy output values that can enhance classification, and good generalization for use with multiple images. The disadvantages of the method are slow training time, inconsistent results due to random initial weights, and the requirement of obscure initialization values (e.g., learning rate and hidden layer size). Possible techniques for ameliorating these problems are discussed. It is concluded that, although the neural network method has several unique capabilities, it will become a useful tool in remote sensing only if it is made faster, more predictable, and easier to use.

Paola, Justin D.↗