Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerics”

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 1,369 records · Page 76
All sourcesNASA NTRSDOE OSTI
NASA NTRS2010

Adaptive Numerical Algorithms in Space Weather Modeling

Space weather describes the various processes in the Sun-Earth system that present danger to human health and technology. The goal of space weather forecasting is to provide an opportunity to mitigate these negative effects. Physics-based space weather modeling is characterized by disparate temporal and spatial scales as well as by different physics in different domains. A multi-physics system can be modeled by a software framework comprising of several components. Each component corresponds to a physics domain, and each component is represented by one or more numerical models. The publicly available Space Weather Modeling Framework (SWMF) can execute and couple together several components distributed over a parallel machine in a flexible and efficient manner. The framework also allows resolving disparate spatial and temporal scales with independent spatial and temporal discretizations in the various models. Several of the computationally most expensive domains of the framework are modeled by the Block-Adaptive Tree Solar wind Roe Upwind Scheme (BATS-R-US) code that can solve various forms of the magnetohydrodynamics (MHD) equations, including Hall, semi-relativistic, multi-species and multi-fluid MHD, anisotropic pressure, radiative transport and heat conduction. Modeling disparate scales within BATS-R-US is achieved by a block-adaptive mesh both in Cartesian and generalized coordinates. Most recently we have created a new core for BATS-R-US: the Block-Adaptive Tree Library (BATL) that provides a general toolkit for creating, load balancing and message passing in a 1, 2 or 3 dimensional block-adaptive grid. We describe the algorithms of BATL and demonstrate its efficiency and scaling properties for various problems. BATS-R-US uses several time-integration schemes to address multiple time-scales: explicit time stepping with fixed or local time steps, partially steady-state evolution, point-implicit, semi-implicit, explicit/implicit, and fully implicit numerical schemes. Depending on the application, we find that different time stepping methods are optimal. Several of the time integration schemes exploit the block-based granularity of the grid structure. The framework and the adaptive algorithms enable physics based space weather modeling and even forecasting.

Toth, Gabor↗
NASA NTRS2010

Middle Atmosphere Dynamics with Gravity Wave Interactions in the Numerical Spectral Model: Tides and Planetary Waves

As Lindzen (1981) had shown, small-scale gravity waves (GW) produce the observed reversals of the zonal-mean circulation and temperature variations in the upper mesosphere. The waves also play a major role in modulating and amplifying the diurnal tides (DT) (e.g., Waltersheid, 1981; Fritts and Vincent, 1987; Fritts, 1995a). We summarize here the modeling studies with the mechanistic numerical spectral model (NSM) with Doppler spread parameterization for GW (Hines, 1997a, b), which describes in the middle atmosphere: (a) migrating and non-migrating DT, (b) planetary waves (PW), and (c) global-scale inertio gravity waves. Numerical experiments are discussed that illuminate the influence of GW filtering and nonlinear interactions between DT, PW, and zonal mean variations. Keywords: Theoretical modeling, Middle atmosphere dynamics, Gravity wave interactions, Migrating and non-migrating tides, Planetary waves, Global-scale inertio gravity waves.

Mayr, Hans G.↗
NASA NTRS2010

Application of Raytracing Through the High Resolution Numerical Weather Model HIRLAM for the Analysis of European VLBI

An important limitation for the precision in the results obtained by space geodetic techniques like VLBI and GPS are tropospheric delays caused by the neutral atmosphere, see e.g. [1]. In recent years numerical weather models (NWM) have been applied to improve mapping functions which are used for tropospheric delay modeling in VLBI and GPS data analyses. In this manuscript we use raytracing to calculate slant delays and apply these to the analysis of Europe VLBI data. The raytracing is performed through the limited area numerical weather prediction (NWP) model HIRLAM. The advantages of this model are high spatial (0.2 deg. x 0.2 deg.) and high temporal resolution (in prediction mode three hours).

Garcia-Espada, Susana↗
NASA NTRS2003

Numerical Study of Pyrolysis of Biomass in Fluidized Beds

