SEARCH · Search NASA

Results for “parabolic 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 379 records · Page 21

All sources NASA NTRS DOE OSTI

NASA NTRS1990

Collisionless tearing instability in magnetotail plasmas

The problem of the linear stability of collisionless tearing modes in the earth's magnetotail is revisited. It is found that the collisionless tearing mode is linearly unstable with wavelengths of the order of 10 R(E). It is shown that an important feature neglected in earlier theories is a nonzero equilibrium B(y)-field. The physics of the instability is elucidated in the context of a simple slab model and a sheared parabolic model which is representative of the magnetotail in which all three components of the magnetic field are nonzero. The dispersion equation for the instability is obtained by a boundary-layer analysis. The implications of the theory for recent observations on current disruption and diversion during substorms is discussed.

Wang, Xiaogang↗

NASA NTRS1986

Parameter estimation problems for distributed systems using a multigrid method

The problem of estimating spatially varying coefficients of partial differential equations is considered from observation of the solution and of the right hand side of the equation. It is assumed that the observations are distributed in the domain and that enough observations are given. A method of discretization and an efficient multigrid method for solving the resulting discrete systems are described. Numerical results are presented for estimation of coefficients in an elliptic and a parabolic partial differential equation.

Taasan, S.↗

NASA NTRS1988

Finite-dimensional approximation for optimal fixed-order compensation of distributed parameter systems

In controlling distributed parameter systems it is often desirable to obtain low-order, finite-dimensional controllers in order to minimize real-time computational requirements. Standard approaches to this problem employ model/controller reduction techniques in conjunction with LQG theory. In this paper we consider the finite-dimensional approximation of the infinite-dimensional Bernstein/Hyland optimal projection theory. This approach yields fixed-finite-order controllers which are optimal with respect to high-order, approximating, finite-dimensional plant models. The technique is illustrated by computing a sequence of first-order controllers for one-dimensional, single-input/single-output, parabolic (heat/diffusion) and hereditary systems using spline-based, Ritz-Galerkin, finite element approximation. Numerical studies indicate convergence of the feedback gains with less than 2 percent performance degradation over full-order LQG controllers for the parabolic system and 10 percent degradation for the hereditary system.

Bernstein, Dennis S.↗

NASA NTRS1990

Parameter estimation problems for distributed systems using a multigrid method

Ta'asan, Shlomo↗

NASA NTRS1991

Spatial Direct Numerical Simulation of Boundary-Layer Transition Mechanisms: Validation of PSE Theory

A study of instabilities in incompressible boundary-layer flow on a flat plate is conducted by spatial direct numerical simulation (DNS) of the Navier-Stokes equations. Here, the DNS results are used to critically evaluate the results obtained using parabolized stability equations (PSE) theory and to study mechanisms associated with breakdown from laminar to turbulent flow. Three test cases are considered: two-dimensional Tollmien-Schlichting wave propagation, subharmonic instability breakdown, and oblique-wave break-down. The instability modes predicted by PSE theory are in good quantitative agreement with the DNS results, except a small discrepancy is evident in the mean-flow distortion component of the 2-D test problem. This discrepancy is attributed to far-field boundary- condition differences. Both DNS and PSE theory results show several modal discrepancies when compared with the experiments of subharmonic breakdown. Computations that allow for a small adverse pressure gradient in the basic flow and a variation of the disturbance frequency result in better agreement with the experiments.

Joslin, R. D.↗

NASA NTRS1987

Design of Optimum Ducts Using an Efficient 3-D Viscous Computational Flow Analysis

Design of fluid dynamically efficient ducts is addressed through the combination of an optimization analysis with a three-dimensional viscous fluid dynamic analysis code. For efficiency, a parabolic fluid dynamic analysis was used. Since each function evaluation in an optimization analysis is a full three-dimensional viscous flow analysis requiring 200,000 grid points, it is important to use both an efficient fluid dynamic analysis and an efficient optimization technique. Three optimization techniques are evaluated on a series of test functions. The Quasi-Newton (BFGS, eta = .9) technique was selected as the preferred technique. A series of basic duct design problems are performed. On a two-parameter optimization problem, the BFGS technique is demonstrated to require half as many function evaluations as a steepest descent technique.

Madabhushi, Ravi K.↗

NASA NTRS1994

Parabolized stability equations

