SEARCH · Engineering Papers

Results for “numerical computation”

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 343 records · Page 19

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Massively parallel modeling and inversion of electrical resistivity tomography data using PFLOTRAN

Abstract. Electrical resistivity tomography (ERT) is a broadly accepted geophysical method for subsurface investigations. Interpretation of field ERT data usually requires the application of computationally intensive forward modeling and inversion algorithms. For large-scale ERT data, the efficiency of these algorithms depends on the robustness, accuracy, and scalability on high-performance computing resources. In this regard, we present a robust and highly scalable implementation of forward modeling and inversion algorithms for ERT data. The implementation is publicly available and developed within the framework of PFLOTRAN, an open-source, state-of-the-art massively parallel subsurface flow and transport simulation code. The forward modeling is based on a finite-volume discretization of the governing differential equations, and the inversion uses a Gauss–Newton optimization scheme. To evaluate the accuracy of the forward modeling, two examples are first presented by considering layered (1D) and 3D earth conductivity models. The computed numerical results show good agreement with the analytical solutions for the layered earth model and results from a well-established code for the 3D model. Inversion of ERT data, simulated for a 3D model, is then performed to demonstrate the inversion capability by recovering the conductivity of the model. To demonstrate the parallel performance of PFLOTRAN's ERT process model and inversion capabilities, large-scale scalability tests are performed by using up to 131 072 processes on a leadership class supercomputer. These tests are performed for the two most computationally intensive steps of the ERT inversion: forward modeling and Jacobian computation. For the forward modeling, we consider models with up to 122 ×106 degrees of freedom (DOFs) in the resulting system of linear equations and demonstrate that the code exhibits almost linear scalability on up to 10 000 DOFs per process. On the other hand, the code shows superlinear scalability for the Jacobian computation, mainly because all computations are fairly evenly distributed over each process with no parallel communication.

58 GEOSCIENCES↗

NASA NTRS1985

Three-dimensional, transonic rotor flow field reconstructed from holographic interferogram data

Holographic interferometry and computer-assisted tomography (CAT) are used to determine the transonic flow field of a model rotor blade in hover. A pulsed ruby laser records 40 interferograms with a 61 cm-diam view field near the model rotor-blade tip operating at a tip Mach number of 0.90. After digitizing the interferograms and extracting fringe-order functions, the data are transferred to a CAT code. The CAT code then calculates pressure coefficients in several planes above the blade surface. The values from the holography-CAT method compare favorably with previously obtained numerical computations and laser velocimeter measurements at most locations near the blade tip. The results demonstrate the technique's potential for three-dimensional transonic rotor flow studies.

Kittleson, J. K.↗

NASA NTRS1985

Reconstruction of a 3-dimensional transonic rotor flow field from holographic interferogram data

Holographic interferometry and computer-assisted tomography (CAT) are used to determine the transonic velocity field of a model rotor blade in hover. A pulsed ruby laser recorded 40 interferograms with a 2-ft-diam view field near the model rotor-blade tip operating at a tip Mach number of 0.90. After digitizing the interferograms and extracting fringe-order functions, the data are transferred to a CAT code. The CAT code then calculates the perturbation velocity in seeral planes above the blade surface. The values from the holography-CAT method compare favorably with previously obtained numerical computations in most locations near the blade tip. The results demonstrate the technique's potential for three-dimensional transonic rotor flow studies.

Yu, Y. H.↗

NASA NTRS1985

Reconstruction of a three-dimensional, transonic rotor flow field from holographic interferogram data

Holographic interferometry and computer-assisted tomography (CAT) are used to determine the transonic velocity field of a model rotor blade in hover. A pulsed ruby laser recorded 40 interferograms with a 2-ft-diam view field near the model rotor-blade tip operating at a tip Mach number of 0.90. After digitizing the interferograms and extracting fringe-order functions, the data are transferred to a CAT code. The CAT code then calculates the perturbation velocity in several planes above the blade surface. The value from the holography-CAT method compare favorably with previously obtained numerical computations in most locations near the blade tip. The results demonstrate the technique's potential for three-dimensional transonic rotor flow studies.