A report presents a numerical-simulation study of pyrolysis of biomass in fluidized-bed reactors, performed by use of the mathematical model described in Model of Fluidized Bed Containing Reacting Solids and Gases (NPO-30163), which appears elsewhere in this issue of NASA Tech Briefs. The purpose of the study was to investigate the effect of various operating conditions on the efficiency of production of condensable tar from biomass. The numerical results indicate that for a fixed particle size, the fluidizing-gas temperature is the foremost parameter that affects the tar yield. For the range of fluidizing-gas temperatures investigated, and under the assumption that the pyrolysis rate exceeds the feed rate, the optimum steady-state tar collection was found to occur at 750 K. In cases in which the assumption was not valid, the optimum temperature for tar collection was found to be only slightly higher. Scaling up of the reactor was found to exert a small negative effect on tar collection at the optimal operating temperature. It is also found that slightly better scaling is obtained by use of shallower fluidized beds with greater fluidization velocities.

Bellan, Josette↗
NASA NTRS2007

Numerical Modeling of Propellant Boil-Off in a Cryogenic Storage Tank

A numerical model to predict boil-off of stored propellant in large spherical cryogenic tanks has been developed. Accurate prediction of tank boil-off rates for different thermal insulation systems was the goal of this collaboration effort. The Generalized Fluid System Simulation Program, integrating flow analysis and conjugate heat transfer for solving complex fluid system problems, was used to create the model. Calculation of tank boil-off rate requires simultaneous simulation of heat transfer processes among liquid propellant, vapor ullage space, and tank structure. The reference tank for the boil-off model was the 850,000 gallon liquid hydrogen tank at Launch Complex 39B (LC- 39B) at Kennedy Space Center, which is under study for future infrastructure improvements to support the Constellation program. The methodology employed in the numerical model was validated using a sub-scale model and tank. Experimental test data from a 1/15th scale version of the LC-39B tank using both liquid hydrogen and liquid nitrogen were used to anchor the analytical predictions of the sub-scale model. Favorable correlations between sub-scale model and experimental test data have provided confidence in full-scale tank boil-off predictions. These methods are now being used in the preliminary design for other cases including future launch vehicles

Majumdar, A. K.↗
NASA NTRS2012

Improved Impact of Atmospheric Infrared Sounder (AIRS) Radiance Assimilation in Numerical Weather Prediction

Improvements to global and regional numerical weather prediction (NWP) have been demonstrated through assimilation of data from NASA s Atmospheric Infrared Sounder (AIRS). Current operational data assimilation systems use AIRS radiances, but impact on regional forecasts has been much smaller than for global forecasts. Retrieved profiles from AIRS contain much of the information that is contained in the radiances and may be able to reveal reasons for this reduced impact. Assimilating AIRS retrieved profiles in an identical analysis configuration to the radiances, tracking the quantity and quality of the assimilated data in each technique, and examining analysis increments and forecast impact from each data type can yield clues as to the reasons for the reduced impact. By doing this with regional scale models individual synoptic features (and the impact of AIRS on these features) can be more easily tracked. This project examines the assimilation of hyperspectral sounder data used in operational numerical weather prediction by comparing operational techniques used for AIRS radiances and research techniques used for AIRS retrieved profiles. Parallel versions of a configuration of the Weather Research and Forecasting (WRF) model with Gridpoint Statistical Interpolation (GSI) that mimics the analysis methodology, domain, and observational datasets for the regional North American Mesoscale (NAM) model run at the National Centers for Environmental Prediction (NCEP)/Environmental Modeling Center (EMC) are run to examine the impact of each type of AIRS data set. The first configuration will assimilate the AIRS radiance data along with other conventional and satellite data using techniques implemented within the operational system; the second configuration will assimilate AIRS retrieved profiles instead of AIRS radiances in the same manner. Preliminary results of this study will be presented and focus on the analysis impact of the radiances and profiles for selected cases.

Zavodsky, Bradley↗
NASA NTRS2012

An Experimental and Numerical Study of Roughness-Induced Instabilities in a Mach 3.5 Boundary Layer

