Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “flow 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.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Constraint energy minimizing generalized multiscale finite element method for multi-continuum Richards equations

In fluid flow simulation, the multi-continuum model is a useful strategy. When the heterogeneity and contrast of coefficients are high, the system becomes multiscale, and some kinds of reduced order methods are demanded. Combining these techniques with nonlinearity, we will consider in this paper a dual-continuum model which is generalized as a multi-continuum model for a coupled system of nonlinear Richards equations as unsaturated flows, in complex heterogeneous fractured porous media; and we will solve it by a novel multiscale approach utilizing the constraint energy minimizing generalized multiscale finite element method (CEM-GMsFEM). In particular, such a nonlinear system will be discretized in time and then linearized by Picard iteration (whose global convergence is proved theoretically). Subsequently, we tackle the resulting linearized equations by the CEM-GMsFEM and obtain proper offline multiscale basis functions to span the multiscale space (which contains the pressure solution). More specifically, we first introduce two new sources of samples, and the GMsFEM is used over each coarse block to build local auxiliary multiscale basis functions via solving local spectral problems, that are crucial for detecting high-contrast channels. Second, per oversampled coarse region, local multiscale basis functions are created through the CEM as constrainedly minimizing an energy functional. Various numerical tests for our approach reveal that the error converges with the coarse-grid size and that only few oversampling layers as well as basis functions are needed.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1991

Pressure Fluctuations In Simulated Turbulent Channel Flow

Report presents study of structures of fluctuations in pressure in numerically simulated incompressible, turbulent channel flow. Based on statistical analyses of pressure and velocity fields, part of continuing effort to understand relationships among pressures, velocities, vorticities, kinetic energies, and other features of turbulent flows in vicinities of walls.

Kim, John↗
DOE OSTI2022

PeleC: An adaptive mesh refinement solver for compressible reacting flows

Reacting flow simulations for combustion applications require extensive computing capabilities. Leveraging the AMReX library, the Pele suite of combustion simulation tools targets the largest supercomputers available and future exascale machines. We introduce PeleC, the compressible solver in the Pele suite, and detail its capabilities, including complex geometry representation, chemistry integration, and discretization. We present a comparison of development efforts using both OpenACC and AMReX’s C++ performance portability framework for execution on multiple GPU architectures. We discuss relevant details that have allowed PeleC to achieve high performance and scalability. PeleC’s performance characteristics are measured through relevant simulations on multiple supercomputers. The success of PeleC’s design for exascale is exhibited through demonstration of a 160 billion cell simulation and weak scaling onto 100% of Summit, an NVIDIA-based GPU supercomputer at Oak Ridge National Laboratory. Our results provide confidence that PeleC will enable future combustion science simulations with unprecedented fidelity.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1989

Thermocapillary convection in zone-melting crystal growth - An open-boat physical simulation

Thermocapillary convection in a molten zone of NaNO3 contained in a boat with a free horizontal surface, that is heated from above by a centered wire heater, was studied to simulate flow in zone-melting crystal growth. Using a laser-light-cut technique and fine SiO powder as a tracer, convection in the melt zone was visualized in two different cases. In the first case, the entire melt surface was free, while in the second the melt surface was free only in the immediate vicinity of one vertical wall and was covered elsewhere, this wall being to simulate the melt/crystal interface during crystal growth. It was observed that thermocapillary convection near this wall prevailed in the first case, but was reduced significantly in the second. Since thermocapillary rather than natural convection dominated in the melt, the effect of the partial covering of the melt surface on thermocapillary convection in the melt observed in this study is expected to be similar under microgravity.

Kim, Y. J.↗
NASA NTRS1995

Tracking a turbulent spot in an immersive environment

We describe an interactive, immersive 3D system called Tracktur, which allows a viewer to track the development of a turbulent flow. Tracktur displays time-varying vortex structures extracted from a numerical flow simulation. The user navigates the space and probes the data within a windy 3D landscape. In order to sustain a constant frame rate, we enforce a fixed polygon budget on the geometry. We supplement the interactive system with interactive, hypertext documentation available on the World Wide Web. In actual use by a fluid dynamicist, the system has yielded new insights into the transition to turbulence in a laminar flow.

