SEARCH · Engineering Papers

Results for “Direct numerical simulation”

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 1,387 records · Page 77

All sources NASA NTRS DOE OSTI

NASA NTRS1973

Numerical calculation of the three dimensional transonic flow over a yawed wing.

Results are presented of calculations of the three dimensional steady transonic flow over a finite yawed wing. The full potential flow equation is solved in a transformed coordinate system which permits the boundary conditions to be satisfied exactly. The correct differential properties are enforced by rotating the difference scheme to conform with the flow direction, and fast convergence is assured by simulating a time dependent equation designed to settle quickly to a steady state. Computed lift drag ratios are consistent with the results of wind tunnel tests of a yawed wing conducted by R. T. Jones (1972).-

Jameson, A.↗

NASA NTRS1985

Evaluation of a pulse control law for flexible spacecraft

The following analytical and experimental studies were conducted: (1) A simple algorithm was developed to suppress the structural vibrations of 3-dimensional distributed parameter systems, subjected to interface motion and/or directly applied forces. The algorithm is designed to cope with structural oscillations superposed on top of rigid-body motion: a situation identical to that encountered by the SCOLE components. A significant feature of the method is that only local measurements of the structural displacements and velocities relative to the moving frame of reference are needed. (2) A numerical simulation study was conducted on a simple linear finite element model of a cantilevered plate which was subjected to test excitations consisting of impulsive base motion and of nonstationary wide-band random excitation applied at its root. In each situation, the aim was to suppress the vibrations of the plate relative to the moving base. (3) A small mechanical model resembling an aircraft wing was designed and fabricated to investigate the control algorithm under realistic laboratory conditions.

Source record↗

NASA NTRS1990

A study of the noise mechanisms of transonic blade-vortex interactions

Transonic blade-vortex interactions (BVI) are simulated numerically and the noise mechanisms are investigated. The two-dimensional high frequency transonic small disturbance equation is solved numerically (VTRAN2 code). An ADI scheme with monotone switches is used; viscous effects are included on the boundary, and the vortex is simulated by the cloud in cell method. The Kirchhoff method is used for the extension of the numerical two-dimensional near-field aerodynamic results to the linear acoustic three dimensional far field. The viscous effects (shock/boundary layer interactions) on BVI is investigated. The different types of shock motion are identified and compared. Two important disturbances with different directivity exist in the pressure signal and are believed to be related to the fluctuating lift and drag forces. Noise directivity for different cases is shown. The maximum radiation occurs at an angle between 60 and 90 degrees below the horizontal for an airfoil-fixed coordinate system and depends on the details of the airfoil shape. Different airfoil shapes are studied and classified according to the BVI noise produced.

Lyrintzis, Anastasios S.↗

NASA NTRS1993

Multigrid methods for numerical simulation of laminar diffusion flames

This paper documents the result of a computational study of multigrid methods for numerical simulation of 2D diffusion flames. The focus is on a simplified combustion model, which is assumed to be a single step, infinitely fast and irreversible chemical reaction with five species (C3H8, O2, N2, CO2 and H2O). A fully-implicit second-order hybrid scheme is developed on a staggered grid, which is stretched in the streamwise coordinate direction. A full approximation multigrid scheme (FAS) based on line distributive relaxation is developed as a fast solver for the algebraic equations arising at each time step. Convergence of the process for the simplified model problem is more than two-orders of magnitude faster than other iterative methods, and the computational results show good grid convergence, with second-order accuracy, as well as qualitatively agreement with the results of other researchers.

Liu, C.↗

NASA NTRS1997

A Thermal and Electrical Analysis of Power Semiconductor Devices

The state-of-art power semiconductor devices require a thorough understanding of the thermal behavior for these devices. Traditional thermal analysis have (1) failed to account for the thermo-electrical interaction which is significant for power semiconductor devices operating at high temperature, and (2) failed to account for the thermal interactions among all the levels involved in, from the entire device to the gate micro-structure. Furthermore there is a lack of quantitative studies of the thermal breakdown phenomenon which is one of the major failure mechanisms for power electronics. This research work is directed towards addressing. Using a coupled thermal and electrical simulation, in which the drift-diffusion equations for the semiconductor and the energy equation for temperature are solved simultaneously, the thermo-electrical interactions at the micron scale of various junction structures are thoroughly investigated. The optimization of gate structure designs and doping designs is then addressed. An iterative numerical procedure which incorporates the thermal analysis at the device, chip and junction levels of the power device is proposed for the first time and utilized in a BJT power semiconductor device. In this procedure, interactions of different levels are fully considered. The thermal stability issue is studied both analytically and numerically in this research work in order to understand the mechanism for thermal breakdown.

