Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical Simulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,081 records · Page 60
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Numerical simulation of parametric resonance in point absorbers using a simplified model

Abstract Parametric resonance is a non‐linear phenomenon in which a system can oscillate at a frequency different from its exciting frequency. Some wave energy converters are prone to this phenomenon, which is usually detrimental to their performance. Here, a computationally efficient way of simulating parametric resonance in point absorbers is presented. The model is based on linear potential theory, so the wave forces are evaluated at the mean position of the body. However, the first‐order variation of the body's centres of gravity and buoyancy is taken into account. This gives essentially the same result as a more rigorous approach of keeping terms in the equation of motion up to second order in the body motions. The only difference from a linear model is the presence of non‐zero off‐diagonal elements in the mass matrix. The model is benchmarked against state‐of‐the‐art non‐linear Froude–Krylov and computational fluid dynamics models for free decay, regular wave, and focused wave group cases. It is shown that the simplified model is able to simulate parametric resonance in pitch to a reasonable accuracy even though no non‐linear wave forces are included. The simulation speed on a standard computer is up to two orders of magnitude faster than real time.

16 TIDAL AND WAVE POWER↗
NASA NTRS2009

2D and 3D Numerical Simulations of Flux Cancellation

Cancellation of magnetic flux in the solar photosphere and chromosphere has been linked observationally and theoretically to a broad range of solar activity, from filament channel formation to CME initiation. Because this phenomenon is typically measured at only a single layer in the atmosphere, in the radial (line of sight) component of the magnetic field, the actual processes behind this observational signature are ambiguous. It is clear that reconnection is involved in some way, but the location of the reconnection sites and associated connectivity changes remain uncertain in most cases. We are using numerical modeling to demystify flux cancellation, beginning with the simplest possible configuration: a subphotospheric Lundquist flux tube surrounded by a potential field, immersed in a gravitationally stratified atmosphere, spanning many orders of magnitude in plasma beta. In this system, cancellation is driven slowly by a 2-cell circulation pattern imposed in the convection zone, such that the tops of the cells are located around the beta=1 level (i.e., the photosphere) and the flows converge and form a downdraft at the polarity inversion line; note however that no flow is imposed along the neutral line. We will present the results of 2D and 3D MHD-AMR simulations of flux cancellation, in which the flux at the photosphere begins in either an unsheared or sheared state. In all cases, a low-lying flux rope is formed by reconnection at the polarity inversion line within a few thousand seconds. The flux rope remains stable and does not rise, however, in contrast to models which do not include the presence of significant mass loading.

Karpen, Judith T.↗
NASA NTRS2009

20 and 3D Numerical Simulations of Flux Cancellation

Cancellation of magnetic flux in the solar photosphere and chromosphere has been linked observationally and theoretically to a broad range of solar activity, from filament channel formation to CME initiation. Because this phenomenon is typically measured at only a single layer in the atmosphere, in the radial (line of sight) component of the magnetic field, the actual processes behind this observational signature are ambiguous. It is clear that reconnection is involved in some way, but the location of the reconnection sites and associated connectivity changes remain uncertain in most cases. We are using numerical modeling to demystify flux cancellation, beginning with the simplest possible configuration: a subphotospheric Lundquist flux tube surrounded by a potential field, immersed in a gravitationally stratified atmosphere, spanning many orders of magnitude in plasma beta. In this system, cancellation is driven slowly by a 2-cell circulation pattern imposed in the convection zone, such that the tops of the cells are located around the beta= 1 level (Le., the photosphere) and the flows converge and form a downdraft at the polarity inversion line; note however that no flow is imposed along the neutral line. We will present the results of 2D and 3D MHD-AMR simulations of flux cancellation, in which the flux at the photosphere begins in either an unsheared or sheared state. In all cases, a lOW-lying flux rope is formed by reconnection at the polarity inversion line within a few thousand seconds. The flux rope remains stable and does not rise, however, in contrast to models which do not include the presence of significant mass loading.

Karpen, Judith T.↗
DOE OSTI2023

Core-scale numerical simulation and comparison of breakdown of shale and resulting fractures using sc-CO 2 and water as injectants