Banks, David C.↗
NASA NTRSDate not supplied

Verification of Anisotropic Mesh Adaptation for Complex Aerospace Applications

The stabilized finite element solver, FUN3D-SFE, along with the grid mechanics package refine are verified for aerospace applications of laminar and turbulent flow simulations. The current verification exercise represents an extension of previous research using FUN3D-SFEwith adjoint-based mesh adaptation to generate highly anisotropic adapted meshes for inviscid problems. Adaptations are performed using a solution-based approach that controls the Lpnorm of Mach number interpolation error and an adjoint-based approach that controls the error in some output functional. Adaptive results are shown for laminar subsonic flow over a delta wing, laminar subsonic flow over ONERA M6 wing, inviscid supersonic flow over a sonic boom test case, and a turbulent flow over a high lift configuration (JAXA Standard Model). Mesh convergence results are also compared with results from FUN3D-FV (Finite Volume) whenever available. For all test cases considered, FUN3D-SFE gives significantly better accurate results than FUN3D-FV on coarse meshes.

Aravind Balan↗
DOE OSTI2023

Interinstitutional Study of the New EURO-GANEX Process Resistance by Gamma Irradiation Test Loops

As part of the homogeneous actinides recycling strategy, the EURO-GANEX process is one of the most promising options to achieve the goal of minor actinides recovery. However, EURO-GANEX also has various drawbacks that are being addressed. Improvements made to EURO-GANEX system have resulted in the emergence of the so-called New EURO-GANEX process, where the composition of the solvent has been modified by replacing TODGA and DMDOHEMA with cis-mTDDGA in the organic phase and SO3-Ph-BTP with PyTri-Diol in the aqueous phase in order to resolve important issues. The objective of this work is twofold: evaluate the gamma radiolysis resistance of the new EURO-GANEX process by dynamic irradiation conditions simulating the three main steps of the process and validate the design of CIEMAT Náyade, CEA Marcel and INL irradiation loop devices since each of them mimics different aspects of the real process. Náyade and the INL loops could irradiate together the organic and aqueous phases, whereas in the CEA loop, the irradiated solvent is recycled continuously inside a platform with several stages of mixer-settlers containing aqueous flows simulating the three main steps of the process. The extraction performances and changes in the composition of the solvent have been analysed during the irradiation experiment by different techniques: gamma spectrometry and ICP-MS/OES for cations or tracers extraction, and HPLC-MS to quantify the degradation and evaluate the degradation compounds. Additionally, in spite of some differences between the three-irradiation facilities, this inter-institutional study shows that they are three comparative tools, which provide results on the stability towards radiolysis of a liquid-liquid extraction system. Favourable extraction results for the different steps are obtained according to the static irradiation studies found in literature. However, the degradation of cis-mTDDGA is appreciable leading to degradation compounds, some of which form precipitates and produce important changes in viscosity, important aspects that must be addressed prior to the successful industrial application of the new EURO-GANEX process.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
NASA NTRS2001

Computation of Coupled Thermal-Fluid Problems in Distributed Memory Environment

The thermal-fluid coupling problems are very important to aerospace and engineering applications. Instead of analyzing heat transfer and fluid flow separately, this study merged two well-accepted engineering solution methods, SINDA for thermal analysis and FDNS for fluid flow simulation, into a unified multi-disciplinary thermal fluid prediction method. A fully conservative patched grid interface algorithm for arbitrary two-dimensional and three-dimensional geometry has been developed. The state-of-the-art parallel computing concept was used to couple SINDA and FDNS for the communication of boundary conditions through PVM (Parallel Virtual Machine) libraries. Therefore, the thermal analysis performed by SINDA and the fluid flow calculated by FDNS are fully coupled to obtain steady state or transient solutions. The natural convection between two thick-walled eccentric tubes was calculated and the predicted results match the experiment data perfectly. A 3-D rocket engine model and a real 3-D SSME geometry were used to test the current model, and the reasonable temperature field was obtained.

