SEARCH · Engineering Papers

Results for “implicit methods”

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 361 records · Page 20

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Impact of the numerical solution approach of a plant hydrodynamic model (v0.1) on vegetation dynamics

Abstract. Numerous plant hydrodynamic models have started to be implemented in vegetation dynamics models, reflecting the central role of plant hydraulic traits in driving water, energy, and carbon cycles, as well as plant adaptation to climate change. Different numerical approximations of the governing equations of the hydrodynamic models have been documented, but the numerical accuracy of these models and its subsequent effects on the simulated vegetation function and dynamics have rarely been evaluated. Using different numerical solution methods (including implicit and explicit approaches) and vertical discrete grid resolutions, we evaluated the numerical performance of a plant hydrodynamic module in the Functionally Assembled Terrestrial Ecosystem Simulator (FATES-HYDRO version 0.1) based on single-point and global simulations. Our simulation results showed that when near-surface vertical grid spacing is coarsened (grid size >10 cm), the model significantly overestimates aboveground biomass (AGB) in most of the temperate forest locations and underestimates AGB in the boreal forest locations, as compared to a simulation with finer vertical grid spacing. Grid coarsening has a small effect on AGB in the tropical zones of Asia and South America. In particular, coarse surface grid resolution should not be used when there are large and prolonged water content differences among soil layers at depths due to long dry-season duration and/or well-drained soil or when soil evaporation is a dominant fraction of evapotranspiration. Similarly, coarse surface grid resolution should not be used when there is lithologic discontinuity along the soil depth. This information is useful for uncertainty quantification, sensitivity analysis, or the training of surrogate models to design the simulations when computational cost limits the use of ensemble simulations.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1971

Finite difference solution to the three-dimensional, incompressible thermal energy boundary-layer equation

An implicit numerical method has been adopted for the solution to the three-dimensional, energy boundary-layer equation. The energy equation is written in terms of a dimensionless temperature function, the relative stagnation-enthalpy difference, and transformed by the introduction of a Blasius-type transformation of coordinates as well as dimensionless stream functions. The method is applied to the problem consisting of an infinite cylinder joined with its axis perpendicular to a thin, flat, heated plate. A Prandtl number equal to one is simply considered.

Fillo, J. A.↗

NASA NTRS1975

Viscous flow solutions with a cubic spline approximation

A cubic spline approximation is used for the solution of several problems in fluid mechanics. This procedure provides a high degree of accuracy even with a nonuniform mesh, and leads to a more accurate treatment of derivative boundary conditions. The truncation errors and stability limitations of several typical integration schemes are presented. For two-dimensional flows a spline-alternating-direction-implicit (SADI) method is evaluated. The spline procedure is assessed and results are presented for the one-dimensional nonlinear Burgers' equation, as well as the two-dimensional diffusion equation and the vorticity-stream function system describing the viscous flow in a driven cavity. Comparisons are made with analytic solutions for the first two problems and with finite-difference calculations for the cavity flow.

Rubin, S. G.↗

NASA NTRS1975

Numerical study of supersonic turbulent flow over small protuberances

Supersonic turbulent boundary layers over two-dimensional protuberances are investigated, using the numerical finite difference alternating direction implicit (ADI) method. The turbulence is modeled mathematically. The turbulence is represented here by the eddy viscosity approach. The turbulent boundary layer structure as well as an interest in thick boundary layers and much larger protuberance heights than in the laminar case lead to new difficulties. The problems encountered and the means to remove them are discussed.

Polak, A.↗

NASA NTRS1976

Condensation of water vapor in rarefaction waves. I - Homogeneous nucleation

A detailed theoretical investigation has been made of the condensation of water vapor/carrier gas mixtures in the nonstationary rarefaction wave generated in a shock tube. It is assumed that condensation takes place by homogeneous nucleation. The equations of motion together with the nucleation rate and the droplet growth equations were solved numerically by the method of characteristics and Lax's method of implicit artificial viscosity. It is found that, for the case considered, the condensation wave formed by the collapse of the metastable nonequilibrium state is followed by a shock wave generated by the intersection of characteristics of the same family. The expansion is practically isentropic up to the onset of condensation. The condensation front accelerates in the x,t plane. The results of the computations for a chosen case of water vapor/nitrogen mixture are presented by plotting variations of pressure, nucleation rate, number density of critical clusters, and condensate mass-fraction along three particle paths. Some consideration is given to homogeneous condensation experiments conducted in a shock tube. Although a direct comparison of the present theoretical work and these experiments is not possible, several worthwhile interpretative features have resulted nevertheless.

Sislian, J. P.↗

NASA NTRS1978

Large eddy simulation of incompressible turbulent channel flow

