Search NASA⌕ Search

SEARCH · Search NASA

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,297 records · Page 72
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

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↗
NASA NTRS1994

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.↗
NASA NTRS1994

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↗
NASA NTRS1994

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↗
NASA NTRS1996

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↗
DOE OSTI2023

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↗
DOE OSTI2016

Predicting Ice Sheet and Climate Evolution at Extreme Scales

A main research objectives of PISCEES is the development of formal methods for quantifying uncertainties in ice sheet modeling. Uncertainties in simulating and projecting mass loss from the polar ice sheets arise primarily from initial conditions, surface and basal boundary conditions, and model parameters. In general terms, two main chains of uncertainty propagation may be identified: 1. inverse propagation of observation and/or prior onto posterior control variable uncertainties; 2. forward propagation of prior or posterior control variable uncertainties onto those of target output quantities of interest (e.g., climate indices or ice sheet mass loss). A related goal is the development of computationally efficient methods for producing initial conditions for an ice sheet that are close to available present-day observations and essentially free of artificial model drift, which is required in order to be useful for model projections (“initialization problem”). To be of maximum value, such optimal initial states should be accompanied by “useful” uncertainty estimates that account for the different sources of uncerainties, as well as the degree to which the optimum state is constrained by available observations. The PISCEES proposal outlined two approaches for quantifying uncertainties. The first targets the full exploration of the uncertainty in model projections with sampling-based methods and a workflow managed by DAKOTA (the main delivery vehicle for software developed under QUEST). This is feasible for low-dimensional problems, e.g., those with a handful of global parameters to be inferred. This approach can benefit from derivative/adjoint information, but it is not necessary, which is why it often referred to as “non-intrusive”. The second approach makes heavy use of derivative information from model adjoints to address quantifying uncertainty in high-dimensions (e.g., basal boundary conditions in ice sheet models). The use of local gradient, or Hessian information (i.e., second derivatives of the cost function), requires additional code development and implementation, and is thus often referred to as an “intrusive” approach. Within PISCEES, MIT has been tasked to develop methods for derivative-based UQ, the ”intrusive” approach discussed above. These methods rely on the availability of first (adjoint) and second (Hessian) derivative code, developed through intrusive methods such as algorithmic differentiation (AD). While representing a significant burden in terms of code development, derivative-baesd UQ is able to cope with very high-dimensional uncertainty spaces. That is, unlike sampling methods (all variations of Monte Carlo), calculational burden is independent of the dimension of the uncertainty space. This is a significant advantage for spatially distributed uncertainty fields, such as threedimensional initial conditions, three-dimensional parameter fields, or two-dimensional surface and basal boundary conditions. Importantly, uncertainty fields for ice sheet models generally fall into this category.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Inertial impedance of coalescence during collision of liquid drops

The fluid dynamics of the collision and coalescence of liquid drops has intrigued scientists and engineers for more than a century owing to its ubiquitousness in nature, e.g. raindrop coalescence, and industrial applications, e.g. breaking of emulsions in the oil and gas industry. The complexity of the underlying dynamics, which includes occurrence of hydrodynamic singularities, has required study of the problem at different scales – macroscopic, mesoscopic and molecular – using stochastic and deterministic methods. In this work, a multi-scale, deterministic method is adopted to simulate the approach, collision, and eventual coalescence of two drops where the drops as well as the ambient fluid are incompressible, Newtonian fluids. The free boundary problem governing the dynamics consists of the Navier–Stokes system and associated initial and boundary conditions that have been augmented to account for the effects of disjoining pressure as the separation between the drops becomes of the order of a few hundred nanometres. This free boundary problem is solved by a Galerkin finite element-based algorithm. The interplay of inertial, viscous, capillary and van der Waals forces on the coalescence dynamics is investigated. It is shown that, in certain situations, because of inertia two drops that are driven together can first bounce before ultimately coalescing. This bounce delays coalescence and can result in the computed value of the film drainage time departing significantly from that predicted from existing scaling theories.

Mechanics↗
NASA NTRS1972

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.↗
NASA NTRS1985

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.↗
NASA NTRS1991

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.↗
NASA NTRS1993

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.↗
NASA NTRS1942

New method of extrapolation of the resistance of a model planing boat to full size

The previously employed method of extrapolating the total resistance to full size with lambda(exp 3) (model scale) and thereby foregoing a separate appraisal of the frictional resistance, was permissible for large models and floats of normal size. But faced with the ever increasing size of aircraft a reexamination of the problem of extrapolation to full size is called for. A method is described by means of which, on the basis of an analysis of tests on planing surfaces, the variation of the wetted surface over the take-off range is analytically obtained. The friction coefficients are read from Prandtl's curve for turbulent boundary layer with laminar approach. With these two values a correction for friction is obtainable.

Sottorf, W↗
DOE OSTI2014

Asymptotic theory of double layer and shielding of electric field at the edge of illuminated plasma