Yu, Y. H.↗

NASA NTRS1993

Normal modes of the world's oceans: A numerical investigation using Proudman functions

The numerical modeling of the normal modes of the global oceans is addressed. The results of such modeling could be expected to serve as a guide in the analysis of observations and measurements intended to detect these modes. The numerical computation of normal modes of the global oceans is a field in which several investigations have obtained results during the past 15 years. The results seem to be model-dependent to an unsatisfactory extent. Some modeling areas, such as higher resolution of the bathymetry, inclusion of self-attraction and loading, the role of the Arctic Ocean, and systematic testing by means of diagnostic models are addressed. The results show that the present state of the art is such that a final solution to the normal mode problem still lies in the future. The numerical experiments show where some of the difficulties are and give some insight as to how to proceed in the future.

Sanchez, Braulio V.↗

DOE OSTI2024

Correlation function distributions for O ( N ) lattice field theories in the disordered phase

Numerical computations in strongly interacting quantum field theories are often performed using Monte Carlo sampling methods. A key task in these calculations is to estimate the value of a given physical quantity from the distribution of stochastic samples that are generated using the Monte Carlo method. Typically, the sample mean and sample variance are used to define the expectation values and uncertainties of computed quantities. However, the Monte Carlo sample distribution contains more information than these basic properties, and it is useful to investigate it more generally. In this work, the exact form of the probability distributions of two-point correlation functions at zero momentum in O ( N ) lattice field theories in the disordered phase and in infinite volume are determined. These distributions allow for a robust investigation of the efficacy of the Monte Carlo sampling procedure and are shown also to allow for improved estimators of the target physical quantity to be constructed. The theoretical expectations are shown to agree with numerical calculations in the O ( 2 ) model. Published by the American Physical Society 2024

Astronomy & Astrophysics↗

NASA NTRS1987

Transonic airfoil computation using the integral equation with and without embedded Euler domains

Two transonic computational schemes which are based on the Integral Equation Formulation of the full potential equation were presented. The first scheme is a Shock Capturing-Shock Fitting (SCSF) scheme which uses the full potential equation throughout with the exception of the shock wave where the Rankine-Hugoniot relations are used to cross and fit the shock. The second scheme is an Integral Equation with Embedded Euler (IEEE) scheme which uses the full potential equation with an embedded region where the Euler equations are used. The two schemes are applied to several transonic airfoil flows and the results were compared with numerous computational results and experimental domains with fine grids. The SCSF-scheme is restricted to flows with weak shock, while the IEEE-scheme can handle strong shocks. Currently, the IEEE scheme is applied to other transonic flows with strong shocks as well as to unsteady pitching oscillations.

Kandil, Osama A.↗

NASA NTRS1988

The development of flux-split algorithms for flows with non-equilibrium thermodynamics and chemical reactions

A finite-volume method for the numerical computation of flows with nonequilibrium thermodynamics and chemistry is presented. A thermodynamic model is described which simplifies the coupling between the chemistry and thermodynamics and also results in the retention of the homogeneity property of the Euler equations (including all the species continuity and vibrational energy conservation equations). Flux-splitting procedures are developed for the fully coupled equations involving fluid dynamics, chemical production and thermodynamic relaxation processes. New forms of flux-vector split and flux-difference split algorithms are embodied in a fully coupled, implicit, large-block structure, including all the species conservation and energy production equations. Several numerical examples are presented, including high-temperature shock tube and nozzle flows. The methodology is compared to other existing techniques, including spectral and central-differenced procedures, and favorable comparisons are shown regarding accuracy, shock-capturing and convergence rates.

Grossman, B.↗

NASA NTRS1987

Computation of turbulent flows using an extended k-epsilon turbulence closure model

