Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “grid codes”

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 109 records · Page 6

Computational investigation of circular-to-rectangular transition ducts

A study demonstrating that flows inside transition ducts with unusual geometry can be analyzed by employing proper selections of a Navier-Stokes code, grid topology, and turbulence modeling is presented. Based on comparison between existing experimental data and the computed results for the same configurations, reasonable agreements were obtained for wall static pressure in the transition duct. Static pressure comparisons in the supersonic nozzle section were excellent, as well as agreement between computed and measured mass flow and thrust performance.

Pao, S. P.↗

Evaluation of CFD Turbulent Heating Prediction Techniques and Comparison With Hypersonic Experimental Data

Results from a study to assess the accuracy of turbulent heating and skin friction prediction techniques for hypersonic applications are presented. The study uses the original and a modified Baldwin-Lomax turbulence model with a space marching code. Grid converged turbulent predictions using the wall damping formulation (original model) and local damping formulation (modified model) are compared with experimental data for several flat plates. The wall damping and local damping results are similar for hot wall conditions, but differ significantly for cold walls, i.e., T(sub w) / T(sub t) < 0.3, with the wall damping heating and skin friction 10-30% above the local damping results. Furthermore, the local damping predictions have reasonable or good agreement with the experimental heating data for all cases. The impact of the two formulations on the van Driest damping function and the turbulent eddy viscosity distribution for a cold wall case indicate the importance of including temperature gradient effects. Grid requirements for accurate turbulent heating predictions are also studied. These results indicate that a cell Reynolds number of 1 is required for grid converged heating predictions, but coarser grids with a y(sup +) less than 2 are adequate for design of hypersonic vehicles. Based on the results of this study, it is recommended that the local damping formulation be used with the Baldwin-Lomax and Cebeci-Smith turbulence models in design and analysis of Hyper-X and future hypersonic vehicles.

Dilley, Arthur D.↗

Giant Molecular Cloud Structure and Evolution

Bodenheimer and Burkert extended earlier calculations of cloud core models to study collapse and fragmentation. The initial condition for an SPH collapse calculation is the density distribution of a Bonnor-Ebert sphere, with near balance between turbulent plus thermal energy and gravitational energy. The main parameter is the turbulent Mach number. For each Mach number several runs are made, each with a different random realization of the initial turbulent velocity field. The turbulence decays on a dynamical time scale, leading the cloud into collapse. The collapse proceeds isothermally until the density has increased to about 10(exp 13) g cm(exp -3). Then heating is included in the dense regions. The nature of the fragmentation is investigated. About 15 different runs have been performed with Mach numbers ranging from 0.3 to 3.5 (the typical value observed in molecular cloud cores is 0.7). The results show a definite trend of increasing multiplicity with increasing Mach number (M), with the number of fragments approximately proportional to (1 + M). In general, this result agrees with that of Fisher, Klein, and McKee who published three cases with an AMR grid code. However our results show that there is a large spread about this curve. For example, for M=0.3 one case resulted in no fragmentation while a second produced three fragments. Thus it is not only the value of M but also the details of the superposition of the various velocity modes that play a critical role in the formation of binaries. Also, the simulations produce a wide range of separations (10-1000 AU) for the multiple systems, in rough agreement with observations. These results are discussed in two conference proceedings.

Hollenbach, David↗

NPARC Code Upgraded with Two-Equation Turbulence Models

The National PARC (NPARC) Alliance was established by the NASA Lewis Research Center and the Air Force Arnold Engineering Development Center to provide the U.S. aeropropulsion community with a reliable Navier-Stokes code for simulating the nonrotating components of propulsion systems. Recent improvements to the turbulence model capabilities of the NPARC code have significantly improved its capability to simulate turbulent flows. Specifically, the Chien k-epsilon and Wilcox k-omega turbulence models were implemented at Lewis. Lewis researchers installed the Chien k-epsilon model into NPARC to improve the code's ability to calculate turbulent flows with attached wall boundary layers and free shear layers. Calculations with NPARC have demonstrated that the Chien k-epsilon model provides more accurate calculations than those obtained with algebraic models previously available in the code. Grid sensitivity investigations have shown that computational grids must be packed against the solid walls such that the first point off of the wall is placed in the laminar sublayer. In addition, matching the boundary layer and momentum thicknesses entering mixing regions is necessary for an accurate prediction of the free shear-layer growth.

