Search NASA⌕ Search

SEARCH · Search NASA

Results for “3-D flow”

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 19 records

A generalized hyperbolic marching method for chemically reacting 3-D supersonic flow using a splitting technique

A generalized hyperbolic marching method employing a nonorthogonal coordinate system and using a split differencing scheme for calculating steady supersonic flow over aerodynamic shapes is presented. It is a second-order-accurate mixed explicit-implicit procedure that solves the inviscid adiabatic and nondiffusive equations for chemically reacting flow in integral conservation-law form. The relationship between the integral and differential forms of the equations are examined and the relative merits of each discussed. The method admits initial Cauchy data situated on any arbitrary surface and integrates them forward along a general curvilinear coordinate, distorting and deforming the surface as it advances. The chemical kinetics term is split from the convective terms which are themselves dimensionally split, thereby freeing the fluid operators from the restricted step-size imposed by the chemical reactions and increasing the computational efficiency. The accuracy of this splitting technique is analyzed, a sufficient stability criterion is established, and comparisons are made with another method.

Rizzi, A. W.↗

Quasi-three-dimensional flow solution by meridional plane analysis

A computer program has been developed to obtain subsonic or shockfree transonic, nonviscous flow analysis on the hub-shroud mid-channel flow surface of a turbomachine. The analysis may be for any annular passage, with or without blades. The blades may be fixed or rotating and may be twisted and leaned. The flow may be axial, radial or mixed. Blade surface velocities over the entire blade are approximated based on the rate of change of angular momentum. This gives a 3-D flow picture based on a 2-D analysis. The paper discusses the method used for the program and shows examples of the type of passages and blade rows which can be analyzed. Also, some numerical examples are given to show how the program can be used for practical assistance in design of blading, annular passages, and annular diffusers.

Katsanis, T.↗

Nearly spherical constant-power detonation waves as driven by focused radiation.

Analysis is made of the flow within a 3-D spark created in a gas absorbing energy from a steady conical beam of radiation (typically obtained from an array of lasers) with nearly spherical symmetry. First-order deviations are represented by series of Legendre polynomials of the polar angle. Self-similarity is demonstrated. For harmonics 1 to 5 radial profiles of velocity components, density and pressure are presented. The linear perturbation scheme exhibits a singularity at the focus. This reflects the existence of a 'forbidden zone', as is shown by a study of the full non-linear flow equations.

George, Y. H.↗

Preparation of the Ice Certification of the Dornier 328 Regional Airliner By Numerical Simulation and By Ground Test

The Dornier 328, a new regional airliner, is to be qualified to FAR25/JAR25 requirements for operation into known icing conditions. All lifting surfaces are pneumatically deiced. Impingement limits were determined using a validated computational fluid dynamic (CFD) code. Wind tunnel tests with a model of the whole aircraft and the isolated empennage equipped with artificial ice shapes showed a degradation of handling characteristics due to ice accretion. Extensive two and three dimensional CFD calculations evaluated ice accretion at the deflection elevator horn. Icing tunnel tests were performed with a model of the horizontal tail with icing similitude fulfilled. The Dornier Do328 powerplant air induction system has successfully completed a program of icing tests in support of aircraft certification for operation into known icing conditions. Potential flight icing conditions were extensively analyzed using a CFD potential flow code and 3-D particle tracking routine to identify critical flight icing conditions. The tunnel test program verified the analytically predicted critical ice accretion surfaces of the air induction system, the adequacy of the ice protection provisions and demonstrated compliance with the applicable (JAR25) regulatory material.

D Welte↗

Effect of Surface Roughness on Dynamic Stall in Pitching Motion