An extended kappa-epsilon turbulence model is proposed and tested with successful results. An improved transport equation for the rate of dissipation of the turbulent kinetic energy, epsilon, is proposed. The proposed model gives more effective response to the energy production rate than does the standard kappa-epsilon turbulence model. An extra time scale of the production range is included in the dissipation rate equation. This enables the present model to perform equally well for several turbulent flows with different characteristics, e.g., plane and axisymmetric jets, turbulent boundary layer flow, turbulent flow over a backward-facing step, and a confined turbulent swirling flow. A second-order accurate finite difference boundary layer code and a nearly second-order accurate finite difference elliptic flow solver are used for the present numerical computations.

Chen, Y.-S.↗

NASA NTRS2006

Stability of the Wave Bearing on an Elastic Support

Numerical computation predicts that an elastic support can substantially improve the stability of the wave bearing if the dynamic stiffness and damping of this support are in a specific range of values. To experimentally validate this prediction, the housing of a gas bearing was mounted on elastic O-rings and the threshold of sub-synchronous whirl motion was experimentally observed when the bearing runs unloaded with a rotating speed up to 30,000 RPM. The O-ring system was also dynamically characterized by measuring its stiffness and damping at various frequencies up to 500 Hz. Good correlation exists between the experimental data and numerical prediction.

Dimofte, Florin↗

NASA NTRS1951

General method and thermodynamic tables for computation of equilibrium composition and temperature of chemical reactions

A rapidly convergent successive approximation process is described that simultaneously determines both composition and temperature resulting from a chemical reaction. This method is suitable for use with any set of reactants over the complete range of mixture ratios as long as the products of reaction are ideal gases. An approximate treatment of limited amounts of liquids and solids is also included. This method is particularly suited to problems having a large number of products of reaction and to problems that require determination of such properties as specific heat or velocity of sound of a dissociating mixture. The method presented is applicable to a wide variety of problems that include (1) combustion at constant pressure or volume; and (2) isentropic expansion to an assigned pressure, temperature, or Mach number. Tables of thermodynamic functions needed with this method are included for 42 substances for convenience in numerical computations.

Huff, Vearl N↗

NASA NTRS1993

On acoustic radiation from a vibrating panel

An experimental and numerical study of the response and radiation from an acoustically loaded aircraft panel is presented. In the experiment, the panel is excited by a normally incident, harmonic wave. Measurements of the panel response and the resulting transmitted pressure are made in both the near- and acoustic far-fields. The numerical computations model the experiment and in particular account for the full coupling between the panel and the surrounding three dimensional acoustic fluid. The results demonstrate that for a sufficiently high excitation level, the panel response becomes nonlinear. The nonlinearity is characterized by the appearance of harmonics and subharmonics in the power spectral densities of the panel motion and consequently in the resulting acoustic radiation. The primary characteristic of the far-field acoustic radiation is the increase in harmonic content relative to the fundamental and subharmonics with increasing distance from the panel. This is shown to be a result of linear and weakly nonlinear wave propagation effects. The numerical results show that the radiated far-field pressure is strongly dependent on the position of the measurement point with respect to the panel center. The experimental and numerical results are in good qualitative agreement.

Frendi, Abdelkader↗

NASA NTRS2018

Instability Wave–Streak Interactions in a High Mach Number Boundary Layer at Flight Conditions

