Engineering PapersSearch

SEARCH · Engineering Papers

Results for “Computational Physics and Simulations”

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 901 records · Page 50

Transonic Navier-Stokes wing solutions using a zonal approach. Part 2: High angle-of-attack simulation

A computer code is under development whereby the thin-layer Reynolds-averaged Navier-Stokes equations are to be applied to realistic fighter-aircraft configurations. This transonic Navier-Stokes code (TNS) utilizes a zonal approach in order to treat complex geometries and satisfy in-core computer memory constraints. The zonal approach has been applied to isolated wing geometries in order to facilitate code development. Part 1 of this paper addresses the TNS finite-difference algorithm, zonal methodology, and code validation with experimental data. Part 2 of this paper addresses some numerical issues such as code robustness, efficiency, and accuracy at high angles of attack. Special free-stream-preserving metrics proved an effective way to treat H-mesh singularities over a large range of severe flow conditions, including strong leading-edge flow gradients, massive shock-induced separation, and stall. Furthermore, lift and drag coefficients have been computed for a wing up through CLmax. Numerical oil flow patterns and particle trajectories are presented both for subcritical and transonic flow. These flow simulations are rich with complex separated flow physics and demonstrate the efficiency and robustness of the zonal approach.

Chaderjian, N. M.

Transonic Navier-Stokes wing solutions using a zonal approach. Part 2: High angle-of-attack simulation

A computer code is under development whereby the thin-layer Reynolds-averaged Navier-Stokes equations are to be applied to realistic fighter aircraft configurations. This transonic Navier-Stokes code (TNS) utilizes a zonal approach in order to treat complex geometries and satisfy in-core computer memory constraints. The zonal approach was applied to isolated wing geometries in order to facilitate code development. The TNS finite difference algorithm, zonal methodology, and code validation with experimental data is addressed. Also addressed are some numerical issues such as code robustness, efficiency, and accuracy at high angles of attack. Special free-stream-preserving metrics proved an effective way to treat H-mesh singularities over a large range of severe flow conditions, including strong leading edge flow gradients, massive shock induced separation, and stall. Furthermore, lift and drag coefficients were computed for a wing up through CLmax. Numerical oil flow patterns and particle trajectories are presented both for subcritical and transonic flow. These flow simulations are rich with complex separated flow physics and demonstrate the efficiency and robustness of the zonal approach.

Chaderjian, Neal M.

Identification of structural interface characteristics using component mode synthesis

The inability to adequately model connections has limited the ability to predict overall system dynamic response. Connections between structural components are often mechanically complex and difficult to accurately model analytically. Improved analytical models for connections are needed to improve system dynamic predictions. This study explores combining Component Mode synthesis methods for coupling structural components with Parameter Identification procedures for improving the analytical modeling of the connections. Improvements in the connection properties are computed in terms of physical parameters so the physical characteristics of the connections can be better understood, in addition to providing improved input for the system model. Two sample problems, one utilizing simulated data, the other using experimental data from a rotor dynamic test rig are presented.

Huckelbridge, A. A.

Identification of structural interface characteristics using component mode synthesis

The inability to adequately model connections has limited the ability to predict overall system dynamic response. Connections between structural components are often mechanically complex and difficult to accurately model analytically. Improved analytical models for connections are needed to improve system dynamic predictions. This study explores combining Component Mode synthesis methods for coupling structural components with Parameter Identification procedures for improving the analytical modeling of the connections. Improvements in the connection properties are computed in terms of physical parameters so the physical characteristics of the connections can be better understood, in addition to providing improved input for the system model. Two sample problems, one utilizing simulated data, the other using experimental data from a rotor dynamic test rig, are presented.

Huckelbridge, A. A.

Numerical tests of a modified full implicit continuous Eulerian (FICE) scheme with projected normal characteristic boundary conditions for MHD flows

A numerical method has been developed based on a modified full implicit continuous Eulerian (FICE) scheme and projected normal characteristic boundary conditions for simulating MHD flows which undergo a long process of evolution. An astrophysical flow is chosen for illustration of this procedure, and numerical tests are made to verify the computational stability and physically realistic solution. Three computational tests have been accomplished; they are tests of solving methods, characteristic boundary condition, and time steps. The tests show that the program from the modified FICE scheme with proper boundary conditions and time steps can be made numerically stable for a time long enough to obtain physically plausible solutions.

Wu, S. T.

Identification of structural interface characteristics using component mode synthesis