Vafai, Kambiz↗

NASA NTRS1997

Methods for the Organogenesis of Skeletal Muscle in Tissue Culture

Skeletal muscle structure is regulated by many factors, including nutrition, hormones, electrical activity, and tension. The muscle cells are subjected to both passive and active mechanical forces at all stages of development and these forces play important but poorly understood roles in regulating muscle organogenesis and growth. For example, during embryogenesis, the rapidly growing skeleton places large passive mechanical forces on the attached muscle tissue. These forces not only help to organize the proliferating mononucleated myoblasts into the oriented, multinucleated myofibers of a functional muscle but also tightly couple the growth rate of muscle to that of bone. Postnatally, the actively contracting, innervated muscle fibers are subjected to different patterns of active and passive tensions which regulate longitudinal and cross sectional myofiber growth. These mechanically-induced organogenic processes have been difficult to study under normal tissue culture conditions, resulting in the development of numerous methods and specialized equipment to simulate the in vivo mechanical environment.These techniques have led to the "engineering" of bioartificial muscles (organoids) which display many of the characteristics of in vivo muscle including parallel arrays of postmitotic fibers organized into fascicle-like structures with tendon-like ends. They are contractile, express adult isoforms of contractile proteins, perform directed work, and can be maintained in culture for long periods. The in vivo-like characteristics and durability of these muscle organoids make them useful for long term in vitro studies on mechanotransduction mechanisms and on muscle atrophy induced by decreased tension. In this report, we described a simple method for generating muscle organoids from either primary embrionic avain or neonatal rodent myoblasts.

Vandenburgh, Herman↗

NASA NTRS1999

El Nino and the Global Ocean Observing System

Until a decade ago, an often-quoted expression in oceanography is that very few observations are recorded throughout the ocean. Now, the sentiment is no longer valid in the uppermost 10% of the tropical Pacific Ocean nor at the surface of the global ocean. One of the remarkable legacies of the 1985-1994 Tropical Oceans Global Atmosphere (TOGA) Program is an in situ marine meteorological and upper oceanographic measurement array throughout the equatorial Pacific to monitor the development and maintenance of El Nino episodes. The TOGA Observing System, which initially consisted of moored- and drifting-buoy arrays, a network of commercial ships, and coastal and island stations, now includes a constellation of satellites and data-assimilating models to simulate subsurface oceanographic conditions. The El Nino and La Nina tropical Pacific Ocean observing system represents the initial phase of an integrated global ocean observing system. Remarkable improvements have been made in ocean model simulation of subsurface currents, but some problems persist. For example, the simulation of the South Equatorial Current (SEC) remains an important challenge in the 2S-2N Pacific equatorial wave guide. During El Nino the SEC at the equator is reduced and sometimes the direction is reversed, becoming eastward. Both conditions allow warm water stored in the western Pacific to invade the eastern region, creating an El Nino episode. Assimilation of data is a tenet of faith to correct simulation errors caused by deficiencies in surface fluxes (especially wind stress) and parameterizations of subgrid-scale physical processes. In the first of two numerical experiments, the Pacific SEC was simulated with and without assimilation of subsurface temperature data. Along the equator, a very weak SEC occurred throughout the eastern Pacific, independent of assimilation of data. However, as displayed in the diagram, in the western Pacific there was no satisfactory agreement between the two simulations. To help determine reliability of the simulated SEC in the western Pacific, current measurements recorded during the 9-19 October 1994 voyage of the French research vessel L'Atalante are also shown in the diagram. With data assimilation, the simulated SEC was in much better agreement with L'Atalante observations. The simulated SEC with data assimilation was far from perfect, in part because of the sparsity of subsurface temperature observations. In the next experiment, TOPEX/POSEIDON sea surface height data in combination with subsurface temperatures will be assimilated to assess further improvement of the simulation of the SEC.

Halpern, David↗

NASA NTRS1989

A perspective on future directions in aerospace propulsion system simulation

The design and development of aircraft engines is a lengthy and costly process using today's methodology. This is due, in large measure, to the fact that present methods rely heavily on experimental testing to verify the operability, performance, and structural integrity of components and systems. The potential exists for achieving significant speedups in the propulsion development process through increased use of computational techniques for simulation, analysis, and optimization. This paper outlines the concept and technology requirements for a Numerical Propulsion Simulation System (NPSS) that would provide capabilities to do interactive, multidisciplinary simulations of complete propulsion systems. By combining high performance computing hardware and software with state-of-the-art propulsion system models, the NPSS will permit the rapid calculation, assessment, and optimization of subcomponent, component, and system performance, durability, reliability and weight-before committing to building hardware.

