Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computer experiment”

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 1,513 records · Page 84
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

NeuroCoreX: Brain-Inspired Computing from Code to Circuit

NeuroCoreX is an open-source codebase that enables the implementation of brain-inspired, energy-efficient neuromorphic computing models on FPGA hardware. Designed to support real-time learning, all-to-all neural connectivity, and flexible network architectures, NeuroCoreX offers a hands-on, accessible platform for exploring biologically inspired models of neural computation. It empowers researchers, students, and developers to implement and experiment with adaptive systems—bringing the power of neuromorphic computing to a broader community through a low-cost, scalable, and reconfigurable framework.

Gautam, Ashish [Oak Ridge National Laboratory (ORN↗
DOE OSTI2021

Simulations of Quantum Circuits with Approximate Noise using qsim and Cirq

We introduce multinode quantum trajectory simulations with qsim, an open source high performance simulator of quantum circuits. qsim can be used as a backend of Cirq, a Python software library for writing quantum circuits. We present a novel delayed inner product algorithm for quantum trajectories which can result in an order of magnitude speedup for low noise simulation. We also provide tools to use this framework in Google Cloud Platform, with high performance virtual machines in a single mode or multinode setting. Multinode configurations are well suited to simulate noisy quantum circuits with quantum trajectories. Finally, we introduce an approximate noise model for Google's experimental quantum computing platform and compare the results of noisy simulations with experiments for several quantum algorithms on Google's Quantum Computing Service.

Isakov, Sergei V.↗
NASA NTRS1988

Unsteady transonic aerodynamics and aeroelastic calculations at low-supersonic freestreams

A computational procedure is presented to simulate transonic unsteady flows and corresponding aeroelasticity of wings at low-supersonic freestreams. The flow is modeled by using the transonic small-perturbation theory. The structural equations of motions are modeled using modal equations of motion directly coupled with aerodynamics. Supersonic freestreams are simulated by properly accounting for the boundary conditions based on pressure waves along the flow characteristics in streamwise planes. The flow equations are solved using the time-accurate, alternating-direction implicit finite-difference scheme. The coupled aeroelastic equations of motion are solved by an integration procedure based on the time-accurate, linear-acceleration method. The flow modeling is verified by comparing calculations with experiments for both steady and unsteady flows at supersonic freestreams. The unsteady computations are made for oscillating wings. Comparisons of computed results with experiments show good agreement. Aeroelastic responses are computed for a rectangular wing at Mach numbers ranging from subtransonic to upper-transonic (supersonic) freestreams. The extension of the transonic dip into the upper transonic regime is illustrated.

Guruswamy, Guru P.↗
NASA NTRS1994

Numerical simulation of the edge tone phenomenon

Time accurate Navier-Stokes computations were performed to study a class 2 (acoustic) whistle, the edge tone, and to gain knowledge of the vortex-acoustic coupling mechanisms driving production of these tones. Results were obtained by solving the full Navier-Stokes equations for laminar compressible air flow of a two dimensional jet issuing from a slit interacting with a wedge. Cases considered were determined by varying the distance from the slit to the wedge. Flow speed was kept constant at 1,750 cm/s as was the slit thickness of 0.1 cm, corresponding to conditions in the experiments of Brown. The analytical computations revealed edge tones to be present in four harmonic stages of jet flow instability over the wedge as the jet length was varied from 0.3 to 1.6 cm. Excellent agreement was obtained in all four edge tone stage cases between the present computational results and the experimentally obtained frequencies and flow visualization results of Brown. Specific edge tone generation phenomena and further confirmation of certain theories and empirical formulas concerning these phenomena were brought to light in this analytical simulation of edge tones.

Dougherty, N. S.↗
NASA NTRS2007

Potential of VIIRS Data for Regional Monitoring of Gypsy Moth Defoliation: Implications for Forest Threat Early Warning System

A NASA RPC (Rapid Prototyping Capability) experiment was conducted to assess the potential of VIIRS (Visible/Infrared Imager/Radiometer Suite) data for monitoring non-native gypsy moth (Lymantria dispar) defoliation of forests. This experiment compares defoliation detection products computed from simulated VIIRS and from MODIS (Moderate Resolution Imaging Spectroradiometer) time series products as potential inputs to a forest threat EWS (Early Warning System) being developed for the USFS (USDA Forest Service). Gypsy moth causes extensive defoliation of broadleaved forests in the United States and is specifically identified in the Healthy Forest Restoration Act (HFRA) of 2003. The HFRA mandates development of a national forest threat EWS. This system is being built by the USFS and NASA is aiding integration of needed satellite data products into this system, including MODIS products. This RPC experiment enabled the MODIS follow-on, VIIRS, to be evaluated as a data source for EWS forest monitoring products. The experiment included 1) assessment of MODIS-simulated VIIRS NDVI products, and 2) evaluation of gypsy moth defoliation mapping products from MODIS-simulated VIIRS and from MODIS NDVI time series data. This experiment employed MODIS data collected over the approximately 15 million acre mid-Appalachian Highlands during the annual peak defoliation time frame (approximately June 10 through July 27) during 2000-2006. NASA Stennis Application Research Toolbox software was used to produce MODIS-simulated VIIRS data and NASA Stennis Time Series Product Tool software was employed to process MODIS and MODIS-simulated VIIRS time series data scaled to planetary reflectance. MODIS-simulated VIIRS data was assessed through comparison to Hyperion-simulated VIIRS data using data collected during gypsy moth defoliation. Hyperion-simulated MODIS data showed a high correlation with actual MODIS data (NDVI R2 of 0.877 and RMSE of 0.023). MODIS-simulated VIIRS data for the same date showed moderately high correlation with Hyperion-simulated VIIRS data (NDVI R2 of 0.62 and RMSE of 0.035), even though the datasets were collected about a half an hour apart during changing weather conditions. MODIS products (MOD02, MOD09, and MOD13) and MOD02-simulated VIIRS time series data were used to generate defoliation mapping products based on image classification and image differencing change detection techniques. Accuracy of final defoliation mapping products was assessed by image interpreting over 170 randomly sampled locations found on Landsat and ASTER data in conjunction with defoliation map data from the USFS. The MOD02-simulated VIIRS 400-meter NDVI classification produced a similar overall accuracy (87.28 percent with 0.72 Kappa) to the MOD02 250-meter NDVI classification (86.71 percent with 0.71 Kappa). In addition, the VIIRS 400-meter NDVI, MOD02 250-meter NDVI, and MOD02 500-meter NDVI showed good user and producer accuracies for the defoliated forest class (70 percent) and acceptable Kappa values (0.66). MOD02 and MOD02-simulated VIIRS data both showed promise as data sources for regional monitoring of forest disturbance due to insect defoliation.

