Search NASA⌕ Search

SEARCH · Search NASA

Results for “Simulation techniques”

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 811 records · Page 45

On the Direct Construction of the Steady Traveling Solution to High Explosive Sandwich, Cylinder and Aquarium Tests via a Streamline Finite Volume Approximation

The cylinder, sandwich and aquarium tests are routine experiments used to investigate the performance of high explosives. Here, we present a solution technique which simulates these experiments as two dimensional problems in steady traveling coordinates. Using a coordinate system aligned with the streamlines of the flow, an efficient numerical scheme is developed which uses an approximate Riemann solver to model the change in state across streamlines. The technique can also accommodate arbitrary equations of state for the materials. Two verification exercises show that the algorithm converges at second order. Validation studies using the aquarium experiment agreed well with simulations, though the omission of a material strength model leads to some discrepancies between simulation and experiments for the cylinder test. This computational technique is able to predict the behavior of these common high explosive experiments on the order of minutes on a single processor with very fine spatial resolution, on the order of 10 μm.

36 MATERIALS SCIENCE↗

High-fidelity retrieval from instantaneous line-of-sight returns of nacelle-mounted lidar including supervised machine learning

Abstract. Wind turbine applications that leverage nacelle-mounted Doppler lidar are hampered by several sources of uncertainty in the lidar measurement, affecting both bias and random errors. Two problems encountered especially for nacelle-mounted lidar are solid interference due to intersection of the line of sight with solid objects behind, within, or in front of the measurement volume and spectral noise due primarily to limited photon capture. These two uncertainties, especially that due to solid interference, can be reduced with high-fidelity retrieval techniques (i.e., including both quality assurance/quality control and subsequent parameter estimation). Our work compares three such techniques, including conventional thresholding, advanced filtering, and a novel application of supervised machine learning with ensemble neural networks, based on their ability to reduce uncertainty introduced by the two observed nonideal spectral features while keeping data availability high. The approach leverages data from a field experiment involving a continuous-wave (CW) SpinnerLidar from the Technical University of Denmark (DTU) that provided scans of a wide range of flows both unwaked and waked by a field turbine. Independent measurements from an adjacent meteorological tower within the sampling volume permit experimental validation of the instantaneous velocity uncertainty remaining after retrieval that stems from solid interference and strong spectral noise, which is a validation that has not been performed previously. All three methods perform similarly for non-interfered returns, but the advanced filtering and machine learning techniques perform better when solid interference is present, which allows them to produce overall standard deviations of error between 0.2 and 0.3 m s−1, or a 1 %–22 % improvement versus the conventional thresholding technique, over the rotor height for the unwaked cases. Between the two improved techniques, the advanced filtering produces 3.5 % higher overall data availability, while the machine learning offers a faster runtime (i.e., ∼ 1 s to evaluate) that is therefore more commensurate with the requirements of real-time turbine control. The retrieval techniques are described in terms of application to CW lidar, though they are also relevant to pulsed lidar. Previous work by the authors (Brown and Herges, 2020) explored a novel attempt to quantify uncertainty in the output of a high-fidelity lidar retrieval technique using simulated lidar returns; this article provides true uncertainty quantification versus independent measurement and does so for three techniques rather than one.

47 OTHER INSTRUMENTATION↗

G and C boost and abort study summary, exhibit B

A six degree of freedom simulation of rigid vehicles was developed to study space shuttle vehicle boost-abort guidance and control techniques. The simulation was described in detail as an all digital program and as a hybrid program. Only the digital simulation was implemented. The equations verified in the digital simulation were adapted for use in the hybrid simulation. Study results were obtained from four abort cases using the digital program.

Backman, H. D.↗

Multiple Scale Simulation for Transitional and Turbulent Flow

A new concept, Multiple Scale Simulation (MSS), is presented in this paper. The basic idea is that the flow is decomposed into several component groups according to spatial and temporal length scales. Each group has its own subdomain, governing system, mesh size, and discretization method. The simulation is then performed groupwise. This approach has been successfully applied in combination with the intergrid dissipation technique for simulation of transitional and turbulent flow in 3-D boundary layers, and it is feasible for 3-D airfoils and other more complex configurations. MSS should prove to ameliorate the scale problems associated with conventional direct numerical simulation.

Liu, Chaoqun↗

Space Software for Automotive Design

