SEARCH · Engineering Papers

Results for “Numerical Simulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 667 records · Page 37

All sources NASA NTRS DOE OSTI

NASA NTRS1990

Numerical simulation of three-dimensional transonic flows

The three-dimensional flow over a projectile has been computed using an implicit, approximately factored, partially flux-split algorithm. A simple composite grid scheme has been developed in which a single grid is partitioned into a series of smaller grids for applications which require an external large memory device such as the SSD of the CRAY X-MP/48 or multi-tasking. The accuracy and stability of the composite grid scheme have been tested by numerically simulating the flow over an ellipsoid at an angle of attack and comparing the solution with a single-grid solution. The flow field over a projectile at M = 0.96 and 1.1, and 4-deg angle of attack has been computed using a fine grid and compared with experiment.

Sahu, Jubaraj↗

DOE OSTI2015

Numerical simulations of turbulent molecular clouds regulated by reprocessed radiation feedback from nascent super star clusters

Radiation feedback from young star clusters embedded in giant molecular clouds (GMCs) is believed to be important to the control of star formation. For the most massive and dense clouds, including those in which super star clusters (SSCs) are born, pressure from reprocessed radiation exerted on dust grains may disperse a significant portion of the cloud mass back into the interstellar medium. Using our radiation hydrodynamics code, Hyperion, we conduct a series of numerical simulations to test this idea. Our models follow the evolution of self-gravitating, strongly turbulent clouds in which collapsing regions are replaced by radiating sink particles representing stellar clusters. We evaluate the dependence of the star formation efficiency (SFE) on the size and mass of the cloud and κ, the opacity of the gas to infrared (IR) radiation. We find that the single most important parameter determining the evolutionary outcome is κ, with κ≳15 cm{sup 2} g{sup −1} needed to disrupt clouds. For κ=20−40 cm{sup 2} g{sup −1}, the resulting SFE = 50%−70% is similar to empirical estimates for some SSC-forming clouds. The opacities required for GMC disruption likely apply only in dust-enriched environments. We find that the subgrid model approach of boosting the direct radiation force L/c by a “trapping factor” equal to a cloud’s mean IR optical depth can overestimate the true radiation force by factors of ∼4−5. We conclude that feedback from reprocessed IR radiation alone is unlikely to significantly reduce star formation within GMCs unless their dust abundances or cluster light-to-mass ratios are enhanced.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Influence of Blood Rheology and Turbulence Models in the Numerical Simulation of Aneurysms

An aneurysm is a vascular malformation that can be classified according to its location (cerebral, aortic) or shape (saccular, fusiform, and mycotic). Recently, the study of blood flow interaction with aneurysms has gained attention from physicians and engineers. Shear stresses, oscillatory shear index (OSI), gradient oscillatory number (GON), and residence time have been used as variables to describe the hemodynamics as well as the origin and evolution of aneurysms. However, the causes and hemodynamic conditions that promote their growth are still under debate. The present work presents numerical simulations of three types of aneurysms: two aortic and one cerebral. Simulation results showed that the blood rheology is not relevant for aortic aneurysms. However, for the cerebral aneurysm case, blood rheology could play a relevant role in the hemodynamics. The evaluated turbulence models showed equivalent results in both cases. Lastly, a simulation considering the fluid–structure interaction (FSI) showed that this phenomenon is the dominant factor for aneurysm simulation.

59 BASIC BIOLOGICAL SCIENCES↗

NASA NTRS1983

Numerical simulation of separated flows

A new numerical method, based on the Vortex Method, for the simulation of two-dimensional separated flows, was developed and tested on a wide range of gases. The fluid is incompressible and the Reynolds number is high. A rigorous analytical basis for the representation of the Navier-Stokes equation in terms of the vorticity is used. An equation for the control of circulation around each body is included. An inviscid outer flow (computed by the Vortex Method) was coupled with a viscous boundary layer flow (computed by an Eulerian method). This version of the Vortex Method treats bodies of arbitrary shape, and accurately computes the pressure and shear stress at the solid boundary. These two quantities reflect the structure of the boundary layer. Several versions of the method are presented and applied to various problems, most of which have massive separation. Comparison of its results with other results, generally experimental, demonstrates the reliability and the general accuracy of the new method, with little dependence on empirical parameters. Many of the complex features of the flow past a circular cylinder, over a wide range of Reynolds numbers, are correctly reproduced.

Spalart, P. R.↗