The method of matched asymptotic expansions is applied to the problem of a collisionless plasma generated by UV illumination localized in a central part of the plasma in the limiting case of small Debye length λ{sub D}. A second-approximation asymptotic solution is found for the double layer positioned at the boundary of the illuminated region and for the un-illuminated plasma for the plane geometry. Numerical calculations for different values of λ{sub D} are reported and found to confirm the asymptotic results. The net integral space charge of the double layer is asymptotically small, although in the plane geometry it is just sufficient to shield the ambipolar electric field existing in the illuminated region and thus to prevent it from penetrating into the un-illuminated region. The double layer has the same mathematical nature as the intermediate transition layer separating an active plasma and a collisionless sheath, and the underlying physics is also the same. In essence, the two layers represent the same physical object: a transonic layer.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

fPINNs: Fractional Physics-Informed Neural Networks

Physics-informed neural networks (PINNs), introduced, are effective in solving integer-order partial differential equations (PDEs) based on scattered and noisy data. PINNs employ standard feedforward neural networks (NNs) with the PDEs explicitly encoded into the NN using automatic differentiation, while the sum of the mean-squared PDE residuals and the mean-squared error in initial-boundary conditions is minimized with respect to the NN parameters. Here we extend PINNs to fractional PINNs (fPINNs) to solve space-time fractional advection-diffusion equations (fractional ADEs), and we study systematically their convergence, hence explaining both fPINNs and PINNs for the first time. Specifically, we demonstrate their accuracy and effectiveness in solving multidimensional forward and inverse problems with forcing terms whose values are only known at randomly scattered spatio-temporal coordinates (black-box (BB) forcing terms). A novel element of the fPINNs is the hybrid approach that we introduce for constructing the residual in the loss function using both automatic differentiation for the integer-order operators and numerical discretization for the fractional operators. This approach bypasses the difficulties stemming from the fact that automatic differentiation is not applicable to fractional operators because the standard chain rule in integer calculus is not valid in fractional calculus. To discretize the fractional operators, we employ the Grünwald--Letnikov (GL) formula in one-dimensional fractional ADEs and the vector GL formula in conjunction with the directional fractional Laplacian in two- and three-dimensional fractional ADEs. We first consider the one-dimensional fractional Poisson equation and compare the convergence of the fPINNs against the finite difference method (FDM). We present the solution convergence using both the mean L 2 error as well as the standard deviation due to sensitivity to NN parameter initializations. Using different GL formulas we observe first-, second-, and third-order convergence rates for small size training sets but the error saturates for larger training sets. We explain these results by analyzing the four sources of numerical errors due to discretization, sampling, NN approximation, and optimization. The total error decays monotonically (below 10 –5 for a third-order GL formula) but it saturates beyond that point due to the optimization error. We also analyze the relative balance between discretization and sampling errors and observe that the sampling size and the number of discretization points (auxiliary points) should be comparable to achieve the highest accuracy. As we increase the depth of the NN up to certain value, the mean error decreases and the standard deviation increases, whereas the width has essentially no effect unless its value is either too small or too large. We next consider time-dependent fractional ADEs and compare white-box (WB) and BB forcing. We observe that for the WB forcing, our results are similar to the aforementioned cases; however, for the BB forcing fPINNs outperform FDM. Subsequently, we consider multidimensional time-, space-, and space-time-fractional ADEs using the directional fractional Laplacian and we observe relative errors of 10 –3 ~ 10 –4 . Lastly, we solve several inverse problems in one, two, and three dimensions to identify the fractional orders, diffusion coefficients, and transport velocities and obtain accurate results given proper initializations even in the presence of significant noise.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1974

A curvature-corrected Kirchhoff formulation for radar sea-return from the near vertical

A new theoretical treatment of the problem of electromagnetic wave scattering from a randomly rough surface is given. A high frequency correction to the Kirchhoff approximation is derived from a field integral equation for a perfectly conducting surface. The correction, which accounts for the effect of local surface curvature, is seen to be identical with an asymptotic form found by Fock (1945) for diffraction by a paraboloid. The corrected boundary values are substituted into the far field Stratton-Chu integral, and average backscattered powers are computed assuming the scattering surface is a homogeneous Gaussian process. Preliminary calculations for K(-4) ocean wave spectrum indicate a resonable modelling of polarization effects near the vertical, theta 45 deg. Correspondence with the results of small perturbation theory is shown.

Jackson, F. C.↗
NASA NTRS1975

A study of reacting free and ducted hydrogen/air jets

The mixing and reaction of a supersonic jet of hydrogen in coaxial free and ducted high temperature test gases were investigated. The importance of chemical kinetics on computed results, and the utilization of free-jet theoretical approaches to compute enclosed flow fields were studied. Measured pitot pressure profiles were correlated by use of a parabolic mixing analysis employing an eddy viscosity model. All computations, including free, ducted, reacting, and nonreacting cases, use the same value of the empirical constant in the viscosity model. Equilibrium and finite rate chemistry models were utilized. The finite rate assumption allowed prediction of observed ignition delay, but the equilibrium model gave the best correlations downstream from the ignition location. Ducted calculations were made with finite rate chemistry; correlations were, in general, as good as the free-jet results until problems with the boundary conditions were encountered.

Beach, H. L., Jr.↗