SEARCH · Search NASA

Results for “Numerical methods and modeling”

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 307 records · Page 17

All sources NASA NTRS DOE OSTI

NASA NTRS1980

A dynamic analysis of the motion of a low-wing general aviation aircraft about its calculated equilibrium flat spin mode

The properties of the flat spin mode of a general aviation configuration have been studied through analysis of rotary balance data, numerical simulation, and analytical study of the equilibrium state. The equilibrium state is predicted well from rotary balance data. The variations of yawing moment and pitching moment as functions of sideslip have been shown to be of great importance in obtaining accurate modeling. These dependencies are not presently available with sufficient accuracy from previous tests or theories. The stability of the flat spin mode has been examined extensively using numerical linearization, classical perturbation methods, and reduced order modeling. The stability exhibited by the time histories and the eigenvalue analyses is shown to be strongly dependent on certain static cross derivatives and more so on the dynamic derivatives. Explicit stability criteria are obtained from the reduced order models.

Tischler, M. B.↗

DOE OSTI2025

Water Mass Transformation Budgets in Finite‐Volume Generalized Vertical Coordinate Ocean Models

Water Mass Transformation (WMT) theory provides conceptual tools that in principle enable innovative analyses of numerical ocean models; in practice, however, these methods can be challenging to implement and interpret, and therefore remain under-utilized. Our aim is to demonstrate the feasibility of diagnosing all terms in the water mass budget and to exemplify their usefulness for scientific inquiry and model development by quantitatively relating water mass changes, overturning circulations, boundary fluxes, and interior mixing. We begin with a pedagogical derivation of key results of classical WMT theory. We then describe best practices for diagnosing each of the water mass budget terms from the output of Finite-Volume Generalized Vertical Coordinate (FV-GVC) ocean models, including the identification of a non-negligible remainder term as the spurious numerical mixing due to advection scheme discretization errors. We illustrate key aspects of the methodology through the analysis of a polygonal region of the Greater Baltic Sea in a regional demonstration simulation using the Modular Ocean Model v6 (MOM6). We verify the convergence of our WMT diagnostics by brute-force, comparing time-averaged (“offline”) diagnostics on various vertical grids to timestep-averaged (“online”) diagnostics on the native model grid. Finally, we briefly describe a stack of xarray-enabled Python packages for evaluating WMT budgets in FV-GVC models (culminating in the new xwmb package), which is intended to be model-agnostic and available for community use and development.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1992

Models of convection-driven tectonic plates - A comparison of methods and results

Recent numerical studies of convection in the earth's mantle have included various features of plate tectonics. This paper describes three methods of modeling plates: through material properties, through force balance, and through a thin power-law sheet approximation. The results obtained are compared using each method on a series of simple calculations. From these results, scaling relations between the different parameterizations are developed. While each method produces different degrees of deformation within the surface plate, the surface heat flux and average plate velocity agree to within a few percent. The main results are not dependent upon the plate modeling method and herefore are representative of the physical system modeled.

King, Scott D.↗

NASA NTRS1996

Computation of Static Shapes and Voltages for Micromachined Deformable Mirrors with Nonlinear Electrostatic Actuators

In modeling micromachined deformable mirrors with electrostatic actuators whose gap spacings are of the same order of magnitude as those of the surface deformations, it is necessary to use nonlinear models for the actuators. In this paper, we consider micromachined deformable mirrors modeled by a membrane or plate equation with nonlinear electrostatic actuator characteristics. Numerical methods for computing the mirror deformation due to given actuator voltages and the actuator voltages required for producing the desired deformations at the actuator locations are presented. The application of the proposed methods to circular deformable mirrors whose surfaces are modeled by elastic membranes is discussed in detail. Numerical results are obtained for a typical circular micromachined mirror with electrostatic actuators.

Wang, P. K. C.↗

NASA NTRS1995

Applications of direct numerical simulation of turbulence in second order closures

