Search NASA⌕ Search

SEARCH · Search NASA

Results for “numerical models”

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.

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

At least 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
NASA NTRS2009

Three-Dimensional Numerical Modeling of Magnetohydrodynamic Augmented Propulsion Experiment

Over the past several years, NASA Marshall Space Flight Center has engaged in the design and development of an experimental research facility to investigate the use of diagonalized crossed-field magnetohydrodynamic (MHD) accelerators as a possible thrust augmentation device for thermal propulsion systems. In support of this effort, a three-dimensional numerical MHD model has been developed for the purpose of analyzing and optimizing accelerator performance and to aid in understanding critical underlying physical processes and nonideal effects. This Technical Memorandum fully summarizes model development efforts and presents the results of pretest performance optimization analyses. These results indicate that the MHD accelerator should utilize a 45deg diagonalization angle with the applied current evenly distributed over the first five inlet electrode pairs. When powered at 100 A, this configuration is expected to yield a 50% global efficiency with an 80% increase in axial velocity and a 50% increase in centerline total pressure.

Turner, M. W.↗
DOE OSTI2025

Microstructure and CPO evolution of dynamically recrystallized olivine during complex deformation conditions: a full-field numerical modeling approach

The rheological properties of mantle rocks are strongly dependent on their crystallographic preferred orientation (CPO). Olivine CPO, defined by the orientation of seismically fast [100] axes parallel to flow direction, is also thought to be a dominant contributor to seismic anisotropy in the Earth's upper mantle. However, the amount of deformation needed to overprint a new CPO on a pre-existing fabric and the impact of the inherited CPOs on the transient microstructure evolution, remain unknown. This study employs a full-field numerical approach (VPFFT-ELLE) to explore the dynamic recrystallization and microstructural evolution of olivine polycrystalline aggregates under complex deformation conditions. We test four combinations of successive pure shear and simple shear boundary conditions. Findings indicate that inherited CPOs influence subsequent deformation in a manner dependent on the kinematic relationship between successive stages. In all cases, a minor strain increment (ε ∼0.3–0.6) is sufficient to erase the previous microstructure and CPO. However, when deformation conditions change dramatically (e.g., stretching direction changes orthogonally), the intensity of the new CPO developed is significantly lower and strain distribution are specially altered. During a transient strain stage, pre-existing microstructures undergo extensive reworking, especially when deformation conditions are changed dramatically, such as switching from simple shear to a pure shear condition with a parallel shortening direction relative to the stretching direction. We estimate the significance of these results in interpreting observations of seismic velocity anisotropy, concluding that P-wave seismic anisotropy is significantly and positively correlated with the evolution of olivine CPO with deformation history. This research underscores the transient nature of microstructural rearrangement in olivine aggregates and the necessity for caution in interpreting seismic anisotropy in regions with complex deformation histories, as inherited CPOs can influence current fabric development and induce deviation to the present deformation conditions.

36 MATERIALS SCIENCE↗
NASA NTRS1985

Advances in parameterization of cumulus convection in numerical models

Observational studies over the last decade have revealed new sources of difficulty in the already complex problem of parameterization of precipitating cumulus convection in large scale numerical weather prediction models. The structure of mesoscale updrafts and downdrafts and the difficulties they present for cumulus parameterization is reviewed. An approach is proposed which follows Kuo (1974) but includes cumulus and mesoscale downdrafts in the vertical distributions of heating and moistening. The development of a realistic life cycle behavior in the numerical prediction of a mesoscale convective complex will be shown.

Molinari, J.↗
NASA NTRS1992

A numerical model including PID control of a multizone crystal growth furnace

This paper presents a 2D axisymmetric combined conduction and radiation model of a multizone crystal growth furnace. The model is based on a programmable multizone furnace (PMZF) designed and built at NASA Lewis Research Center for growing high quality semiconductor crystals. A novel feature of this model is a control algorithm which automatically adjusts the power in any number of independently controlled heaters to establish the desired crystal temperatures in the furnace model. The control algorithm eliminates the need for numerous trial and error runs previously required to obtain the same results. The finite element code, FIDAP, used to develop the furnace model, was modified to directly incorporate the control algorithm. This algorithm, which presently uses PID control, and the associated heat transfer model are briefly discussed. Together, they have been used to predict the heater power distributions for a variety of furnace configurations and desired temperature profiles. Examples are included to demonstrate the effectiveness of the PID controlled model in establishing isothermal, Bridgman, and other complicated temperature profies in the sample. Finally, an example is given to show how the algorithm can be used to change the desired profile with time according to a prescribed temperature-time evolution.

