Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simplified Simulator”

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 181 records · Page 10

ULF waves at Comets Halley and Giacobini-Zinner - Comparison with theory

Because the Doppler shift occurring in ion-cyclotron resonance is equal and opposite to that associated with the flow of the solar wind plasma past the observer, one should be able to monitor the amplitude of waves due to the cyclotron resonant instability by examining the amplitude of waves near the gyrofrequency of the respective ions as measured in the spacecraft frame. This expectation has been checked by using one-dimensional electromagnetic hybrid simulations of two ion pickup. This simplifies the comparison of theoretical predictions of wave growth with the observations. Distant from both comets the wave amplitudes are very weak. As the comet is approached the wave amplitude grows. These observations agree qualitatively with the simulations which predict amplitudes that depend on the properties of the injected beams and the local injection rate.

Le, G.↗

Slat Cove Noise Modeling: A Posteriori Analysis of Unsteady RANS Simulations

A companion paper by Khorrami et al demonstrates the feasibility of simulating the (nominally) self-sustained, large-scale unsteadiness within the leading-edge slat-cove region of multi-element airfoils using unsteady Reynolds-Averaged Navier-Stokes (URANS) equations, provided that the turbulence production term in the underlying two-equation turbulence model is switched off within the cove region. In conjunction with a FfowesWilliams-Hawkings solver, the URANS computations were shown to capture the dominant portion of the acoustic spectrum attributed to slat noise, as well as reproducing the increased intensity of slat cove motions (and, correspondingly, far-field noise as well) at the lower angles of attack. This paper examines that simulation database, augmented by additional simulations, with the objective of transitioning this apparent success to aeroacoustic predictions in an engineering context. As a first step towards this goal, the simulated flow and acoustic fields are compared with experiment and simplified analytical model. Rather intense near-field fluctuations in the simulated flow are found to be associated with unsteady separation along the slat bottom surface, relatively close to the slat cusp. Accuracy of the laminar-cove simulations in this near-wall region is raised to be an open issue. The adjoint Green's function approach is also explored in an attempt to identify the most efficient noise source locations.

Choudhari, Meelan↗

Prediction of high temperature metal matrix composite ply properties