The three-dimensional, time-dependent primitive equations of motion were numerically integrated for the case of turbulent channel flow. A partially implicit numerical method was developed. An important feature of this scheme is that the equation of continuity is solved directly. The residual field motions were simulated through an eddy viscosity model, while the large-scale field was obtained directly from the solution of the governing equations. An important portion of the initial velocity field was obtained from the solution of the linearized Navier-Stokes equations. The pseudospectral method was used for numerical differentiation in the horizontal directions, and second-order finite-difference schemes were used in the direction normal to the walls. The large eddy simulation technique is capable of reproducing some of the important features of wall-bounded turbulent flows. The resolvable portions of the root-mean square wall pressure fluctuations, pressure velocity-gradient correlations, and velocity pressure-gradient correlations are documented.

Moin, P.↗

NASA NTRS1982

A composite velocity procedure for the compressible Navier-Stokes equations

A new boundary-layer relaxation procedure is presented. In the spirit of the theory of matched asymptotic expansions, a multiplicative composite of the appropriate velocity representations for the inviscid and viscous regions is prescribed. The resulting equations are structured so that far from the surface of the body the momentum equations lead to the Bernoulli relation for the pressure, while the continuity equation reduces to the familiar compressible potential equation. Close to the body surface, the governing equations and solution techniques are characteristic of those describing interacting boundary-layers; although, the full Navier-Stokes equations are considered here. Laminar flow calculations for the subsonic flow over an axisymmetric boattail simulator geometry are presented for a variety of Reynolds and Mach numbers. A strongly implicit solution method is applied for the coupled velocity components.

Khosla, P. K.↗

NASA NTRS1982

Numerical simulation of the flow and fuel-air mixing in an axisymmetric piston-cylinder arrangement

The implicit factored method of Beam and Warming was employed to describe the flow and the fuel-air mixing in an axisymmetric piston-cylinder configuration during the intake and compression strokes. The governing equations were established on the basis of laminar flow. The increased mixing due to turbulence was simulated by appropriately chosen effective transport properties. Calculations were performed for single-component gases and for two-component gases and for two-component gas mixtures. The flow field was calculated as functions of time and position for different geometries, piston speeds, intake-charge-to-residual-gas-pressure ratios, and species mass fractions of the intake charge. Results are presented in graphical form which show the formation, growth, and break-up of those vortices which form during the intake stroke and the mixing of fuel and air throughout the intake and compression strokes. It is shown that at bore-to-stroke ratio of less than unity, the vortices may break-up during the intake stroke. It is also shown that vortices which do not break-up during the intake stroke coalesce during the compression stroke. The results generated were compared to existing numerical solutions and to available experimental data.

Shih, T. I. P.↗

NASA NTRS1982

Vortex motion in axisymmetric piston-cylinder configurations

By using the Beam and Warming implicit-factored method of solution of the Navier-Stokes equations, velocities were calculated inside axisymmetric piston cylinder configurations during the intake and compression strokes. Results are presented in graphical form which show the formation, growth and breakup of those vortices which form during the intake stroke by the jet issuing from the valve. It is shown that at bore-to-stroke ratio of less than unity, the vortices may breakup during the intake stroke. It is also shown that vortices which do not breakup during the intake stroke coalesce during the compression stroke.

Shih, T. I. P.↗

NASA NTRS1983

Unconditionally stable implicit-explicit algorithms for coupled thermal stress waves

An unconditionally stable implicit-explicit method is proposed for the analysis of transient coupled thermal stress waves, and the computer-implementation aspects of the method are discussed. In the method proposed here, the mechanical displacement and temperature are used as independent variables. The resulting coupled finite element matrix equations are symmetric.

Liu, W. K.↗

NASA NTRS1983

An interacting boundary layer model for cascades

A laminar, incompressible interacting boundary layer model is developed for two-dimensional cascades. In the limit of large cascade spacing these equations reduce to the interacting boundary layer equations for a single body immersed in an infinite stream. A fully implicit numerical method is used to solve the governing equations, and is found to be at least as efficient as the same technique applied to the single body problem. Solutions are then presented for a cascade of finite flat plates and a cascade of finite sine-waves, with cusped leading and trailing edges.

Davis, R. T.↗

NASA NTRS1983

XTRAN2L: A program for solving the general-frequency unsteady transonic small disturbance equation

A program, XTRAN2L, for solving the general-frequency unsteady transonic small disturbance potential equation was developed. It is a modification of the LTRAN2-NLR code. The alternating-direction-implicit (ADI) method of Rizzetta and Chin is used to advance solutions of the potential equation in time Engquist-Osher monotone spatial differencing is used in the ADI solution algorithm. As a result, the XTRAN2L code is more robust and more efficient than similar codes that use Murman-Cole type-dependent spatial differencing. Nonreflecting boundary conditions that are consistent with the general-frequency equation have been developed and implemented at the far-field boundaries. Use of those conditions allow the computational boundaries to be moved closer to the airfoil with no loss of accuracy. This makes the XTRAN2L code more economical to use.