Panzarella, Charles H.↗
NASA NTRS2019

Air Pollution Inputs to the Mojave Desert By Fusing Surface Mobile and Airborne in Situ and Airborne and Satellite Remote Sensing: A Case Study of Interbasin Transport With Numerical Model Validation

Deserts are fragile and highly sensitive ecosystems that increasingly are affected by upwind urban areas and industrial activities. The Los Angeles Basin (LAB) contributes to poor air quality in downwind deserts including the Mojave Desert. Few studies have investigated potential air pollution inputs to the Mojave, whose fragile ecosystem includes endangered plant and animal species. Data were collected on 19 August 2015 by a mobile air quality laboratory, AMOG (AutoMObile trace Gas) Surveyor, that observed inputs can arise from the LAB as well as the San Joaquin Valley (SJV), California. The campaign used a strong methane (CH4) plume as a tracer for the downwind fate of emissions from Bakersfield area petroleum production and also measured ozone (O3). Additional in situ concurrent airborne GHG and O3 data were collected by AJAX - Alpha Jet Atmospheric eXperiment. Both AMOG and AJAX measure winds. Mojave Desert air quality was very poor (visibility ~4 km). Based on the winds, an additional source was inferred beyond the LAB and SJV Basins. Numerical transport modeling and analysis of aerosol lidar data collected the same day by the Cloud Profiling LiDAR onboard the Earth Research-2 stratospheric airplane demonstrated that fires in Northern California were responsible, with prevailing winds transporting air southwards along the eastern Sierra Nevada Range (Bishop Valley) to the Mojave. Whereas the southern and eastern Mojave are impacted by SJV and LAB outflow, the north Mojave generally avoids these inputs. This study shows it can be affected by even distant wildfires, which likely will increase in occurrence and intensity from climate change. Thus, regulatory efforts to reduce air quality impacts on the endangered Mojave ecosystem must include wildfires and also account for the significant differences between different regions of the Mojave. Currently, there is a paucity of studies, highlighting the critical need for field research.

Ira Leifer↗
DOE OSTI2026

Numerical Modeling & Size Optimization of Thermal Energy Storage for Iron & Steel Production

Iron and steel production are responsible for 90 million MtCO2 per year in the United States. Hydrogen direct reduction of iron (H2DRI) is a promising pathway for a more sustainable iron production than commercially deployed technologies which rely on natural gas. The H2DRI process requires hydrogen at a temperature of up to 950 degrees C fed into a reduction furnace to produce pellets or briquettes that are used in the downstream iron and steelmaking process. In this work, we propose to use an electrical thermal energy storage (ETES) system, that can use renewable electricity to store high-temperature heat and dispatch it upon demand. Such a system can buffer the H2DRI plant from the variability of electricity prices by charging during curtailment and running the plant from storage during times of peak electricity price. We have developed heat transfer models for two different ETES systems that can be used to heat up hydrogen to the required temperatures: a particle-based ETES and a firebrick ETES. These models are used to evaluate the performance of such a system and support the sizing and preliminary cost estimation. The preliminary results using both models show that designing ETES systems for an industrial-scale H2DRI furnace is feasible. The firebrick ETES system has limited operational duration, which might limit the price buffering effect unless significantly oversized. The particle ETES system heat exchanger has industry-feasible dimensions, but its storage capacity would be decided upon the number of particle storage silos.

25 ENERGY STORAGE↗
NASA NTRS1992

Numerical modeling of fluid and electromagnetic phenomena in an arcjet

An explicit numerical technique is used to solve the axisymmetric reduced electromagnetic field equation. The effect of an electrical arc on a viscous, axisymmetric flow is approximated using an implicit thin layer Navier-Stokes solver with additional electromagnetic source terms in conjunction with the explicit finite difference code.

Flowe, Anita C.↗
NASA NTRS1991

Two-dimensional numerical modeling and solution of convection heat transfer in turbulent He II