Source record↗

Computational Analysis of a Chevron Nozzle Uniquely Tailored for Propulsion Airframe Aeroacoustics

A computational flow field and predicted jet noise source analysis is presented for asymmetrical fan chevrons on a modern separate flow nozzle at take off conditions. The propulsion airframe aeroacoustic asymmetric fan nozzle is designed with an azimuthally varying chevron pattern with longer chevrons close to the pylon. A baseline round nozzle without chevrons and a reference nozzle with azimuthally uniform chevrons are also studied. The intent of the asymmetric fan chevron nozzle was to improve the noise reduction potential by creating a favorable propulsion airframe aeroacoustic interaction effect between the pylon and chevron nozzle. This favorable interaction and improved noise reduction was observed in model scale tests and flight test data and has been reported in other studies. The goal of this study was to identify the fundamental flow and noise source mechanisms. The flow simulation uses the asymptotically steady, compressible Reynolds averaged Navier-Stokes equations on a structured grid. Flow computations are performed using the parallel, multi-block, structured grid code PAB3D. Local noise sources were mapped and integrated computationally using the Jet3D code based upon the Lighthill Acoustic Analogy with anisotropic Reynolds stress modeling. In this study, trends of noise reduction were correctly predicted. Jet3D was also utilized to produce noise source maps that were then correlated to local flow features. The flow studies show that asymmetry of the longer fan chevrons near the pylon work to reduce the strength of the secondary flow induced by the pylon itself, such that the asymmetric merging of the fan and core shear layers is significantly delayed. The effect is to reduce the peak turbulence kinetic energy and shift it downstream, reducing overall noise production. This combined flow and noise prediction approach has yielded considerable understanding of the physics of a fan chevron nozzle designed to include propulsion airframe aeroacoustic interaction effects.

Massey, Steven J.↗

Numerical Methodology for Coupled Time-Accurate Simulations of Primary and Secondary Flowpaths in Gas Turbines

Detailed information of the flow-fields in the secondary flowpaths and their interaction with the primary flows in gas turbine engines is necessary for successful designs with optimized secondary flow streams. Present work is focused on the development of a simulation methodology for coupled time-accurate solutions of the two flowpaths. The secondary flowstream is treated using SCISEAL, an unstructured adaptive Cartesian grid code developed for secondary flows and seals, while the mainpath flow is solved using TURBO, a density based code with capability of resolving rotor-stator interaction in multi-stage machines. An interface is being tested that links the two codes at the rim seal to allow data exchange between the two codes for parallel, coupled execution. A description of the coupling methodology and the current status of the interface development is presented. Representative steady-state solutions of the secondary flow in the UTRC HP Rig disc cavity are also presented.

Przekwas, A. J.↗

Application of the FUN3D Unstructured-Grid Navier-Stokes Solver to the 4th AIAA Drag Prediction Workshop Cases

FUN3D Navier-Stokes solutions were computed for the 4th AIAA Drag Prediction Workshop grid convergence study, downwash study, and Reynolds number study on a set of node-based mixed-element grids. All of the baseline tetrahedral grids were generated with the VGRID (developmental) advancing-layer and advancing-front grid generation software package following the gridding guidelines developed for the workshop. With maximum grid sizes exceeding 100 million nodes, the grid convergence study was particularly challenging for the node-based unstructured grid generators and flow solvers. At the time of the workshop, the super-fine grid with 105 million nodes and 600 million elements was the largest grid known to have been generated using VGRID. FUN3D Version 11.0 has a completely new pre- and post-processing paradigm that has been incorporated directly into the solver and functions entirely in a parallel, distributed memory environment. This feature allowed for practical pre-processing and solution times on the largest unstructured-grid size requested for the workshop. For the constant-lift grid convergence case, the convergence of total drag is approximately second-order on the finest three grids. The variation in total drag between the finest two grids is only 2 counts. At the finest grid levels, only small variations in wing and tail pressure distributions are seen with grid refinement. Similarly, a small wing side-of-body separation also shows little variation at the finest grid levels. Overall, the FUN3D results compare well with the structured-grid code CFL3D. The FUN3D downwash study and Reynolds number study results compare well with the range of results shown in the workshop presentations.