Dynamic stall plays a critical role in determining the performance and stability of a wide range of fluid-dynamics systems in various engineering applications. This unsteady aerodynamic phenomenon is particularly significant for maneuvering aircraft wings, jet aircraft subjected to gust encounters, helicopter rotor blades, and wind turbine blades [1–3]. The prediction of the dynamic stall vortex (DSV) is challenging due to factors such as unsteady aerodynamics, three-dimensional (3-D) effects, turbulence and flow separation, incoming gust, and surface impact effects [4–6]. Hence, advanced computational techniques and modeling approaches in computational fluid dynamics (CFD) would be required to enhance the accuracy and reliability of DSV predictions in dynamic stall scenarios. In the past, Batther and Lee [7] employed delayed detached eddy simulations (DDES) to understand the flow physics associated with the onset of dynamic stall. Their approach demonstrated that DDES achieves results comparable to those obtained from large-eddy simulations at a reduced computational cost. In another study, Khalifa et al. [8] examined the 3-D aspects of dynamic stall on a NACA 0012 airfoil using DES solvers. In conclusion, the findings underlined the superiority of 3-D simulations over two-dimensional approaches, particularly in predicting the lift coefficient values and capturing dynamic stall stages more precisely.

97 MATHEMATICS AND COMPUTING↗

Comparison of interference-free numerical results with sample experimental data for the AEDC wall-interference model at transonic and subsonic flow conditions

Numerical results obtained from two computer programs recently developed with NASA support and now available for use by others are compared with some sample experimental data taken on a rectangular-wing configuration in the AEDC 16-Foot Transonic Tunnel at transonic and subsonic flow conditions. This data was used in an AEDC investigation as reference data to deduce the tunnel-wall interference effects for corresponding data taken in a smaller tunnel. The comparisons were originally intended to see how well a current state-of-the-art transonic flow calculation for a simple 3-D wing agreed with data which was felt by experimentalists to be relatively interference-free. As a result of the discrepancies between the experimental data and computational results at the quoted angle of attack, it was then deduced from an approximate stress analysis that the sting had deflected appreciably. Thus, the comparisons themselves are not so meaningful, since the calculations must be repeated at the proper angle of attack. Of more importance, however, is a demonstration of the utility of currently available computational tools in the analysis and correlation of transonic experimental data.

Newman, P. A.↗

Evaluating Chemical Kinetics Predictions for Propane Using 3-D and 0-D Models in a Boosted Spark-Ignited Engine

Propane has been shown to be a promising alternative fuel to reduce emissions while simultaneously achieving high efficiencies in medium- and heavy-duty engines. These high-power density applications require boosted engines which, combined with high compression ratio, can lead to auto-ignition and knock. While three-dimensional (3-D) computational fluid dynamics (CFD) models are often used for resolving the complex fluid flow in engines, these models can become computationally expensive when simulating detailed chemical kinetics. Likewise, zero-dimensional (0-D) models are computationally concise enough for kinetics development, but lack any flow-field information which governs the flame propagation processes in spark ignition (SI) engines. This work presents a comprehensive comparison between 3-D and 0-D closed cycle simulations at knocking conditions in a high compression ratio high stroke-to-bore ratio propane engine. In order to initialize the flow-field for the 3-D closed cycle (intake valve closing, (IVC) to exhaust valve opening, (EVO)) simulation, a motored multi-cycle 3-D model was run using Converge to create a map at IVC, reducing the computational time. The map allowed a non-homogeneous 3-D closed cycle simulation to be satisfactorily validated against experiments, while a homogeneous case using only the turbulence field mapping was also simulated, mimicking 0-D modeling. The 3-D simulations were used to prescribe the initial conditions (e.g., IVC thermodynamics, speciation, burn-rate profile) for a 2-zone 0-D SI engine model in Chemkin Pro for both cases. It was found that 2-zone 0-D modeling underpredicted the knock onset timing, likely due to the lack of thermal stratification in the unburned gas region. Future work will carry multi-zone 0-D modeling to capture the fuel auto-ignition in the unburned region.