John Thousand of Wolverine Western Corp. put his aerospace group to work on an unfamiliar job, designing a brake drum using computer design techniques. Computer design involves creation of a mathematical model of a product and analyzing its effectiveness in simulated operation. Technique enables study of performance and structural behavior of a number of different designs before settling on a final configuration. Wolverine employees attacked a traditional brake drum problem, the sudden buildup of heat during fast and repeated braking. Part of brake drum not confined tends to change its shape under combination of heat, physical pressure and rotational forces, a condition known as bellmouthing. Since bellmouthing is a major factor in braking effectiveness, a solution of problem would be a major advance in automotive engineering. A former NASA employee, now a Wolverine employee, knew of a series of NASA computer programs ideally suited to confronting bellmouthing. Originally developed as aids to rocket engine nozzle design, it's capable of analyzing problems generated in a rocket engine or automotive brake drum by heat, expansion, pressure and rotational forces. Use of these computer programs led to new brake drum concept featuring a more durable axle, and heat transfer ribs, or fins, on hub of drum.

Source record↗

Three-dimensional Navier-Stokes calculations of multiple interacting vortex rings

Results from a finite-difference Navier-Stokes code for three-dimensional, unsteady, vortical flows in unbounded domains are presented and analyzed in this paper. The vortical flows presented are representative of vortex rings and other closed vortical tubes or structures in fluid mechanics. Such structures are important elements in fluid flows such as jets, atmospheric turbulence, and the far-field wakes of aircraft, and studies of their interaction may aid in an understanding of complex fluid flows. The paper demonstrates that computational methods can be used as a viable alternative or supplement to experimental techniques for studying the physics of vortex flows. The separate visualization of vortex stretching, convection, and diffusion is presented in this paper for a single elliptical vortex ring.The calculations employ a truncated series expansion technique to simulate the unbounded nature of the fluid flow with a finite computational domain, which is a more accurate technique than the conventional freestream boundary specification. The numerical divergence of the three-dimensional vorticity field is considered as a useful estimate of truncation error, and the use of a kinetic energy decay law as a calculation check is demonstrated. Results from the Navier-Stokes code are presented for the unsteady motion of two and four vortex rings along parallel axes, and the results agree qualitatively with experimental flow visualization.

Chamberlain, J. P.↗

Direct numerical simulations of a reacting turbulent mixing layer by a pseudospectral-spectral element method

An account is given of the implementation of the spectral-element technique for simulating a chemically reacting, spatially developing turbulent mixing layer. Attention is given to experimental and numerical studies that have investigated the development, evolution, and mixing characteristics of shear flows. A mathematical formulation is presented of the physical configuration of the spatially developing reacting mixing layer, in conjunction with a detailed representation of the spectral-element method's application to the numerical simulation of mixing layers. Results from 2D and 3D calculations of chemically reacting mixing layers are given.

Mcmurtry, Patrick A.↗

Visualization techniques for the gyrokinetic tokamak simulation code

Gyrokinetic simulations of plasma microturbulence in tokamaks are challenging to visualize because the compute grid follows the magnetic field lines that spiral around the torus. We have overcome this challenge by developing three new approaches that improve visualization of gyrokinetics. Our techniques work directly with the topology of magnetic flux surfaces where the simulation stores variables in concentric rings on poloidal planes (vertical cross sections of the torus). Our visualization preview step triangulates each consecutive pair of rings to display the data on a poloidal plane. The second visualization technique follows spiral field lines around the torus and constructs polygons to visualize a flux surface. Third, the poloidal triangles are connected between planes to form prisms that compose a 3-D model of the entire torus. The visualization workflow produces detailed geometry that matches the high resolution, irregular compute grid for every time step. The surface and solid models are displayed in scientific visualization programs to effectively explore and communicate the results, including fluctuation of electron density, ion temperature, and electrostatic potential. Highly detailed renderings verify plasma behavior along magnetic field lines over time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Computer simulation of carbonization and graphitization of coal

Abstract This study describes computer simulations of carbonization and graphite formation, including the effects of hydrogen, nitrogen, oxygen, and sulfur. We introduce a novel technique to simulate carbonization, ‘Simulation of Thermal Emission of Atoms and Molecules (STEAM)’, designed to elucidate volatile outgassing and density variations in the intermediate material during carbonization. The investigation analyzes the functional groups that endure through high-temperature carbonization and examines the graphitization processes in carbon-rich materials containing non-carbon impurity elements. The physical, vibrational, and electronic attributes of impure amorphous graphite are analyzed, and the impact of nitrogen on electronic conduction is investigated.

Ugwumadu, C. (ORCID:0000000199207594)↗