This paper discusses two methods of developing models for the rapid pressure-strain correlation term in the Reynolds stress transport equation using direct numerical simulation (DNS) data. One is a perturbation about isotropic turbulence, the other is a perturbation about two-component turbulence -- an extremely anisotropic turbulence. A model based on the latter method is proposed and is found to be very promising when compared with DNS data and other models.

Shih, Tsan-Hsing↗

DOE OSTI2025

Machine Learned Empirical Numerical Integrator from Simulated Data

Recently, a number of state-of-the-art surrogate machine learning (ML) models have been designed for global weather and climate prediction, which have been trained using reanalysis data products. Reanalysis data products are constructed using numerical model simulations that combine numerical integration of partial differential equations and parameterization schemes. These products are typically only archived and made available using coarsened spatial and temporal resolutions. This study explores the impact of the numerical generation methods used to produce the training datasets and the temporal resolution of those datasets on machine learning surrogate models. Using the nonlinear vector autoregression (NVAR) machine as an explainable ML technique, simple dynamical systems are emulated with ML models trained on data produced by three classical numerical integration schemes. NVAR is validated as a skillful ML method, capable of producing accurate predictions and, more importantly, reconstructing both the underlying dynamics and the numerical integration scheme used to generate the training data. However, the machine fails to generalize predictions on unseen test data generated by different numerical integration schemes, despite the underlying dynamical system being the same. This result provides a word of caution for the growing field of machine learning emulation of weather and climate dynamics. Furthermore, we illustrate using NVAR that training on temporally coarsened data may increase the required complexity of ML models and potentially introduce new numerical challenges. Finally, we discover that empirical integration schemes with arbitrary time-stepping sizes can be constructed directly from the data, which implies a potential for the development of empirical numerical integration schemes.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2016

Implications of Projected Climate Change for Groundwater Recharge in the Western United States

Existing studies on the impacts of climate change on groundwater recharge are either global or basin/ location-specific. The global studies lack the specificity to inform decision making, while the local studies do little to clarify potential changes over large regions (major river basins, states, or groups of states), a scale often important in the development of water policy. An analysis of the potential impact of climate change on groundwater recharge across the western United States (west of 100 degrees longitude) is presented synthesizing existing studies and applying current knowledge of recharge processes and amounts. Eight representative aquifers located across the region were evaluated. For each aquifer published recharge budget components were converted into four standard recharge mechanisms: diffuse, focused, irrigation, and mountain-systems recharge. Future changes in individual recharge mechanisms and total recharge were then estimated for each aquifer. Model-based studies of projected climate-change effects on recharge were available and utilized for half of the aquifers. For the remainder, forecasted changes in temperature and precipitation were logically propagated through each recharge mechanism producing qualitative estimates of direction of changes in recharge only (not magnitude). Several key patterns emerge from the analysis. First, the available estimates indicate average declines of 10-20% in total recharge across the southern aquifers, but with a wide range of uncertainty that includes no change. Second, the northern set of aquifers will likely incur little change to slight increases in total recharge. Third, mountain system recharge is expected to decline across much of the region due to decreased snowpack, with that impact lessening with higher elevation and latitude. Factors contributing the greatest uncertainty in the estimates include: (1) limited studies quantitatively coupling climate projections to recharge estimation methods using detailed, process-based numerical models; (2) a generally poor understanding of hydrologic flowpaths and processes in mountain systems; (3) difficulty predicting the response of focused recharge to potential changes in the frequency and intensity of extreme precipitation events; and (4) unconstrained feedbacks between climate, irrigation practices, and recharge in highly developed aquifer systems.

longitude↗

NASA NTRS1975

A numerical method for solving the Navier-Stokes equations with application to shock-boundary layer interactions

A numerical method for solving the compressible form of the unsteady Navier-Stokes equations is described. This method was originally presented in 1970 and has since been modified during the development of computer programs at Ames for implementing models that account for the effects of turbulence in shock-induced separated flows. Although this paper does not describe the turbulence models themselves, a complete description of the basic numerical method is given with emphasis on the choice of a computational mesh for high Reynolds number flows, finite-difference approximations for mixed partial derivatives, extension of the Courant-Friedrichs-Lewy stability condition for viscous flows, mesh boundary conditions, and numerical smoothing for strong shock-wave calculations.

