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.

SourceResearch connectionFrom yearThrough year

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

At least 1,153 records · Page 64

All sources NASA NTRS DOE OSTI

NASA NTRS1979

A new approach to the correlation of boundary layer mass transfer rates with thermal diffusion and/or variable properties

A rational approach to the correlation of boundary layer mass transport rates, applicable to many commonly encountered laminar flow conditions with thermal diffusion and/or variable properties, is outlined. The correlation scheme builds upon already available constant property blowing/suction solutions by introducing appropriate correction factors to account for the additional ('pseudo' blowing and source) effects identified with variable properties and thermal diffusion. Applications of the scheme to the particular laminar boundary layer mass transfer problems considered herein (alkali and transition metal compound vapor transport) indicates satisfactory accuracy up to effective blowing factors equivalent to about one third of the 'blow off' value. As a useful by-product of the variable property correlation, we extend the heat-mass transfer analogy, for a wide range of Lewis numbers, to include variable property effects.

Srivastava, R.↗

NASA NTRS1982

A numerical method for the dynamics of non-spherical cavitation bubbles

A boundary integral numerical method for the dynamics of nonspherical cavitation bubbles in inviscid incompressible liquids is described. Only surface values of the velocity potential and its first derivatives are involved. The problem of solving the Laplace equation in the entire domain occupied by the liquid is thus avoided. The collapse of a bubble in the vicinity of a solid wall and the collapse of three bubbles with collinear centers are considered.

Lucca, G.↗

NASA NTRS1981

Eigenvalues of the Rayleigh-Benard and Marangoni problems

The eigenvalues of the linear Benard-Marangoni stability problem are discussed. Pearson and Nield boundary conditions, which correspond to a rigid, isothermal lower boundary and a stress-free conducting upper boundary are considered. It is shown that although a critical value of the Marangoni number can be determined, the number is not, strictly speaking, an eigenvalue and cannot be used as an eigenvalue parameter for the determination of an eigenvector set.

Rosenblat, S.↗

NASA NTRS2013

Automatic Data Filter Customization Using a Genetic Algorithm

This work predicts whether a retrieval algorithm will usefully determine CO2 concentration from an input spectrum of GOSAT (Greenhouse Gases Observing Satellite). This was done to eliminate needless runtime on atmospheric soundings that would never yield useful results. A space of 50 dimensions was examined for predictive power on the final CO2 results. Retrieval algorithms are frequently expensive to run, and wasted effort defeats requirements and expends needless resources. This algorithm could be used to help predict and filter unneeded runs in any computationally expensive regime. Traditional methods such as the Fischer discriminant analysis and decision trees can attempt to predict whether a sounding will be properly processed. However, this work sought to detect a subsection of the dimensional space that can be simply filtered out to eliminate unwanted runs. LDAs (linear discriminant analyses) and other systems examine the entire data and judge a "best fit," giving equal weight to complex and problematic regions as well as simple, clear-cut regions. In this implementation, a genetic space of "left" and "right" thresholds outside of which all data are rejected was defined. These left/right pairs are created for each of the 50 input dimensions. A genetic algorithm then runs through countless potential filter settings using a JPL computer cluster, optimizing the tossed-out data s yield (proper vs. improper run removal) and number of points tossed. This solution is robust to an arbitrary decision boundary within the data and avoids the global optimization problem of whole-dataset fitting using LDA or decision trees. It filters out runs that would not have produced useful CO2 values to save needless computation. This would be an algorithmic preprocessing improvement to any computationally expensive system.

Mandrake, Lukas↗

DOE OSTI2015

Generalized Multiscale Finite-Element Method (GMsFEM) for elastic wave propagation in heterogeneous, anisotropic media