Whitlow, W., Jr.↗

NASA NTRS1985

Development of a boundary-layer-type solver based on simultaneous iteration technique for axisymmetric separated flows

A boundary-layer-type solver is developed for the numerical solution of axisymmetric separated flows. A new fully implicit coupling scheme for the viscous and inviscid regions is demonstrated. This fully implicit coupling technique is similar to the work of Carter, Veldman, and an extension of an earlier work of Halim and Hafez. A comparison is made for the convergence rate using this new fully implicit coupling technique and the semiimplicit coupling of Halim and Hafez. Numerical results using the fully implicit coupling are obtained for laminar incompressible separated flows, including a boattail and a series of trough geometries. Also, the near-wake flow problem is considered using the present formulation. A clear conclusion of this investigation is that the present scheme using the fully implicit coupling method converges at a faster rate than the semiimplicit coupling and the partially parabolized Navier-Stokes (PPNS) procedures.

Halim, A. A. M.↗

NASA NTRS1985

Computing with high-resolution upwind schemes for hyperbolic equations

Computational aspects of modern high-resolution upwind finite-difference schemes for hyperbolic systems of conservation laws are examined. An operational unification is demonstrated for constructing a wide class of flux-difference-split and flux-split schemes based on the design principles underlying total variation diminishing (TVD) schemes. Consideration is also given to TVD scheme design by preprocessing, the extension of preprocessing and postprocessing approaches to general control volumes, the removal of expansion shocks and 'glitches', relaxation methods for implicit TVD schemes, and a new family of high-accuracy TVD schemes.

Chakravarthy, S. R.↗

NASA NTRS1985

Theoretical and experimental comparison of vapor cavitation in dynamically loaded journal bearings

Vapor cavitation for a submerged journal bearing under dynamically loaded conditions was investigated. The observation of vapor cavitation in the laboratory was done by high-speed photography. It was found that vapor cavitation occurs when the tensile stress applied to the oil exceeded the tensile strength of the oil or the binding of the oil to the surface. The theoretical solution to the Reynolds equation is determined numerically using a moving boundary algorithm. This algorithm conserves mass throughout the computational domain including the region of cavitation and its boundaries. An alternating direction implicit (MDI) method is used to effect the time march. A rotor undergoing circular whirl was studied. Predicted cavitation behavior was analyzed by three-dimensional computer graphic movies. The formation, growth, and collapse of the bubble in response to the dynamic conditions is shown. For the same conditions of dynamic loading, the cavitation bubble was studied in the laboratory using high-speed photography.

Brewe, D. E.↗

NASA NTRS1986

A two-dimensional numerical study of the flow inside the combustion chambers of a motored rotary engine

A numerical study was performed to investigate the unsteady, multidimensional flow inside the combustion chambers of an idealized, two-dimensional, rotary engine under motored conditions. The numerical study was based on the time-dependent, two-dimensional, density-weighted, ensemble-averaged conservation equations of mass, species, momentum, and total energy valid for two-component ideal gas mixtures. The ensemble-averaged conservation equations were closed by a K-epsilon model of turbulence. This K-epsilon model of turbulence was modified to account for some of the effects of compressibility, streamline curvature, low-Reynolds number, and preferential stress dissipation. Numerical solutions to the conservation equations were obtained by the highly efficient implicit-factored method of Beam and Warming. The grid system needed to obtain solutions were generated by an algebraic grid generation technique based on transfinite interpolation. Results of the numerical study are presented in graphical form illustrating the flow patterns during intake, compression, gaseous fuel injection, expansion, and exhaust.

Shih, T. I. P.↗

NASA NTRS1986

Fluid flow and fuel-air mixing in a motored two-dimensional Wankel rotary engine

The implicit-factored method of Beam and Warming was employed to obtain numerical solutions to the conservation equations of mass, species, momentum, and energy to study the unsteady, multidimensional flow and mixing of fuel and air inside the combustion chambers of a two-dimensional Wankel rotary engine under motored conditions. The effects of the following engine design and operating parameters on fluid flow and fuel-air mixing during the intake and compression cycles were studied: engine speed, angle of gaseous fuel injection during compression cycle, and speed of the fuel leaving fuel injector.

Shih, T. I.-P.↗

NASA NTRS1986

A two-dimensional numerical study of the flow inside the combustion chamber of a motored rotary engine

Shih, T. I-P.↗