Maccormack, R. W.↗

NASA NTRS1978

Large-eddy simulation of a turbulent mixing layer

The three dimensional, time dependent (incompressible) vorticity equations were used to simulate numerically the decay of isotropic box turbulence and time developing mixing layers. The vorticity equations were spatially filtered to define the large scale turbulence field, and the subgrid scale turbulence was modeled. A general method was developed to show numerical conservation of momentum, vorticity, and energy. The terms that arise from filtering the equations were treated (for both periodic boundary conditions and no stress boundary conditions) in a fast and accurate way by using fast Fourier transforms. Use of vorticity as the principal variable is shown to produce results equivalent to those obtained by use of the primitive variable equations.

Mansour, N. N.↗

NASA NTRS1992

Assessment of higher order turbulence models for complex two- and three-dimensional flowfields

A numerical method is presented to solve the three-dimensional Navier-Stokes equations in combination with a full Reynolds-stress turbulence model. Computations will be shown for three complex flowfields. The results of the Reynolds-stress model will be compared with those predicted by two different versions of the k-omega model. It will be shown that an improved version of the k-omega model gives as accurate results as the Reynolds-stress model.

Menter, Florian R.↗

DOE OSTI2026

Closure models for the feedback of energetic particles on plasma turbulence

Energetic particles interact with the plasma surrounding them, resonating with certain types of plasma waves to stabilize them while destabilizing others, and changing the character of the background turbulence in ways that have not been fully quantified or understood. Interaction with the turbulent background plasma is key to the acceleration of many types of energetic particles including high-energy cosmic rays, solar energetic particles, and pick-up ions. The acceleration of particles is a process that would ideally be described by a kinetic model, a type of model that follows a probability distribution function (PDF) for all particles in 7-dimensional (x, y, z, v x , v y , v z , t) space. Because of the high dimensionality of a kinetic model, simulations that solve kinetic equations use the largest computational resources currently available, and are yet unable to simulate a realistic number of particles, reach the large scales necessary for astrophysical problems, and use high-precision numerical methods. Two available alternatives to kinetic plasma models have been explored for this problem, with limited success. One is a multi-fluid model produced by a cumulant discarding closure, which evolves coupled equations for the velocity, magnetic field, and internal energy for both the background plasma and the fluid of energetic particles. However, simulations that solve multi-fluid magnetohydrodynamic (MHD) equations are able to include the interaction with energetic particles only in crude ways, typically as an add-on pressure term. The second alternative is to use a hybrid method to couple a fluid description of the background plasma to a kinetic model or a Fokker–Planck model for the energetic particles. These methods are hampered by the physical modeling of the coupling. In this work, we develop a new model, which follows the PDF for all particles; this can be viewed as a step toward physical realism above a multi-fluid MHD model, while also being more computationally efficient than a kinetic model. The equations we develop model both the background plasma and the energetic particles self-consistently. Over the last decade, similar PDF methods have been developed to a high level of sophistication to model reactive flows and turbulent combustion for engineering applications. For treatment of the feedback of the energetic particles on a background plasma, a PDF closure approach should evaluate the mean characteristics, including the density, with better statistical quality than will particle-sampling procedures.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2003

A Discrete Constraint for Entropy Conservation and Sound Waves in Cloud-Resolving Modeling

Ideal cloud-resolving models contain little-accumulative errors. When their domain is so large that synoptic large-scale circulations are accommodated, they can be used for the simulation of the interaction between convective clouds and the large-scale circulations. This paper sets up a framework for the models, using moist entropy as a prognostic variable and employing conservative numerical schemes. The models possess no accumulative errors of thermodynamic variables when they comply with a discrete constraint on entropy conservation and sound waves. Alternatively speaking, the discrete constraint is related to the correct representation of the large-scale convergence and advection of moist entropy. Since air density is involved in entropy conservation and sound waves, the challenge is how to compute sound waves efficiently under the constraint. To address the challenge, a compensation method is introduced on the basis of a reference isothermal atmosphere whose governing equations are solved analytically. Stability analysis and numerical experiments show that the method allows the models to integrate efficiently with a large time step.