Spruce, Joseph P.↗
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 NTRS1991

High-temperature properties of gases

Computational chemistry is being applied at NASA Ames Research Center to a variety of problems in chemistry, physics, and materials sciences. Supercomputers and modern computational chemistry techniques have provided a powerful new tool to help fill NASA's continuing need for information about the properties of gases and materials as well as their interaction. For example, radiative transition probabilities, spectroscopic constants, bond dissociation energies, transport properties, and chemical reaction rates for molecules can be determined computationally just as reliably as by experiment. Recent results on the radiative transition probabilities of the strongly radiating systems of air species, based on state-of-the-art computational chemistry calculations, are presented and compared with experimental data.

Cooper, David M.↗
NASA NTRS2000

Simulation of Hypervelocity Impact on Aluminum-Nextel-Kevlar Orbital Debris Shields

An improved hybrid particle-finite element method has been developed for hypervelocity impact simulation. The method combines the general contact-impact capabilities of particle codes with the true Lagrangian kinematics of large strain finite element formulations. Unlike some alternative schemes which couple Lagrangian finite element models with smooth particle hydrodynamics, the present formulation makes no use of slidelines or penalty forces. The method has been implemented in a parallel, three dimensional computer code. Simulations of three dimensional orbital debris impact problems using this parallel hybrid particle-finite element code, show good agreement with experiment and good speedup in parallel computation. The simulations included single and multi-plate shields as well as aluminum and composite shielding materials. at an impact velocity of eleven kilometers per second.

Fahrenthold, Eric P.↗
NASA NTRS2018

A Numerical and Experimental Investigation of Flow Separation Control over a Wall-Mounted Hump Model

Numerical simulations and wind tunnel experiments were performed to investigate flow separation and its control over a wall-mounted hump model. Three-dimensional, unsteady fluid dynamic simulations were used to capture the relevant flow physics of unsteady flow separation and its control via active and passive flow control methods. Surface pressure measurements and surface oilflow visualization were obtained in the experiments. Current measurements documented a considerably thinner incoming boundary layer, but similar pressure distribution, compared to the reference data reported in the CFD (Computational Fluid Dynamics) Validation (CFDVAL 2004) Experiments. Consistent with the data available in the literature, numerical simulations predicted a longer separation bubble, which manifested itself as a shift in the pressure distributions. Streamwise vortices generated by passive vortex generators increased the suction pressure and provided substantial pressure recovery. Numerical simulations of vortex generators qualitatively predicted the effects of control, and revealed patterns of attached/separated flow regions. On the other hand, active flow control using steady discrete jets reduced, but was not able to eliminate, flow separation for the tested configuration. Although the major disagreement between the numerical and experimental results was in the pressure recovery region due to the inaccurate prediction of separation bubbles, the relative pressure recoveries from the baselines were comparable. Surface oilflow visualization and simulated surface streamlines were qualitatively in good agreement revealing the key flow features.