NASA NTRS1998

Studying Turbulence Using Numerical Simulation Databases: Proceedings of the Summer Program - No. 7

The Seventh Summer Program of the Center for Turbulence Research took place in the four-week period, July 5 to July 31, 1998. This was the largest CTR Summer Program to date, involving thirty-six participants from the U. S. and nine other countries. Thirty-one Stanford and NASA-Ames staff members facilitated and contributed to most of the Summer projects. A new feature, and perhaps a preview of the future programs, was that many of the projects were executed on non-NASA computers. These included supercomputers located in Europe as well as those operated by the Departments of Defense and Energy in the United States. In addition, several simulation programs developed by the visiting participants at their home institutions were used. Another new feature was the prevalence of lap-top personal computers which were used by several participants to carry out some of the work that in the past were performed on desk-top workstations. We expect these trends to continue as computing power is enhanced and as more researchers (many of whom CTR alumni) use numerical simulations to study turbulent flows. CTR's main role continues to be in providing a forum for the study of turbulence for engineering analysis and in facilitating intellectual exchange among the leading researchers in the field. Once again the combustion group was the largest. Turbulent combustion has enjoyed remarkable progress in using simulations to address increasingly complex and practically more relevant questions. The combustion group's studies included such challenging topics as fuel evaporation, soot chemistry, and thermonuclear reactions. The latter study was one of three projects related to the Department of Energy's ASCI Program (www.llnl.gov/asci); the other two (rocket propulsion and fire safety) were carried out in the turbulence modeling group. The flow control and acoustics group demonstrated a successful application of the so-called evolution algorithms which actually led to a previously unknown forcing strategy for jets yielding increased spreading rate. A very efficient algorithm for flow in complex geometries with moving boundaries based on the immersed boundary forcing technique was tested with very encouraging results. Also a new strategy for the destruction of aircraft trailing vortices was introduced and tested. The Reynolds Averaged Modeling (RANS) group demonstrated that the elliptic relaxation concept for RANS calculations is also applicable to transonic flows with shocks; however, prediction of laminar/turbulent transition remains an important pacing item. A large fraction of the LES effort was devoted to the development and testing of a new algorithmic procedure (as opposed to phenomenological model) for subgrid scale modeling based on regularized de-filtering of the flow variables. This appears to be a very promising approach, and a significant effort is currently underway to assess its robustness in high Reynolds number flows and in conjunction with numerical methods for complex flows. As part of the Summer Program two review tutorials were given on Turbulent structures in hydrocarbon pool fires (Sheldon Tieszen), and Turbulent combustion modeling: from RANS to LES via DNS (Luc Vervisch); and two seminars entitled Assessment of turbulence models for engineering applications (Paul Durbin) and Subgrid-scale modeling for non-premixed, turbulent reacting flows (James Riley) were presented. A number of colleagues from universities, government agencies, and industry attended the final presentations of the participants on July 31 and participated in the discussions. There are twenty-six papers in this volume grouped in five areas. Each group is preceded with an overview by its coordinator.

Source record↗

DOE OSTI2017

Numerical Simulation and Experimental Study of the Gas–Solid Flow Behavior Inside a Full-Loop Circulating Fluidized Bed: Evaluation of Different Drag Models

Both experimental and computational studies of the fluidization of high-density polyethylene (HDPE) particles in a small-scale full-loop circulating fluidized bed (CFB) are conducted. Experimental measurements of pressure drop are taken at various locations along the bed. The solids circulation rate is measured with an advanced particle image velocimetry (PIV) technique. Bed height of the quasi-static region in the standpipe is also measured. Comparative numerical simulations are performed with a computational fluid dynamics solver utilizing a discrete element method (CFD-DEM). Here, this paper examines the effect of different drag laws used in the CFD simulations through a detailed and direct comparison with experimental data from a small-scale, full-loop circulating fluidized bed. The Hill–Koch–Ladd drag correlation was shown to have good agreement with respect to system component pressure drop and inventory height in the standpipe.

42 ENGINEERING↗

DOE OSTI2015

Numerical simulation of the environmental impact of hydraulic fracturing of tight/shale gas reservoirs on near‐surface groundwater: Background, base cases, shallow reservoirs, short‐term gas, and water transport