Miller, Brent A.↗

NASA NTRS1996

Multigrid Methods for Fully Implicit Oil Reservoir Simulation

In this paper we consider the simultaneous flow of oil and water in reservoir rock. This displacement process is modeled by two basic equations: the material balance or continuity equations and the equation of motion (Darcy's law). For the numerical solution of this system of nonlinear partial differential equations there are two approaches: the fully implicit or simultaneous solution method and the sequential solution method. In the sequential solution method the system of partial differential equations is manipulated to give an elliptic pressure equation and a hyperbolic (or parabolic) saturation equation. In the IMPES approach the pressure equation is first solved, using values for the saturation from the previous time level. Next the saturations are updated by some explicit time stepping method; this implies that the method is only conditionally stable. For the numerical solution of the linear, elliptic pressure equation multigrid methods have become an accepted technique. On the other hand, the fully implicit method is unconditionally stable, but it has the disadvantage that in every time step a large system of nonlinear algebraic equations has to be solved. The most time-consuming part of any fully implicit reservoir simulator is the solution of this large system of equations. Usually this is done by Newton's method. The resulting systems of linear equations are then either solved by a direct method or by some conjugate gradient type method. In this paper we consider the possibility of applying multigrid methods for the iterative solution of the systems of nonlinear equations. There are two ways of using multigrid for this job: either we use a nonlinear multigrid method or we use a linear multigrid method to deal with the linear systems that arise in Newton's method. So far only a few authors have reported on the use of multigrid methods for fully implicit simulations. Two-level FAS algorithm is presented for the black-oil equations, and linear multigrid for two-phase flow problems with strong heterogeneities and anisotropies is studied. Here we consider both possibilities. Moreover we present a novel way for constructing the coarse grid correction operator in linear multigrid algorithms. This approach has the advantage in that it preserves the sparsity pattern of the fine grid matrix and it can be extended to systems of equations in a straightforward manner. We compare the linear and nonlinear multigrid algorithms by means of a numerical experiment.

Molenaar, J.↗

NASA NTRS1973

Analytical and Radio-Histo-Chemical Experiments of Plants and Tissue Culture Cells Treated with Lunar and Terrestrial Materials

The nature and mechanisms of the apparent simulation of growth originally observed in plants growing in contact with lunar soil during the Apollo project quarantine are examined. Preliminary experiments employing neutron activated lunar soil indicate uptake of a few elements by plants. It was found that while the preliminary neutron activation technique allowed demonstration of uptake of minerals it presented numerous disadvantages for use in critical experiments directed at elucidating possible mechanisms of stimulation.

Halliwell, R. S.↗

NASA NTRS1989

The dynamics of charged particles in the near wake of a very negatively charged body - Laboratory experiment and numerical simulation

A numerical simulation that is cylindrical in configuration space and three-dimensional in velocity space has been initiated to test a model for the near-wake dynamics of a very negatively charged body, with reference to the plasma environment around spacecraft. The simulation parameters were closely matched to those of a laboratory experiment so that the results can be compared directly. The laboratory study showed that the electrons and ions can display different temporal features in the filling-in of the wake; and that they can both be found within one body diameter of an object with a highly negative body potential. It was also found that the temperature of the electrons in the very near wake could be somewhat colder than the ambient value, suggesting the possibility of a filtering mechanism being operative there. The simulation results to date largely corroborate the density findings.

Morgan, M. Alvin↗

NASA NTRS2001

A Highly Accurate Technique for the Treatment of Flow Equations at the Polar Axis in Cylindrical Coordinates using Series Expansions

Numerical methods for solving the flow equations in cylindrical or spherical coordinates should be able to capture the behavior of the exact solution near the regions where the particular form of the governing equations is singular. In this work we focus on the treatment of these numerical singularities for finite-differences methods by reinterpreting the regularity conditions developed in the context of pseudo-spectral methods. A generally applicable numerical method for treating the singularities present at the polar axis, when nonaxisymmetric flows are solved in cylindrical, coordinates using highly accurate finite differences schemes (e.g., Pade schemes) on non-staggered grids, is presented. Governing equations for the flow at the polar axis are derived using series expansions near r=0. The only information needed to calculate the coefficients in these equations are the values of the flow variables and their radial derivatives at the previous iteration (or time) level. These derivatives, which are multi-valued at the polar axis, are calculated without dropping the accuracy of the numerical method using a mapping of the flow domain from (0,R)*(0,2pi) to (-R,R)*(0,pi), where R is the radius of the computational domain. This allows the radial derivatives to be evaluated using high-order differencing schemes (e.g., compact schemes) at points located on the polar axis. The proposed technique is illustrated by results from simulations of laminar-forced jets and turbulent compressible jets using large eddy simulation (LES) methods. In term of the general robustness of the numerical method and smoothness of the solution close to the polar axis, the present results compare very favorably to similar calculations in which the equations are solved in Cartesian coordinates at the polar axis, or in which the singularity is removed by employing a staggered mesh in the radial direction without a mesh point at r=0, following the method proposed recently by Mohseni and Colonius (1). Extension of the method described here for incompressible flows or for any other set of equations that are solved on a non-staggered mesh in cylindrical or spherical coordinates with finite-differences schemes of various level of accuracy is immediate.

Constantinescu, George S.↗

NASA NTRS2022

Wind-Optimal Lateral Trajectories for a Multirotor Aircraft in Urban Air Mobility

The primary motivation for this paper is to quantify the operational benefits (energy consumption and flight duration) of flying wind-optimal lateral trajectories for short flights (less than 60 miles) anticipated in the urban environment. The optimal control model presented includes a wind model for quantifying the effect of wind on the lateral trajectory. The optimal control problem is numerically solved using the direct collocation method. Energy consumption and flight duration flying wind-optimal lateral trajectories are compared with corresponding values obtained flying great-circle paths between the same origin and destination pairs to determine the operational benefits of wind-optimal routing for short flights. The flight duration results for different scenarios are validated using a simulation tool designed and developed at NASA for exploring advanced air traffic management concepts. This research study suggests that for short flights in an urban environment, operational benefits of the wind-optimal lateral trajectories over the corresponding great-circle trajectories in terms of energy consumption and flight duration per flight are dependent on: i) wind field’s spatial variability, ii) wind magnitude, iii) the direction of route relative to the wind field, and iv) cruise segment length. The operational benefits observed in realistic flyable wind scenarios are less than 2.5%; these could be translated to an equivalent of a maximum of 2 min of cruise flight duration savings in the urban air mobility environment. As expected, headwinds and tailwinds along the flight route most significantly impact energy consumption and flight duration.