Koklu, Mehti↗
DOE OSTI2023

Bentonite buffer under high temperature: laboratory experiments and coupled process modeling

Abstract. Bentonite buffer in a geological repository will be simultaneously heated from decaying radioactive waste and hydrated from the surrounding host rock, triggering complex and coupled THMC (thermal–hydrological–mechanical–chemical) processes. Understanding the THMC behavior of bentonite-based engineered barrier system (EBS) is key to the evaluation and prediction of its long-term performance. Studies on the THMC process have been focused on conditions under 100 ∘C, as most design concepts impose a thermal limit of 100 ∘C in bentonite. Recently, studies under high-temperature conditions have been conducted to evaluate the possibility of raising the thermal limit and expanding the data/knowledge base to increase the confidence level. In this abstract, we present a series of bench-scale laboratory experiments at high temperatures (up to 200 ∘C) and the corresponding modeling work. Two sets of column tests were conducted, and each set consisted of two test columns: a control column undergoing only hydration (non-heated) and an experiment column experiencing both heating and hydration (heated). During the experiment, frequent X-ray computed tomography (CT) images were collected to provide a 3D visualization of the density distribution and present the spatiotemporal evolution of (1) hydration/dehydration, (2) clay swelling/shrinkage, (3) displacement, and (4) mineral precipitation. The two sets of tests differ with respect to several experimental conditions, such as bentonite type, compacted density and water content, water chemistry, and hydration pressure, but the important difference is that the first set used bentonite powder with a dry density of 1.28 g cm−3, whereas the second set used granulated bentonite (mixture of pellets and powder) with a dry density of 1.45–1.5 g cm−3. In both sets of experiments, a comprehensive post-dismantling characterization of bentonite samples was carried out after the column tests had been running for 1.5 years. Comparing non-heated and heated columns, the temperature gradient led to lower degree of homogenization of bentonite after bentonite became fully saturated; comparing the first and second sets, granulated and powdered bentonite exhibited drastically different hydration behavior. A THM model with a 2D axisymmetric grid system was used to interpret the data from the first set of tests. The model considers the combined impact of saturation, fluid pressure, and porosity change due to swelling/compression on the spatiotemporal distribution of bulk density and movement of the thermocouple modules. Observations from the tests help us understand the early perturbation of bentonite buffer under high temperature, and data from these tests improve the calibration of key constitutive hydrological and mechanical models and, therefore, enhance the modeling capability with respect to calculating the long-term evolution of bentonite buffer.

Zheng, Liange↗
NASA NTRS2002

The NASA Low-Pressure Turbine Flow Physics Program: A Review