Zeng, Xi-Ping↗

NASA NTRS1976

Numerical simulation of turbulent jet noise, part 2

Results on the numerical simulation of jet flow fields were used to study the radiated sound field, and in addition, to extend and test the capabilities of the turbulent jet simulation codes. The principal result of the investigation was the computation of the radiated sound field from a turbulent jet. In addition, the computer codes were extended to account for the effects of compressibility and eddy viscosity, and the treatment of the nonlinear terms of the Navier-Stokes equations was modified so that they can be computed in a semi-implicit way. A summary of the flow model and a description of the numerical methods used for its solution are presented. Calculations of the radiated sound field are reported. In addition, the extensions that were made to the fundamental dynamical codes are described. Finally, the current state-of-the-art for computer simulation of turbulent jet noise is summarized.

Metcalfe, R. W.↗

NASA NTRS1978

On least squares approximations to indefinite problems of the mixed type

A least squares method is presented for computing approximate solutions of indefinite partial differential equations of the mixed type such as those that arise in connection with transonic flutter analysis. The method retains the advantages of finite difference schemes namely simplicity and sparsity of the resulting matrix system. However, it offers some great advantages over finite difference schemes. First, the method is insensitive to the value of the forcing frequency, i.e., the resulting matrix system is always symmetric and positive definite. As a result, iterative methods may be successfully employed to solve the matrix system, thus taking full advantage of the sparsity. Furthermore, the method is insensitive to the type of the partial differential equation, i.e., the computational algorithm is the same in elliptic and hyperbolic regions. In this work the method is formulated and numerical results for model problems are presented. Some theoretical aspects of least squares approximations are also discussed.

Fix, G. J.↗

NASA NTRS1980

Mountain winds (revisited)

The prediction of extremely high wind speeds, at ground level on the downstream side of a mountain range, is possible by solving the initial value problem for a two-layered nonlinear shallow water model of the atmosphere. Three different numerical methods are described to find the solutions which may involve shocks: (1) the vonNeumann-Richtmyer artificial viscosity method, (2) a filtering scheme, and (3) a hybrid method.

Isaacson, E.↗

NASA NTRS1981

Numerical studies of laminar and turbulent drag reduction

Two-dimensional incompressible flow over wavy surfaces is studied numerically by spectral methods. Turbulence effects are modeled. Results for symmetric and asymmetric wave forms are presented. Effect of propagating surface waves on drag reduction is studied. Comparisons between computer simulations and experimental results are made.

Balasubramanian, R.↗

NASA NTRS1982

Fluid flow in crystallization processes

Investigations of the interactions of heat, mass, and momentum transport in crystal growth from the melt are described. Studies of the mall-scale floating zone process and on a prototype of the vertical Bridgman growth system were emphasized. In both systems detailed numerical calculations are used to dissect the interplay between fluid convection and dopant segregation. These calculations are based on finite element techniques that make feasible the complete solution of solidification problems which include convection. Other numerical methods were developed for solving thermal models of crystal growth processes with melt/gas menisci and for simulating the microscale instabilities in solidification interfaces.

Brown, R. A.↗

NASA NTRS1986

Star counts and galactic structure

Star counts are discussed as an aid in studies of galactic structure. Iterative computations using a method-of-color-bins to match numerical models with observational data are described. The star counts provide bounds for the luminosity function and the color-magnitude diagram and allow spheroid normalization in a copycat Galaxy model which posits a galaxy with a Population I disk and a Population II spheroid. Good agreement has been obtained between star counts predicted by the model and observed distributions, including color and magnitude ranges and direction. Modifications being considered for the model include consideration of the halo, a galactic bulge and a thick disk. Several galactic structural features which have yet to be modeled and which may be tractable to the copycat model corrected with star counts are identified.

Bahcall, John N.↗