SEARCH · Engineering Papers

Results for “direct 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 1,261 records · Page 70

All sources NASA NTRS DOE OSTI

NASA NTRS1990

Efficient Monte Carlo simulation of rarefied flow in a small nozzle

The flow of carbon dioxide through a small conical nozzle is computed using two different numerical approaches. The first of these is the direct simulation Monte Carlo method in which the physical phenomena are modeled at the molecular level. The second technique solves the Navier-Stokes equations of fluid flow. The purpose of the paper is to compare the solutions generated by these two methods for the conditions prevailing in a small, low density nozzle. The comparison between the two sets of solutions reveals only small differences. However, the flow properties in the thick viscous layer close to the nozzle at the exit plane show important deviations which have implications for the subsequent determination of the plume backflow.

Boyd, Iain D.↗

NASA NTRS1989

Computational fluid dynamics at the Lewis Research Center: An overview

Lewis is a multidisciplinary Center with strong research and development programs in aeronautical and space propulsion, power, space communications, space experiments and materials. Computational fluid dynamics (CFD) is playing an important and growing role in most of these areas. Described here is how CFD is integrated into these programs and highlights elements of the CFD activities. Examples are presented of codes developed to predict flow fields in advanced propulsion systems and several of the code validation experiments are described. The CFD effort at Lewis ranges from basic research on new and improved algorithms through code development to the application of these codes to specific engineering problems. Because of the substantial improvement in CFD's predictive capability, its use at Lewis is on a steep growth path, spreading rapidly into new areas which had not traditionally taken advantage of the techniques of numerical simulation. Multidisciplinary codes and the future direction of CFD at Lewis are discussed.

Stubbs, Robert M.↗

NASA NTRS1989

Direct simulation of rarefied hypersonic flows