It is important to develop fast yet accurate numerical methods for seismic wave propagation to characterize complex geological structures and oil and gas reservoirs. However, the computational cost of conventional numerical modeling methods, such as finite-difference method and finite-element method, becomes prohibitively expensive when applied to very large models. We propose a Generalized Multiscale Finite-Element Method (GMsFEM) for elastic wave propagation in heterogeneous, anisotropic media, where we construct basis functions from multiple local problems for both the boundaries and interior of a coarse node support or coarse element. The application of multiscale basis functions can capture the fine scale medium property variations, and allows us to greatly reduce the degrees of freedom that are required to implement the modeling compared with conventional finite-element method for wave equation, while restricting the error to low values. We formulate the continuous Galerkin and discontinuous Galerkin formulation of the multiscale method, both of which have pros and cons. Applications of the multiscale method to three heterogeneous models show that our multiscale method can effectively model the elastic wave propagation in anisotropic media with a significant reduction in the degrees of freedom in the modeling system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Physics-Informed Neural Network Method for Forward and Backward Advection-Dispersion Equations

Advection-dispersion equations (ADEs) are commonly used to describe transport phenomena in porous media. Even though mature discretization-based numerical methods for ADEs exist, some challenges still remain, especially when it comes to solving advection-dominated forward ADEs and diffusion-dominated backward ADEs. The latter problem usually arises in the source identification context and leads to numerically unstable grid-based solutions that require a form of regularization or should be treated as an inverse problem that is computationally more expensive because it requires solving the forward problem multiple times. In this study, we propose a discretization-free approach based on the physics-informed neural network (PINN) method for solving coupled ADE and Darcy flow equations with space-dependent hydraulic conductivity. In this approach, the hydraulic conductivity, hydraulic head, and concentration fields are approximated with deep neural networks (DNNs). We assume that the conductivity field is given by its values on a grid, and we use these values to train the conductivity DNN. The head and concentration DNNs are trained by minimizing the residuals of the flow equation and ADE and using the initial and boundary conditions as additional constraints. The PINN method is applied to one- and two-dimensional forward ADE problems, where its performance for various P\'{e}clet numbers ($Pe$) is compared with the analytical and numerical solutions. We find that the PINN method is accurate with errors of less than 1\% and outperforms some conventional discretization-based methods for $Pe$ larger that 100. Next, we demonstrate that the PINN method remains accurate for the backward ADEs, with the relative errors in most cases staying under 5\% compared to the reference concentration field. Finally, we show that when available, the concentration measurements can be easily incorporated in the PINN method and significantly improve (by more than 50\% in the considered cases) the accuracy of the PINN solution of the backward ADE.

He, Qizhi↗

DOE OSTI2020

Bi-cross validation of spectral clustering hyperparameters

One challenge impeding the analysis of terabyte scale X-ray scattering data from the Linac Coherent Light Source (LCLS) is determining the number of clusters required for the execution of traditional clustering algorithms. Here, we demonstrate that the previous work using bi-cross validation to determine the number of singular vectors directly maps to the spectral clustering problem of estimating both the number of clusters and hyperparameter values. Here by applying this method to LCLS X-ray scattering data enables the identification of dropped shots without manually setting boundaries on detector fluence and provides a path toward identifying rare and anomalous events.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1990

Numerical simulation of the incompressible internal flow through a tilting disk valve

A numerical simulation of the incompressible viscous flow through a prosthetic tilting disk heart valve is presented in order to demonstrate the current capability to model unsteady flows with moving boundaries. Both steady and unsteady flow calculations are performed by solving the incompressible Navier-Stokes equations in three-dimensional generalized curvilinear coordinates. In order to handle the moving boundary problems, the chimera grid embedding scheme which decomposes a complex computational domain into several simple subdomains is used. An algebraic turbulence model for internal flows is incorporated to reach the physiological values of Reynolds number. Good agreement is obtained between the numerical results and experimental measurements. It is found that the tilting disk valve causes large regions of separated flow, and regions of high shear.

Chang, I-Dee↗

NASA NTRS1991

Simulation of blood flow through an artificial heart