Progress on a joint experimental and numerical study of laminar-to-turbulent transition induced by an isolated roughness element in a high-speed laminar boundary layer is reported in this paper. The numerical analysis suggests that transition is driven by the instability of high- and low-speed streaks embedded in the wake of the isolated roughness element. In addition, spatial stability analysis revealed that the wake flow supports multiple modes (even and odd) of convective instabilities that experience strong enough growth to cause transition. The experimental measurements, which included hot-wire and pitot-probe surveys, confirmed the existence of embedded high- and low-speed streaks in the roughness wake. Furthermore, the measurements indicate the presence of both even and odd modes of instability, although their relative magnitude depends on the specifics of the roughness geometry and flow conditions (e.g., the value of Re(sub kk) or k/delta. For the two test cases considered in the measurements (Re(sub kk) values of 462 and 319), the even mode and the odd mode were respectively dominant and appear to play a primary role in the transition process.

Kegerise, Michael A.↗
NASA NTRS2012

SO2 over Central China: Measurements, Numerical Simulations and the Tropospheric Sulfur Budget

SO2 in central China was measured in situ from an aircraft and remotely using the Ozone Monitoring Instrument (OMI) from the Aura satellite; results were used to develop a numerical tool for evaluating the tropospheric sulfur budget - sources, sinks, transformation and transport. In April 2008, measured ambient SO2 concentrations decreased from approx.7 ppbv near the surface to approx. 1 ppbv at 1800 m altitude (an effective scale height of approx.800 m), but distinct SO2 plumes were observed between 1800 and 4500 m, the aircraft's ceiling. These free tropospheric plumes play a major role in the export of SO2 and in the accuracy of OMI retrievals. The mean SO2 column contents from aircraft measurements (0.73 DU, Dobson Units) and operational OMI SO2 products (0.63+/-0.26 DU) were close. The OMI retrievals were well correlated with in situ measurements (r = 0.84), but showed low bias (slope = 0.54). A new OMI retrieval algorithm was tested and showed improved agreement and bias (r = 0.87, slope = 0.86). The Community Multiscale Air Quality (CMAQ) model was used to simulate sulfur chemistry, exhibiting reasonable agreement (r = 0.62, slope = 1.33) with in situ SO2 columns. The mean CMAQ SO2 loading over central and eastern China was 54 kT, approx.30% more than the estimate from OMI SO2 products, 42 kT. These numerical simulations, constrained by observations, indicate that ",50% (35 to 61 %) of the anthropogenic sulfur emissions were transported downwind, and the overall lifetime of tropospheric SO2 was 38+/-7 h.

He, Hao↗
NASA NTRS2012

Numerical Relativity, Black Hole Mergers, and Gravitational Waves: Part III

This series of 3 lectures will present recent developments in numerical relativity, and their applications to simulating black hole mergers and computing the resulting gravitational waveforms. In this third and final lecture, we present applications of the results of numerical relativity simulations to gravitational wave detection and astrophysics.

Centrella, Joan↗
NASA NTRS2012

Numerical Relativity, Black Hole Mergers, and Gravitational Waves: Part I

This series of 3 lectures will present recent developments in numerical relativity, and their applications to simulating black hole mergers and computing the resulting gravitational waveforms. In this first lecture, we introduce the basic ideas of numerical relativity, highlighting the challenges that arise in simulating gravitational wave sources on a computer.

Centrella, Joan↗
NASA NTRS2011

Numerical Simulations on Origin of Galilean Moons' Magnetic Anomalies

Galileo mission detected the magnetic anomalies originated from Galilean moons. These anomalies are likely generated in the moons interiors, under the influence of a strong ambient Jovian field. Among various possible generation mechanisms of the anomalies, we focus on magneto-convection and dynamos in the interiors via numerical simulation. To mimic the electromagnetic environment of the moons, we introduce in our numerical model an external uniform magnetic field B(sub 0) with a fixed orientation but varying field strength. Our results show that a finite B(sub 0) can substantially alter the dynamo processes inside the core. When the ambient field strength B(sub 0) increases to approximately 40% of the field generated by the pure dynamo action, the convective state in the core changes significantly: the convective flow decreases by 80% in magnitude, but the differential rotation becomes stronger in much of the fluid layer, leading to a stronger field generated in the core. The field morphologies inside the core tend to align with the ambient field, while the flow patterns show the symmetry-breaking effect under the influence of B(sub 0). Furthermore, the generated field tends to be temporally more stable.

Jiao, LiQuo↗
NASA NTRS2011

A Fully Nonlinear, Dynamically Consistent Numerical Model for Solid-Body Ship Motion. I. Ship Motion with Fixed Heading

In this paper, we describe the details of our numerical model for simulating ship solidbody motion in a given environment. In this model, the fully nonlinear dynamical equations governing the time-varying solid-body ship motion under the forces arising from ship wave interactions are solved with given initial conditions. The net force and moment (torque) on the ship body are directly calculated via integration of the hydrodynamic pressure over the wetted surface and the buoyancy effect from the underwater volume of the actual ship hull with a hybrid finite-difference/finite-element method. Neither empirical nor free parametrization is introduced in this model, i.e. no a priori experimental data are needed for modelling. This model is benchmarked with many experiments of various ship hulls for heave, roll and pitch motion. In addition to the benchmark cases, numerical experiments are also carried out for strongly nonlinear ship motion with a fixed heading. These new cases demonstrate clearly the importance of nonlinearities in ship motion modelling.

Lin, Ray-Quing↗
NASA NTRS2014

Numerical Uncertainty Analysis for Computational Fluid Dynamics using Student T Distribution -- Application of CFD Uncertainty Analysis Compared to Exact Analytical Solution

Computational Fluid Dynamics (CFD) is the standard numerical tool used by Fluid Dynamists to estimate solutions to many problems in academia, government, and industry. CFD is known to have errors and uncertainties and there is no universally adopted method to estimate such quantities. This paper describes an approach to estimate CFD uncertainties strictly numerically using inputs and the Student-T distribution. The approach is compared to an exact analytical solution of fully developed, laminar flow between infinite, stationary plates. It is shown that treating all CFD input parameters as oscillatory uncertainty terms coupled with the Student-T distribution can encompass the exact solution.

Groves, Curtis E.↗
NASA NTRS2013

On the Theory and Numerical Simulation of Cohesive Crack Propagation with Application to Fiber-Reinforced Composites

The phenomenon of crack propagation is among the predominant modes of failure in many natural and engineering structures, often leading to severe loss of structural integrity and catastrophic failure. Thus, the ability to understand and a priori simulate the evolution of this failure mode has been one of the cornerstones of applied mechanics and structural engineering and is broadly referred to as "fracture mechanics." The work reported herein focuses on extending this understanding, in the context of through-thickness crack propagation in cohesive materials, through the development of a continuum-level multiscale numerical framework, which represents cracks as displacement discontinuities across a surface of zero measure. This report presents the relevant theory, mathematical framework, numerical modeling, and experimental investigations of through-thickness crack propagation in fiber-reinforced composites using the Variational Multiscale Cohesive Method (VMCM) developed by the authors.

Rudraraju, Siva Shankar↗
NASA NTRS2013

A Numerical Combination of Extended Boundary Condition Method and Invariant Imbedding Method Applied to Light Scattering by Large Spheroids and Cylinders

The extended boundary condition method (EBCM) and invariant imbedding method (IIM) are two fundamentally different T-matrix methods for the solution of light scattering by nonspherical particles. The standard EBCM is very efficient but encounters a loss of precision when the particle size is large, the maximum size being sensitive to the particle aspect ratio. The IIM can be applied to particles in a relatively large size parameter range but requires extensive computational time due to the number of spherical layers in the particle volume discretization. A numerical combination of the EBCM and the IIM (hereafter, the EBCM+IIM) is proposed to overcome the aforementioned disadvantages of each method. Even though the EBCM can fail to obtain the T-matrix of a considered particle, it is valuable for decreasing the computational domain (i.e., the number of spherical layers) of the IIM by providing the initial T-matrix associated with an iterative procedure in the IIM. The EBCM+IIM is demonstrated to be more efficient than the IIM in obtaining the optical properties of large size parameter particles beyond the convergence limit of the EBCM. The numerical performance of the EBCM+IIM is illustrated through representative calculations in spheroidal and cylindrical particle cases.

invariant imbeddings↗
NASA NTRS2014

Numerical Simulations For the F-16XL Aircraft Configuration

Numerical simulations of flow around the F-16XL are presented as a contribution to the Cranked Arrow Wing Aerodynamic Project International II (CAWAPI-II). The NASA Tetrahedral Unstructured Software System (TetrUSS) is used to perform numerical simulations. This CFD suite, developed and maintained by NASA Langley Research Center, includes an unstructured grid generation program called VGRID, a postprocessor named POSTGRID, and the flow solver USM3D. The CRISP CFD package is utilized to provide error estimates and grid adaption for verification of USM3D results. A subsonic high angle-of-attack case flight condition (FC) 25 is computed and analyzed. Three turbulence models are used in the calculations: the one-equation Spalart-Allmaras (SA), the two-equation shear stress transport (SST) and the kε turbulence models. Computational results, and surface static pressure profiles are presented and compared with flight data. Solution verification is performed using formal grid refinement studies, the solution of Error Transport Equations, and adaptive mesh refinement. The current study shows that the USM3D solver coupled with CRISP CFD can be used in an engineering environment in predicting vortex-flow physics on a complex configuration at flight Reynolds numbers.

Elmiligui, Alaa A.↗
NASA NTRS2011

Direct Numerical Simulations of Transitional/Turbulent Wakes

The interest in transitional/turbulent wakes spans the spectrum from an intellectual pursuit to understand the complex underlying physics to a critical need in aeronautical engineering and other disciplines to predict component/system performance and reliability. Cylinder wakes have been studied extensively over several decades to gain a better understanding of the basic flow phenomena that are encountered in such flows. Experimental, computational and theoretical means have been employed in this effort. While much has been accomplished there are many important issues that need to be resolved. The physics of the very near wake of the cylinder (less than three diameters downstream) is perhaps the most challenging of them all. This region comprises the two detached shear layers, the recirculation region and wake flow. The interaction amongst these three components is to some extent still a matter of conjecture. Experimental techniques have generated a large percentage of the data that have provided us with the current state of understanding of the subject. More recently computational techniques have been used to simulate cylinder wakes, and the data from such simulations are being used to both refine our understanding of such flows as well as provide new insights. A few large eddy and direct numerical simulations (LES and DNS) of cylinder wakes have appeared in the literature in the recent past. These investigations focus on the low Reynolds number range where the cylinder boundary layer is laminar (sub-critical range). However, from an engineering point of view, there is considerable interest in the situation where the upper and/or lower boundary layer of an airfoil is turbulent, and these turbulent boundary layers separate from the airfoil to contribute to the formation of the wake downstream. In the case of cylinders, this only occurs at relatively large unit Reynolds numbers. However, in the case of airfoils, the boundary layer has the opportunity to transition to turbulence on the airfoil surface at a relatively lower unit Reynolds number because the characteristic length of the airfoil is typically one to two orders of magnitude larger than the trailing edge diameter. This transition to turbulence would occur unless there is a strong favorable pressure gradient that results in the boundary layer remaining laminar or transitional over the surface of the airfoil. This presentation will focus on two direct numerical simulations that have been performed at NASA ARC. The first is of a cylinder wake with laminar separating boundary layers. The second is the wake of a flat plate with a circular trailing edge. The upper and lower plate surface boundary layers are both turbulent and statistically identical. Thus the computed wake is symmetric in a statistical sense. This flow is more representative of airfoil wakes than cylinder wakes. Results from the two simulations including flow visualization and turbulence statistics in the near wake will be presented at the seminar.

wakes↗
NASA NTRS2014

Numerical Simulation of a High-Lift Configuration with Embedded Fluidic Actuators

Numerical simulations have been performed for a vertical tail configuration with deflected rudder. The suction surface of the main element of this configuration is embedded with an array of 32 fluidic actuators that produce oscillating sweeping jets. Such oscillating jets have been found to be very effective for flow control applications in the past. In the current paper, a high-fidelity computational fluid dynamics (CFD) code known as the PowerFLOW(Registered TradeMark) code is used to simulate the entire flow field associated with this configuration, including the flow inside the actuators. The computed results for the surface pressure and integrated forces compare favorably with measured data. In addition, numerical solutions predict the correct trends in forces with active flow control compared to the no control case. Effect of varying yaw and rudder deflection angles are also presented. In addition, computations have been performed at a higher Reynolds number to assess the performance of fluidic actuators at flight conditions.

Vatsa, Veer N.↗