The inability to adequately model connections has limited the ability to predict overall system dynamic response. Connections between structural components are often mechanically complex and difficult to accurately model analytically. Improved analytical models for connections are needed to improve system dynamic predictions. This study explores combining Component Mode synthesis methods for coupling structural components with Parameter Identification procedures for improving the analytical modeling of the connections. Improvements in the connection properties are computed in terms of physical parameters so the physical characteristics of the connections can be better understood, in addition to providing improved input for the system model. Two sample problems, one utilizing simulated data, the other using experimental data from a rotor dynamic test rig, are presented.

Huckelbridge, A. A.

Local-Area-Network Simulator

Local Area Network Extensible Simulator (LANES) computer program provides method for simulating performance of high-speed local-area-network (LAN) technology. Developed as design and analysis software tool for networking computers on board proposed Space Station. Load, network, link, and physical layers of layered network architecture all modeled. Mathematically models according to different lower-layer protocols: Fiber Distributed Data Interface (FDDI) and Star*Bus. Written in FORTRAN 77.

Gibson, Jim

SRS Computer Animation and Drive Train System

The spinning rocket simulator (SRS) is an ongoing project at Oral Roberts University. The goal of the SRS is to gather crucial data concerning a spinning rocket under thrust for the purpose of analysis and correction of the coning motion experienced by this type of spacecraft maneuver. The computer animation simulates a virtual, scale model of the component of the SRS that represents the spacecraft itself. This component is known as the (VSM), or virtual spacecraft model. During actual physical simulation, this component of the SRS will experience a coning. The goal of the animation is to cone the VSM within that range to accurately represent the motion of the actual simulator. The drive system of the SRS is the apparatus that turns the actual simulator. It consists of a drive motor, motor mount and chain to power the simulator into motion. The motor mount is adjustable and rigid for high torque application. A digital stepper motor controller actuates the main drive motor for linear acceleration. The chain transfers power from the motor to the simulator via sprockets on both ends.

Arthun, Daniel

Multiagent Work Practice Simulation: Progress and Challenges

Modeling and simulating complex human-system interactions requires going beyond formal procedures and information flows to analyze how people interact with each other. Such work practices include conversations, modes of communication, informal assistance, impromptu meetings, workarounds, and so on. To make these social processes visible, we have developed a multiagent simulation tool, called Brahms, for modeling the activities of people belonging to multiple groups, situated in a physical environment (geographic regions, buildings, transport vehicles, etc.) consisting of tools, documents, and a computer system. We are finding many useful applications of Brahms for system requirements analysis, instruction, implementing software agents, and as a workbench for relating cognitive and social theories of human behavior. Many challenges remain for representing work practices, including modeling: memory over multiple days, scheduled activities combining physical objects, groups, and locations on a timeline (such as a Space Shuttle mission), habitat vehicles with trajectories (such as the Shuttle), agent movement in 3D space (e.g., inside the International Space Station), agent posture and line of sight, coupled movements (such as carrying objects), and learning (mimicry, forming habits, detecting repetition, etc.).

Clancey, William J.

Predictive Models for Semiconductor Device Design and Processing

The device feature size continues to be on a downward trend with a simultaneous upward trend in wafer size to 300 mm. Predictive models are needed more than ever before for this reason. At NASA Ames, a Device and Process Modeling effort has been initiated recently with a view to address these issues. Our activities cover sub-micron device physics, process and equipment modeling, computational chemistry and material science. This talk would outline these efforts and emphasize the interaction among various components. The device physics component is largely based on integrating quantum effects into device simulators. We have two parallel efforts, one based on a quantum mechanics approach and the second, a semiclassical hydrodynamics approach with quantum correction terms. Under the first approach, three different quantum simulators are being developed and compared: a nonequlibrium Green's function (NEGF) approach, Wigner function approach, and a density matrix approach. In this talk, results using various codes will be presented. Our process modeling work focuses primarily on epitaxy and etching using first-principles models coupling reactor level and wafer level features. For the latter, we are using a novel approach based on Level Set theory. Sample results from this effort will also be presented.

Meyyappan, Meyya

Fluid Structure Interaction of Parachutes in Supersonic Planetary Entry

A research program to provide physical insight into disk-gap-band parachute operation in the supersonic regime on Mars was conducted. The program included supersonic wind tunnel tests, computational fluid dynamics and fluid structure interaction simulations. Specifically, the nature and cause of the "area oscillation" phenomenon were investigated to determine the scale, aerodynamic, and aero-elastic dependence of the supersonic parachute collapse and re-inflation event. A variety of non-intrusive, temporally resolved, and high resolution diagnostic techniques were used to interrogate the flow and generate validation datasets. The results of flow visualization, particle image velocimetry, load measurements, and photogrammetric reconstruction will be presented. Implications to parachute design, use, and verification will also be discussed.