Wei, H.↗
NASA NTRS2004

Preparing and Analyzing Iced Airfoils

SmaggIce version 1.2 is a computer program for preparing and analyzing iced airfoils. It includes interactive tools for (1) measuring ice-shape characteristics, (2) controlled smoothing of ice shapes, (3) curve discretization, (4) generation of artificial ice shapes, and (5) detection and correction of input errors. Measurements of ice shapes are essential for establishing relationships between characteristics of ice and effects of ice on airfoil performance. The shape-smoothing tool helps prepare ice shapes for use with already available grid-generation and computational-fluid-dynamics software for studying the aerodynamic effects of smoothed ice on airfoils. The artificial ice-shape generation tool supports parametric studies since ice-shape parameters can easily be controlled with the artificial ice. In such studies, artificial shapes generated by this program can supplement simulated ice obtained from icing research tunnels and real ice obtained from flight test under icing weather condition. SmaggIce also automatically detects geometry errors such as tangles or duplicate points in the boundary which may be introduced by digitization and provides tools to correct these. By use of interactive tools included in SmaggIce version 1.2, one can easily characterize ice shapes and prepare iced airfoils for grid generation and flow simulations.

Vickerman, Mary B.↗
DOE OSTI2021

DOE-GENIORS Dynamic Irradiation Testing

As part of the homogeneous actinides recycling strategy, the EURO-GANEX process is one of the most promising options to achieve the goal of minor actinides recovery. Improvements made to EURO-GANEX system have resulted in the emergence of the so-called New EURO-GANEX process, where the composition of the solvent has been modified by replacing TODGA and DMDOHEMA with cis-mTDDGA in the organic phase and SO3-Ph-BTP with PyTri-Diol in the aqueous phase in order to resolve important issues. The objective of this work is twofold: evaluate the gamma radiolytic resistance of the new EURO-GANEX process by dynamic irradiation conditions simulating the three main steps of the process, and validate the design of CIEMAT Náyade, CEA Marcel and INL irradiation loop devices since each of them mimics different aspects of the real process. Náyade and the INL loops could irradiate the organic and aqueous phases together, whereas in the CEA loop, the irradiated solvent is recycled continuously inside a platform with several stages of mixer-settlers containing aqueous flows simulating the three main steps of the process. The extraction performances and changes in the composition of the solvent have been analyzed during the irradiation experiment by different techniques: gamma spectrometry and ICP-MS/OES for cations or radioactive tracer extraction, and HPLC-MS to identify and quantify the degradation compounds. Despite some differences between the three-irradiation facilities, this inter-institutional study shows that these three comparative tools provide similar trends in the radiolytic stability of a liquid-liquid extraction system. Favourable extraction results for the different steps are obtained according to the static irradiation studies found in literature. However, the degradation of cis-mTDDGA is appreciable leading to degradation compounds, some of which form precipitates and produce important changes in viscosity, important aspects that must be addressed prior to the successful industrial application of the new EURO-GANEX process.

38 - RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCL↗
DOE OSTI2017

Regulation-Structured Dynamic Metabolic Model Provides a Potential Mechanism for Delayed Enzyme Response in Denitrification Process

In a recent study of denitrification dynamics in hyporheic zone sediments, we observed a significant time lag (up to several days) in enzymatic response to the changes in substrate concentration. To explore an underlying mechanism and understand the interactive dynamics between enzymes and nutrients, we developed a trait-based model that associates a community’s traits with functional enzymes, instead of typically used species guilds (or functional guilds). This enzyme-based formulation allows to collectively describe biogeochemical functions of microbial communities without directly parameterizing the dynamics of species guilds, therefore being scalable to complex communities. As a key component of modeling, we accounted for microbial regulation occurring through transcriptional and translational processes, the dynamics of which was parameterized based on the temporal profiles of enzyme concentrations measured using a new signature peptide-based method. The simulation results using the resulting model showed several days of a time lag in enzymatic responses as observed in experiments. Further, the model showed that the delayed enzymatic reactions could be primarily controlled by transcriptional responses and that the dynamics of transcripts and enzymes are closely correlated. The developed model can serve as a useful tool for predicting biogeochemical processes in natural environments, either independently or through integration with hydrologic flow simulators.