The parabolized stability equations (PSE) are a new approach to analyze the streamwise evolution of single or interacting Fourier modes in weakly nonparallel flows such as boundary layers. The concept rests on the decomposition of every mode into a slowly varying amplitude function and a wave function with slowly varying wave number. The neglect of the small second derivatives of the slowly varying functions with respect to the streamwise variable leads to an initial boundary-value problem that can be solved by numerical marching procedures. The PSE approach is valid in convectively unstable flows. The equations for a single mode are closely related to those of the traditional eigenvalue problems for linear stability analysis. However, the PSE approach does not exploit the homogeneity of the problem and, therefore, can be utilized to analyze forced modes and the nonlinear growth and interaction of an initial disturbance field. In contrast to the traditional patching of local solutions, the PSE provide the spatial evolution of modes with proper account for their history. The PSE approach allows studies of secondary instabilities without the constraints of the Floquet analysis and reproduces the established experimental, theoretical, and computational benchmark results on transition up to the breakdown stage. The method matches or exceeds the demonstrated capabilities of current spatial Navier-Stokes solvers at a small fraction of their computational cost. Recent applications include studies on localized or distributed receptivity and prediction of transition in model environments for realistic engineering problems. This report describes the basis, intricacies, and some applications of the PSE methodology.

Herbert, Thorwald↗

NASA NTRS1995

Japan's research on particle clouds and sprays

Most of energy used by us is generated by combustion of liquid and solid fuels. These fuels are burned in combustors mainly as liquid sprays and pulverized solids, respectively. A knowledge of the combustion processes in combustors is needed to achieve proper designs that have stable operation, high efficiency, and low emission levels. However, current understanding of liquid and solid particle cloud combustion is far from complete. If combustion experiments for these fuels are performed under a normal gravity field, some experimental difficulties are encountered. These difficulties encountered include, that since the particles fall by the force of gravity it is impossible to stop the particles in the air, the falling speeds of particles are different from each other, and are depend on the particle size, the flame is lifted up and deformed by the buoyancy force, and natural convection makes the flow field more complex. Since these experimental difficulties are attributable to the gravity force, a microgravity field can eliminate the above problems. This means that the flame propagation experiments in static homogeneous liquid and solid particle clouds can be carried out under a microgravity field. This will provide much information for the basic questions related to combustion processes of particle clouds and sprays. In Japan, flame propagation processes in the combustible liquid and solid particle clouds have been studied experimentally by using a microgravity field generated by a 4.5 s dropshaft, a 10 s dropshaft, and by parabolic flight. Described in this presentation are the recent results of flame propagations studies in a homogeneous liquid particle cloud, in a mixture of liquid particles/gas fuel/air, in a PMMA particle cloud, and in a pulverized coal particle cloud.

Sato, Jun'ichi↗

NASA NTRS1971

Optimum transfer between hyperbolic asymptotes with turn angle less than the maximum natural turn angle.

The impulsive transfer of minimum characteristic velocity between two given hyperbolic asymptotes associated with a real planet is considered. The constraint of remaining above the surface of the planet introduces a maximum natural turn angle. This divides the problem naively into two categories. The first is the case when the turn angle required by the given asymptotes is less than or equal to the maximum natural turn angle. This is the more reasonable situation for a swingby maneuver, and is the case considered in this paper. The resulting optimal transfers are of five basic types. In boundary conditions space, by far the largest area of transfers are either single impulse transfers or transfers through the parabolic level requiring two finite impulses. In no case has an optimal transfer been found which uses more than two finite impulses and three infinitesimal impulses. Complete results are presented for the case in which the given turn angle is less than or equal to the maximum natural turn angle.

Walton, J. M.↗

NASA NTRS1988

An algorithm for unsteady flows with strong convection

An implicit ADI numerical method for the calculation of 2-D unsteady flows with strong convection effects is described. The method is based on the conventional Crank-Nicholson approach for parabolic equations but an upwind-downwind differencing is used for the first order spatial derivatives associated with convection. The differencing is carried out in the current and previous time plane in such a way that the algorithm is second order accurate in both space and time. The difference equations are factored into sequential operators, one in each independent spatial variable; the solution at each time step may then be computed as a sequence of tridiagonal matrix problems. The method may be used in a noniterative manner although iteration at each time step is recommended in situations where the effects of convection are strong.

Peridier, Vallorie↗

NASA NTRS1992

Parabolized reduced Navier-Stokes computational techniques

A review is presented of methods in which composite or reduced Navier-Stokes (RNS) equations are treated with a pressure-gradient-based flux-vector splitting. The methods are similar to large Re asymptotic formulations, and streamwise diffusion terms are ignored in favor of an explicit deferred corrector based on higher-order diffusion terms. The methods can be used for 2D and 3D supersonic flows in which the effects of real gas are incorporated. Several examples of the procedure are given, and subsonic and supersonic flows are handled well with relaxation procedures that incorporate multigrid acceleration. Effective solutions are described for problems ranging from incompressible flows and supersonic flows to sharp shocks and reverse-flow capturing.