A numerical simulation of the incompressible viscous flow through a prosthetic tilting disk heart valve is presented in order to demonstrate the current capability to model unsteady flows with moving boundaries. Both steady state and unsteady flow calculations are done by solving the incompressible Navier-Stokes equations in 3-D generalized curvilinear coordinates. In order to handle the moving boundary problems, the chimera grid embedding scheme which decomposes a complex computational domain into several simple subdomains is used. An algebraic turbulence model for internal flows is incorporated to reach the physiological values of Reynolds number. Good agreement is obtained between the numerical results and experimental measurements. It is found that the tilting disk valve causes large regions of separated flow, and regions of high shear.

Kiris, Cetin↗

DOE OSTI2016

A thermodynamic study of shear banding in polymer solutions

Although shear banding is a ubiquitous phenomenon observed in soft materials, the mechanisms that give rise to shear-band formation are not always the same. In this work, we develop a new two-fluid model for semi-dilute entangled polymer solutions using the generalized bracket approach of nonequilibrium thermodynamics. The model is based on the hypothesis that the direct coupling between polymer stress and concentration is the driving mechanism of steady shear-band formation. To obtain smooth banded profiles in the two-fluid framework, a new stress-diffusive term is added to the time evolution equation for the conformation tensor. The advantage of the new model is that the differential velocity is treated as a state variable. This allows a straightforward implementation of the additional boundary conditions arising from the derivative diffusive terms with respect to this new state variable. To capture the overshoot of the shear stress during the start of a simple shear flow, we utilize a nonlinear Giesekus relaxation. Moreover, we include an additional relaxation term that resembles the term used in the Rouse linear entangled polymer model to account for convective constraint release and chain stretch to generate the upturn of the flow curve at large shear rates. Numerical calculations performed for cylindrical Couette flow confirm the independency of the solution from the deformation history and initial conditions. Furthermore, we find that stress-induced migration is the responsible diffusive term for steady-state shear banding. Because of its simplicity, the new model is an ideal candidate for the use in the simulation of more complex flows.

42 ENGINEERING↗

NASA NTRS1974

A finite element formulation for supersonic flows around complex configurations

The problem of small perturbation potential supersonic flow around complex configurations is considered. This problem requires the solution of an integral equation relating the values of the potential on the surface of the body to the values of the normal derivative, which is known from the small perturbation boundary conditions. The surface of the body is divided into small (hyperboloidal quadrilateral) surface elements which are described in terms of the Cartesian components of the four corner points. The values of the potential (and its normal derivative) within each element are assumed to be constant and equal to its value at the centroid of the element. This yields a set of linear algebraic equations whose coefficients are given by source and doublet integrals over the surface elements. Closed form evaluations of the integrals are presented.

Morino, L.↗

NASA NTRS1974

A finite-element analysis for steady and oscillatory supersonic flows around complex configurations

The problem of small perturbation potential supersonic flow around complex configurations is considered. This problem requires the solution of an integral equation relating the values of the potential on the surface of the body to the values of the normal derivative, which is known from the small perturbation boundary conditions. The surface of the body is divided into small (hyperboloidal quadrilateral) surface elements, sigma sub i, which are described in terms of the Cartesian components of the four corner points. The values of the potential (and its normal derivative) within each element is assumed to be constant and equal to its value at the centroid of the element, and this yields a set of linear algebraic equations. The coefficients of the equation are given by source and doublet integrals over the surface elements, sigma sub i. The results obtained using the above formulation are compared with existing analytical and experimental results.

Morino, L.↗

NASA NTRS1989

Thrust vectoring effect on time-optimal 90 degrees angle of attack pitch up maneuvers of a high alpha fighter aircraft

Using thrust vectoring, the problem of pitching up a high alpha fighter aircraft to 90 degrees angle of attack in minimum time is considered. Pontryagin's maximum principle is used together with a two-point boundary value numerical algorithm to derive open-loop controls for various parameterized limits on thrust vectoring angles. Without thrust vectoring, a high alpha fighter with initial conditions 0.6 Mach number and 15,000 feet requires 7.0 seconds to pitch-up to 90 degrees. However, with 20-degree angle thrust vectoring, the pitch-up takes only 1.5 seconds with a final Mach number of 0.47. Also, the gain in altitude is less than 100 feet. More important, additional thrust vectoring yields litle benefit beyond the 20-degree angles. This paper describes control solutions in detail for thrust vectoring angles between 0 and 90 degrees.