Wind-Optimal Trajectory, Lateral Path, Urban Air M↗

NASA NTRS1982

An explanation for the systematic flow of plasma in the solar transition region

Using numerical simulations, it is shown that the systematic flow of plasma along a coronal magnetic flux tube is easily produced by a change in the spatial dependence of the heating rate from a symmetric deposition which supports a stationary equilibrium to a time-independent asymmetric deposition. The velocity of the flow is roughly proportional to the heating symmetry and is directed to the side of the loop away from the bulk of the energy deposition.

Boris, J. P.↗

NASA NTRS2005

Spaceborne Doppler Radar Measurements of Rainfall: Correction of Errors Induced by Pointing Uncertainties

In this paper a sea surface radar echo spectral analysis technique to correct for the rainfall velocity error caused by radar-pointing uncertainty is presented. The correction procedure is quite straightforward when the radar is observing a homogeneous rainfall field. When nonuniform beam filling (NUBF) occurs and attenuating frequencies are used, however, additional steps are necessary in order to correctly estimate the antenna-pointing direction. This new technique relies on the application of the combined frequency-time (CFT) algorithm to correct for uneven attenuation effects on the observed sea surface Doppler spectrum. The performance of this correction technique was evaluated by a Monte Carlo simulation of the Doppler precipitation radar backscatter from high-resolution 3D rain fields (either generated by a cloud resolving numerical model or retrieved from airborne radar measurements). The results show that the antenna-pointing-induced error can, indeed, be reduced by the proposed technique in order to achieve 1 m s(exp -1) accuracy on rainfall vertical velocity estimates.

remote sensing↗

NASA NTRSDate not supplied

Progress Toward Evaluating Roughness-Resolved Aerothermodynamic Simulations for the Mars Sample Return Earth Entry System