The application of the finite element method (superelement technique) in conjunction with basic concepts from mechanics of materials theory is demonstrated to predict the thermomechanical behavior of high temperature metal matrix composites (HTMMC). The simulated behavior is used as a basis to establish characteristic properties of a unidirectional composite idealized an as equivalent homogeneous material. The ply properties predicted include: thermal properties (thermal conductivities and thermal expansion coefficients) and mechanical properties (moduli and Poisson's ratio). These properties are compared with those predicted by a simplified, analytical composite micromechanics model. The predictive capabilities of the finite element method and the simplified model are illustrated through the simulation of the thermomechanical behavior of a P100-graphite/copper unidirectional composite at room temperature and near matrix melting temperature. The advantage of the finite element analysis approach is its ability to more precisely represent the composite local geometry and hence capture the subtle effects that are dependent on this. The closed form micromechanics model does a good job at representing the average behavior of the constituents to predict composite behavior.

Caruso, J. J.↗

Modeling of the Ablation of Fibrous Materials in the Knudsen Regime

During atmospheric entry of planetary probes, the thermal protection system (TPS) of the probe is exposed to high temperatures under low pressures. In these conditions, carbonous TPS materials undergo gasification in the Knudsen regime leading to mass loss and wall recession called ablation. This work aims to improve the understanding of materiaVenvironment interactions through a study of the coupling between carbon dioxide transport in the Knudsen regime, heterogeneous oxidation of carbon, and sutface recession. A 3D Monte-Carlo simulation tool is used for this study. The fibrous architecture of the materiils, consisting of high porosity random array of carbon fibers, is numerically reproduced on a 3D Cartesian grid. Mass transport in the Knudsen regime from the boundary layer to the surface, and inside this porous material is simulated by random walk. A reaction probability is used to simulate the heterogeneous oxidation reaction. The surface recession is followed by front tracking using a simplified marching cube approach. The output data of the simulations are ablation velocity and dynamic evolution of the material porosity. A parametric study is carried out to analyze the material behavior as a function of Knudsen number for the porous media (length of the mean free path compared to the mean pore diameter) and the intrinsic reactivity of the carbon fibers. The results enable extrapolation of laboratory experimental data to actual entry conditions.

Lachaud, J.↗

CFD and Experimental Data Comparisons for Conventional and AFC-Enabled CRM High-Lift Configurations

Numerical simulations have been performed for a simplified high-lift (SHL) version of the Common Research Model (CRM) configuration, where the Fowler flaps of a representative conventional high-lift (CRM-HL) configuration are replaced by a set of simple hinged flaps. These hinged flaps are equipped with integrated modular active flow control (AFC) cartridges on the suction surface of the flap shoulder, and the resulting geometry is known as the CRM-SHL-AFC configuration. The main objective is to make use of AFC devices on the CRM-SHL-Configuration to produce the aerodynamic performance (lift) comparable to that of the CRM-HL configuration over a large angle of attack range encompassing maximum lift conditions. For comparison purposes, computations are also per-formed for the CRM-HL configuration. In the current paper, PowerFLOWR©, a CFD code based on the Lattice Boltzmann method (LBM), is used to simulate the entire flow field associated with the CRM-SHL-AFC configuration equipped with a combination of two different types of AFC actuators. The transonic version of the PowerFLOWR© code that has been validated for high-speed flows is used to simulate the flow field generated by the high-momentum actuators required to mitigate reversed flow regions on the suction surfaces of the main wing and the flap. This study is focused on the AFC systems and actuator arrangements that had emerged based on the parametric studies conducted previously at the nominal landing condition. Comparisons of the numerical solutions for lift and surface pressures are presented here with the experimental data, demonstrating the usefulness of CFD for predicting the flow field and lift characteristics of AFC-enabled high-lift configurations over a broad angle of attack range. The numerical solutions predict the expected trends in aerodynamic forces with angle of attack variation.

Veer N Vatsa↗

Fiber Composite Sandwich Thermostructural Behavior: Computational Simulation

Several computational levels of progressive sophistication/simplification are described to computationally simulate composite sandwich hygral, thermal, and structural behavior. The computational levels of sophistication include: (1) three-dimensional detailed finite element modeling of the honeycomb, the adhesive and the composite faces; (2) three-dimensional finite element modeling of the honeycomb assumed to be an equivalent continuous, homogeneous medium, the adhesive and the composite faces; (3) laminate theory simulation where the honeycomb (metal or composite) is assumed to consist of plies with equivalent properties; and (4) derivations of approximate, simplified equations for thermal and mechanical properties by simulating the honeycomb as an equivalent homogeneous medium. The approximate equations are combined with composite hygrothermomechanical and laminate theories to provide a simple and effective computational procedure for simulating the thermomechanical/thermostructural behavior of fiber composite sandwich structures.

Chamis, C. C.↗

Propagation of Electrical Excitation in a Ring of Cardiac Cells: A Computer Simulation Study

The propagation of electrical excitation in a ring of cells described by the Noble, Beeler-Reuter (BR), Luo-Rudy I (LR I), and third-order simplified (TOS) mathematical models is studied using computer simulation. For each of the models it is shown that after transition from steady-state circulation to quasi-periodicity achieved by shortening the ring length (RL), the action potential duration (APD) restitution curve becomes a double-valued function and is located below the original ( that of an isolated cell) APD restitution curve. The distributions of APD and diastolic interval (DI) along a ring for the entire range of RL corresponding to quasi-periodic oscillations remain periodic with the period slightly different from two RLs. The 'S' shape of the original APD restitution curve determines the appearance of the second steady-state circulation region for short RLs. For all the models and the wide variety of their original APD restitution curves, no transition from quasi-periodicity to chaos was observed.

Kogan, B. Y.↗

Interannual Variability of Martian Global Dust Storms: Simulations with a Low-Order Model of the General Circulation

We present simulations of the interannual variability of martian global dust storms (GDSs) with a simplified low-order model (LOM) of the general circulation. The simplified model allows one to conduct computationally fast long-term simulations of the martian climate system. The LOM is constructed by Galerkin projection of a 2D (zonally averaged) general circulation model (GCM) onto a truncated set of basis functions. The resulting LOM consists of 12 coupled nonlinear ordinary differential equations describing atmospheric dynamics and dust transport within the Hadley cell. The forcing of the model is described by simplified physics based on Newtonian cooling and Rayleigh friction. The atmosphere and surface are coupled: atmospheric heating depends on the dustiness of the atmosphere, and the surface dust source depends on the strength of the atmospheric winds. Parameters of the model are tuned to fit the output of the NASA AMES GCM and the fit is generally very good. Interannual variability of GDSs is possible in the IBM, but only when stochastic forcing is added to the model. The stochastic forcing could be provided by transient weather systems or some surface process such as redistribution of the sand particles in storm generating zones on the surface. The results are sensitive to the value of the saltation threshold, which hints at a possible feedback between saltation threshold and dust storm activity. According to this hypothesis, erodable material builds up its a result of a local process, whose effect is to lower the saltation threshold until a GDS occurs. The saltation threshold adjusts its value so that dust storms are barely able to occur.

Pankine, A. A.↗

Coniferous Canopy BRF Simulation Based on 3-D Realistic Scene

It is difficulties for the computer simulation method to study radiation regime at large-scale. Simplified coniferous model was investigate d in the present study. It makes the computer simulation methods such as L-systems and radiosity-graphics combined method (RGM) more powerf ul in remote sensing of heterogeneous coniferous forests over a large -scale region. L-systems is applied to render 3-D coniferous forest scenarios: and RGM model was used to calculate BRF (bidirectional refle ctance factor) in visible and near-infrared regions. Results in this study show that in most cases both agreed well. Meanwhiie at a tree and forest level. the results are also good.

Wang, Xin-yun↗

Analysis of the Two-Level NO PLIF Model for Low-Temperature High-Speed Flow Applications

The current work compares experimentally obtained nitric oxide (NO) laser-induced fluorescence (LIF) spectra with the equivalent spectra obtained analytically. The experimental spectra are computed from captured images of fluorescence in a gas cell and from a laser sheet passing through the fuel-air mixing flowfield produced by a high-speed fuel injector. The fuel injector is a slender strut that is currently being studied as a part of the Enhanced Injection and Mixing Project (EIMP) at the NASA Langley Research Center. This injector is placed downstream of a Mach 6 facility nozzle, which simulates the high Mach number airflow at the entrance of a scramjet combustor, and injects helium, which is used as an inert substitute for hydrogen fuel. Experimental planar (P) LIF is obtained by using a UV laser to excite fluorescence from the NO molecules that are present in either a gas cell or the facility air used for the EIMP experiments. The experimental data are obtained for several segments of the NO fluorescence spectrum. The selected segments encompass LIF lines with rotational quantum numbers appropriate for low-to-moderate temperature flows similar to those corresponding to the nominal experimental flow conditions. The experimental LIF spectra are then evaluated from the data and compared with those obtained from the theoretical models. The theoretical spectra are obtained from LIFBASE and LINUS software, and from a simplified version of the two-level fluorescence model. The equivalent analytic PLIF images are also obtained by applying only the simplified model to the results of the Reynolds-averaged simulations (RAS) of the mixing flowfield. Good agreement between the experimental and theoretical results provides increased confidence in both the simplified LIF modeling and CFD simulations for further investigations of high-speed injector performance using this approach.

Drozda, Tomasz G.↗

Analysis of the Two-Level NO PLIF Model for Low-Temperature High-Speed Flow Applications

The current work compares experimentally obtained nitric oxide (NO) laser-induced fluorescence (LIF) spectra with the equivalent spectra obtained analytically. The experimental spectra are computed from captured images of fluorescence in a gas cell and from a laser sheet passing through the fuel-air mixing flowfield produced by a high-speed fuel injector. The fuel injector is a slender strut that is currently being studied as a part of the Enhanced Injection and Mixing Project (EIMP) at the NASA Langley Research Center. This injector is placed downstream of a Mach 6 facility nozzle, which simulates the high Mach number airflow at the entrance of a scramjet combustor, and injects helium, which is used as an inert substitute for hydrogen fuel. Experimental planar (P) LIF is obtained by using a UV laser to excite fluorescence from the NO molecules that are present in either a gas cell or the facility air used for the EIMP experiments. The experimental data are obtained for several segments of the NO fluorescence spectrum. The selected segments encompass LIF lines with rotational quantum numbers appropriate for low-to-moderate temperature flows similar to those corresponding to the nominal experimental flow conditions. The experimental LIF spectra are then evaluated from the data and compared with those obtained from the theoretical models. The theoretical spectra are obtained from LIFBASE and LINUS software, and from a simplified version of the two-level fluorescence model. The equivalent analytic PLIF images are also obtained by applying only the simplified model to the results of the Reynolds-averaged simulations (RAS) of the mixing flowfield. Good agreement between the experimental and theoretical results provides increased confidence in both the simplified LIF modeling and CFD simulations for further investigations of high-speed injector performance using this approach.

Tomasz G Drozda↗

Comparison between Computational and Experimental Wall-Shear Stress and Pressure Measurements in a Supersonic Wind Tunnel

We study the performance of two selected RANS models (SST-$k$-$\omega$ and Baldwin-Lomax) to simulate the mean wall-shear and pressure data obtained from a distributed-roughness-related experimental campaign conducted at Ahmic Aerospace, LLC, located in Ohio, US. A structured compressible Navier-Stokes in-house solver, DPLR, is used to conduct the RANS simulations. As the final objective is to perform the studies under high-enthalpy flow conditions, such as inside the AHF facility at NASA-Ames Research Center (ARC), the current investigation focuses on the low-enthalpy cases using Ahmic's high-Reynolds-number supersonic wind tunnel as the initial step to simplify the flow physics. Furthermore, the wind tunnel simulations are conducted with perfectly-smooth walls to further simplify the problem and perform the initial testing of different solver-based RANS models. The Baldwin-Lomax model is observed to be more accurate than the SST model. However, when these RANS models and their roughness augmentation models are utilized to study their accuracy with respect to the experimental datasets, where the roughness effects are activated, both models perform well for a low roughness height case. However, as the roughness height increases, the SST model more accurately predicts the wall data (shear and pressure) than the Baldwin-Lomax model. Therefore, this present work has built a solid foundation for selecting the best RANS model to study the wall data under high-enthalpy flow conditions at the AHF facility soon.

Prakash Shrestha↗

The combustion program at CTR

Understanding and modeling of turbulent combustion are key problems in the computation of numerous practical systems. Because of the lack of analytical theories in this field and of the difficulty of performing precise experiments, direct numerical simulation (DNS) appears to be one of the most attractive tools to use in addressing this problem. The general objective of DNS of reacting flows is to improve our knowledge of turbulent combustion but also to use this information for turbulent combustion models. For the foreseeable future, numerical simulation of the full three-dimensional governing partial differential equations with variable density and transport properties as well as complex chemistry will remain intractable; thus, various levels of simplification will remain necessary. On one hand, the requirement to simplify is not necessarily a handicap: numerical simulations allow the researcher a degree of control in isolating specific physical phenomena that is inaccessible in experiments. CTR has pursued an intensive research program in the field of DNS for turbulent reacting flows since 1987. DNS of reacting flows is quite different from DNS of non-reacting flows: without reaction, the equations to solve are clearly the five conservation equations of the Navier Stokes system for compressible situations (four for incompressible cases), and the limitation of the approach is the Reynolds number (or in other words the number of points in the computation). For reacting flows, the choice of the equations, the species (each species will require one additional conservation equation), the chemical scheme, and the configuration itself is more complex.

Poinsot, Thierry J.↗

Integrated Process-Structure-Property Simulations for Additive Manufacturing Using the Open-Source Materialite Package

The microstructure and properties of additively manufactured (AM) metals are strongly dependent on process conditions. Therefore, process-structure-property (PSP) simulations are a useful tool for exploring process parameter space, studying process variations, and quantifying uncertainty in material properties. However, integrating process-structure and structure-property simulations often involves connecting multiple software packages. Each package may use unique data structures and require substantial domain knowledge. This presentation demonstrates PSP simulation capabilities of Materialite, an open-source package developed at NASA Langley Research Center. Materialite simplifies model linkages by using a common data structure and model interface, enabling straightforward simulation across a PSP model chain. Physics-based models, including kinetic Monte Carlo and crystal plasticity, are implemented within the package. The model interface is also intended to simplify implementation of new models and enable integration with external simulation tools. Example use cases include uncertainty quantification with PSP models and GPU-accelerated powder bed fusion AM process models.

additive manufacturing↗

Performance verification and system parameter identification of spacecraft tape recorder control servo

Design adequacy of the lead-lag compensator of the frequency loop, accuracy checking of the analytical expression for the electrical motor transfer function, and performance evaluation of the speed control servo of the digital tape recorder used on-board the 1976 Viking Mars Orbiters and Voyager 1977 Jupiter-Saturn flyby spacecraft are analyzed. The transfer functions of the most important parts of a simplified frequency loop used for test simulation are described and ten simulation cases are reported. The first four of these cases illustrate the method of selecting the most suitable transfer function for the hysteresis synchronous motor, while the rest verify and determine the servo performance parameters and alternative servo compensation schemes. It is concluded that the linear methods provide a starting point for the final verification/refinement of servo design by nonlinear time response simulation and that the variation of the parameters of the static/dynamic Coulomb friction is as expected in a long-life space mission environment.

Mukhopadhyay, A. K.↗

A simplified model for the viscous crossflow in a slotted test section

A simplified physical model was constructed which simulates the viscous crossflow in a fluid layer near the slots at a fixed streamwise location in a slotted wind tunnel. For low to moderate Reynolds numbers, numerical solutions of the two dimensional, incompressible Navier-Stokes equations in stream function and vorticity, which govern the model flow, were obtained. Fairly general slot geometry was incorporated by means of the Thompson-Thames-Mastin transformation. An approximate factorization scheme with cyclic acceleration parameters was employed to solve a finite difference analog of the stream function equation. The vorticity equation was numerically solved with a modified version of the classical alternating direction implicit scheme. Although no quantitative assessment of solution accuracy can be made, numerical results for variations in incremental wall pressure around the slat are at least qualitatively similar to some experimental results.

Cooke, C. H.↗

A simplified model for the viscous crossflow in a slotted test section

A simplified physical model is constructed which simulates the viscous crossflow in a fluid layer near the slots at a fixed streamwise location in a slotted wind tunnel. For low to moderate Reynolds numbers, numerical solutions of the two-dimensional, incompressible Navier-Stokes equations in stream function and vorticity, which govern the model flow, are obtained. Fairly general slot geometry is incorporated by means of the Thompson-Thames-Mastin transformation. An approximate factorization scheme with cyclic acceleration parameters is employed to solve a finite difference analog of the stream function equation. The vorticity equation is numerically solved with a modified version of the classical alternating direction implicit (ADI) scheme. Although no quantitative assessment of solution accuracy can be made, numerical results for variations in incremental wall pressure around the slot are at least qualitatively similar to some experimental results of Berndt and Sorenson.

Cooke, C. H.↗

Snowmelt and Infiltration Deficiencies of SSiB and Their Resolution with a New Snow-Physics Scheme

A two-year 1987-1988 integration of SSiB forced with ISLSCP Initiative I surface data (as part of the Global Soil Wetness Project, GSWP, evaluation and intercomparison) produced generally realistic land surface fluxes and hydrology. Nevertheless, the evaluation also helped to identify some of the deficiencies of the current version of the Simplified Simple Biosphere (SSiB) model. The simulated snowmelt was delayed in most regions, along with excessive runoff and lack of an spring soil moisture recharge. The SSIB model had previously been noted to have a problem producing accurate soil moisture as compared to observations in the Russian snowmelt region. Similarly, various GSWP implementations of SSIB found deficiencies in this region of the simulated soil moisture and runoff as compared to other non-SSiB land-surface models (LSMs). The origin of these deficiencies was: 1) excessive cooling of the snow and ground, and 2) deep frozen soil disallowing snowmelt infiltration. The problem was most severe in regions that experience very cold winters. In SSiB, snow was treated as a unified layer with the first soil layer, causing soil and snow to cool together in the winter months, as opposed to snow cover acting as an insulator. In the spring season, a large amount of heat was required to thaw a hard frozen snow plus deep soil layers, delaying snowmelt and causing meltwater to become runoff over the frozen soil rather than infiltrate into it.

Sud, Y. C.↗