SEARCH · Search NASA

Results for “advection equation”

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 217 records · Page 12

All sources NASA NTRS DOE OSTI

NASA NTRS1984

A spectral collocation method for the Navier-Stokes equations

A Fourier-Chebyshev spectral method for the incompressible Navier-Stokes equations is described. It is applicable to a variety of problems including some with fluid properties which vary strongly both in the normal direction and in time. In this fully spectral algorithm, a preconditioned iterative technique is used for solving the implicit equations arising from semi-implicit treatment of pressure, mean advection and vertical diffusion terms. The algorithm is tested by applying it to hydrodynamic stability problems in channel flow and in external boundary layers with both constant and variable viscosity.

Malik, M. R.↗

NASA NTRS1985

A spectral collocation method for the Navier-Stokes equations

Malik, M. R.↗

NASA NTRS1991

Simplification of the Kalman filter for meteorological data assimilation

The paper proposes a new statistical method of data assimilation that is based on a simplification of the Kalman filter equations. The forecast error covariance evolution is approximated simply by advecting the mass-error covariance field, deriving the remaining covariances geostrophically, and accounting for external model-error forcing only at the end of each forecast cycle. This greatly reduces the cost of computation of the forecast error covariance. In simulations with a linear, one-dimensional shallow-water model and data generated artificially, the performance of the simplified filter is compared with that of the Kalman filter and the optimal interpolation (OI) method. The simplified filter produces analyses that are nearly optimal, and represents a significant improvement over OI.

Dee, Dick P.↗

NASA NTRS1992

The diagnosis of a preblocking explosively-developing extratropical cyclone system

This paper presents the diagnosis of an extratropical cyclone that developed explosively from 18 to 19 January 1979 over the North Atlantic Ocean. The diagnosis applies data obtained from the Goddard Laboratory for Atmospheres FGGE SOP-1 level III-b global analyses on a 4 deg latitude by 5 deg longitude grid to the extended height tendency and the Zwack-Okossi development equations. The cyclone developed initially in response to cyclonic vorticity advection downstream from an upper air through, warm air advection in a strongly baroclinic region, and latent heat release in the cyclone domain. As development continued, thermal advection and latent heat release increased their role in forcing height falls in the cyclone domain, while the influence of vorticity advection decreased. Finally, development ceased when anticyclonic vorticity advection below 700 mb and adiabatic cooling in the ascending air combined to neutralize the influence of warm air advection and latent heat release.

Uhl, Mary A.↗

NASA NTRS1992

Three-dimensional analytical model for the mixed layer depth

A three-dimensional (3D) semianalytical mixed-layer thermal inertia model for the mixed layer depth is developed, which includes the horizontal advection and diffusion processes modeled previously by Yan et al. (1991) in a one-dimensional (1D) model. The 3D model, which is based on the 3D thermal energy conservation equation, can be used as a supplement to the 1D model in the area where advection is important and for a finer resolution grid prediction. The model can be easily forced by remotely sensed data and is much simpler to compute than numerical 3D mixed-layer turbulence closure models. The new model was verified using data from the AVHRR/multichannel sea surface temperature data set and the Cooperative Ocean Atmosphere Data Set for the North Pacific Ocean, showing good agreement.

Yan, Xiao-Hai↗

NASA NTRS1997

A Three-Wave Model of the Stratosphere with Coupled Dynamics, Radiation and Photochemistry

A zonal mean chemistry transport model (2-D CTM) coupled with a semi-spectral dynamical model is used to simulate the distributions of trace gases in the present day atmosphere. The zonal-mean and eddy equations for the velocity and the geopotential height are solved in the semi-spectral dynamical model. The residual mean circulation is derived from these dynamical variables and used to advect the chemical species in the 2- D CTM. Based on a linearized wave transport equation, the eddy diffusion coefficients for chemical tracers are expressed in terms of the amplitude, frequency and growth rate of dynamical waves; local chemical loss rates; and a time constant parameterizing small scale mixing. The contributions to eddy flux are from the time varying wave amplitude (transient eddy), chemical reactions (chemical eddy) and small scale mixing. In spite of the high truncation in the dynamical module (only three longest waves are resolved), the model has simulated many observed characteristics of stratospheric dynamics and distribution of chemical species including ozone. Compared with the values commonly used in 2-D CTMs, the eddy diffusion coefficients for chemical species calculated in this model are smaller, especially in the subtropics. It is also found that the chemical eddy diffusion has only a small effects in determining the distribution of most slow species, including ozone in the stratosphere.