Rubin, Stanley G.↗

DOE OSTI2017

Low- and high-order accurate boundary conditions: From Stokes to Darcy porous flow modeled with standard and improved Brinkman lattice Boltzmann schemes

The present contribution focuses on the accuracy of reflection-type boundary conditions in the Stokes–Brinkman–Darcy modeling of porous flows solved with the lattice Boltzmann method (LBM), which we operate with the two-relaxation-time (TRT) collision and the Brinkman-force based scheme (BF), called BF-TRT scheme. In parallel, we compare it with the Stokes–Brinkman–Darcy linear finite element method (FEM) where the Dirichlet boundary conditions are enforced on grid vertices. In bulk, both BF-TRT and FEM share the same defect: in their discretization a correction to the modeled Brinkman equation appears, given by the discrete Laplacian of the velocity-proportional resistance force. This correction modifies the effective Brinkman viscosity, playing a crucial role in the triggering of spurious oscillations in the bulk solution. While the exact form of this defect is available in lattice-aligned, straight or diagonal, flows; in arbitrary flow/lattice orientations its approximation is constructed. At boundaries, we verify that such a Brinkman viscosity correction has an even more harmful impact. Already at the first order, it shifts the location of the no-slip wall condition supported by traditional LBM boundary schemes, such as the bounce-back rule. For that reason, this work develops a new class of boundary schemes to prescribe the Dirichlet velocity condition at an arbitrary wall/boundary-node distance and that supports a higher order accuracy in the accommodation of the TRT-Brinkman solutions. For their modeling, we consider the standard BF scheme and its improved version, called IBF; this latter is generalized in this work to suppress or to reduce the viscosity correction in arbitrarily oriented flows. Our framework extends the one- and two-point families of linear and parabolic link-wise boundary schemes, respectively called B-LI and B-MLI, which avoid the interference of the Brinkman viscosity correction in their closure relations. The performance of LBM and FEM is thoroughly evaluated in three benchmark tests, which are run throughout three distinctive permeability regimes. The first configuration is a horizontal porous channel, studied with a symbolic approach, where we construct the exact solutions of FEM and BF/IBF with different boundary schemes. The second problem refers to an inclined porous channel flow, which brings in as new challenge the formation of spurious boundary layers in LBM; that is, numerical artefacts that arise due to a deficient accommodation of the bulk solution by the low-accurate boundary scheme. The third problem considers a porous flow past a periodic square array of solid cylinders, which intensifies the previous two tests with the simulation of a more complex flow pattern. The ensemble of numerical tests provides guidelines on the effect of grid resolution and the TRT free collision parameter over the accuracy and the quality of the velocity field, spanning from Stokes to Darcy permeability regimes. It is shown that, with the use of the high-order accurate boundary schemes, the simple, uniform-mesh-based TRT-LBM formulation can even surpass the accuracy of FEM employing hardworking body-fitted meshes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS1985

Long-time asymptotics of a system for plasma diffusion

A system of parabolic nonlinear equations that describe the diffusion of a fully collisional plasma across a strong magnetic field is discussed. It is demonstrated that the solution to this system tends to a time asymptotic state which is of space-time separable form, theta(t)f(x). Furthermore, f(x) is independent of the initial conditions and theta(t) depends slightly on the initial conditions. The rate of decay of the temporal part is governed by a nonlinear eigenvalue problem. Since the equations are considered in a bounded domain we are able to analyze the effect of boundary conditions on the evolution of the system. Additional effects on radiation, heating, and particle injection can also be accounted for. Essential differences between the behavior of a fully-coupled system and a scalar equation are observed.

Rosenau, P.↗

NASA NTRS1993

Microgravity combustion of dust suspensions

Unlike the combustion of homogeneous gas mixtures, there are practically no reliable fundamental data (i.e., laminar burning velocity, flammability limits, quenching distance, minimum ignition energy) for the combustion of heterogeneous dust suspensions. Even the equilibrium thermodynamic data such as the constant pressure volume combustion pressure and the constant pressure adiabatic flame temperature are not accurately known for dust mixtures. This is mainly due to the problem of gravity sedimentation. In normal gravity, turbulence, convective flow, electric and acoustic fields are required to maintain a dust in suspension. These external influences have a dominating effect on the combustion processes. Microgravity offers a unique environment where a quiescent dust cloud can in principle be maintained for a sufficiently long duration for almost all combustion experiments (dust suspensions are inherently unstable due to Brownian motion and particle aggregation). Thus, the microgravity duration provided by drop towers, parabolic flights, and the space shuttle, can all be exploited for different kinds of dust combustion experiments. The present paper describes some recent studies on microgravity combustion of dust suspension carried out on the KC-135 and the Caravelle aircraft. The results reported are obtained from three parabolic flight campaigns.