Abstract Hydrocarbon production from unconventional resources and the use of reservoir stimulation techniques, such as hydraulic fracturing, has grown explosively over the last decade. However, concerns have arisen that reservoir stimulation creates significant environmental threats through the creation of permeable pathways connecting the stimulated reservoir with shallower freshwater aquifers, thus resulting in the contamination of potable groundwater by escaping hydrocarbons or other reservoir fluids. This study investigates, by numerical simulation, gas and water transport between a shallow tight‐gas reservoir and a shallower overlying freshwater aquifer following hydraulic fracturing operations, if such a connecting pathway has been created. We focus on two general failure scenarios: (1) communication between the reservoir and aquifer via a connecting fracture or fault and (2) communication via a deteriorated, preexisting nearby well. We conclude that the key factors driving short‐term transport of gas include high permeability for the connecting pathway and the overall volume of the connecting feature. Production from the reservoir is likely to mitigate release through reduction of available free gas and lowering of reservoir pressure, and not producing may increase the potential for release. We also find that hydrostatic tight‐gas reservoirs are unlikely to act as a continuing source of migrating gas, as gas contained within the newly formed hydraulic fracture is the primary source for potential contamination. Such incidents of gas escape are likely to be limited in duration and scope for hydrostatic reservoirs. Reliable field and laboratory data must be acquired to constrain the factors and determine the likelihood of these outcomes.

03 NATURAL GAS↗

DOE OSTI2017

Mineralogical and transport controls on the evolution of porous media texture using direct numerical simulation

The evolution of porous media due to mineral dissolution and precipitation can change the bulk properties of subsurface materials. The pore-scale structure of the media, including its physical and mineralogical heterogeneity, exerts controls on porous media evolution via transport limitations to reactive surfaces and mineral accessibility. Here in this paper we explore how these controls affect the evolution of the texture in porous media at the pore scale. For this purpose, a pore-scale flow and reactive transport model is developed that explicitly tracks mineral surfaces as they evolve using a direct numerical simulation approach. Simulations of dissolution in single-mineral domains provide insights into the transport controls at the pore scale, while the simulation of a fracture surface composed of bands of faster-dissolving calcite and slower-dissolving dolomite provides insights into the mineralogical controls on evolution. Transport-limited conditions at the grain-pack scale may result in unstable evolution, a situation in which dissolution is focused in a fast-flowing, fast-dissolving path. Due to increasing velocities, the evolution in these regions is like that observed under conditions closer to strict surface control at the pore scale. That is, grains evolve to have oblong shapes with their long dimensions aligning with the local flow directions. Another example of an evolving reactive transport regime that affects local rates is seen in the evolution of the fracture surface. As calcite dissolves, the diffusive length between the fracture flow path and the receding calcite surfaces increases. Thus, the calcite dissolution reaction becomes increasingly limited by diffusion.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Numerical simulation of gas production from natural gas hydrate deposits with multi-branch wells: Influence of reservoir properties

Vast amounts of natural gas hydrate are buried in subseafloor sediments without impermeable boundaries, which is recognized as an essential energy source for the future. Previously, the multi-branch well was proposed to enhance the recovery efficiency of natural gas hydrate, and the gas production rate has been dramatically improved comparing with the vertical well. However, the multi-branch well shows a terrible performance in gas production duration. Here, as a continuation of the previous study, numerical simulations were conducted to investigate the influence of hydrate reservoir properties on the gas production potential. Results indicate that it is hard to extract hydrate commercially for hydrate accumulations without impermeable boundaries. A high initial hydrate saturation leads to a long gas production duration but a low gas production rate. An increase in the intrinsic permeability of isotropic reservoirs would shorten the gas production duration and result in a low gas recovery ratio. Permeability anisotropy shows a noticeable effect on enhancing the gas recovery ratio and the gas production duration due to the improved pressure propagation pattern. Therefore, in the upcoming field tests, reservoir reconstructions that enhance permeability anisotropy are strongly suggested to obtain better outcomes.

03 NATURAL GAS↗

DOE OSTI2015

Numerical simulations of Mach stem formation via intersecting bow shocks

Hubble Space Telescope observations show bright knots of Ha emission within outflowing young stellar jets. Velocity variations in the flow create secondary bow shocks that may intersect and lead to enhanced emission. When the bow shocks intersect at or above a certain critical angle, a planar shock called a Mach stem is formed. These shocks could produce brighter Ha emission since the incoming flow to the Mach stem is parallel to the shock normal. In this paper we report first results of a study using 2-D numerical simulations designed to explore Mach stem formation at the intersection of bow shocks formed by hypersonic “bullets” or “clumps”. Our 2-D simulations show how the bow shock shapes and intersection angles change as the adiabatic index g changes. We show that the formation or lack of a Mach stem in our simulations is consistent with the steady-state Mach stem formation theory. In conclusion, our ultimate goal, which is part of an ongoing research effort, is to characterize the physical and observational conse- quences of bow shock intersections including the formation of Mach stems.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2017