Shia, Run-Lie↗

NASA NTRS2021

Evidence that Horizontal Moisture Advection Regulates the Ubiquitous Amplification of Rainfall Variability over Tropical Oceans

The column moist static energy (MSE) budget equation approximates the processes associated with column moistening and drying in the tropics, and is therefore predictive of precipitation amplification and decay. We use ERA-Interim (ERA-I) and TRMM 3B42 data to investigate day-to-day convective variability and distinguish the roles of horizontal MSE (or moisture) advection versus vertical advection, sources, and sinks. Over tropical convergence zones, results suggest that horizontal moisture advection is a primary driver of day-to-day precipitation fluctuations; when drying via horizontal moisture advection is smaller (greater) than Chikira’s “column process,” precipitation tends to amplify (decay). In the absence of horizontal moisture advection, precipitation tends to increase spontaneously almost universally through a positive column process feedback. This bulk positive feedback is characterized by negative effective gross moist stability (GMS), which is maintained throughout the tropical convergence zones. How this positive feedback is achieved varies geographically, depending on the shape of vertical velocity (omega) profiles. In regions where omega profiles are top-heavy, the effective GMS is negative primarily owing to strong feedbacks between convection and diabatic MSE sources (radiative and surface fluxes). In these regions, vertical MSE advection stabilizes the atmosphere (positive vertical GMS). Where omega profiles are bottom-heavy, by contrast, a positive feedback is primarily driven by import of MSE through a shallow circulation (negative vertical GMS). The diabatic feedback and vertical GMS are in a seesaw balance, offsetting one another. Our results suggest that ubiquitous convective variability is amplified by the same mechanism as moisture-mode instability.

horizontal moisture advection↗

NASA NTRS1992

Simulations of surface winds at the Viking Lander sites using a one-level model

The one-level model developed by Mass and Dempsey for use in predicting surface flows in regions of complex terrain was adapted to simulate surface flows at the Viking lander sites on Mars. In the one-level model, prediction equations for surface winds and temperatures are formulated and solved. Surface temperatures change with time in response to diabatic heating, horizontal advection, adiabatic heating and cooling effects, and horizontal diffusion. Surface winds can change in response to horizontal advection, pressure gradient forces, Coriolis forces, surface drag, and horizontal diffusion. Surface pressures are determined by integration of the hydrostatic equation from the surface to some reference level. The model has successfully simulated surface flows under a variety of conditions in complex-terrain regions on Earth.

Bridger, Alison F. C.↗

NASA NTRS1991

Synoptic/planetary-scale interactions and blocking over the North Atlantic Ocean

One segment of work in the past year focused on the diagnosis of a major blocking anticyclone and its interacting synoptic scale circulations that occurred during January 1979 over the North Atlantic Ocean. Another segment focused on the diagnosis of a second explosive cyclone development that occurred over the southeastern United States during the time of block formation. The diagnoses were accomplished using the diagnostic relationship known as the Zwack-Okossi (Z-O) development equation. Results indicate that in both cyclone cases the development occurred as a result of the favorable influence of positive vorticity advection, warm air advection, and latent heat release and ceased when one or more of these influences diminished. The advantages of the Z-O equation are described.

Smith, Phillip J.↗

NASA NTRS1998

Complete Galilean-Invariant Lattice BGK Models for the Navier-Stokes Equation

Galilean invariance has been an important issue in lattice-based hydrodynamics models. Previous models concentrated on the nonlinear advection term. In this paper, we take into account the nonlinear response effect in a systematic way. Using the Chapman-Enskog expansion up to second order, complete Galilean invariant lattice BGK models in one dimension (theta = 3) and two dimensions (theta = 1) for the Navier-Stokes equation have been obtained.

Qian, Yue-Hong↗

NASA NTRS2013

An Assessment of the Skill of GEOS-5 Seasonal Forecasts