Numerical schemes are employed to investigate heat transfer in the turbulent flow of He II. FEM is used to solve a set of equations governing the heat transfer and hydrodynamics of He II in the turbulent regime. Numerical results are compared with available experimental data and interpreted in terms of conventional heat transfer parameters such as the Prandtl number, the Peclet number, and the Nusselt number. Within the prescribed Reynolds number domain, the Gorter-Mellink thermal counterflow mechanism becomes less significant, and He II acts like an ordinary fluid. The convection heat transfer characteristics of He II in the highly turbulent regime can be successfully described by using the conventional turbulence and heat transfer theories.

Zhang, Burt X.↗
DOE OSTI2025

Sensitivity analysis of numerical modeling input parameters on floating offshore wind turbine loads in extreme idling conditions

Abstract. Floating offshore wind turbine (FOWT) systems are subject to complex environmental loads, with significant potential for damage in extreme storm conditions. Design simulations in these conditions are required to assess the survivability of the device with some level of confidence. Aero-hydro-servo-elastic engineering tools can be used with a reasonable balance of accuracy and computational efficiency. The models require many input parameters to describe the air and water conditions, the system properties, and the load calculations. Each of these parameters has some possible range, due to either statistical uncertainty or variations with time. Variation in the input parameters can have important effects on the uncertainty in the resulting loads, but it is not practical to perform detailed assessments of the impact of this uncertainty for every input parameter. This work demonstrates a method to identify the input parameters that have the most impact on the loads to focus further inspection. The process is done specifically for extreme storm load cases defined in the International Electrotechnical Commission design requirements for floating offshore wind turbines. The analysis was performed using the International Energy Agency Wind 15 MW offshore reference wind turbine atop the University of Maine VolturnUS-S reference platform in two US offshore wind regions, the Gulf of Maine and Humboldt Bay. It was found that the direction of incident waves and current, yaw misalignment, and the length of mooring line sections were among the primary sensitivities.

17 WIND ENERGY↗
DOE OSTI2024

Efficient general method for numerically modeling laser pulse propagation, overlap, and lifetime effects in amplifiers

An efficient numerical time-dependent general method is developed to address incoherent pulse overlap and lifetime effects in laser amplifiers. The alternating propagation-population laser energetics method (APPLE) has been validated against a semi-discrete coupled rate equation numerical method (SDRE) and analytic formalisms in bounding cases. APPLE is based on decoupled rates applied to a time-dependent framework where both space-time-dependent populations and pulse energetics are consistently updated in each time step. A significant advantage of APPLE lies in its conceptual simplicity, ease of implementation, and relatively small computational cost. SDRE tracks the populations through coupled rates and uses the method of lines to discretize the hyperbolic partial differential transport equations allowing for use of ordinary differential equation solvers. With reasonably sized mesh, we report both energetic and power pulse shape relative differences on the order of one percent between the models over a large range of initial conditions.

47 OTHER INSTRUMENTATION↗
NASA NTRS2017

Numerical Modeling of Saturated Boiling in a Heated Tube

This paper describes a mathematical formulation and numerical solution of boiling in a heated tube. The mathematical formulation involves a discretization of the tube into a flow network consisting of fluid nodes and branches and a thermal network consisting of solid nodes and conductors. In the fluid network, the mass, momentum and energy conservation equations are solved and in the thermal network, the energy conservation equation of solids is solved. A pressure-based, finite-volume formulation has been used to solve the equations in the fluid network. The system of equations is solved by a hybrid numerical scheme which solves the mass and momentum conservation equations by a simultaneous Newton-Raphson method and the energy conservation equation by a successive substitution method. The fluid network and thermal network are coupled through heat transfer between the solid and fluid nodes which is computed by Chen's correlation of saturated boiling heat transfer. The computer model is developed using the Generalized Fluid System Simulation Program and the numerical predictions are compared with test data.

Majumdar, Alok↗
NASA NTRS1988

Numerical models of thermally and mechanically coupled two-layer convection of highly viscous fluids

Thermal convection in the earth mantle is investigated by means of numerical simulations. The mantle models comprise two horizontal layers of viscous incompressible fluid (with identical or differing properties) separated by a fixed horizontal interface; an isothermal, fixed-heat-flux, or insulating lower boundary; and heat supplied either from below or internally (in the lower layer only, equally in both layers, or primarily in the upper layer). The mathematical formulation of the models is explained, and results for linear stability and finite-amplitude convection are presented in extensive tables and graphs and discussed in detail. Particular attention is given to the presence of thermal coupling without interface distortion in many cases, the predominance of long-wavelength cells in the finite-amplitude models, and the large (50 percent) temperature difference across the interface in all cases.