trait-based modeling↗
NASA NTRS2020

A strongly changing accretion morphology during the outburst decay of the neutron star x-ray binary 4U 1608-52

It is commonly assumed that the properties and geometry of the accretion flow in transient low-mass X-ray binaries (LMXBs) significantly change when the X-ray luminosity decays below ∼10−2 of the Eddington limit (LEdd). However, there are few observational cases where the evolution of the accretion flow is tracked in a single X-ray binary over a wide dynamic range. In this work, we use NuSTAR and NICER observations obtained during the 2018 accretion outburst of the neutron star LMXB 4U 1608−52, to study changes in the reflection spectrum. We find that the broad Fe–Kα line and Compton hump, clearly seen during the peak of the outburst when the X-ray luminosity is ∼1037 erg s−1 (∼0.05 LEdd), disappear during the decay of the outburst when the source luminosity drops to ∼4.5 × 1035 erg s−1 (∼0.002 LEdd). We show that this non-detection of the reflection features cannot be explained by the lower signal-to-noise ratio at lower flux, but is instead caused by physical changes in the accretion flow. Simulating synthetic NuSTAR observations on a grid of inner disc radius, disc ionization, and reflection fraction, we find that the disappearance of the reflection features can be explained by either increased disc ionization (log ξ  4.1) or a much decreased reflection fraction. A changing disc truncation alone, however, cannot account for the lack of reprocessed Fe–Kα emission. The required increase in ionization parameter could occur if the inner accretion flow evaporates from a thin disc into a geometrically thicker flow, such as the commonly assumed formation of a radiatively inefficient accretion flow at lower mass accretion rates.

J. van den Eijnden↗
NASA NTRS1979

Computational aerodynamics development and outlook /Dryden Lecture in Research for 1979/

Some past developments and current examples of computational aerodynamics are briefly reviewed. An assessment is made of the requirements on future computer memory and speed imposed by advanced numerical simulations, giving emphasis to the Reynolds averaged Navier-Stokes equations and to turbulent eddy simulations. Experimental scales of turbulence structure are used to determine the mesh spacings required to adequately resolve turbulent energy and shear. Assessment also is made of the changing market environment for developing future large computers, and of the projections of micro-electronics memory and logic technology that affect future computer capability. From the two assessments, estimates are formed of the future time scale in which various advanced types of aerodynamic flow simulations could become feasible. Areas of research judged especially relevant to future developments are noted.

Chapman, D. R.↗
NASA NTRS1980

Numerical method for predicting flow characteristics and performance of nonaxisymmetric nozzles. Part 2: Applications

A computer implemented numerical method for predicting the flow in and about an isolated three dimensional jet exhaust nozzle is summarized. The approach is based on an implicit numerical method to solve the unsteady Navier-Stokes equations in a boundary conforming curvilinear coordinate system. Recent improvements to the original numerical algorithm are summarized. Equations are given for evaluating nozzle thrust and discharge coefficient in terms of computed flowfield data. The final formulation of models that are used to simulate flow turbulence effect is presented. Results are presented from numerical experiments to explore the effect of various quantities on the rate of convergence to steady state and on the final flowfield solution. Detailed flowfield predictions for several two and three dimensional nozzle configurations are presented and compared with wind tunnel experimental data.

Thomas, P. D.↗
NASA NTRS1984

Elements of computational engine-airframe integration

Computerized flow simulation methods are being developed to refine the traditional engine-airframe integration process. These analyses reduce design risk, minimize wind tunnel and flight test requirements, and provide detailed design information not available by any other means. This paper describes computational tools which have matured sufficiently to permit project applications. Several applications are combined with research studies to illustrate the prediction of high-speed nonlinear phenomena associated with engine inlet flow rate, nacelle shape and position, afterbody shape and integration, propeller slipstream interactions, airframe-induced inlet flow, and pylon integration. Both fighter and transport configurations are used to highlight current capabilities and future requirements.