The seasonal forecast skill of the NASA Global Modeling and Assimilation Office coupled global climate model (CGCM) is evaluated based on an ensemble of 9-month lead forecasts for the period 1993 to 2010. The results from the current version (V2) of the CGCM consisting of the GEOS-5 AGM coupled to the MOM4 ocean model are compared with those from an earlier version (V1) in which the AGCM (the NSIPP model) was coupled to the Poseidon Ocean Model. It was found that the correlation skill of the Sea Surface Temperature (SST) forecasts is generally better in V2, especially over the sub-tropical and tropical central and eastern Pacific, Atlantic, and Indian Ocean. Furthermore, the improvement in skill in V2 mainly comes from better forecasts of the developing phase of ENSO from boreal spring to summer. The skill of ENSO forecasts initiated during the boreal winter season, however, shows no improvement in terms of correlation skill, and is in fact slightly worse in terms of root mean square error (RMSE). The degradation of skill is found to be due to an excessive ENSO amplitude. For V1, the ENSO amplitude is too strong in forecasts starting in boreal spring and summer, which causes large RMSE in the forecast. For V2, the ENSO amplitude is slightly stronger than that in observations and V1 for forecasts starting in boreal winter season. An analysis of the terms in the SST tendency equation, shows that this is mainly due to an excessive zonal advective feedback. In addition, V2 forecasts that are initiated during boreal winter season, exhibit a slower phase transition of El Nino, which is consistent with larger amplitude of ENSO after the ENSO peak season. It is found that this is due to weak discharge of equatorial Warm Water Volume (WWV). In both observations and V1, the discharge of equatorial WWV leads the equatorial geostrophic easterly current so as to damp the El Nino starting in January. This process is delayed by about 2 months in V2 due to the slower phase transition of the equatorial zonal current from westerly to easterly.

GMAO↗

NASA NTRS1999

Physically-Derived Dynamical Cores in Atmospheric General Circulation Models

The algorithm chosen to represent the advection in atmospheric models is often used as the primary attribute to classify the model. Meteorological models are generally classified as spectral or grid point, with the term grid point implying discretization using finite differences. These traditional approaches have a number of shortcomings that render them non-physical. That is, they provide approximate solutions to the conservation equations that do not obey the fundamental laws of physics. The most commonly discussed shortcomings are overshoots and undershoots which manifest themselves most overtly in the constituent continuity equation. For this reason many climate models have special algorithms to model water vapor advection. This talk focuses on the development of an atmospheric general circulation model which uses a consistent physically-based advection algorithm in all aspects of the model formulation. The shallow-water model is generalized to three dimensions and combined with the physics parameterizations of NCAR's Community Climate Model. The scientific motivation for the development is to increase the integrity of the underlying fluid dynamics so that the physics terms can be more effectively isolated, examined, and improved. The expected benefits of the new model are discussed and results from the initial integrations will be presented.

Rood, Richard B.↗

NASA NTRS1999

Physically-Derived Dynamical Cores in Atmospheric General Circulation Models

The algorithm chosen to represent the advection in atmospheric models is often used as the primary attribute to classify the model. Meteorological models are generally classified as spectral or grid point, with the term grid point implying discretization using finite differences. These traditional approaches have a number of shortcomings that render them non-physical. That is, they provide approximate solutions to the conservation equations that do not obey the fundamental laws of physics. The most commonly discussed shortcomings are overshoots and undershoots which manifest themselves most overtly in the constituent continuity equation. For this reason many climate models have special algorithms to model water vapor advection. This talk focuses on the development of an atmospheric general circulation model which uses a consistent physically-based advection algorithm in all aspects of the model formulation. The shallow-water model of Lin and Rood (QJRMS, 1997) is generalized to three dimensions and combined with the physics parameterizations of NCAR's Community Climate Model. The scientific motivation for the development is to increase the integrity of the underlying fluid dynamics so that the physics terms can be more effectively isolated, examined, and improved. The expected benefits of the new model are discussed and results from the initial integrations will be presented.

Rood, Richard B.↗

DOE OSTI2020

A buoyancy–shear–drag-based turbulence model for Rayleigh–Taylor, reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz mixing