This work pursues the efficient Navier-Stokes simulation, using NASA’s standard LAURA and DPLR simulation tools, of the flowfield and convective heating resulting from the patterned-roughness surface of the Mars Sample Return (MSR) Earth Entry System (EES). This type of simulation has been infeasible for past NASA missions. However, for EES, the combination of a patterned roughness heatshield approximated with a sinusoidal pattern, the zero degrees angle of attack nominal flow, and the supersonic boundary layer edge Mach number, all combine to make these roughness resolved grid simulations feasible. The patterned roughness and zero degrees angle of attack nominal flow allow the simulation domain to be reduced to a thin axial flowfield slice that captures a minimum of single period of the roughness pattern. With a periodic boundary condition applied along the boundaries of this slice, simulation on the reduced domain is identical to the full domain but with significant cost saving. The supersonic boundary layer edge Mach number prevents the upstream influence, which allows the smooth stagnation region to be frozen. This freezing of the stagnation region flow is required because of the numerical instability resulting from the tight spacing of the grid in the axial direction as the stagnation point is approached. Applying the developed analysis approach to a peak heating EES case results in a roughness augmentation factor for the convective heating of 1.32, which agrees well with correlations that predict values between 1.30 and 1.37. Analysis approach developed in this work provides the framework for future studies of turbulence modeling impacts and the influence of an ablating surface.

Earth Entry System↗

NASA NTRSDate not supplied

Statistical Approach on Utilizing Ground-based Experiments to Model Break-up Events

Two laboratory-based impact tests have been conducted to develop and extend the capability of NASA satellite breakup models. The first experiment, the Satellite Orbital Debris Characterization Impact Test (SOCIT), was conducted by the U.S. Department of Defense and NASA in 1992. It employed a fully functional U.S. Navy Transit spacecraft, fabricated from materials commonly used in the 1960s. SOCIT fragmentation data formed the basis of the current NASA Standard Satellite Breakup Model (SSBM). A second experiment, DebriSat, was conducted in 2014 by the DebriSat consortium: the NASA Orbital Debris Program Office; the United States Space Force Space Systems Command, formerly the Air Force Space and Missile Systems Center, the Air Force Arnold Engineering Development Complex, and the University of Florida. This impact test was performed on a high-fidelity mock-up satellite assembled from modern components. Data from both experiments are expected to contribute to the next-generation models for on-orbit breakup analyses, long-term environment predictions, and debris risk assessment. This paper uses the direct statistical sampling of the SOCIT and DebriSat data ensembles to model an on-orbit breakup event, rather than the analytic expressions estimated from the samples as with as with the NASA SSBM and its predecessor. This direct method involves drawing fragments (i.e., sampling with replacement) from each fragmentation data-subset containing mass, shape category, material density, characteristic length, mass, and cross-sectional area. As part of the sampling process, the properties of sampled fragment data (e.g., area-to-mass ratios) are numerically checked to ensure they do not contain unrealistic quantities. The process to simulate a breakup cloud composed of fragments from SOCIT and DebriSat datasets is discussed. The methodology for ensuring that conservation of overall mass of the sampled fragment cloud under this sampling approach is highlighted. Finally, the results are compared with simulated clouds generated from SSBM for specific historical breakup events.

Andrew Vavrin↗

NASA NTRS2023

Statistical Approach on Utilizing Ground-Based Experiments to Model Break-Up Events

Two laboratory-based impact tests have been conducted to develop and extend the capability of NASA satellite breakup models. The first experiment, the Satellite Orbital Debris Characterization Impact Test (SOCIT), was conducted by the U.S. Department of Defense and NASA in 1992. It employed a fully functional U.S. Navy Transit spacecraft, fabricated from materials commonly used in the 1960s. SOCIT fragmentation data formed the basis of the current NASA Standard Satellite Breakup Model (SSBM). A second experiment, DebriSat, was conducted in 2014 by the DebriSat consortium: the NASA Orbital Debris Program Office, the United States Space Force Space Systems Command, formerly the Air Force Space and Missile Systems Center, the Air Force Arnold Engineering Development Complex, and the University of Florida. This impact test was performed on a high-fidelity mock-up satellite assembled from modern components. Data from both experiments are expected to contribute to the next-generation models for on-orbit breakup analyses, long-term environment predictions, and debris risk assessment. This paper uses the direct statistical sampling of the SOCIT and DebriSat data ensembles to model an on-orbit breakup event, rather than the analytic expressions estimated from the samples in the SSBM and its predecessor. This direct method involves drawing fragments (i.e., sampling with replacement) from each fragmentation data-subset containing mass, shape category, material density, characteristic length, mass, and cross-sectional area. As part of the sampling process, the properties of sampled fragment data (e.g., area-to-mass ratios) are numerically checked to ensure they do not contain unrealistic quantities. The process to simulate a breakup cloud composed of fragments from SOCIT and DebriSat datasets is discussed. The methodology for ensuring that conservation of overall mass of the sampled fragment cloud under this sampling approach is highlighted. Finally, the results are compared with simulated clouds generated from SSBM for specific historical breakup events.

Andrew Vavrin↗