Stalford, Harold↗

NASA NTRS1971

The problem of uniqueness in thermosphere dynamics

The unique solution of tidal wave propagation within the thermosphere depends on the boundary conditions of the model. It is shown that the radiation condition leads to such a unique solution. Any other boundary values give rise to deviations from the physical solution with unrealistic physical parameters in the environment of the boundaries. The thickness of these boundary layers with unrealistic solutions is a few scale heights below the upper boundary of the model and a few tens of scale heights above the lower boundary of the model.

Volland, H.↗

DOE OSTI2018

Iterated Discontinuity Factors Based on Reference Outgoing Currents for Discontinuous Finite Element Diffusion

We present a novel formulation for iteratively computing Discontinuity Factors (DF) designed to reproduce reference partial outgoing currents at the interface between macro-regions using Discontinuous Galerkin Finite Elements to solve the multigroup neutron diffusion equation. Similarly, an expression for the boundary coefficients (BCf) is proposed to match the partial outgoing current at vacuum boundaries. Upon convergence, macro-region averaged reaction and leakage rates are preserved. In addition, we consider a source-free steady-state problem depending on the DFs. The existence of a non-trivial solution is due to the non-linearity of the problem. This form greatly simplifies the implementation within Rattlesnake and efficiently converges by updating the DFs and BCfs at every linear iteration of the Preconditioned Jacobian-Free Newton-Krylov solver. The main challenge then lies in determining an appropriate scaling factor to ensure that the initial values for the DFs and BCfs are not too far from their converged value. While this factor does not change the final solution, it has the potential to significantly affect convergence, for purely numerical reasons. An expression based on the reference fission source integral is proposed and satisfying results are presented for problems and compared to Superhomogénéisation (SPH) results with both reflecting and vacuum boundary conditions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2018

Iterated Discontinuity Factors Based on Reference Outgoing Currents for Discontinuous Finite Element Diffusion

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2019

Self-Induction of a Fine Cylindrical Metallic Wire Versus a Mechanism of Electron Surface Scattering

The self-induction of a thin cylindrical metallic wire has been calculated. A general case when the ratio of the electron free path to the wire radius may take an arbitrary value has been considered. The dependence of the regular reflectance on surface defect density and angle of incidence of electrons on the wire’s inner surface has been taken for boundary conditions of the problem.

36 MATERIALS SCIENCE↗

NASA NTRS1993

Improved Finite Element Modeling of the Turbofan Engine Inlet Radiation Problem

Improvements have been made in the finite element model of the acoustic radiated field from a turbofan engine inlet in the presence of a mean flow. The problem of acoustic radiation from a turbofan engine inlet is difficult to model numerically because of the large domain and high frequencies involved. A numerical model with conventional finite elements in the near field and wave envelope elements in the far field has been constructed. By employing an irrotational mean flow assumption, both the mean flow and the acoustic perturbation problem have been posed in an axisymmetric formulation in terms of the velocity potential; thereby minimizing computer storage and time requirements. The finite element mesh has been altered in search of an improved solution. The mean flow problem has been reformulated with new boundary conditions to make it theoretically rigorous. The sound source at the fan face has been modeled as a combination of positive and negative propagating duct eigenfunctions. Therefore, a finite element duct eigenvalue problem has been solved on the fan face and the resulting modal matrix has been used to implement a source boundary condition on the fan face in the acoustic radiation problem. In the post processing of the solution, the acoustic pressure has been evaluated at Gauss points inside the elements and the nodal pressure values have been interpolated from them. This has significantly improved the results. The effect of the geometric position of the transition circle between conventional finite elements and wave envelope elements has been studied and it has been found that the transition can be made nearer to the inlet than previously assumed.

Roy, Indranil Danda↗