Lee-Rausch, Elizabeth M.↗

Application of a Full Reynolds Stress Model to High Lift Flows

A recently developed second-moment Reynolds stress model was applied to two challenging high-lift flows: (1) transonic flow over the ONERA M6 wing, and (2) subsonic flow over the DLR-F11 wing-body configuration from the second AIAA High Lift Prediction Workshop. In this study, the Reynolds stress model results were contrasted with those obtained from one- and two{equation turbulence models, and were found to be competitive in terms of the prediction of shock location and separation. For an ONERA M6 case, results from multiple codes, grids, and models were compared, with the Reynolds stress model tending to yield a slightly smaller shock-induced separation bubble near the wing tip than the simpler models, but all models were fairly close to the limited experimental surface pressure data. For a series of high-lift DLR{F11 cases, the range of results was more limited, but there was indication that the Reynolds stress model yielded less-separated results than the one-equation model near maximum lift. These less-separated results were similar to results from the one-equation model with a quadratic constitutive relation. Additional computations need to be performed before a more definitive assessment of the Reynolds stress model can be made.

Lee-Rausch, E. M.↗

Hardware-in-the-loop Testing of Signal Temporal Logic Frequency Control in a Microgrid

This paper presents a hardware-in-the-loop (HITL) testing and lessons learned for a signal temporal logic (STL) frequency control in a microgrid setting. In contrast to most, if not all, traditional control mechanisms, the STL control allows for including temporal constraints in the control formulation to handle sophisticated specifications from grid codes and protection relays. This feature provides richer descriptions of control specifications to address both magnitude and time simultaneously, hence a valuable tool for power system operation and control. An ultra-low voltage (24V) HITL microgrid platform consisting of two type-3 wind turbine generators and one diesel generator is used to test the proposed STL frequency control. In experimental case studies, we demonstrate that the HITL frequency responses closely match the simulated ones and strictly satisfy the designed STL control specifications.

Ferrari Maglia, Max↗

Fault Response of Distributed Energy Resources Considering the Requirements of IEEE 1547-2018

Inverter-based DER responses to faults on the electric power system are different than those of conventional generators and are often poorly understood. Inverter fault responses are largely software-defined, within physics-based constraints. DERs are regulated by standards and grid codes that also constrain their responses. Though fault response of DERs will vary widely, for any DER following any particular standard, a pattern of fault response can be synthesized by observing that particular standard. In this paper, the recently published IEEE 1547-2018 is examined to explore the general fault response of any DER that follows IEEE 1547-2018. Additionally, an IEEE 1547-2018 compliant inverter model is developed and tested in simulation to find the fault response, and that response is compared with the fault response of a commercial off-the-shelf inverter.

41 EE - Solar Energy Technologies Office (EE-4S)↗

US 2023/0182605 A1 Network constraint energy management system for electric vehicle depot charging and scheduling

Network constraint energy management system for electric vehicle (EV) depot charging and scheduling. In an embodiment, a power schedule is received from an economic dispatch application for a charging depot comprising EV charging station(s) and distributed energy resource(s). The power schedule may be simulated on a distribution network model of the charging depot, according to load flow analysis, to determine whether any grid-code violations occur. In response to the detection of violation(s), a constraint may be generated for each violating node, and the economic dispatch application may be re-executed with the constraint(s) to produce a new power schedule, until no violations are detected. When not all load demand can be satisfied by the power schedule, a charging schedule may be adjusted to ensure that critical energy requirements are satisfied. The final power and charging schedules may be used to schedule and control power generation and charging in the charging depot.

Hafiz, Faeza↗

Charging location selection of an incoming electric vehicle for large EV depot infrastructure

As the use of renewable energy resources increases within charging infrastructures, the risk that charging an additional load within the charging infrastructure will create a violation of the applicable grid code also increases. In an embodiment, a charging location is selected for an incoming load based on a load flow analysis, to thereby mitigate the risk of a violation. This selection process may also account for the charging schedule of the load and the capacity of available charging stations (e.g., fast charging vs. slow charging).

Hafiz, Faeza↗

MAG3D and its application to internal flowfield analysis

MAG3D (multiblock adaptive grid, 3D) is a 3D solution-adaptive grid generation code which redistributes grid points to improve the accuracy of a flow solution without increasing the number of grid points. The code is applicable to structured grids with a multiblock topology. It is independent of the original grid generator and the flow solver. The code uses the coordinates of an initial grid and the flow solution interpolated onto the new grid. MAG3D uses a numerical mapping and potential theory to modify the grid distribution based on properties of the flow solution on the initial grid. The adaptation technique is discussed, and the capability of MAG3D is demonstrated with several internal flow examples. Advantages of using solution-adaptive grids are also shown by comparing flow solutions on adaptive grids with those on initial grids.

Lee, K. D.↗

Euler Technology Assessment program for preliminary aircraft design employing SPLITFLOW code with Cartesian unstructured grid method

This report documents results from the Euler Technology Assessment program. The objective was to evaluate the efficacy of Euler computational fluid dynamics (CFD) codes for use in preliminary aircraft design. Both the accuracy of the predictions and the rapidity of calculations were to be assessed. This portion of the study was conducted by Lockheed Fort Worth Company, using a recently developed in-house Cartesian-grid code called SPLITFLOW. The Cartesian grid technique offers several advantages for this study, including ease of volume grid generation and reduced number of cells compared to other grid schemes. SPLITFLOW also includes grid adaptation of the volume grid during the solution convergence to resolve high-gradient flow regions. This proved beneficial in resolving the large vortical structures in the flow for several configurations examined in the present study. The SPLITFLOW code predictions of the configuration forces and moments are shown to be adequate for preliminary design analysis, including predictions of sideslip effects and the effects of geometry variations at low and high angles of attack. The time required to generate the results from initial surface definition is on the order of several hours, including grid generation, which is compatible with the needs of the design environment.

Finley, Dennis B.↗

Marshall Space Flight Center surface modeling and grid generation applications

The Solid Rocket Motors (SRM) used by NASA to propel the Space Shuttle employ gimballing nozzles as a means for vehicular guidance during launch and ascent. Gimballing a nozzle renders the pressure field of the exhaust gases nonaxisymmetric. This has two effects: (1) it exerts a torque and side load on the nozzle; and (2) the exhaust gases flow circumferentially in the aft-dome region, thermally loading the flexible boot, case-to-nozzle joint, and casing insulation. The use of CFD models to simulate such flows is imperative in order to assess SRM design. The grids for these problems were constructed by obtaining information from drawings and tabulated coordinates. The 2D axisymmetric grids were designed and generated using the EZ-Surf and GEN2D surface and grid generation codes. These 2D grids were solved using codes such as FDNS, GASP, and MINT. These axisymmetric grids were rotated around the center-line to form 3D nongimballed grids. These were then gimballed around the pivot point and the gaps or overlaps resurfaced to obtain the final domains, which contained approximately 366,000 grid points. The 2D solutions were then rotated and manipulated as appropriate for geometry and used as initial guesses in the final solution. The analyses were used in answering questions about flight criteria.

Williams, Robert W.↗

TPSAS-NF1676L-10980-DND

We plan to perform the following sets of computations on unadapted (fixed) grids: 1) Structured RANS set 1 (Code: CFL3D, Grid: Str-OnetoOne-A-v1 (supplied by HiLiftPW-1 committee), Turbulence model: Spalart-Allmaras), 2) Structured RANS set 2 (Code: CFL3D, Grid: Str-OnetoOne-A-v1 (supplied by HiLiftPW-1 committee), Turbulence model: Menter SST), 3) Structured RANS set 3 (time permitting) (Code: CFL3D, Grid: Str-OnetoOne-B-v1 (supplied by HiLiftPW-1 committee), Turbulence model: Menter SST), 4) Unstructured RANS set 1 (Code: FUN3D, Grid: Unst-Mixed-FromTet-Nodecentered-A-v1 (supplied by HiLiftPW-1 committee), Turbulence model: Spalart-Allmaras), and 5) Unstructured RANS set 2 (time permitting) (Code: FUN3D, Grid: Unst-Hex-FromOnetoOne-A-v1 (supplied by HiLiftPW-1 committee), Turbulence model: Spalart-Allmaras),. Optional case 3 is not being computed. CFL3D is a structured upwind-biased cell-centered RANS code,1 and FUN3D is an unstructured upwind-biased node-centered RANS code