parachute

Verification and Validation of Elastodynamic Simulation Software for Aerospace Research

Physics-based simulation of nondestructive evaluation (NDE) inspection can help to advance the inspectability and reliability of mechanical systems. However, NDE simulations applicable to non-idealized mechanical components often require large compute domains and long run times. This has prompted development of custom NDE simulation software tailored to high performance computing (HPC) hardware. Verification and validation (V&V) is an integral part of developing this software to ensure implementations are robust and applicable to inspection problems, producing tools and simulations suitable for computational NDE research. This presentation addresses factors common to V&V of several elastodynamic simulation codes applicable to ultrasonic NDE. Examples are drawn from in-house simulation software at NASA Langley Research Center, ranging from ensuring reliability in a 1D heterogeneous media wave equation solver to the V&V needs of 3D cluster-parallel elastodynamic software. Factors specific to a research environment are addressed, where individual simulation results can be as relevant as the software product itself. Distinct facets of V&V are discussed including testing to establish software reliability, employing systematic approaches for consistency with fundamental conservation laws, establishing the numerical stability of algorithms, and demonstrating concurrence with empirical data. This talk also addresses V&V practices for small groups of researchers. This includes establishing resources (e.g. time and personnel) for V&V during project planning to mitigate and control the risk of setbacks. Similarly, we identify ways for individual researchers to use V&V during simulation software development itself to both speed up the development process and reduce incurred technical debt.

NDE

Charged-particle transport in high energy density plasmas

This Special Topic Collection grew out of two gatherings of researchers active in the high energy density (HED) physics community: a mini-conference on charged-particle transport in HED plasma held during the 64th annual meeting of the American Physical Society's Division of Plasma Physics (Spokane, WA, November 2022) and a dedicated charged-particle transport coefficient code comparison workshop (Livermore, CA, July 2023). These gatherings provided opportunities for theoretical, computational, and experimental researchers to discuss the state of the field, including current capabilities and methods, needs of hydrodynamic simulations, and frontiers for future research. Finally, this special issue collects a total of 13 research and review articles on charged-particle transport in HED plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY

Effect of artificial viscosity on shocked particle-laden flows for staggered grid Lagrangian methods

Abstract Shocked particle-laden flows are important to many natural and industrial processes. When simulating these systems, artificial viscosity is often required to prevent numerical artifacts, such as ringing, from arising in the pressure and density fields. The linear and quadratic coefficients of the artificial viscosity determine the amount of smoothing that occurs in these fields. For particle-laden flows, however, many of the fluid–particle interaction forces, for example, the pressure gradient force and unsteady forces, depend on gradients in the fluid fields. Furthermore, while the shock passes over a particle, these forces can be more dominant than drag. This means that the artificial viscosity coefficients affect how a particle and fluid interact when simulating shocked particle systems. Here this effect is investigated for isolated particles and for a particle curtain using a staggered grid Lagrangian approach. The artificial viscosity coefficients have a significant impact on the maximum force that a fluid imparts to a particle, which is important for determining whether a particle will break up in response to the shock. Furthermore, it is found that the density ratio between the particle and the fluid is important in determining whether the artificial viscosity coefficients have a significant impact on the particle’s motion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC

Enabling in-time Prognostics with Surrogate Modeling through Physics-enhanced Dynamic Mode Decomposition Method

Computational models provide essential quantitative tools for assessing and predicting the health and performance of physical systems. However, high-fidelity models are rarely used in real-time operations or large optimization loops, due to their time-intensive nature. A common approach to improving computational efficiency of prognosis is to employ surrogate models. Such models can significantly decrease computation time for some accuracy loss. In this context, use of Dynamic Mode Decomposition (DMD) is proposed to generate surrogate models for lithium-ion (Li-ion) battery discharge. DMD has been suggested and used successfully in the area of fluid dynamics for over a decade, but it has not been applied to the PHM domain, where far-ahead prediction of nonlinear behavior is crucial to propagate faults or predict Remaining Useful Life (RUL). For Li-ion battery health management, the standard application of DMD using only the observable quantities of interest was unable to capture the nonlinear discharge of batteries exhibited in lab testing. The Koopman theory, however, provides a mechanism to tradeoff low dimensional nonlinear models with high-dimensional linear ones in a DMD framework, by augmenting nonlinear state variables into the system representation. In this way, DMD allows for configurable simulation accuracy dependent on the dimensionality of the Koopman operator. For battery health management, we augmented the observable variables with the hidden states of a higher-fidelity physics model to build the DMD surrogate. In comparison to a high-fidelity model, the surrogate improved computational efficiency with only a minimal loss of accuracy, and enabled long-term prognostics horizons. A generalized method for this was implemented in the prog models python package.