Unbiased simulation of near-Clifford quantum circuits

Modeling and simulation are essential for predicting and verifying the behavior of fabricated quantum circuits, but existing simulation methods are either impractically costly or require an unrealistic simplification of error processes. In this paper, we present a method of simulating noisy Clifford circuits that is both accurate and practical in experimentally relevant regimes. In particular, the cost is weakly exponential in the size and the degree of non-Cliffordness of the circuit. Our approach is based on the construction of exact representations of quantum channels as quasiprobability distributions over stabilizer operations, which are then sampled, simulated, and weighted to yield unbiased statistical estimates of circuit outputs and other observables. As a demonstration of these techniques, we simulate a Steane [[7,1,3]]-encoded logical operation with non-Clifford errors and compute its fault tolerance error threshold. Finally, we expect that the method presented here will enable studies of much larger and more realistic quantum circuits than was previously possible.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Thematic mapper design parameter investigation

This study simulated the multispectral data sets to be expected from three different Thematic Mapper configurations, and the ground processing of these data sets by three different resampling techniques. The simulated data sets were then evaluated by processing them for multispectral classification, and the Thematic Mapper configuration, and resampling technique which provided the best classification accuracy were identified.

Colby, C. P., Jr.↗

A New Method for Accurate Treatment of Flow Equations in Cylindrical Coordinates Using Series Expansions

The motivation of this work is the ongoing effort at the Center for Turbulence Research (CTR) to use large eddy simulation (LES) techniques to calculate the noise radiated by jet engines. The focus on engine exhaust noise reduction is motivated by the fact that a significant reduction has been achieved over the last decade on the other main sources of acoustic emissions of jet engines, such as the fan and turbomachinery noise, which gives increased priority to jet noise. To be able to propose methods to reduce the jet noise based on results of numerical simulations, one first has to be able to accurately predict the spatio-temporal distribution of the noise sources in the jet. Though a great deal of understanding of the fundamental turbulence mechanisms in high-speed jets was obtained from direct numerical simulations (DNS) at low Reynolds numbers, LES seems to be the only realistic available tool to obtain the necessary near-field information that is required to estimate the acoustic radiation of the turbulent compressible engine exhaust jets. The quality of jet-noise predictions is determined by the accuracy of the numerical method that has to capture the wide range of pressure fluctuations associated with the turbulence in the jet and with the resulting radiated noise, and by the boundary condition treatment and the quality of the mesh. Higher Reynolds numbers and coarser grids put in turn a higher burden on the robustness and accuracy of the numerical method used in this kind of jet LES simulations. As these calculations are often done in cylindrical coordinates, one of the most important requirements for the numerical method is to provide a flow solution that is not contaminated by numerical artifacts. The coordinate singularity is known to be a source of such artifacts. In the present work we use 6th order Pade schemes in the non-periodic directions to discretize the full compressible flow equations. It turns out that the quality of jet-noise predictions using these schemes is especially sensitive to the type of equation treatment at the singularity axis. The objective of this work is to develop a generally applicable numerical method for treating the singularities present at the polar axis, which is particularly suitable for highly accurate finite-differences schemes (e.g., Pade schemes) on non-staggered grids. The main idea is to reinterpret the regularity conditions developed in the context of pseudo-spectral methods. A set of exact equations at the singularity axis is derived using the appropriate series expansions for the variables in the original set of equations. The present treatment of the equations preserves the same level of accuracy as for the interior scheme. We also want to point out the wider utility of the method, proposed here in the context of compressible flow equations, as its extension for incompressible flows or for any other set of equations that are solved on a non-staggered mesh in cylindrical coordinates with finite-differences schemes of various level of accuracy is straightforward. The robustness and accuracy of the proposed technique is assessed by comparing results from simulations of laminar forced-jets and turbulent compressible jets using LES with similar calculations in which the equations are solved in Cartesian coordinates at the polar axis, or in which the singularity is removed by employing a staggered mesh in the radial direction without a mesh point at r = 0.

Constantinescu, G.S.↗

Modeling of Stick-Slip Behavior in Sheared Granular Fault Gouge Using the Combined Finite-Discrete Element Method