Ellsworth, Kirk↗
NASA NTRS2011

Numerical Modeling of Cometary Meteoroid Streams Encountering Mars and Venus

We have simulated numerically the existence of meteoroid streams that encounter the orbits of Mars and Venus, potentially producing meteor showers at those planets. We find that 17 known comets can produce such showers, the intensity of which can be determined through observations. Six of these streams contain dense dust trails capable of producing meteor outbursts.

Christou, A. A.↗
DOE OSTI2024

Numerical modeling of impurity powder injection in W7-X

In this work, we present numerical simulation results of powder injection in W7-X using the EMC3-EIRENE and DIS codes. First, we model powder injection experiments performed in W7-X with the Probe Mounted Powder Injector. The simulation results qualitatively agree with visible imaging measurements. Secondly, we perform predictive simulations to guide the installation of an Impurity Powder Dropper in W7-X, allowing to choose in between several available non-vertical ports to maximize the amount of powder penetrating into the plasma, as well as the verticality of the port, to minimize sticking of the powders in the in-vessel stainless steel guiding tube. Port AEM41 is selected as the best candidate for IPD installation. The robustness of the simulation results has been verified for different plasma densities, powder materials and sizes, powder friction coefficient and changes in the plasma flow.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1981

Regional eddy vorticity transport and the equilibrium vorticity budgets of a numerical model ocean circulation

A mean vorticity budget analysis is presented of Holland's (1978) numerical ocean general circulation experiment. The stable budgets are compared with classical circulation theory to emphasize the ways in which the mesoscale motions of the model alter (or leave unaltered) classical vorticity balances. The basinwide meridional transports of vorticity by the mean flow and by the mesoscale flow in the mean are evaluated to establish the role(s) of the mesoscale in the larger scale equilibrium vorticity transports. The vorticity equation for this model fluid system is presented and the budget analysis method is described. Vorticity budgets over the selected regions and on a larger scale are given, and a summary of budget results is provided along with remarks about the utility of this type of analysis.

Harrison, D. E.↗
NASA NTRS1977

Comparisons among second-order, fourth-order, and pseudospectral advection techniques in a mesoscale numerical model

The reported investigation is concerned with mesoscale squall-line simulation utilizing numerical techniques. Particular attention is given to fourth-order results and their implications for the dynamical evolution of mesoscale squall-line systems which contain severe local storms. An adiabatic inviscid set of prognostic equations is utilized in a z coordinate system for the fundamental experiment. The fourth-order advection is employed for all time-dependent equations. The Euler-backward time marching scheme is utilized with a 60 second time step. The horizontal mesh length is 42 km. Horizontal diffusion is accomplished by utilizing a smoother-desmoother. Lateral boundary conditions are designed to reduce the development of strong gradients in dependent variables near the boundaries by bringing interior values from the grid to the boundaries.

Kaplan, M.↗
NASA NTRS1988

Numerical modeling of chemical vapor deposition (CVD) in a horizontal reactor

In the present numerical prediction of the deposition rate of silicon from silane in a CVD process, the conservation equations for mass, momentum, energy, and chemical species are solved on a staggered grid using the SIMPLE algorithm, while the rate of chemical reactions in the gas phase and on the susceptor surface is obtained from an Arrhenius rate equation. Predicted deposition rates as a function of position along the susceptor with and without the gas phase chemical reaction are compared with the available experimental and numerical data; agreement is excellent except at the leading edge of the susceptor, where the deposition rate is overpredicted.

Sheikholeslami, M. Z.↗
NASA NTRS1999

Numerical Modeling of Three-Dimensional Fluid Flow with Phase Change

We present a numerical method to compute phase change dynamics of three-dimensional deformable bubbles. The full Navier-Stokes and energy equations are solved for both phases by a front tracking/finite difference technique. The fluid boundary is explicitly tracked by discrete points that are connected by triangular elements to form a front that is used to keep the stratification of material properties sharp and to calculate the interfacial source terms. Two simulations are presented to show robustness of the method in handling complex phase boundaries. In the first case, growth of a vapor bubble in zero gravity is studied where large volume increase of the bubble is managed by adaptively increasing the front resolution. In the second case, growth of a bubble under high gravity is studied where indentation at the rear of the bubble results in a region of large curvature which challenges the front tracking in three dimensions.

Esmaeeli, Asghar↗