Supercritical carbon dioxide (sc-CO 2 ) is an alternative to water for stimulation of low permeability systems such as shale gas and geothermal resources. Previously core-scale experimental studies have compared the behavior of CO 2 to water injection for sample breakdown. Due to differences in experimental setup and core sample preparation, inconsistent or even apparently contradictory conclusions have resulted. To reconcile this contradiction, a phase-field numerical model is applied to understand hydraulic fracturing experiments using Green River shale found in the literature. The finite element numerical model incorporates a rate-dependent phase-field fracture model developed separately to describe fracture initiation and growth. We investigate the impact of various material and fluid properties on the resulting fractures. Most importantly, we study the effect of fluid properties and boundary conditions on the breakdown pressure, including the direction of the resulting fracture plane. Model results predict that (1) sc-CO 2 injection in the laboratory may result in greater breakdown pressure than that of water under no-flow boundary conditions because lower viscosity sc-CO 2 may result in pressure build up at the core boundary that opposes fracture initiation and (2) lower viscosity sc-CO 2 also produces fast-propagating fractures that are less influenced by the bedding plane on their resulting fracture topology. Here our model offers a straightforward explanation and reconciliation of existing experimental observations, as well as a means to extrapolate to new conditions. Exploration of field-scale conditions suggests less pronounced or no elevation in breakdown pressure when sc-CO 2 is injected because the pressure build up effect at the system boundary is significantly less or absent at field length scales.

42 ENGINEERING↗
NASA NTRS1990

Numerical simulations of supersonic flow through oscillating cascade sections

A finite difference code was developed for modeling inviscid, unsteady supersonic flow by solution of the compressible Euler equations. The code uses a deforming grid technique to capture the motion of the airfoils and can model oscillating cascades with any arbitrary interblade phase angle. A flat plate cascade is analyzed, and results are compared with results from a small perturbation theory. The results show very good agreement for both the unsteady pressure distributions and the integrated force predictions. The reason for using the numerical Euler code over a small perturbation theory is the ability to model real airfoils that have thickness and camber. Sample predictions are presented for a cascade of loaded airfoils and show appreciable differences in the unsteady surface pressure distributions when compared with the flat plate results.

Huff, Dennis L.↗
NASA NTRS1990

Numerical simulations of supersonic flow through oscillating cascade sections

A finite difference code has been developed for modeling inviscid, unsteady supersonic flow by solution of the compressible Euler equations. The code uses a deforming grid technique to capture the motion of the airfoils and can model oscillating cascades with any arbitrary interblade phase angle. A flat plate cascade is analyzed, and results are compared with results from a small-perturbation theory. The results show very good agreement for both the unsteady pressure distributions and the integrated force predictions. The reason for using the numerical Euler code over a small-perturbation theory is the ability to model 'real' airfoils that have thickness and camber. Sample predictions are presented for a cascade of loaded airfoils and show appreciable differences in the unsteady surface pressure distributions when compared with the flat plate results.

Huff, Dennis L.↗
NASA NTRS1998

Numerical Simulations of Steady and Unsteady Transition in Low-Pressure Turbine Blade Rows

Transition plays an important role in the prediction of losses and performance in low-pressure turbines. The transition location on a turbine blade may vary significantly because of the wakes from upstream blade rows, and intermittent flow separation can also affect the transition process in an unsteady flow environment. In the present investigation, an unsteady Navier-Stokes analysis is used to predict transition in a low-pressure turbine cascade and a low-pressure turbine stage. The numerical flow analysis is third-order spatially accurate and second-order temporally accurate, and the equations of motion are integrated using an implicit time-marching procedure. The Baldwin-Lomax and k-epsilon turbulence models, in conjunction with several algebraic transition models, have been used to predict the location of transition. Predicted results include unsteady blade loadings, time-histories of the pressure, transition locations and boundary layer quantities, as well as performance quantities and comparisons with the available experimental/design data.

Dorney, Daniel↗
NASA NTRS1980

Numerical simulation of advection fog formation on multi-disperse aerosols due to combustion-related pollutants

The effects of multi-disperse distribution of the aerosol population are presented. Single component and multi-component aerosol species on the condensation/nucleation processes which affect the reduction in visibility are described. The aerosol population with a high particle concentration provided more favorable conditions for the formation of a denser fog than the aerosol population with a greater particle size distribution when the value of the mass concentration of the aerosols was kept constant. The results were used as numerical predictions of fog formation. Two dimensional observations in horizontal and vertical coordinates, together with time-dependent measurements were needed as initial values for the following physical parameters: (1)wind profiles; (2) temperature profiles; (3) humidity profiles; (4) mass concentration of aerosol particles; (5) particle size distribution of aerosols; and (6) chemical composition of aerosols. Formation and dissipation of advection fog, thus, can be forecasted numerically by introducing initial values obtained from the observations.