Sheared granular layers undergoing stick-slip behavior are broadly employed to study the physics and dynamics of earthquakes. In this paper, a two-dimensional implementation of the combined finite-discrete element method (FDEM), which merges the finite element method (FEM) and the discrete element method (DEM), is used to explicitly simulate a sheared granular fault system including both gouge and plate, and to investigate the influence of different normal loads on seismic moment, macroscopic friction coefficient, kinetic energy, gouge layer thickness, and recurrence time between slips. In the FDEM model, the deformation of plates and particles is simulated using the FEM formulation while particle-particle and particle-plate interactions are modeled using DEM-derived techniques. The simulated seismic moment distributions are generally consistent with those obtained from the laboratory experiments. In addition, the simulation results demonstrate that with increasing normal load, (i) the kinetic energy of the granular fault system increases; (ii) the gouge layer thickness shows a decreasing trend; and (iii) the macroscopic friction coefficient does not experience much change. Analyses of the slip events reveal that, as the normal load increases, more slip events with large kinetic energy release and longer recurrence time occur, and the magnitude of gouge layer thickness decrease also tends to be larger; while the macroscopic friction coefficient drop decreases. Finally, the simulations not only reveal the influence of normal loads on the dynamics of sheared granular fault gouge, but also demonstrate the capabilities of FDEM for studying stick-slip dynamic behavior of granular fault systems.

58 GEOSCIENCES↗

Comparative evaluation of workload estimation techniques in piloting tasks

Techniques to measure operator workload in a wide range of situations and tasks were examined. The sensitivity and intrusion of a wide variety of workload assessment techniques in simulated piloting tasks were investigated. Four different piloting tasks, psychomotor, perceptual, mediational, and communication aspects of piloting behavior were selected. Techniques to determine relative sensitivity and intrusion were applied. Sensitivity is the relative ability of a workload estimation technique to discriminate statistically significant differences in operator loading. High sensitivity requires discriminable changes in score means as a function of load level and low variation of the scores about the means. Intrusion is an undesirable change in the task for which workload is measured, resulting from the introduction of the workload estimation technique or apparatus.

Wierwille, W. W.↗

Direct numerical simulation of laminar breakdown in high-speed, axisymmetric boundary layers

The compressible Navier-Stokes equations are solved using spectral collocation and high-order compact-difference techniques to simulate the laminar breakdown in high-speed, axisymmetric boundary-layer flow. Mach 4.5 flow along a hollow cylinder and Mach 6.8 flow along a sharp cone are considered. Data obtained replicate two previously unexplained phenomena, namely, the appearance of so-called 'rope-like waves' and 'the precursor transition effect', in which transitional flow originates near the critical layer well upstream of the transition location at the wall. The numerical data also reveal that neither of these effects can be explained, even qualitatively, by linear stability theory alone. It is shown that rope-like appearance arises from secondary instability. Certain features of the precursor transition effect also emerge from secondary instability but its nature is revealed to be fundamentally nonlinear.

Pruett, C. D.↗

Accurate low-cost methods for performance evaluation of cache memory systems

Methods of simulation based on statistical techniques are proposed to decrease the need for large trace measurements and for predicting true program behavior. Sampling techniques are applied while the address trace is collected from a workload. This drastically reduces the space and time needed to collect the trace. Simulation techniques are developed to use the sampled data not only to predict the mean miss rate of the cache, but also to provide an empirical estimate of its actual distribution. Finally, a concept of primed cache is introduced to simulate large caches by the sampling-based method.

Laha, Subhasis↗

Coupled fuel performance calculations in VERA and demonstration on Watts Bar unit 1, cycle 1

As the core simulator capabilities in the Virtual Environment for Reactor Applications (VERA) have become more mature and stable, increased attention has been focused on coupling Bison to provide fuel performance simulations. This technique has been a very important driver for the pellet-clad interaction challenge problem being addressed by the Consortium for Advanced Simulation of Light Water Reactors (CASL). Here, two coupling approaches are demonstrated on quarter core problems based on Watts Bar Nuclear Plant Unit 1, Cycle 1: (1) an inline approach in which a one-way coupling is used between neutronics/thermal hydraulics (through MPACT/CTF) and fuel performance (through Bison) but no fuel temperature information is passed back to MPACT/CTF, and (2) a two-way approach (coupled) in which the fuel temperature is passed from Bison to MPACT/CTF. In both approaches, power and temperature distributions from MPACT/CTF are used to inform the Bison simulations for each rod in the core.The demonstrations presented here are the first integrated fuel performance simulations in VERA, which opens many possibilities for future work, including applications to accident-tolerant fuel efforts and transient simulations, which are of critical importance to CASL. These demonstrations also highlight the potential to move away from the current Bison-informed fuel temperature lookup table approach, which is the default in MPACT/CTF simulations, if performance improvements are made in the near future.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