Lee, John H. S.↗

DOE OSTI2016

Observability of localized magnetoplasmons in quantum dots: Scrutinizing the eligibility of far-infrared, Raman, and electron-energy-loss spectroscopies

We investigate a one-component, quasi-zero dimensional, quantum plasma exposed to a parabolic potential and an applied magnetic field in the symmetric gauge. If the size of such a system as can be realized in the semiconducting quantum dots is on the order of the de-Broglie wavelength, the electronic and optical properties become highly tunable. Then the quantum size effects challenge the observation of many-particle phenomena such as the magneto-optical absorption, Raman intensity, and electron energy-loss spectrum. An exact analytical solution of the problem leads us to infer that these many-particle phenomena are, in fact, dictated by the generalized Kohn’s theorem (GKT) in the long-wavelength limit. Maneuvering the confinement and/or the magnetic field furnishes the resonance energies capable of being explored with the FIR, Raman, and/or electron-energy-loss spectroscopy. This implies that either of these probes is competent in observing the localized magnetoplasmons in the system. As an application of the rigorous analytical diagnosis of the system, we have presented various pertinent single-particle, such as Fock-Darwin spectrum, Fermi energy, zigzag excitation spectrum, and magneto-optical transitions, and the many-particle phenomena, such as magneto-optical absorption, Raman intensity, and electron energy-loss probability. In the latter, the energy position of the resonance peaks is observed to be independent of the electron-electron interactions and hence of the number of electrons in the quantum dot in compliance with the GKT. It is found that both confinement potential and magnetic field play a decisive role in influencing the aforementioned many-particle phenomena. Specifically, increasing (decreasing) the strength of the confining potential is found to be analogous to shrinking (expanding) the size of the quantum dots and results into a blue (red) shift in the respective spectra. Intensifying the magnetic field has two-fold effects in the resonance spectra associated with all three spectroscopies considered here: the lower resonance peak observes a red shift, whereas the higher one experiences a blue shift. This is a unique and intriguing behavior observed in the quantum dots with complete confinement. A deeper insight into the physics of the quantum dots is paving the way for their implementation in such diverse fields as quantum computing and medical imaging.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2023

Theory of an Automatic Seepage Meter and Ramifications for Applications

A new approach for measuring fluxes across surface water—groundwater interfaces was recently proposed. The Automatic Seepage Meter (ASM) is equipped with a precise water level sensor and digital memory that analyzes water level time series in a vertical tube inserted into a streambed. The ability to infer flux values with high temporal resolution relies on an accurate interpretation of water level dynamics inside the tube. Here, we reduce the three-dimensional hydrodynamic problem that describes the ASM water level in a variety of field conditions to a single ordinary differential equation. This novel general analytical solution for estimating ASM responses is more comprehensive and flexible than previous approaches and is applicable to the entire range of field conditions, including steady or transient stream stages, evaporation, rainfall, and noise. For example, our analysis determines the timing of the nonmonotonic ASM response to a monotonic linear stream stage variation and explains previously used empirical parabolic approximation for estimating fluxes. We present algorithms for simultaneous inference of vertical interface flux and hydraulic conductivity values together with an example code. We quantify how the accuracy of parameter estimation depends on test duration and noise amplitude and propose how our analysis can be used to optimize field test protocols. On this basis, changing the ASM geometry by increasing the radius and decreasing tube insertion depth may enable ASM field test protocols that estimate interface flux and hydraulic conductivity faster while maintaining desired accuracy. Potential applications of joint parameter estimation are suggested.

47 OTHER INSTRUMENTATION↗

NASA NTRS1984

Tidal disruption and the origin of the Moon

The dynamic problem of the tidal disruption of a rocky planetismal was solved by a direct integration of the fully three-dimensional, nonlinear equations of motion. The hypothesis that any object that passes within the Roche limit is disrupted was disproven. A time dependent solution was performed numerically, treating the planetismal as a fluid with a Murnaghan equation of state in the solid regions and zero pressure otherwise. Calculations show that a rocky body which passes by the Earth on a parabolic orbit with a perigee within the Roche limit is not tidally disrupted. Objects on hyperbolic orbits would experience even less tidal disruption. The results herein do not apply to bodies with very low viscosity. It is shown, however, that tidal disruption can be ruled out as a mechanism for reducing planetismal masses. Mechanisms for forming the Moon which rely upon tidal disruption are unlikely to be correct.

Mizuno, H.↗

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