Numerical simulation of star formation by the bow shock of the centaurus a jet

Recent Hubble Space Telescope (HST) observations of the extragalactic radio source Centaurus A (Cen A) display a young stellar population around the southwest tip of the inner filament 8.5 kpc from the galactic center of Cen A, with ages in the range 1–3 Myr. Crockett et al. have argued that the transverse bow shock of the Cen A jet triggered this star formation as it impacted dense molecular cores of clouds in the filament. To test this hypothesis, we perform three-dimensional numerical simulations of star formation induced by the jet bow shock in the inner filament of Cen A, using a positivity-preserving, weighted, essentially non-oscillatory method to solve the equations of gas dynamics with radiative cooling. We find that star clusters form inside a bow-shocked molecular cloud when the maximum initial density of the cloud is ⩾40 H{sub 2} molecules cm{sup −3}. In a typical molecular cloud of mass 10{sup 6} M{sub ⊙} and diameter 200 pc, approximately 20 star clusters of mass 10{sup 4} M{sub ⊙} are formed, matching the HST images.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1994

Numerical simulation of the interaction of transport, diffusion and chemical reactions in an urban plume

A model system is presented that takes into account the main physical and chemical processes on the regional scale here in an area of 100x100 sq km. The horizontal gridsize used is 2x2 sq km. For a case study, it is demonstrated how the model system can be used to separate the contributions of the processes advection, turbulent diffusion, and chemical reactions to the diurnal cycle of ozone. In this way, typical features which are visible in observations and are reproduced by the numerical simulations can be interpreted.

Vogel, Bernhard↗

DOE OSTI2015

Numerical Simulation of Big Area Additive Manufacturing (3D Printing) of a Full Size Car

Distortion and buildup of residual stresses are common design problems in structures manufactured using additive manufacturing techniques, popularly known as 3D printing. These problems lead to poor product quality, which often is improved by trial and error, where the sensitivity of different manufacturing configurations to product quality are experimentally determined. With the introduction of Big Area Additive Manufacturing (BAAM), where significant amount of material is deposited per unit time, trial and error philosophy reduces the economical efficiency of the process. In this work, a numerical simulation methodology based on Finite Element Method (FEM), multi-scale damage mechanics and fracture mechanics is introduced to simulate the BAAM process for determining product quality in terms of distortions, material damage and interface fracture due to manufacturing. In addition, this tool can be used to obtain the sensitivities of the manufacturing configurations with respect to product quality in an economically efficient manner compared to the trial and error philosophy. Computational tools providing a fully coupled thermo-mechanical solution of polymer additive manufacturing with chopped fiber reinforced plastics are presented.

42 ENGINEERING↗

DOE OSTI2019

Numerical Simulation of Metal Melt Flow in a One-Strand Tundish Regarding Active Filtration and Reactive Cleaning

In the paper, two different cleaning strategies for nonmetallic inclusions in steel melts, active filtration and reactive cleaning, are examined in a prototype tundish configuration. In active filtration, nonmetallic inclusions are deposited at the filter surfaces. In reactive cleaning, nonmetallic inclusions stick to the filter surfaces, too. In addition, they are floated by the action of carbon monoxide bubbles, which are generated by reaction between carbon and oxygen in the steel melt. In order to compare the performance of both strategies, numerical simulations of the two-phase flows of steel melt and dispersed nonmetallic inclusions are performed. Turbulence is resolved with implicit large eddy simulation. If necessary, species transports of dissolved carbon in the melt and reaction with oxygen are employed. Cleaning efficiencies are deduced from the simulations which demonstrate that reactive cleaning is much more efficient than active filtration.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Numerical Simulation Studies of Ultrasonic De-Icing for Heating, Ventilation, Air Conditioning, and Refrigeration Structures