Boppe, C. W.↗
NASA NTRS2006

Single- and Multiple-Objective Optimization with Differential Evolution and Neural Networks

Genetic and evolutionary algorithms have been applied to solve numerous problems in engineering design where they have been used primarily as optimization procedures. These methods have an advantage over conventional gradient-based search procedures became they are capable of finding global optima of multi-modal functions and searching design spaces with disjoint feasible regions. They are also robust in the presence of noisy data. Another desirable feature of these methods is that they can efficiently use distributed and parallel computing resources since multiple function evaluations (flow simulations in aerodynamics design) can be performed simultaneously and independently on ultiple processors. For these reasons genetic and evolutionary algorithms are being used more frequently in design optimization. Examples include airfoil and wing design and compressor and turbine airfoil design. They are also finding increasing use in multiple-objective and multidisciplinary optimization. This lecture will focus on an evolutionary method that is a relatively new member to the general class of evolutionary methods called differential evolution (DE). This method is easy to use and program and it requires relatively few user-specified constants. These constants are easily determined for a wide class of problems. Fine-tuning the constants will off course yield the solution to the optimization problem at hand more rapidly. DE can be efficiently implemented on parallel computers and can be used for continuous, discrete and mixed discrete/continuous optimization problems. It does not require the objective function to be continuous and is noise tolerant. DE and applications to single and multiple-objective optimization will be included in the presentation and lecture notes. A method for aerodynamic design optimization that is based on neural networks will also be included as a part of this lecture. The method offers advantages over traditional optimization methods. It is more flexible than other methods in dealing with design in the context of both steady and unsteady flows, partial and complete data sets, combined experimental and numerical data, inclusion of various constraints and rules of thumb, and other issues that characterize the aerodynamic design process. Neural networks provide a natural framework within which a succession of numerical solutions of increasing fidelity, incorporating more realistic flow physics, can be represented and utilized for optimization. Neural networks also offer an excellent framework for multiple-objective and multi-disciplinary design optimization. Simulation tools from various disciplines can be integrated within this framework and rapid trade-off studies involving one or many disciplines can be performed. The prospect of combining neural network based optimization methods and evolutionary algorithms to obtain a hybrid method with the best properties of both methods will be included in this presentation. Achieving solution diversity and accurate convergence to the exact Pareto front in multiple objective optimization usually requires a significant computational effort with evolutionary algorithms. In this lecture we will also explore the possibility of using neural networks to obtain estimates of the Pareto optimal front using non-dominated solutions generated by DE as training data. Neural network estimators have the potential advantage of reducing the number of function evaluations required to obtain solution accuracy and diversity, thus reducing cost to design.

Rai, Man Mohan↗
DOE OSTI2023

Quasi‐dynamic breathing model of the lung incorporating viscoelasticity of the lung tissue

Abstract We advanced a novel model to calculate viscoelastic lung compliance and airflow resistance in presence of mucus, accounting for the quasi‐linear viscoelastic stress–strain response of the parenchyma (alveoli) tissue. We adapted a continuum‐based numerical modeling approach for the lung, integrating the fluid mechanics of the airflow within individual generations of the bronchi and alveoli. The model accounts for elasticity of the deformable bronchioles, resistance to airflow due to the presence of mucus within the bronchioles, and subsequent mucus flow. Simulated quasi‐dynamic inhalation and expiration cycles were used to characterize the net compliance and resistance of the lung, considering the rheology of the mucus and viscoelastic properties of the parenchyma tissue. The structure and material properties of the lung were identified to have an important contribution to the lung compliance and airflow resistance. The secondary objective of this work was to assess whether a higher frequency and smaller volume of harmonic air flow rate compared to a normal ventilator breathing cycle enhanced mucus outflow. Results predict, lower mucus viscosity and higher excitation frequency of breathing are favorable for the flow of mucus up the bronchi tree, towards the trachea.

60 APPLIED LIFE SCIENCES↗