Elizabeth M Lee-Rausch↗

TPSAS-NF1676L-10454-DND

We plan to perform the following sets of computations on unadapted (fixed) grids: 1) Structured RANS set 1 (Code: CFL3D, Grid: Str-OnetoOne-A-v1 (supplied by HiLiftPW-1 committee), Turbulence model: Spalart-Allmaras), 2) Structured RANS set 2 (Code: CFL3D, Grid: Str-OnetoOne-A-v1 (supplied by HiLiftPW-1 committee), Turbulence model: Menter SST), 3) Structured RANS set 3 (time permitting) (Code: CFL3D, Grid: Str-OnetoOne-B-v1 (supplied by HiLiftPW-1 committee), Turbulence model: Menter SST), 4) Unstructured RANS set 1 (Code: FUN3D, Grid: Unst-Mixed-FromTet-Nodecentered-A-v1 (supplied by HiLiftPW-1 committee), Turbulence model: Spalart-Allmaras), and 5) Unstructured RANS set 2 (time permitting) (Code: FUN3D, Grid: Unst-Hex-FromOnetoOne-A-v1 (supplied by HiLiftPW-1 committee), Turbulence model: Spalart-Allmaras),. Optional case 3 is not being computed. CFL3D is a structured upwind-biased cell-centered RANS code,1 and FUN3D is an unstructured upwind-biased node-centered RANS code