Douvry-Rabjeau, Julien [Oakland University, Roches↗

Aerodynamic heating on 3-D bodies including the effects of entropy-layer swallowing

A relatively simple method was developed previously (authors, 1973) for calculating laminar, transitional, and turbulent heating rates on three-dimensional bodies in hypersonic flows. This method was shown to yield reasonably accurate results for laminar heating on blunted circular and elliptical cones and an earlier version of the space shuttle vehicle. As the boundary layer along the surface grows, more and more of the inviscid-flow mass is entrained into the boundary layer, and the streamlines which passed through the nearly normal portion of the bow shock wave are 'swallowed' by the boundary layer. This phenomenon is often referred to as entropy-layer or streamline swallowing, and it can have a significant effect on the calculated heating rates. An approximate, yet simple, method for including the effects of entropy-layer swallowing in the heating-rate calculations is given.

Dejarnette, F. R.↗

Characterization of Most Promising Sequestration Formations in the Rocky Mountain Region

The project Characterization of Most Promising Sequestration Formations in the Rocky Mountain Region is one of 9 site characterization projects that were implemented as part of ARRA (American Recovery and Reinvestment Act). Data from this project was used to improve resolution of data in NATCARB in the area of study. Data related to this study has already been incorporated in NATCARB Atlas. The Rocky Mountain Carbon Capture and Storage (RMCCS) project investigated multiple geologic formations and characterized a local site on the Colorado Plateau for future CCS opportunities. The RMCCS project focused on the Cretaceous Dakota, Jurassic Entrada, and Pennsylvanian Weber Sandstones, the three largest regional formations. All formations in this project are potential CO2 storage resources for future power plants, natural gas processing plants, cement plants, and oil shale development projects. The area adjacent to Craig, Colorado, (Sand Wash Basin) was the area selected for detailed geologic characterization on the RMCCS project. The basin was selected in part because the geology can be extrapolated to other sites on the Colorado Plateau. Field mapping and seismic surveys were conducted to identify and evaluate the basin's structural configuration. A 9,745-foot deep characterization well was drilled to collect 131 feet of core and a suite of geophysical well log data. Petrophysical tests on samples of core were used to calibrate geophysical log data, which can be used to obtain storage resource estimates and evaluate associated uncertainty as well as simulate the hydrologic behavior of injected CO2. A detailed analysis of the primary formations (Dakota, Entrada and Weber sandstones) yielded a more accurate CO2 storage resource assessment for these formations within the Colorado Plateau; RMCCS estimates indicate a total CO2 storage resource of more than 38,000 million metric tons. The characterization of the Sand Wash Basin (2-D seismic surveys, multiple well logs and lithological, petrophysical and geochemical analyses) allowed for a detailed 3-D model to be constructed. The model served as the framework for analyses ranging from CO2 storage resource, injectivity, and subsurface flow to uncertainty estimates to evaluation of risk.

2-D seismic↗

Space shuttle orbiter reentry flow field and heating analysis

The reentry flow field and thermal environment around the straight wing shuttle orbiter vehicle were determined. Both rarefied and continuum flow fields and associated heating rates on various configurations representative of the orbiter at high angle of attack were calculated. Rarefied flow fields and heating rates were computed by the Monte Carlo direct simulation technique for altitudes above 82.3 km. Continuum inviscid flow fields were calculated by 2-D unsteady and 3-D steady finite difference/artificial viscosity methods and also by a 2-D shock layer analysis technique. Viscous flow fields and heating rates in the continuum regime were computed by a boundary layer integral matrix method for laminar flow and by an aerodynamic surface heating technique for turbulent flow. Shapes considered in the study included flat plates (representing the underside of the orbiter fuselage or the wing MAC), orbiter fuselage cross sections, orbiter wing airfoils, and 3-D orbiter configurations, all at high angle of attack (40 - 60 deg). The theoretical results showed good agreement with measured pressure and heat transfer data.

Rochelle, W. C.↗

Free-spray characteristics and spray-wall interactions of methanol on a gasoline direct injector under flash-boiling and non-flash-boiling conditions

Methanol is considered a promising alternative fuel for internal combustion engines (ICEs) due to its high-octane number, fast laminar flame speed, and elevated latent heat of vaporization, all of which support higher compression ratios and improved thermal efficiency. However, its substantial latent heat of vaporization also poses cold-start challenges, such as misfire and fuel film deposition. This study aims to investigate methanol spray morphology and spray-wall interaction using the Spray M injector from the Engine Combustion Network within a constant-pressure flow vessel. A recently developed unified numerical framework capable of modeling both flash and non-flash boiling sprays is validated against experimental liquid volume fraction data acquired via 3-D computed tomography. Here, the results reveal that flash boiling significantly alters the spray morphology, leading to smaller droplets and spray collapse due to enhanced air-entrainment-induced turbulence. Quantitative agreement between experiments and simulations confirms this behavior. Coupled 0-D equilibrium and 3-D computational fluid dynamics analyses show that flash boiling accelerates evaporation and reduces fuel residence time, while non-flash conditions maintain a persistent liquid core more susceptible to wall wetting. Wall temperature diagnostics reveal that spray collapse alters heat transfer patterns by shifting cooling effects. Mixture fraction analysis indicates that evaporation is primarily governed by shear-layer turbulence, though deviations from adiabatic equilibrium mixing emerge under low-turbulence conditions. Finally, increasing fuel, ambient, and wall temperatures reduces wall wetting and film thickness, mitigating cold-start risks. These findings enhance the understanding of methanol sprays’ behavior and support its adoption as a viable, alternative fuel for ICEs.

Engine Combustion Network↗

Propelling sustainable energy: Multi-omics analysis of pennycress FATTY ACID ELONGATION1 knockout for biofuel production

Abstract The aviation industry’s growing interest in renewable jet fuel has encouraged the exploration of alternative oilseed crops. Replacing traditional fossil fuels with a sustainable, domestically sourced crop can substantially reduce carbon emissions, thus mitigating global climate instability. Pennycress (Thlaspi arvense L.) is an emerging oilseed intermediate crop that can be grown during the offseason between maize (Zea mays) and soybean (Glycine max) to produce renewable biofuel. Pennycress is being domesticated through breeding and mutagenesis, providing opportunities for trait enhancement. Here, we employed metabolic engineering strategies to improve seed oil composition and bolster the plant's economic competitiveness. FATTY ACID ELONGATION1 (FAE1) was targeted using CRISPR-Cas 9 gene editing to eliminate very long chain fatty acids (VLCFAs) from pennycress seed oil, thereby enhancing its cold flow properties. Through an integrated multiomics approach, we investigated the impact of eliminating VLCFAs in developing and mature plant embryos. Our findings revealed improved cold-germination efficiency in fae1, with seedling emergence occurring up to 3 d earlier at 10 °C. However, these alterations led to a tradeoff between storage oil content and composition. Additionally, these shifts in lipid biosynthesis were accompanied by broad metabolic changes, such as the accumulation of glucose and ADP-glucose quantities consistent with increased starch production. Furthermore, shifts to shorter FA chains triggered the upregulation of heat shock proteins, underscoring the importance of VLCFAs in stress signaling pathways. Overall, this research provides crucial insights for optimizing pennycress seed oil while preserving essential traits for biofuel applications.

Rasoul, Amira (ORCID:0009000714290175)↗

Numerical Investigation of Fluid Flow and Space Charge in Liquid Argon Time Projection Chamber (LArTPC) Detectors

Overview This project focused on developing a high-fidelity numerical framework to simulate the multiphysics environment within Liquid Argon Time Projection Chamber (LArTPC) detectors. The primary objective was to characterize the complex interplay between ion transport, background fluid dynamics, and electric field distortions—a critical factor for the calibration and sensitivity of next-generation High Energy Physics experiments, such as DUNE. Technical Achievements The research successfully yielded a hybrid numerical space-charge solver utilizing a Cell-Centered Finite Volume Method (FVM) for ion transport coupled with a Finite Element Method (FEM) for electric potential. Key accomplishments include: • Verification & Validation: The 3-D solver was rigorously verified against 1-D analytical solutions, demonstrating high numerical accuracy in predicting space-charge-induced field deviations. • Field Distortion Analysis: 3D simulations revealed that space charge effects introduce significant non-uniformities in the electric field. Critically, the research identified that background LAr flow velocities, when comparable to ion drift velocities, markedly exacerbate these distortions. • Technology Transfer: The resulting source code and comprehensive user manuals were successfully transferred to collaborators at Fermilab, providing a portable computational tool for the broader scientific community. Challenges and Future Directions While the space-charge solver achieved all performance metrics, the integrated fluid dynamics modeling encountered convergence challenges stemming from the extreme 200-fold disparity in length scales between the detector's 37 mm inlet pipes and the 8-meter global domain. To address this, the project has identified a clear technical pivot toward Hierarchical Geometric Adaptive Mesh Refinement (HG-AMR). By implementing an h-type refinement strategy with hanging nodes, future iterations of this solver will be capable of resolving localized high-gradient inlet flows without the prohibitive computational costs of regular grids. This advancement, combined with data-driven uncertainty quantification based on MicroBooNE-style calibration, will enable the precise modeling of detector responses in large-scale cryogenic environments where direct measurement remains difficult. Impact The computational tools developed under this award provide a foundation for enhancing the energy resolution and spatial reconstruction of noble liquid detectors. By bridging the gap between theoretical fluid dynamics and experimental field calibration, this work supports the DOE’s mission to advance the frontiers of neutrino physics and dark matter detection.

42 ENGINEERING↗

Improving High Resolution Offshore Wind Resource Assessments and Forecasts using Observations in the MA/RI Lease Areas

The third Wind Forecast Improvement Project (WFIP3) sought to improve understanding of the physical phenomena in the atmosphere and ocean that dictate the structure and variability of wind and thermodynamic fields within the Marine Atmospheric Boundary Layer (MABL). WFIP3 conducted a comprehensive 18-month observational study over the Northeast U.S. outer continental shelf, a high use coastal zone, using a 3-D multiscale sensor array to highly resolve the temporal, vertical, and horizontal structure of the coupled atmospheric and oceanic boundary layers. Multiple land-based study sites adjacent to the coastal ocean observed surface meteorology and vertical profiles of atmospheric properties via passive infrared and microwave radiometers, active lidars and radars, and radiosondes. At sea, an array of surface flux buoys and two vertical profiling lidar buoys observed both atmospheric and oceanic properties, augmented by land-based oceanographic radar systems and routine ship-based surveys. Intensive observations of the MABL over the ocean was done from an air-sea interaction flux tower and extended deployments of a large barge platform. WFIP3 focused on mesoscale and sub-mesoscale flows -- including sea breezes, low-level jets, low-level clouds, and coastal storms -- and the ability of advanced numerical model parameterizations to represent them within fully coupled oceanic and atmospheric modeling systems and foundational weather forecast models. Numerous critical forecasting phenomena were observed, however the project was terminated prior to the completion of the field observational period and the analysis period.

54 ENVIRONMENTAL SCIENCES↗

Infrasonic directivity of monopole, dipole and bipole ground-surface reflected sources

Infrasound (acoustic waves below 20 Hz) can be used to detect, locate and quantify activity in the atmosphere such as volcanic eruptions and anthropogenic explosions. Attempts to quantify volcanic eruption parameters such as exit velocity, plume height and mass flow rate using infrasound data depend strongly on assumptions of the acoustic source type. Infrasonic sources may produce omnidirectional or directional wavefields, while propagation effects, such as interaction with topography, can induce further wavefield directivity that is measured by field instrumentation. Limited sampling of these wavefields can hinder our ability to infer the underlying source, and thus our understanding of the eruption characteristics. Equivalent sources are often used to represent acoustic source mechanisms and resultant wavefields. In this study, we review equivalent acoustic sources as they pertain to infrasonic scale and wavelengths commonly encountered in very local (⁠<5 km range) geophysical field deployments. We highlight the equivalent infrasonic bipole source that can be induced by ground-reflection of an elevated monopole; we are not aware of any prior infrasound studies that use the bipole source concept. We use analytical and numerical methods to explore source directivity of monopole, dipole and bipole ground-reflected sources at infrasonic frequencies as well as the additional directivity complications introduced by interactions with topography. We illustrate that for typical volcano-infrasound wavelengths, increasing height above the ground as well as increasing source frequency leads to increased wavefield directivity. Numerical modelling using a simple omnidirectional monopole source embedded in topography further illustrates that both horizontal and vertical infrasound directionality can be induced by topography at the distance scales appropriate for local volcano infrasound monitoring. Information summarized in this analytical and numerical exploration of infrasound directivity may be used to help guide future volcano-infrasound field deployments intended to estimate source parameters or quantify wavefield directivity. Analytic solutions for simple whole-space or half-space atmospheres provide useful formulations for planning or initially analysing geophysical field-scale experimental data; however, especially at very local distances from the source (⁠<5 km), 3-D simulations are necessary to account for complex topography commonly encountered in volcano-infrasound applications.

Infrasound↗

An experimentally informed design process for future inertial confinement fusion facilities

The achievement of ignition in the laboratory has renewed interest in defining the requirements for a future high-gain inertial confinement fusion (ICF) facility. Our best chance of predicting future ICF performance is with 3-D radiation hydrodynamic simulations that have been benchmarked against experimental data, but their high computational cost is prohibitive for use in practical design studies. We introduce a hierarchical approach where 3-D simulations are tuned to match experimental measurements and used to train 3-D degradation models in 1-D simulations allowing for accurate predictions over the entire OMEGA direct-drive database. A genetic algorithm was used in combination with the trained 1-D simulations to search for optimal direct-drive implosion designs at driver energies ranging from 20 kJ to 10 MJ. As the fidelity of 3-D codes improves, this approach will provide a viable experimentally informed tool for defining the next ICF facility.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Multiphysics Analysis of Li Cooled Divertor Substrate During Loss of Coolant Accident (LOCA)

In the ongoing study of potential designs for liquid metal (LM) plasma-facing components (PFCs), so-called “slow” and “fast” Li flow divertor concepts are under investigation. In the previous studies on design and analysis of the slow Li flow divertor and comparison with the fast Li flow divertor, the magnetohydrodynamics (MHD)/heat transfer effects of the Li flowing inside the substrate as a second coolant were comprehensively investigated under the normal steady-state operation conditions. Here, in the present study, the multiphysics analysis is extended to the unsteady abnormal divertor scenario where the Li layer on top of the substrate does not provide full coverage or even totally disappears for a certain period of time, so that the substrate becomes directly exposed to the incident high plasma heat flux. Such an unwanted event may happen regardless of the concept of the divertor and is worth detailed investigations, typically referred to as a loss of coolant accident (LOCA). To address this situation, a simplified scoping analysis is conducted first in 2-D, and then an integrated 3-D modeling is performed using a time-dependent multiphysics model in COMSOL Multiphysics that integrates LM MHD, heat transfer, and solid mechanics. The main goal is to evaluate conditions under which the major material limits, such as the maximum allowable temperature, stress, and displacement of the substrate, can still be met. It was shown that the maximum time over which the substrate of RAFM steel can retain structural integrity during the LOCA is around 0.2 ∼ 0.3 s. Any divertor concept that utilizes RAFM steel as a substrate material and liquid Li as a second coolant should take such a permitted time into consideration.

divertor↗

Physical properties, internal structure, and the three‐dimensional petrography of CI chondrites

physical properties and the nature of their breccation, we investigated nine samples of the Ivuna and Orgueil CI chondrites ranging in size from 1 mm to 4 cm in approximate diameter. The combined mass of unique material investigated in this work is 113 g. For our investigations, we use ideal gas pycnometry, 3-D laser scanning, x-ray computed microtomography (μCT), and accompanying digital data extraction techniques. We found that the bulk density of the samples ranged from 1.61 to 2.10 g cm −3 . Larger samples tend to have a lower bulk density. Grain density (ranging from 2.44 to 2.55 g cm −3 ) is significantly less variable than the bulk density in our samples and the quantity of porosity (ranging from 14.6% to 33.8%) is the dominant factor in determining the bulk density of CI chondrite material. Our μCT results show that the visible porosity across all sizes of our CI chondrite samples is in the form of cracks, but these cracks can account for less than two-thirds of the porosity in the CI chondrites. Other porosity is not visible, even at μCT resolutions of 2.7 μm voxel edge −1 and we conclude that it is sub-micron in nature. It is not clear if the cracks seen in our samples are indigenous to the chondrites or are a result of terrestrial processes. We also find that the CI chondrites are excellent examples of the fractal-like nature of brecciation, where clasts can be observed at all scales we imaged. The breccias are composed of sub-equant-shaped and sub-rounded-textured clasts like melt-free impact breccias on other solar system bodies. From our μCT volume and digital data extraction, we determine that the Ivuna CI chondrite breccia is organized: the mostly sub-equant clasts within our ~2 cm chunk of Ivuna have a mean diameter of 1.33 mm and their aligned longest axes define a lineation structure. We speculate that the lineation was imparted after fragmentation of the clasts by slight shear on the parent asteroid which could be the result of seismic-related granular flow or mild non-axial impact-related compaction. These data will help to place returned asteroidal material from asteroids 162173 Ryugu and 101955 Bennu and the CI chondrites into a mutual geological context.

CI chondrite↗