Ice accumulation on heating, ventilation, air conditioning, and refrigeration (HVACR) structures presents significant operational challenges. These challenges include reduced efficiency, increased energy consumption, and potential damage to equipment. Traditional de-icing methods, such as chemical treatments, mechanical scraping, or heating-based techniques, are often labor-intensive, costly, and environmentally harmful. Here, this study uniquely investigates ultrasonic de-icing as an energy-efficient alternative for HVACR applications, focusing on the specific structural geometries found in these systems. A comprehensive numerical simulation framework was developed using finite element analysis to explore ultrasonic wave propagation across four distinct HVACR structures. Key parameters such as ultrasonic frequency, power levels, and the number and placement of actuators were examined for their impact on ice detachment efficiency. Results from simulations on a plate structure reveal that ultrasonic excitation can propagate effectively across large areas (at least 150 × 150 mm), enhancing the de-icing coverage. Lower frequency (e.g., 30 to 45 kHz) excitation results in greater displacement, improving de-icing performance, while increased actuator numbers with the same total power input also enhance effectiveness. Two actuators seem sufficient for the de-icing of a 300 × 300 mm plate. For tube-and-fin structures, specific high-power ultrasonic frequencies selectively excite the fin plates, demonstrating efficient ice removal when actuated on the tube. However, optimal performance requires careful design of actuator placement and vibration modes to accommodate the irregular shapes of these structures.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2020

3D Numerical Simulation of Kink-driven Rayleigh–Taylor Instability Leading to Fast Magnetic Reconnection

Fast magnetic reconnection involving non-MHD microscale physics is thought to underlie both solar eruptions and laboratory plasma current disruptions. While there is extensive research on both the MHD macroscale physics and the non-MHD microscale physics, the process by which large-scale MHD couples to the microscale physics is not well understood. An MHD instability cascade from a kink to a secondary Rayleigh–Taylor instability in the Caltech astrophysical jet laboratory experiment provides new insights into this coupling and motivates a 3D numerical simulation of this transition from large to small scale. A critical finding from the simulation is that the axial magnetic field inside the current-carrying dense plasma must exceed the field outside. In addition, the simulation verifies a theoretical prediction and experimental observation that, depending on the strength of the effective gravity produced by the primary kink instability, the secondary instability can be Rayleigh–Taylor or mini-kink. Finally, it is shown that the kink-driven Rayleigh–Taylor instability generates a localized electric field sufficiently strong to accelerate electrons to very high energy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Numerical simulation of frost formation and heat transfer on fin-and-tube heat exchangers in turbulent cross-flow

Frost formation in fin-and-tube heat exchangers in turbulent cross-flow presents significant challenges in industrial refrigeration applications, affecting heat transfer efficiency and operational reliability. The purpose of this work is to investigate frost deposition and growth on a staggered bank of a fin-and-tube freezer coil under turbulent forced convection conditions. The focus here is on investigating conditions that closely replicate real-world scenarios in large walk-in industrial freezers. Using a direct numerical simulation approach, we examine the flow dynamics and thermal behaviour in the presence of frost, considering turbulent regimes characterized by a Reynolds number in the range 1050 ≤ R e D , avg ≤ 4800 , with the characteristic length being the outer diameter of the tube and the velocity being the bulk fluid velocity between the plates (fins). Computational fluid dynamics simulations are employed to resolve the interactions between turbulent airflow and the frost layer. Our approach incorporates a modified immersed boundary method and a slow-time acceleration technique to address the complex dynamic interface between the continuously evolving frost layer and the flowing air stream. Our findings indicate that frost forms more on the sides of the finned surfaces (plates) and less on the tubes themselves. This article is part of the theme issue ‘Heat and mass transfer in frost and ice’.

Science & Technology - Other Topics↗

NASA NTRS1990

Numerical simulation of solar wind density fluctuations and their effects on VLF radio interferometry

A mission consisting of an array of small satellites has been proposed to investigate radio sources at frequencies from approximately 2 to 20 MHz, a range unobservable from the ground because of ionospheric absorption and aberrations. Such a mission seems feasible and relatively inexpensive, but further investigation is still necessary. The solar wind is a major concern because it contains turbulence-induced random density variations that cause a fluctuating phase shift in any electromagnetic signal passing through, sometimes making radio source imaging by interferometry impossible. A numerical simulation has been developed which attempts to accurately characterize solar wind density fluctuations. Results have allowed more concrete conclusions to be drawn about observing constraints and have shown that interferometric imaging at radio frequencies below 10 MHz with baselines up to 100 km is possible at solar elongation angles of about 90 deg or more.

Williamson, Robert S., III↗