C L Rumsey↗

Transformation of two and three-dimensional regions by elliptic systems

Several reports are attached to this document which contain the results of our research at the end of this contract period. Three of the reports deal with our work on generating surface grids. One is a preprint of a paper which will appear in the journal Applied Mathematics and Computation. Another is the abstract from a dissertation which has been prepared by Ahmed Khamayseh, a graduate student who has been supported by this grant for the last two years. The last report on surface grids is the extended abstract of a paper to be presented at the 14th IMACS World Congress in July. This report contains results on conformal mappings of surfaces, which are closely related to elliptic methods for surface grid generation. A preliminary report is included on new methods for dealing with block interfaces in multiblock grid systems. The development work is complete and the methods will eventually be incorporated into the National Grid Project (NGP) grid generation code. Thus, the attached report contains only a simple grid system which was used to test the algorithms to prove that the concepts are sound. These developments will greatly aid grid control when using elliptic systems and prevent unwanted grid movement. The last report is a brief summary of some timings that were obtained when the multiblock grid generation code was run on the Intel IPSC/860 hypercube. Since most of the data in a grid code is local to a particular block, only a small fraction of the total data must be passed between processors. The data is also distributed among the processors so that the total size of the grid can be increase along with the number of processors. This work is only in a preliminary stage. However, one of the ERC graduate students has taken an interest in the project and is presently extending these results as a part of his master's thesis.

Mastin, C. Wayne↗