The interaction of stationary streaks undergoing nonmodal growth with modally unstable instability waves in a hypersonic boundary-layer flow is studied using numerical computations. The geometry and flow conditions are selected to match a relevant trajectory location from the ascent phase of the HIFiRE-1 flight experiment; namely, a 7 degree half-angle, circular cone with 2.5 mm nose radius, freestream Mach number equal to 5.30, unit Reynolds number equal to 13.42 m -1 , and wall-to-adiabatic temperature ratio of approximately 0.35 over most of the vehicle. This paper investigates the nonlinear evolution of initially linear optimal disturbances that evolve into finite-amplitude streaks, followed by an analysis of the modal instability characteristics of the perturbed, streaky boundary-layer flow. The investigation is performed with stationary direct numerical simulations (DNS) and plane-marching parabolized stability equations (PSE), in conjunction with partial-differential-equation-based planar eigenvalue analysis. The overall effect of streaks is to reduce the peak amplification factors of instability waves, indicating a possible downstream shift in the onset of laminar-turbulent transition. The present study confirms previous findings that the mean flow distortion of the nonlinear streak perturbation reduces the amplification rates of the Mack-mode instability. More importantly, however, the present results demonstrate that the spanwise varying component of the streak can produce a larger effect on the Mack-mode amplification. The study with selected azimuthal wavenumbers for the stationary streaks reveals that a wavenumber of approximately 1.4 times larger than the optimal wavenumber is more effective in stabilizing the planar Mack-mode instabilities. In the absence of unstable first-mode waves for the present cold-wall condition, transition onset is expected to be delayed until the peak streak amplitude increases to nearly 35 percent of the freestream velocity, when intrinsic instabilities of the boundary-layer streaks begin to dominate the transition process. For streak amplitudes below that limit a significant net stabilization is achieved, yielding a potential transition delay that can exceed 100 percent of the length of the laminar region in the uncontrolled case.

Boundary layer stability↗

NASA NTRSDate not supplied

Zemax Modeling of a Self-aligned Focusing Schlieren System with Simulated Compressible Flow

A self-aligned focusing schlieren system is simulated in the ANSYS Zemax OpticStudio software environment. The model uses a custom scatter-defining dynamic link library on the retroreflective background in order to scatter rays in the direction of the incident ray rather than the specular ray. Simulated compressible flow refractive index fields are imported into Zemax by use of the grid gradient object and a workaround to the maximum point limit is introduced. The Ronchi ruling is built as a high resolution bitmap image in MATLAB and set in Zemax as a slide object allowing the bitmap image to either transmit or terminate rays at each line pair. The Zemax simulations are compared to computational numerical schlieren and experimental conventional, path-integrated schlieren results from a hypersonic stagnation point injection flow configuration and an experimental image of a shock wave from a laser-induced breakdown. In the stagnation point injection simulations, the Zemax schlieren images capture complex flow features such as the Mach disk, shear layers, contact surface, and bow shock well, but contains long, curved features (striations) between the model and the bow shock that are not present in the synthetic schlieren image from the CFD or the experimental image. The image from Zemax of the spark-induced shock wave captures the shape and curvature of the shock but appears blurry due to grid discretization in the computational data.

Nicholas A Mejia↗

DOE OSTI2020

SYSTEMATIC MULTI-INDEX VALIDATIONS OF SIDI ENGINE FLOW FIELD LES COMPUTATIONS USING CRANK ANGLE-RESOLVED PIV MEASUREMENTS

In-cylinder flow fields make significant impacts on the fuel atomization, fuel mixture formation, and combustion process in spark ignition direct injection (SIDI) engines. In recent years, model-based simulation approaches are preferred in regard to investigating the transient in-cylinder flow field characteristics. Most commonly, the simulation models are validated using single representative flow field at a typical crank angle measured by particle image velocimetry (PIV). However, it provides only limited knowledge about the flow field which is highly three-dimensional and of transient nature. In this study, crank angle-resolved Ply measurements are conducted on three distinct planes inside the cylinder to capture the transient process of flow field characteristics which vary with the crank angle. These three planes consist of one tumble plane through the spark plug tip, one tumble plane along two intake ports, and one swirl plane at 30 mm below the cylinder head. Large eddy simulation (LES) is employed for the numerical computations using the CONVERGE codes. On the basis of large datasets for both temporal and spatial domains, a multi-index systematic validation approach is conducted. Crank angle-resolved calculations of global indices and local indices are implemented using the flow fields velocity data obtained from both PIV and LES on select planes. Global indices reveal the trends in similarities of different crank angle degrees and locations, while local indices give the detail comparison results. In summary, with the systematic multi-index validation approach, the crucial crank angle degrees and locations for model verification will be detected. Furthermore, the corresponding critical flow features are analyzed. Practical guideline of flow field validation is proposed.