An overview of the NASA Glenn Low-Pressure Turbine (LPT) Flow Physics Program will be presented. The flow in the LPT is unique for the gas turbine. It is characterized by low Reynolds number and high freestream turbulence intensity and is dominated by interplay of three basic mechanisms: transition, separation and wake interaction. The flow of most interest is on the suction surface, where large losses are generated due to separation. The LPT is a large, multistage, heavy, jet engine component that suffers efficiency degradation between takeoff and cruise conditions due to decrease in Reynolds number with altitude. The performance penalty is around 2 points for large commercial bypass engines and as much as 7 points for small, high cruise altitude, military engines. The gas-turbine industry is very interested in improving the performance of the LPT and in reducing its weight, part count and cost. Many improvements can be accomplished by improved airfoil design, mainly by increasing the airfoil loading that can yield reduction of airfoils and improved performance. In addition, there is a strong interest in reducing the design cycle time and cost. Key enablers of the needed improvements are computational tools that can accurately predict LPT flows. Current CFD tools in use cannot yet satisfactorily predict the unsteady, transitional and separated flow in the LPT. The main reasons are inadequate transition & turbulence models and incomplete understanding of the LPT flow physics. NASA Glenn has established its LPT program to answer these needs. The main goal of the program is to develop and assess models for unsteady CFD of LPT flows. An approach that consists of complementing and augmenting experimental and computational work elements has been adopted. The work is performed in-house and by several academic institutions, in cooperation and interaction with industry. The program was reviewed at the Minnowbrook II meeting in 1997. This review will summarize the progress that was made since and will introduce newly started projects. The LPT program is focused on three areas: acquisition of experimental and numerical databases and on modeling and computation. Priority was initially given to experiments. There are three classes of experiments: simulated LPT passages, linear cascade, both with and without wakes, and low-speed rotating rig. They are being conducted as follows: At NASA GRC on a flat surface with blade pressure distribution, at the US Naval Academy on a curved surface. The addition of wakes is studied at the University of Minnesota in a curved passage with a retractable wake generator, and at Texas A&M University in a linear cascade with continuously running wake generator. The pressure distribution of the Pratt & Whitney blade 'Pak B' is used in all these experiments. Experiments have been performed also in the GEAE Low-Speed Rotating Turbine (LSRT) rig with GE-designed airfoils. Work on numerically generated database is in progress at the University of Kentucky, using the DNS/LES code LESTool developed there. Turbulence/transition model assessment and development is performed also at the University of Kentucky, where a new intermittency transport model was developed and many experimental test cases have been numerically computed. Assessments of models using simulations of multistage LPT experiments were performed at Virginia Commonwealth University using the Corsair code. Work on suction surface separation delay, using passive and active flow-control, has also been initiated. Following the overview, Principal Investigators attending the workshop will present in detail several of the projects supported by NASA.

Ashpis, David E.↗
DOE OSTI2022

A 1D–0D–3D coupled model for simulating blood flow and transport processes in breast tissue

In this work, we present mixed dimensional models for simulating blood flow and transport processes in breast tissue and the vascular tree supplying it. These processes are considered, to start from the aortic inlet to the capillaries and tissue of the breast. Large variations in biophysical properties and flow conditions exist in this system necessitating the use of different flow models for different geometries and flow regimes. In total, we consider four different model types. First, a system of 1D nonlinear hyperbolic partial differential equations (PDEs) is considered to simulate blood flow in larger arteries with highly elastic vessel walls. Second, we assign 1D linearized hyperbolic PDEs to model the smaller arteries with stiffer vessel walls. The third model type consists of ODE systems (0D models). It is used to model the arterioles and peripheral circulation. Finally, homogenized 3D porous media models are considered to simulate flow and transport in capillaries and tissue within the breast volume. Sink terms are used to account for the influence of the venous and lymphatic systems. Combining the four model types, we obtain two different 1D–0D–3D coupled models for simulating blood flow and transport processes: The first model results in a fully coupled 1D–0D–3D model covering the complete path from the aorta to the breast combining a generic arterial network with a patient specific breast network and geometry. The second model is a reduced one based on the separation of the generic and patient specific parts. The information from a calibrated fully coupled model is used as inflow condition for the patient specific sub-model allowing a significant computational cost reduction. Several numerical experiments are conducted to calibrate the generic model parameters and to demonstrate realistic flow simulations compared to existing data on blood flow in the human breast and vascular system. Moreover, we use two different breast vasculature and tissue data sets to illustrate the robustness of our reduced sub-model approach.

multiscale blood flow models↗
NASA NTRS1986

Numerical experiments with a symmetric high-resolution shock-capturing scheme

Characteristic-based explicit and implicit total variation diminishing (TVD) schemes for the two-dimensional compressible Euler equations have recently been developed. This is a generalization of recent work of Roe and Davis to a wider class of symmetric (non-upwind) TVD schemes other than Lax-Wendroff. The Roe and Davis schemes can be viewed as a subset of the class of explicit methods. The main properties of the present class of schemes are that they can be implicit, and, when steady-state calculations are sought, the numerical solution is independent of the time step. In a recent paper, a comparison of a linearized form of the present implicit symmetric TVD scheme with an implicit upwind TVD scheme originally developed by Harten and modified by Yee was given. Results favored the symmetric method. It was found that the latter is just as accurate as the upwind method while requiring less computational effort. Currently, more numerical experiments are being conducted on time-accurate calculations and on the effect of grid topology, numerical boundary condition procedures, and different flow conditions on the behavior of the method for steady-state applications. The purpose here is to report experiences with this type of scheme and give guidelines for its use.

Yee, H. C.↗
NASA NTRS2002

Computational Investigation of the Performance and Back-Pressure Limits of a Hypersonic Inlet