A phenomenological turbulence model for Rayleigh–Taylor, reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz instability-induced mixing is developed using a general buoyancy–shear–drag model. Analytic solutions to the simplified model equations corresponding to each instability are derived separately, which are then used to calibrate the model coefficients to predict specific values of the mixing layer growth parameters and exponents. The buoyancy-shear–drag equations for the bubble and spike mixing layer widths are then solved numerically, and a turbulent diffusivity (or viscosity) is constructed dimensionally as the product of the mixing layer width, h, and its time-derivative, dh/dt. Surrogate turbulent fields are then constructed by multiplying a presumed, approximate self-similar spatial profile by appropriate functions of h and dh/dt. Using several simplifying approximations, the turbulent diffusion equations satisfied by the mean mass fraction and mean shear velocity are solved analytically in a reference frame moving with the mean advection velocity, and it is shown that the mean fields evolve in space and time as expected. The explicit modeling and solution of turbulent transport equations (e.g., for the turbulent kinetic energy and its dissipation rate or a lengthscale) are not required. By using separate equations for the bubble and spike evolution, the model naturally captures the asymmetry induced by an increasing density contrast between the heavy and light fluids. The model also includes molecular dissipation and diffusion, and therefore, can describe transition to fully-developed turbulence if the initial turbulent diffusivity/viscosity is much smaller than the molecular diffusivity/viscosity. This model is applied to constant-acceleration Rayleigh–Taylor, impulsively reshocked Richtmyer–Meshkov, and variable-density Kelvin–Helmholtz mixing layers to demonstrate its utility. It is shown that the numerical solutions of the model calibrated using specific values of the instability growth parameters and exponents: (1) produces mixing layer widths in agreement with the expected self-similar growth power laws; (2) gives spatiotemporal profiles of turbulent fields that are expected and consistent with previous results; (3) predicts the expected power-law growths and decays of the spatially-integrated turbulent fields, and; (4) gives spatiotemporal profiles of the mean fields that are expected and consistent with previous results. This relatively simple model can be implemented numerically, and is expected to be useful for fundamental as well as astrophysical and inertial confinement fusion applications.

42 ENGINEERING↗

NASA NTRS2004

Textbook Multigrid Efficiency for the Steady Euler Equations

A fast multigrid solver for the steady incompressible Euler equations is presented. Unlike time-marching schemes, this approach uses relaxation of the steady equations. Application of this method results in a discretization that correctly distinguishes between the advection and elliptic parts of the operator, allowing efficient smoothers to be constructed. Solvers for both unstructured triangular grids and structured quadrilateral grids have been written. Computations for channel flow and flow over a nonlifting airfoil have computed. Using Gauss-Seidel relaxation ordered in the flow direction, textbook multigrid convergence rates of nearly one order-of-magnitude residual reduction per multigrid cycle are achieved, independent of the grid spacing. This approach also may be applied to the compressible Euler equations and the incompressible Navier-Stokes equations.

Roberts, Thomas W.↗

DOE OSTI2023

An Open‐Source, Physics‐Based, Tropical Cyclone Downscaling Model With Intensity‐Dependent Steering

Abstract An open‐source, physics‐based tropical cyclone (TC) downscaling model is developed, in order to generate a large climatology of TCs. The model is composed of three primary components: (a) a random seeding process that determines genesis, (b) an intensity‐dependent beta‐advection model that determines the track, and (c) a non‐linear differential equation set that determines the intensification rate. The model is entirely forced by the large‐scale environment. Downscaling ERA5 reanalysis data shows that the model is generally able to reproduce observed TC climatology, such as the global seasonal cycle, genesis locations, track density, and lifetime maximum intensity distributions. Inter‐annual variability in TC count and power‐dissipation is also well captured, on both basin‐wide and global scales. Regional TC hazard estimated by this model is also analyzed using return period maps and curves. In particular, the model is able to reasonably capture the observed return period curves of landfall intensity in various sub‐basins around the globe. The incorporation of an intensity‐dependent steering flow is shown to lead to regionally dependent changes in power dissipation and return periods. Advantages and disadvantages of this model, compared to other downscaling models, are also discussed.

Meteorology & Atmospheric Sciences↗

NASA NTRS1991

Simulation of the present-day atmospheric ozone, odd nitrogen, chlorine and other species using a coupled 2-D model in isentropic coordinates