Crank angle-resolved PIV↗

NASA NTRS2006

Beam Heating of Samples: Modeling and Verification

Energy absorbed from the X-ray beam by the sample requires cooling by forced convection (i.e. cryostream) to minimize temperature increase and the damage caused to the sample by the X-ray heating. In this presentation we will first review the current theoretical models and recent studies in the literature, which predict the sample temperature rise for a given set of beam parameters. It should be noted that a common weakness of these previous studies is that none of them provide actual experimental confirmation. This situation is now remedied in our investigation where the problem of x-ray sample heating is taken up once more. We have theoretically investigated, and at the same time, in addition to the numerical computations, performed experiments to validate the predictions. We have modeled, analyzed and experimentally tested the temperature rise of a 1 mm diameter glass sphere (sample surrogate) exposed to an intense synchrotron X-ray beam, while it is being cooled in a uniform flow of nitrogen gas. The heat transfer, including external convection and internal heat conduction was theoretically modeled using CFD to predict the temperature variation in the sphere during cooling and while it was subjected to an undulator (ID sector 19) X-ray beam at the APS. The surface temperature of the sphere during the X-ray beam heating was measured using the infrared camera measurement technique described in a previous talk. The temperatures from the numerical predictions and experimental measurements are compared and discussed. Additional results are reported for the two different sphere sizes and for two different supporting pin orientations.

Kazmierczak, Michael↗

DOE OSTI2021

Assessing data change in scientific datasets

Summary Scientific datasets are growing rapidly and becoming critical to next‐generation scientific discoveries. The validity of scientific results relies on the quality of data used and data are often subject to change, for example, due to observation additions, quality assessments, or processing software updates. The effects of data change are not well understood and difficult to predict. Datasets are often repeatedly updated and recomputing derived data products quickly becomes time consuming and resource intensive and may in some cases not even be necessary, thus delaying scientific advance. Despite its importance, there is a lack of systematic approaches for best comparing data versions to quantify the changes, and ad‐hoc or manual processes are commonly used. In this article, we propose a novel hierarchical approach for analyzing data changes, including real‐time (online) and offline analyses. We employ a variety of fast‐to‐compute numerical analyses, graphical data change representations, and more resource‐intensive recomputations of a subset of the data product. We illustrate the application of our approach using three scientific diverse use cases, namely, satellite, cosmological, and x‐ray data. The results show that a variety of data change metrics should be employed to enable a comprehensive representation and qualitative evaluation of data changes.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Vision-based localization for cooperative robot-CNC hybrid manufacturing

Wire and arc additive manufacturing (WAAM) has shown promise in recent years for producing large-scale parts with higher deposition rates than other additive processes. WAAM is often combined with subtractive machining to form a hybrid manufacturing process. This hybrid process can be realized by retrofitting computer numerical control (CNC) machines with deposition heads, adding spindles and deposition heads to robots, or developing part localization methods to transfer parts from an additive cell to a CNC machine. Here, a novel, robot-CNC hybrid configuration is introduced where a maneuverable robot is placed in front of a CNC machine to deposit material within the machine envelop. Furthermore, this method removes the need for part localization and the extensive machine modifications required for retrofitting; however, the problem of robot localization is also added. In this work, the effects of error in vision-based, contactless robot localization on machining parameters in a robot-machine hybrid process were studied. Performance was characterized on an implementation of this system using classical computer vision techniques. In addition, machining simulations were conducted to evaluate the effects of image-induced error on chip thickness, material removal rate, and machining allowance. Initial tests show that computer vision could adequately locate a robot for the hybrid WAAM process within .5 mm.

Hybrid manufacturing↗