A computational analysis of Mach 6.2 operation of a hypersonic inlet with rectangular-to-elliptical shape transition has been performed. The results of the computations are compared with experimental data for cases with and without a manually imposed back-pressure. While the no-back-pressure numerical solutions match the general trends of the data, certain features observed in the experiments did not appear in the computational solutions. The reasons for these discrepancies are discussed and possible remedies are suggested. Most importantly, however, the computational analysis increased the understanding of the consequences of certain aspects of the inlet design. This will enable the performance of future inlets of this class to be improved. Computational solutions with back-pressure under-estimated the back-pressure limit observed in the experiments, but did supply significant insight into the character of highly back-pressured inlet flows.

Smart, Michael K.↗
NASA NTRS1986

Engineering calculations for communications satellite systems planning

Observed solution times were analyzed for the extended gradient and cyclic coordinate search procedures. The times used in the analysis come from computer runs made during a previously-reported experiment conducted to assess the quality of the solutions to a BSS synthesis problem found by the two search methods. The results of a second experiment with a Fixed Satellite Service (FSS) test problem are also presented. Computational results are summarized for mixed integer programming approaches for solving FSS synthesis problems. A promising heuristic algorithm is described. A synthesis model is discussed for orbital arc allotment optimization. Research plans for the near future are also presented.

Reilly, C. H.↗
NASA NTRS1959

Lift-Drag Ratios for an Arrow Wing With Bodies at Mach Number 3

Force and moment characteristics, including lift-drag ratios, have been measured for bodies of circular and elliptic cross section alone and combined with a warped arrow wing. The test Mach number was 2.94, and the Reynolds number was 3.5 x 10(exp 6) (based on wing mean aerodynamic chord). The experimental results show that for equal volume the use of an elliptical body can result in a noticeably higher maximum lift-drag ratio than that obtained through use of a circular body. Methods for estimating the aerodynamic characteristics have been assessed by comparing computed with experimental results. Because of good agreement of the predictions with experiment, maximum lift-drag ratios have been computed for the arrow wing in combination with bodies of various sizes. These calculations have shown that, for an efficient wing-body combination, little loss in maximum lift-drag ratio results from considerable extension of afterbody length. For example, for a wing-body configuration having a maximum lift-drag ratio of about 7.1, a loss in maximum lift-drag ratio of less than 0.2 results from a 40-percent increase in body volume by extension of afterbody length. It also appears that with body length fixed, maximum lift-drag ratio decreases almost linearly with increase in body diameter. For a wing- body combination employing a body of circular cross section, a decrease in maximum lift-drag ratio from about 9.1 for zero body diameter to about 4.6 for a body diameter of 13.5 percent of the body length was computed.

Jorgensen, Leland H.↗
NASA NTRS2021

Transition Analysis for Isolated Trips on BOLT-II Wind-Tunnel and Flight Configuration

The BOLT-II configuration designed for an upcoming flight experiment exemplifies the complexities of hypersonic boundary layer transition in the presence of multiple and potentially interacting instability mechanisms. The present paper is related to the NASA-led roughness experiment on the secondary side of the BOLT-II configuration. First, computations are performed to aid mean-flow based correlations for the transition measurements acquired during an experimental campaign in the Ludwieg Tube facility at the United States Air Force Academy. To that end, both laminar and turbulent flow computations are performed for selected conditions of the experiment with a 1:3 scale model of the slightly shorter BOLT configuration. The variation in well-known correlation parameters for transition due to isolated roughness elements with respect to element position, height, and flow Reynolds number is delineated. Furthermore, the wake instability theory is used to examine the amplification of unstable disturbances behind the centerline trip at the design conditions of the BOLT-II flight experiment and for its dynamically scaled version on the 1:3 scale model of the wind tunnel experiment. To our knowledge, this study represents the first analysis of the wake instabilities in the presence of a prominent entropy layer outside the boundary layer. At the nominal flight design condition of Re ∞ ) = 5.44 and Re ∞ = 2.5 x 10 6 /m, a diamond planform (“pizza-box”) trip with k/δ = 0.70 and planform-halfwidth-to-height ratio of b/k = 3.0 at X/L = 0.254 induces a strong distortion of the boundary-layer flow in the vicinity of the symmetry plane, such that the wake structure includes a centerline streak produced by the flow retardation behind the trip along with multiple additional, prominent vorticity structures representing the necklace vortex system. Stability results reveal that each of these streaks can support a strong amplification of wake instabilities within a frequency band of 95 kHz–255 kHz, with a peak N-factor of greater than 15 beyond X/L ≈ 0.42. The dynamically scaled wind tunnel configuration with Re ∞ = 20.25 x 10 6 /m is also found to support N > 15 at similar axial stations.

BOLT↗