As the capability of the space transportation vehicles (STV's) expand to meet the requirements for future space exploration and utilization, the effects of rarefied hypersonic flows will play a more significant role in defining the aerodynamic and aerothermodynamic performance of STV's. This is particularly true of the low lift/drag aeroassisted STV's where aerobraking occurs at relatively high altitudes and high velocity. Because of the limitations of the continuum description as expressed by the Navier-Stokes equations and the difficulties of solving the Boltzmann equation, the particle of molecular approach has been developed over the last three decades for modeling rarefied gas effects. The direct simulation Monte Carlo (DSMC) method of Bird is the most used method today for simulating rarefied flows. The DSMC method provides a direct physical simulation as opposed to a numerical solution of a set of model equations. This is accomplished by developing phenomenological models of the relevant physical events. The DSMC method accounts for translational, thermal, chemical, and radiative nonequilibrium effects. The general features of the DSMC method, the numerical requirements for obtaining meaningful results, the modeling used to simulate high temperature gas effects, and applications of the method to calculate the flow about an aeroassist flight experiment vehicle (AFE) are reviewed. The AFE simulates a geosynchronous return while entering the Earth's upper atmosphere at approximately 10 km/s. Results obtained using a general 3-D code are presented for the more rarefied portion of the atmospheric encounter (altitudes of 200 to 100 km) emphasizing surface, flowfield, and aerodynamic characteristics of the AFE. Finally, results obtained using axisymmetric and 1-D versions of the code are presented for lower altitude conditions.

Moss, James N.↗

NASA NTRS1991

Shock wave-boundary layer interactions in rarefied gas flows

A numerical study is presented, using the direct simulation Monte Carlo (DSMC) method, of shock wave-boundary layer interactions in low density supersonic flows. Test cases include two-dimensional, axially-symmetric and three-dimensional flows. The effective displacement angle of the boundary layer is calculated for representative flat plate, wedge, and cone flows. The maximum pressure, shear stress, and heat transfer in the shock formation region is determined in each case. The two-dimensional reflection of an oblique shock wave from a flat plate is studied, as is the three-dimensional interaction of such a wave with a sidewall boundary layer.

Bird, G. A.↗

NASA NTRS2011

Impact of Aerosols on Convective Clouds and Precipitation

Aerosols are a critical factor in the atmospheric hydrological cycle and radiation budget. As a major reason for clouds to form and a significant attenuator of solar radiation, aerosols affect climate in several ways. Current research suggests that aerosol effects on clouds could further extend to precipitation, both through the formation of cloud particles and by exerting persistent radiative forcing on the climate system that disturbs dynamics. However, the various mechanisms behind these effects, in particular the ones connected to precipitation, are not yet well understood. The atmospheric and climate communities have long been working to gain a better grasp of these critical effects and hence to reduce the significant uncertainties in climate prediction resulting from such a lack of adequate knowledge. The central theme of this paper is to review past efforts and summarize our current understanding of the effect of aerosols on precipitation processes from theoretical analysis of microphysics, observational evidence, and a range of numerical model simulations. In addition, the discrepancy between results simulated by models, as well as that between simulations and observations will be presented. Specifically, this paper will address the following topics: (1) fundamental theories of aerosol effects on microphysics and precipitation processes, (2) observational evidence of the effect of aerosols on precipitation processes, (3) signatures of the aerosol impact on precipitation from large-scale analyses, (4) results from cloud-resolving model simulations, and (5) results from large-scale numerical model simulations. Finally, several future research directions on aerosol - precipitation interactions are suggested.

Tao, Wei-Kuo↗

NASA NTRS2012

Impact of Aerosols on Convective Clouds and Precipitation

Aerosols are a critical factor in the atmospheric hydrological cycle and radiation budget. As a major agent for clouds to form and a significant attenuator of solar radiation, aerosols affect climate in several ways. Current research suggests that aerosol effects on clouds could further extend to precipitation, both through the formation of cloud particles and by exerting persistent radiative forcing on the climate system that disturbs dynamics. However, the various mechanisms behind these effects, in particular the ones connected to precipitation, are not yet well understood. The atmospheric and climate communities have long been working to gain a better grasp of these critical effects and hence to reduce the significant uncertainties in climate prediction resulting from such a lack of adequate knowledge. Here we review past efforts and summarize our current understanding of the effect of aerosols on convective precipitation processes from theoretical analysis of microphysics, observational evidence, and a range of numerical model simulations. In addition, the discrepancy between results simulated by models, as well as that between simulations and observations, are presented. Specifically, this paper addresses the following topics: (1) fundamental theories of aerosol effects on microphysics and precipitation processes, (2) observational evidence of the effect of aerosols on precipitation processes, (3) signatures of the aerosol impact on precipitation from largescale analyses, (4) results from cloud-resolving model simulations, and (5) results from large-scale numerical model simulations. Finally, several future research directions for gaining a better understanding of aerosol--cloud-precipitation interactions are suggested.

Tao, Wei-Kuo↗

DOE OSTI2022

Data assimilation for burnup distribution of PWR with three-dimensional variational algorithm and artificial neutral network

In this paper, a data-assimilation method has been proposed and applied for the burnup distribution of PWR. The burnup distribution is significant to the safety and economy of the reactor, as it is essential for the fuel-reloading design and optimization. Due to the burnup distribution cannot be measured directly during the reactor operation, the numerical simulation is widely applied to determine the burnup distribution. However, there is a deviation between the numerical simulation and the actual core due to some unavoidable factors, such as component manufacturing deviation, uneven flow distribution and so on. These differences would induce the errors to the simulation values of power distributions and hence to the burnup distribution. To address this problem, a data-assimilation method for the burnup distribution has been proposed with the application of power-distribution measurements. In our research, the three-dimensional variational (3DVAR) algorithm was applied for burnup-distribution calibration and the artificial neutral network (ANN) was applied to establish the relation between power distribution and corresponding burnup distribution. As engineering verification, the proposed data-assimilation method has been applied to the CNP1000 PWR operated in China. The numerical results indicated that the burnup-distribution errors can be reduced notably, as the maximum value of relative errors for power distribution can be reduced from 5.35% to 3.96%. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

NASA NTRS1992

A domain-specific design architecture for composite material design and aircraft part redesign

Advanced composites have been targeted as a 'leapfrog' technology that would provide a unique global competitive position for U.S. industry. Composites are unique in the requirements for an integrated approach to designing, manufacturing, and marketing of products developed utilizing the new materials of construction. Numerous studies extending across the entire economic spectrum of the United States from aerospace to military to durable goods have identified composites as a 'key' technology. In general there have been two approaches to composite construction: build models of a given composite materials, then determine characteristics of the material via numerical simulation and empirical testing; and experience-directed construction of fabrication plans for building composites with given properties. The first route sets a goal to capture basic understanding of a device (the composite) by use of a rigorous mathematical model; the second attempts to capture the expertise about the process of fabricating a composite (to date) at a surface level typically expressed in a rule based system. From an AI perspective, these two research lines are attacking distinctly different problems, and both tracks have current limitations. The mathematical modeling approach has yielded a wealth of data but a large number of simplifying assumptions are needed to make numerical simulation tractable. Likewise, although surface level expertise about how to build a particular composite may yield important results, recent trends in the KBS area are towards augmenting surface level problem solving with deeper level knowledge. Many of the relative advantages of composites, e.g., the strength:weight ratio, is most prominent when the entire component is designed as a unitary piece. The bottleneck in undertaking such unitary design lies in the difficulty of the re-design task. Designing the fabrication protocols for a complex-shaped, thick section composite are currently very difficult. It is in fact this difficulty that our research will address.

Punch, W. F., III↗

NASA NTRS1984

Numerical simulation of turbulent flows

Computational models of turbulence in incompressible Newtonian fluids governed by the Navier-Stokes equations are reviewed. The governing equations are presented, and both direct and large-eddy-simulations are examined. Resolution requirements and numerical techniques of spatial representation, definition of initial and boundary conditions, and time advancement are considered. Results of simulations of homogeneous turbulence in uniform shear, the evolution of a turbulent mixing layer, and turbulent channel flow are presented graphically and discussed.

Rogallo, R. S.↗

NASA NTRS2019

Investigation of Direct Force Control for Planetary Aerocapture at Neptune

In this work, a direct force control numerical predictor-corrector guidance architecture is developed to enable Neptune aerocapture using flight-heritage blunt body aeroshells. A linear aerodynamics model is formulated for a Mars Science Laboratory-derived aeroshell. The application of calculus of variations shows that the optimal angle of attack and side-slip angle control laws are bang-bang. A closed-loop numerical predictor-corrector direct force control guidance algorithm is developed and numerically simulated using the Program to Optimize Simulated Trajectories II. The Monte Carlo simulated trajectories are demonstrated to be robust to the modeled dispersions in aerodynamics, atmospheric density, and entry state. An aerocapture technology trade study demonstrates that blunt body direct force control aerocapture enables similar performance as slender body bank angle control but halves the peak g-loading.

Deshmukh, Rohan G.↗

NASA NTRS2020

Investigation of Direct Force Control For Aerocapture at Neptune

In this work, a direct force control numerical predictor-corrector guidance architecture is developed to enable Neptune aerocapture using blunt body aeroshells. A linear aerodynamics model is formulated for a Mars Science Laboratory-derived aeroshell. The application of optimal control theory shows that the∆V-minimizing angle of attack and side-slip angle control laws are bang-bang. A closed-loop numerical predictor-corrector direct force control guidance algorithm is developed and numerically simulated using the Program to Optimize Simulated Trajectories II. A series of Monte Carlo simulations are conducted to assess the guidance robustness to uncertainties in vehicle aerodynamics,atmospheric density, and entry state. For the reference set of uncertainties, the direct force control vehicle achieves 99.7% successful science orbit insertion within a 330 m/s total∆V budget for periapsis raise, apoapsis, inclination,and ascending node corrections. Improved atmospheric knowledge and delivery state accuracy are shown to improve the success to 100% and reduce the∆V to 230 m/s. Direct force control is demonstrated to be an enabling technology for blunt body aerocapture at Neptune while providing comparable performance to existing slender body vehicles studied in literature

aerocapture↗

DOE OSTI2022

The ICON-A model for direct QBO simulations on GPUs (version icon-cscs:baf28a514)

Abstract. Classical numerical models for the global atmosphere, as used for numerical weather forecasting or climate research, have been developed for conventional central processing unit (CPU) architectures. This hinders the employment of such models on current top-performing supercomputers, which achieve their computing power with hybrid architectures, mostly using graphics processing units (GPUs). Thus also scientific applications of such models are restricted to the lesser computer power of CPUs. Here we present the development of a GPU-enabled version of the ICON atmosphere model (ICON-A), motivated by a research project on the quasi-biennial oscillation (QBO), a global-scale wind oscillation in the equatorial stratosphere that depends on a broad spectrum of atmospheric waves, which originates from tropical deep convection. Resolving the relevant scales, from a few kilometers to the size of the globe, is a formidable computational problem, which can only be realized now on top-performing supercomputers. This motivated porting ICON-A, in the specific configuration needed for the research project, in a first step to the GPU architecture of the Piz Daint computer at the Swiss National Supercomputing Centre and in a second step to the JUWELS Booster computer at the Forschungszentrum Jülich. On Piz Daint, the ported code achieves a single-node GPU vs. CPU speedup factor of 6.4 and allows for global experiments at a horizontal resolution of 5 km on 1024 computing nodes with 1 GPU per node with a turnover of 48 simulated days per day. On JUWELS Booster, the more modern hardware in combination with an upgraded code base allows for simulations at the same resolution on 128 computing nodes with 4 GPUs per node and a turnover of 133 simulated days per day. Additionally, the code still remains functional on CPUs, as is demonstrated by additional experiments on the Levante compute system at the German Climate Computing Center. While the application shows good weak scaling over the tested 16-fold increase in grid size and node count, making also higher resolved global simulations possible, the strong scaling on GPUs is relatively poor, which limits the options to increase turnover with more nodes. Initial experiments demonstrate that the ICON-A model can simulate downward-propagating QBO jets, which are driven by wave–mean flow interaction.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1999

Numerical Simulations of Crystal Growth of an Alloy Under Microgravity Conditions

The directional solidification of a dilute binary alloy (Bi-1.0 at.%Sn)is investigated. Results are obtained at a gravity level of I pg. Computations are performed in two dimensions with a fixed, non-uniform grid. The simulation involves a solution of the species concentration equation (modified to account for solute rejection at the interface) and energy equation (modified to account for phase-change) for both the solid and liquid phases, in addition to the constitutive equations for describing convective flow in the melt. The effects of conductive heat transfer in the ampoule and in a capillary tube in the sample are included. To gauge the effects of including this growth capillary tube in the apparatus, simulations both with and without the capillary tube are presented and compared. Fully transient simulations have been performed; no simplifying steady-state approximations are used, however, the influence of solute on the melting temperature at the interface is not included. Both thermal and solutal convective cells are seen to form. Convective velocities are significantly damped inside the capillary, causing less segregation due to convection. As solidification proceeds beyond the capillary tube, longitudinal segregation arises as a result of the change in cross-sectional area of solidifying material. The magnitudes of the velocities in this cell increase significantly once the solid/liquid front passes beyond the end of the capillary tube; this causes a corresponding increase in the level of radial solute segregation in the solidified material.

Simpson, James E.↗

DOE OSTI2021

Verification of numerical models for seismic fluid-structure interaction analysis of internal components in liquid-filled advanced reactors

Earthquake shaking of a liquid-filled advanced reactor induces fluid-structure interaction (FSI) between the reactor vessel (tank), its internal components, and the contained liquid. Seismic design, qualification, and risk assessment of such reactors must consider fluid-structure responses, their geometries and support conditions, and three-directional seismic inputs, all of which require numerical simulations. Numerical models used for analysis of safety-related nuclear equipment must be verified and validated. Here, this paper verifies numerical models of submerged components using analytical solutions. Seismic FSI analysis of the numerical models is performed using the Arbitrary Lagrangian-Eulerian (ALE) and Incompressible Computational Fluid Dynamics (ICFD) solvers in LS-DYNA. Prior analytical solutions are reworked in this paper, and calculation errors are identified and corrected. The solutions address frequencies of two concentric cylindrical pipes filled with liquid: the inner pipe is a submerged component, and the outer pipe contains the liquid. A unitless frequency coefficient is defined and presented for pipes of different materials, filled with different liquids, and with a range of dimensions for application to advanced reactors. The numerical models are verified here by comparing the lateral frequencies of submerged components with those calculated using the corrected analytical solutions. Recommendations for verification of numerical models of internal components in advanced reactors for seismic FSI analysis are provided. Although the unitless frequency coefficients and verification procedures are developed for application to advanced nuclear reactors, they are broadly applicable to FSI analysis of submerged components in liquid-filled vessels such as storage tanks, boilers, and steam generators.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2023

Develop a Fast Analysis Solver for Welding Sequence Optimization

During the shipbuilding manufacturing process, materials are exposed to significant stresses, as induced both thermally and mechanically, that alter the intended design and significantly affect the production schedule, labor hours (fitting, welding, rework, etc.), and material structural performance. The type and magnitude of deformation of a given structure depends on many factors such as the material, thickness and quality of components, the process heat input, preheat and inter-pass temperatures, type and size of welds, welding sequence and direction, location, sequence, and degree of fixturing. Numerical simulations using finite element analysis (FEA) have long been used to analyze welding-induced structural distortion. For large assemblies, transient thermal elastic-plastic analysis (TEPA) can take days or weeks to run, and optimization of welding sequence is not feasible. Simplified analysis methods were developed to reduce computational time. However, it is challenging to use these techniques to fully optimize welding sequencing because of their applied simplifications in modeling weld details. A fast analysis solver that could be used by the shipbuilding industry is being developed for optimizing welding sequences by taking full advantage of modern GPU-based HPC hardware and incorporating patented acceleration schemes. The accelerated processing factors are up to 2200 times greater for large, multi-pass welded structures.

Yang, Yu-Ping↗

NASA NTRS1991

An efficient DSMC algorithm applied to a delta wing

A new algorithm for 3D direct simulation Monte Carlo (DSMC) is tested and numerical results are compared with wind tunnel data and results obtained earlier with a more traditional DSMC code. The test case is the flowfield around a delta wing at incidence at Knudsen number of 0.016 and Mach number of 20.2. The results are shown to compare favorably with both experimental and earlier numerical results. The new algorithm is described with special emphasis placed on its distinctive features: Cartesian/unstructured combination grid, special body surface definition, discretization in physical space.

Rault, D. F. G.↗

NASA NTRS1996

Classical problems in computational aero-acoustics

In relation to the expected problems in the development of computational aeroacoustics (CAA), the preliminary applications were to classical problems where the known analytical solutions could be used to validate the numerical results. Such comparisons were used to overcome the numerical problems inherent in these calculations. Comparisons were made between the various numerical approaches to the problems such as direct simulations, acoustic analogies and acoustic/viscous splitting techniques. The aim was to demonstrate the applicability of CAA as a tool in the same class as computational fluid dynamics. The scattering problems that occur are considered and simple sources are discussed.

Hardin, Jay C.↗

NASA NTRS2011

Inlet Flow Control and Prediction Technologies for Embedded Propulsion Systems

Fail-safe, hybrid, flow control (HFC) is a promising technology for meeting high-speed cruise efficiency, low-noise signature, and reduced fuel-burn goals for future, Hybrid-Wing-Body (HWB) aircraft with embedded engines. This report details the development of HFC technology that enables improved inlet performance in HWB vehicles with highly integrated inlets and embedded engines without adversely affecting vehicle performance. In addition, new test techniques for evaluating Boundary-Layer-Ingesting (BLI)-inlet flow-control technologies developed and demonstrated through this program are documented, including the ability to generate a BLI-like inlet-entrance flow in a direct-connect, wind-tunnel facility, as well as, the use of D-optimal, statistically designed experiments to optimize test efficiency and enable interpretation of results. Validated improvements in numerical analysis tools and methods accomplished through this program are also documented, including Reynolds-Averaged Navier-Stokes CFD simulations of steady-state flow physics for baseline, BLI-inlet diffuser flow, as well as, that created by flow-control devices. Finally, numerical methods were employed in a ground-breaking attempt to directly simulate dynamic distortion. The advances in inlet technologies and prediction tools will help to meet and exceed "N+2" project goals for future HWB aircraft.

McMillan, Michelle L.↗