prognostics and health management

Simulation of radar reflectivity and surface measurements of rainfall

Raindrop size distributions (RSDs) are often estimated using surface raindrop sampling devices (e.g., disdrometers) or optical array (2D-PMS) probes. A number of authors have used these measured distributions to compute certain higher-order RSD moments that correspond to radar reflectivity, attenuation, optical extinction, etc. Scatter plots of these RSD moments versus disdrometer-measured rainrates are then used to deduce physical relationships between radar reflectivity, attenuation, etc., which are measured by independent instruments (e.g., radar), and rainrate. In this paper RSDs of the gamma form as well as radar reflectivity (via time series simulation) are simulated to study the correlation structure of radar estimates versus rainrate as opposed to RSD moment estimates versus rainrate. The parameters N0, D0 and m of a gamma distribution are varied over the range normally found in rainfall, as well as varying the device sampling volume. The simulations are used to explain some possible features related to discrepancies which can arise when radar rainfall measurements are compared with surface or aircraft-based sampling devices.

Chandrasekar, V.

Accelerating kinetic plasma simulations with machine-learning-generated initial conditions

Computational models of plasma technologies often solve for the system operating conditions by time-stepping an initial value problem to a quasi-steady solution. However, the strongly nonlinear and multi-timescale nature of plasma dynamics often necessitate millions, or even hundreds of millions, of steps to reach convergence, reducing the effectiveness of these simulations for computer-aided engineering. We consider acceleration of kinetic plasma simulations via data-driven machine-learning-generated initial conditions, which initialize the simulations close to their final quasi-steady-state, thereby reducing the number of steps to reach convergence. Three machine-learning models are developed to predict the density and ion kinetic profiles of capacitively coupled plasma discharges relevant to the microelectronics industry. The models are trained on kinetic simulations over a range of device operating frequencies and pressures. Best performance was observed when simulations were initialized with ion kinetic profiles generated by a convolutional neural network, reducing the mean number of steps to reach convergence by 17.1× when compared to initialization with a zero-dimensional global model. We also outline a workflow for continuous data-driven model improvement and simulation speedup, with the aim of generating sufficient data for full device digital twins.

Artificial neural networks

Numerical simulation of rotorcraft

The objective of the research is to develop and validate accurate, user-oriented viscous CFD codes (with inviscid options) for three-dimensional, unsteady aerodynamic flows about arbitrary rotorcraft configurations. Unsteady, three-dimensional Euler and Navier-Stokes codes are developed, adapted, and extended to rotor-body combinations. Flow solvers are coupled with zonal grid topologies, including rotating and nonrotating blocks. Special grid clustering and wave-fitting techniques were developed to capture low-level radiating acoustic waves. Significant progress was made in computing the propagation of acoustic waves due to the interaction of a concentrated vortex and a helicopter airfoil. The need for higher-order schemes was firmly established in relatively inexpensive two-dimensional calculations. In three dimensions, the number of grid points required to capture the low-level acoustic waves becomes very large, so that large supercomputer memory becomes essential. Good agreement was obtained between the numerical results obtained with a thin-layer Navier-Stokes code and experimental data from a model rotor. In addition, several nonrotating configurations that are sometimes proposed to simulate rotor blade tips in conventional wind tunnels were examined, and the complex flow around the radical tip shape of the world's fastest helicopter is under investigation. These studies demonstrate the flexibility and power of CFD to gain physical insight, study novel ideas, and examine various possibilities that might be difficult or impossible to set up in physical experiments. As a prelude to studies of rotor-body aerodynamic interactions, a preliminary grid topology and moving-interface strategy were developed. A new Euler/Navier-Stokes code using these techniques computes the vortical wake directly, rather than modeling it, as in most previous rotorcraft studies. Several hover cases were run for conventional and advanced-geometry blades. Numerical schemes using multi-zones and/or adaptive grids appear to be necessary to simulate the complex vortical flows in rotor wakes.

Mccroskey, William J.