Hung, R. J.↗
NASA NTRS1995

Numerical Simulation of Flow in a Whirling Annular Seal and Comparison with Experiments

The turbulent flow field in a simulated annular seal with a large clearance/radius ratio (0.015) and a whirling rotor was simulated using an advanced 3D CFD code SCISEAL. A circular whirl orbit with synchronous whirl was imposed on the rotor center. The flow field was rendered quasi-steady by making a transformation to a totaling frame. Standard k-epsilon model with wall functions was used to treat the turbulence. Experimentally measured values of flow parameters were used to specify the seal inlet and exit boundary conditions. The computed flow-field in terms of the velocity and pressure is compared with the experimental measurements inside the seal. The agreement between the numerical results and experimental data with correction is fair to good. The capability of current advanced CFD methodology to analyze this complex flow field is demonstrated. The methodology can also be extended to other whirl frequencies. Half- (or sub-) synchronous (fluid film unstable motion) and synchronous (rotor centrifugal force unbalance) whirls are the most unstable whirl modes in turbomachinery seals, and the flow code capability of simulating the flows in steady as well as whirling seals will prove to be extremely useful in the design, analyses, and performance predictions of annular as well as other types of seals.

Athavale, M. M.↗
NASA NTRS2018

Numerical Simulation of Dynamic Stall Using Near-Body Adaptive Mesh Refinement

Time-dependent Reynolds-averaged Navier-Stokes simulations have been carried out for a NACA 0012 airfoil and a flexible UH-60A rotor undergoing dynamic stall. The OVERFLOW computational fluid dynamics (CFD) code is used to explore the use of near-body adaptive mesh refinement (NB-AMR) for the first time on a flexible helicopter rotor in forward flight. Emphasis is placed on understanding the temporal and spatial convergence of the solutions, their numerical efficiency and stability, and establishing a grid-converged solution. Some of the remaining differences between CFD and flight-test measurements are discussed. Time-dependent flow visualization is used to provide an improved understanding of the physical mechanisms involved with two-dimensional and three-dimensional dynamic stall, and the NB-AMR process.

Near-Body↗
NASA NTRS1998

The Space-Time Conservation Element and Solution Element Method-A New High-Resolution and Genuinely Multidimensional Paradigm for Solving Conservation Laws: Numerical Simulation of Shock Waves and Contact Discontinuities - 2

Without resorting to special treatment for each individual test case, the 1D and 2D CE/SE shock-capturing schemes described previously (in Part I) are used to simulate flows involving phenomena such as shock waves, contact discontinuities, expansion waves and their interactions. Five 1D and six 2D problems are considered to examine the capability and robustness of these schemes. Despite their simple logical structures and low computational cost (for the 2D CE/SE shock-capturing scheme, the CPU time is about 2 micro-secs per mesh point per marching step on a Cray C90 machine), the numerical results, when compared with experimental data, exact solutions or numerical solutions by other methods, indicate that these schemes can accurately resolve shock and contact discontinuities consistently.

Wang, Xiao-Yen↗
NASA NTRS1992

Microgravity experiments and numerical simulations of rotating convection in a hemispherical layer

Hart el al. (1986) has previously described Spacelab experiments on rotating convection in a hemispherical layer with a spherically symmetric body force. A body force is induced on the dielectric fluid (silicon oil) by applying an electrostatic potential across the gap. Various configurations of temperature forcings are applied on the inner and outer spherical boundaries, which along with the rotation rate and the strength of the body force define the external parameters for each experiment. The present paper will present a review of the past experimental results and some new results from a numerical model of previous and planned experiments.

Miller, Timothy L.↗
NASA NTRS1996

Numerical simulation of jet aerodynamics using the three-dimensional Navier-Stokes code PAB3D

This report presents a unified method for subsonic and supersonic jet analysis using the three-dimensional Navier-Stokes code PAB3D. The Navier-Stokes code was used to obtain solutions for axisymmetric jets with on-design operating conditions at Mach numbers ranging from 0.6 to 3.0, supersonic jets containing weak shocks and Mach disks, and supersonic jets with nonaxisymmetric nozzle exit geometries. This report discusses computational methods, code implementation, computed results, and comparisons with available experimental data. Very good agreement is shown between the numerical solutions and available experimental data over a wide range of operating conditions. The Navier-Stokes method using the standard Jones-Launder two-equation kappa-epsilon turbulence model can accurately predict jet flow, and such predictions are made without any modification to the published constants for the turbulence model.