A zonally averaged two-dimensional model for chemical species in the atmosphere (Olaguer et al., 1990), with coupled dynamics, radiation transfer, and chemistry, is used to study the present-day atmospheric ozone, odd nitrogen, chlorine, other species. The model utilizes all zonally averaged physical equations of momentum, energy, and mass and determines self-consistently the advective and diffusive tranasport parameters from the temperature specific to the period of observation. A comparison of the model results with observations showed good agreement.

Yang, H.↗

NASA NTRS1998

Bulk-Flow Analysis of Hybrid Thrust Bearings for Advanced Cryogenic Turbopumps

A bulk-flow analysis and computer program for prediction of the static load performance and dynamic force coefficients of angled injection, orifice-compensated hydrostatic/hydrodynamic thrust bearings have been completed. The product of the research is an efficient computational tool for the design of high-speed thrust bearings for cryogenic fluid turbopumps. The study addresses the needs of a growing technology that requires of reliable fluid film bearings to provide the maximum operating life with optimum controllable rotordynamic characteristics at the lowest cost. The motion of a cryogenic fluid on the thin film lands of a thrust bearing is governed by a set of bulk-flow mass and momentum conservation and energy transport equations. Mass flow conservation and a simple model for momentum transport within the hydrostatic bearing recesses are also accounted for. The bulk-flow model includes flow turbulence with fluid inertia advection, Coriolis and centrifugal acceleration effects on the bearing recesses and film lands. The cryogenic fluid properties are obtained from realistic thermophysical equations of state. Turbulent bulk-flow shear parameters are based on Hirs' model with Moody's friction factor equations allowing a simple simulation for machined bearing surface roughness. A perturbation analysis leads to zeroth-order nonlinear equations governing the fluid flow for the thrust bearing operating at a static equilibrium position, and first-order linear equations describing the perturbed fluid flow for small amplitude shaft motions in the axial direction. Numerical solution to the zeroth-order flow field equations renders the bearing flow rate, thrust load, drag torque and power dissipation. Solution to the first-order equations determines the axial stiffness, damping and inertia force coefficients. The computational method uses well established algorithms and generic subprograms available from prior developments. The Fortran9O computer program hydrothrust runs on a Windows 95/NT personal computer. The program, help files and examples are licensed by Texas A&M University Technology License Office. The study of the static and dynamic performance of two hydrostatic/hydrodynamic bearings demonstrates the importance of centrifugal and advection fluid inertia effects for operation at high rotational speeds. The first example considers a conceptual hydrostatic thrust bearing for an advanced liquid hydrogen turbopump operating at 170,000 rpm. The large axial stiffness and damping coefficients of the bearing should provide accurate control and axial positioning of the turbopump and also allow for unshrouded impellers, therefore increasing the overall pump efficiency. The second bearing uses a refrigerant R134a, and its application in oil-free air conditioning compressors is of great technological importance and commercial value. The computed predictions reveal that the LH2 bearing load capacity and flow rate increase with the recess pressure (i.e. increasing orifice diameters). The bearing axial stiffness has a maximum for a recess pressure rati of approx. 0.55. while the axial damping coefficient decreases as the recess pressure ratio increases. The computer results from three flow models are compared. These models are a) inertialess, b) fluid inertia at recess edges only, and c) full fluid inertia at both recess edges and film lands. The full inertia model shows the lowest flow rates, axial load capacity and stiffness coefficient but on the other hand renders the largest damping coefficients and inertia coefficients. The most important findings are related to the reduction of the outflow through the inner radius and the appearance of subambient pressures. The performance of the refrigerant hybrid thrust bearing is evaluated at two operating speeds and pressure drops. The computed results are presented in dimensionless form to evidence consistent trends in the bearing performance characteristics. As the applied axial load increases, the bearing film thickness and flow rate decrease while the recess pressure increases. The axial stiffness coefficient shows a maximum for a certain intermediate load while the damping coefficient steadily increases. The computed results evidence the paramount of centrifugal fluid inertia at low recess pressures (i.e. low loads), and where there is actually an inflow through the bearing inner diameter, accompanied by subambient pressures just downstream of the bearing recess edge. These results are solely due to centrifugal fluid inertia and advection transport effects. Recommendations include the extension of the computer program to handle flexure pivot tilting pad hybrid bearings and the ability to calculate moment coefficients for shaft angular misalignments.

SanAndres, Luis↗