Pao, S. Paul↗
NASA NTRS1992

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.↗
NASA NTRS1999

User's Manual for LINER: FORTRAN Code for the Numerical Simulation of Plane Wave Propagation in a Lined Two-Dimensional Channel

LINER is a system of Fortran 77 codes which performs a 2D analysis of acoustic wave propagation and noise suppression in a rectangular channel with a continuous liner at the top wall. This new implementation is designed to streamline the usage of the several codes making up LINER, resulting in a useful design tool. Major input parameters are placed in two main data files, input.inc and nurn.prm. Output data appear in the form of ASCII files as well as a choice of GNUPLOT graphs. Section 2 briefly describes the physical model. Section 3 discusses the numerical methods; Section 4 gives a detailed account of program usage, including input formats and graphical options. A sample run is also provided. Finally, Section 5 briefly describes the individual program files.

Reichert, R, S.↗
NASA NTRS1973

A field investigation and numerical simulation of coastal fog

A field investigation of the microphysical and micrometeorological features of fogs occurring near Los Angeles and Vandenberg, California was conducted. Observations of wind speed and direction, temperature, dew point, vertical wind velocity, dew deposition, drop-size distribution, liquid water content, and haze and cloud nucleus concentration were obtained. These observations were initiated in late evening prior to fog formation and continued until the time of dissipation in both advection and radiation fogs. Data were also acquired in one valley fog and several dense haze situations. The behavior of these parameters prior to and during fog are discussed in detail. A two-dimensional numerical model was developed to investigate the formation and dissipation of advection fogs under the influence of horizontal variations in surface temperature. The model predicts the evolution of potential temperature, water vapor content, and liquid water content in a vertical plane as determined by vertical turbulent transfer and horizontal advection. Results are discussed from preliminary numerical experiments on the formation of warm-air advection fog and dissipation by natural and artificial heating from the surface.

Mack, E. J.↗
DOE OSTI2018

Numerical simulations of sand migration during gas production in hydrate-bearing sands interbedded with thin mud layers at site NGHP-02-16

The geomechanical behavior of gas hydrate-bearing sands during gas production is complex, involving changes in hydrate-dependent mechanical characteristics and also redistribution of stresses carried by the hydrate. When an interbedded mud layer is present, the complex response of the hydrate-bearing sediments is more pronounced because hydrate dissociation tends to occur preferentially in the sands adjacent to the hydrate-free mud layer, which can act as a heat source. This leads to non-uniform sediment deformation and potentially to shearing deformation along the sand/mud layers. In addition, since the hydrate-bearing sand permeability becomes far greater than the mud permeability with hydrate dissociation, a greater hydraulic gradient may exist along the interface between the sand/mud layers. These two behaviors may contribute to sand migration and thus a thorough study on the effect of a thin interbedded mud layer on the extent of sand migration is necessary. This paper presents numerical sand migration modeling of sand/mud interbedded hydrate-bearing sediments found at Site NGHP-02-16, one of 25 sites drilled as a part of National Gas Hydrate Program Expedition 02 (NGHP-02) in the offshore of India, subjected to depressurization. For a comparison and better understanding, simulation of homogeneous hydrate-bearing sediments is also conducted. It is found that greater sand migration occurs when the interbedded mud layer is present because of greater shearing deformation and greater pressure gradients at sand-mud interface.

58 GEOSCIENCES↗
DOE OSTI2023

A transfer learning approach for acoustic emission zonal localization on steel plate-like structure using numerical simulation and unsupervised domain adaptation

The detection and localization of damage in metallic structures using acoustic emission (AE) monitoring and artificial intelligence technology such as deep learning has been widely studied. However, a current challenge of this approach is the difficulty of obtaining sufficient labeled historical AE signals for the training process of deep learning models. This problem can be approached through the implementation of transfer learning. The innovation of this paper lies in the development of a transfer learning approach for AE source localization on a stainless-steel structure when no historical labeled AE signals are available for training. A finite element model is developed to generate numerical AE signals for the training. Unsupervised domain adaptation (UDA) technology is utilized to reduce the distribution difference between the numerical and the realistic AE signals and to derive the localization results of the unlabeled realistic AE signals. Finally, the results suggest that the proposed approach is capable of localizing AE signals with high accuracy in the absence of labeled training data.

42 ENGINEERING↗