Decay of weak turbulence
Weak turbulence fields generated by single and multiple stage grids covering Reynolds numbers between 7 and 70 showing decay of energy spectra
Huang, T. T.↗
SEARCH · Search NASA
Results for “weak grid”
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.
Multigrid solution of the Navier-Stokes equations on highly stretched grids with defect correction
Relaxation-based multigrid solvers for the steady incompressible Navier-Stokes equations are examined to determine their computational speed and robustness. Four relaxation methods with a common discretization have been used as smoothers in a single tailored multigrid procedure. The equations are discretized on a staggered grid with first order upwind used for convection in the relaxation process on all grids and defect correction to second order central on the fine grid introduced once per multigrid cycle. A fixed W(1,1) cycle with full weighting of residuals is used in the FAS multigrid process. The resulting solvers have been applied to three 2D flow problems, over a range of Reynolds numbers, on both uniform and highly stretched grids. In all cases the L(sub 2) norm of the velocity changes is reduced to 10(exp -6) in a few 10's of fine grid sweeps. The results from this study are used to draw conclusions on the strengths and weaknesses of the individual relaxation schemes as well as those of the overall multigrid procedure when used as a solver on highly stretched grids.
Sockol, Peter M.↗
Efficacy of the Cell Perturbation Method in Large-Eddy Simulations of Boundary Layer Flow over Complex Terrain
A challenge to simulating turbulent flow in multiscale atmospheric applications is the efficient generation of resolved turbulence motions over an area of interest. One approach is to apply small perturbations to flow variables near the inflow planes of turbulence-resolving simulation domains nested within larger mesoscale domains. While this approach has been examined in numerous idealized and simple terrain cases, its efficacy in complex terrain environments has not yet been fully explored. Here, we examine the benefits of the stochastic cell perturbation method (CPM) over real complex terrain using data from the 2017 Perdigão field campaign, conducted in an approximately 2-km wide valley situated between two nearly parallel ridges. Following a typical configuration for multiscale simulation using nested domains within the Weather Research and Forecasting (WRF) model to downscale from the mesoscale to a large-eddy simulation (LES), we apply the CPM on a domain with horizontal grid spacing of 150 m. At this resolution, spurious coherent structures are often observed under unstable atmospheric conditions with moderate mean wind speeds. Results from such an intermediate resolution grid are often nested down for finer, more detailed LES, where these spurious structures adversely affect the development of turbulence on the subsequent finer grid nest. We therefore examine the impacts of the CPM on the representation of turbulence within the nested LES domain under moderate mean flow conditions in three different stability regimes: weakly convective, strongly convective, and weakly stable. In addition, two different resolutions of the underlying terrain are used to explore the role of the complex topography itself in generating turbulent structures. We demonstrate that the CPM improves the representation of turbulence within the LES domain, relative to the use of high-resolution complex terrain alone. During the convective conditions, the CPM improves the rate at which smaller-scales of turbulence form, while also accelerating the attenuation of the spurious numerically generated roll structures near the inflow boundary. During stable conditions, the coarse mesh spacing of the intermediate LES domain used herein was insufficient to maintain resolved turbulence using CPM as the flow develops downstream, highlighting the need for yet higher resolution under even weakly stable conditions, and the importance of accurate representation of flow on intermediate LES grids.
17 WIND ENERGY↗
A coupled multipoint stress–multipoint flux mixed finite element method for the Biot system of poroelasticity
In this work, we present a mixed finite element method for a five-field formulation of the Biot system of poroelasticity that reduces to a cell-centered pressure–displacement system on simplicial and quadrilateral grids. A mixed stress–displacement–rotation formulation for elasticity with weak stress symmetry is coupled with a mixed velocity–pressure Darcy formulation. The spatial discretization is based on combining the multipoint stress mixed finite element (MSMFE) method for elasticity and the multipoint flux mixed finite element (MFMFE) method for Darcy flow. It uses the lowest order Brezzi–Douglas–Marini mixed finite element spaces for the poroelastic stress and Darcy velocity, piecewise constant displacement and pressure, and continuous piecewise linear or bilinear rotation. A vertex quadrature rule is applied to the velocity, stress, and stress–rotation bilinear forms, which block-diagonalizes the corresponding matrices and allows for local velocity, stress, and rotation elimination. This leads to a cell-centered positive-definite system for pressure and displacement at each time step. We perform error analysis for the semidiscrete and fully discrete formulations, establishing first order convergence for all variables in their natural norms. The numerical tests confirm the theoretical convergence rates and illustrate the locking-free property of the method.
42 ENGINEERING↗
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.↗
Development of Three-Dimensional DRAGON Grid Technology
For a typical three dimensional flow in a practical engineering device, the time spent in grid generation can take 70 percent of the total analysis effort, resulting in a serious bottleneck in the design/analysis cycle. The present research attempts to develop a procedure that can considerably reduce the grid generation effort. The DRAGON grid, as a hybrid grid, is created by means of a Direct Replacement of Arbitrary Grid Overlapping by Nonstructured grid. The DRAGON grid scheme is an adaptation to the Chimera thinking. The Chimera grid is a composite structured grid, composing a set of overlapped structured grids, which are independently generated and body-fitted. The grid is of high quality and amenable for efficient solution schemes. However, the interpolation used in the overlapped region between grids introduces error, especially when a sharp-gradient region is encountered. The DRAGON grid scheme is capable of completely eliminating the interpolation and preserving the conservation property. It maximizes the advantages of the Chimera scheme and adapts the strengths of the unstructured and while at the same time keeping its weaknesses minimal. In the present paper, we describe the progress towards extending the DRAGON grid technology into three dimensions. Essential and programming aspects of the extension, and new challenges for the three-dimensional cases, are addressed.
Zheng, Yao↗
Mass loss in main-sequence A-type stars?
We present new observations of the H-alpha line profile for five main-sequence A-type stars, where very high SNRs were achieved. A search for weak asymmetries was carried out to detect a stellar wind, but the many telluric absorption lines in this range prevent us from taking the full benefit of the high SNRs. This situation was improved by modeling the telluric absorptions to remove them from the observed spectra, but the H-alpha profiles were nevertheless found to be quite symmetric. We calculated the profile of the H-alpha line for a grid of model atmospheres of an A-type star including a weak wind, in order to assess an upper limit on the mass loss rate. The asymmetry is found to be sensitive to the velocity law and to the turbulent velocity of the wind, but its first moment depends much less on them. The upper limits deduced on the mass loss rate are between 1 and 2 x 10 exp -10 solar mass/yr. An improvement up to a factor 10 could be anticipated if new observations could be secured from a much drier location, achieving similar SNRs.
Lanz, T.↗
A nested-grid mesoscale numerical weather prediction model modified for Space Shuttle operational requirements
A nested-grid mesoscale atmospheric simulation system (MASS) is tested over Florida for the case of intense seabreeze-induced convection. The goal of this modeling system is to provide real-time aviation weather support which is designed to fit local terminal operations such as those supporting NASA's STS. Results from a 58 km and a 14.5 km nested-grid simulation show that this version of the MASS is capable of simulating many of the basic characteristics of convective complexes during periods of relatively weak synoptic scale flow regimes. However, it is noted that extensive development work is required with nested-grid cumulus and planetary boundary layer parameterization schemes before many of the meso-beta scale features such as thunderstorm downdraft-produced bubble high pressure centers can be accurately simulated. After these schemes are properly tuned, MASS can be utilized to initialize microscale modeling systems.
Kaplan, M. L.↗
Accelerating nuclear-integrated data center pursuits in the USA: SWOT analysis, power-thermal management strategies and demonstration plan
Here, this study explores the increasing interest in leveraging nuclear power to meet the escalating energy demands of data centers in the United States (U.S.) by focusing on key factors that contribute to accelerated deployment. The study highlights the importance of N+1/N+2 power supplies (where N is the required number of units), outlines research and innovations in nuclear-integrated data center thermal management and demonstration plan. It also provides updates about status and costing of various reactor system designs. A summarized strengths, weaknesses, opportunities, and threats (SWOT) analysis shows the potential options for grid connectivity, reactors, and site selection. Suitable site discussions consider land and water availability, grid access, and optical fiber connectivity, and the study presents graded prospects for Department of Energy (DOE) sites with a specific example. Community engagement and partnerships are emphasized, particularly the roles of local government, federal agencies, utilities, and data center industry partners, which are crucial for accelerating deployment, business outreach, and approvals. The study provides actionable insights for stakeholders to accelerate the deployment of nuclear-powered data centers.
21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗
Enhancements to Program LAURA for computation of three-dimensional hypersonic flow
Changes to Program Laura (Langley Aerothermodynamic Upwind Relaxation Algorithm) are presented which enhance both stability and accuracy of the algorithm. A discussion of iteration/sweeping strategies and their relation to computer architectures is included to best exploit the capabilities of serial, vector, and parallel processor machines. Test cases for Mach 10 perfect gas flow and Mach 32 real gas flow in chemical nonequilibrium over a blunt, raked elliptic cone using the thin-layer Navier-Stokes equations are presented in order to demonstrate the current improved capabilities. Algorithm changes include the use of volume averaging, application of a symmetric total variation diminishing (TVD) scheme, and stronger interaction between the grid/shock alignment routine and the relaxation algorithm. Good comparisons with heat transfer and pitching moment data at three different angles of attack for the Mach 10 tests serve to further validate the present algorithm. Parameters are defined which control the coupling of the specie continuity equations with the solution of the mixture conservation equations. A discussion of the consequences involved in the choice of strong versus weak coupling is presented, and a sample nonequilibrium calculation on a fine grid over a full scale model of the Aeroassist Flight Experiment (AFE) demonstrates current capabilities.
Gnoffo, Peter A.↗
Effects of G-Jitter on Directional Solidification of a Binary Alloy
A study of directional solidification of a weak binary alloy (specifically, Bi - 1 at% Sn) based on the fixed grid single domain approach is being undertaken. The enthalpy method is used to solve for the temperature field over the computational domain including both the solid and liquid phases; latent heat evolution is treated with the aid of an effective specific heat coefficient. A source term accounting for the release of solute into the liquid during solidification has been incorporated into the solute transport equation. The vorticity-stream function formulation is used to describe thermo-solutal convection in the liquid region. In this paper we present a numerical simulation of g-jitter. A background gravity of 1 microgram has been assumed, and new results for the effects of periodic disturbances over a range of amplitudes and frequencies on solute field and segregation have been presented.
Benjapiyaporn, C.↗
Effects of Space Environment on Flow and Concentration During Directional Solidification
A study of directional solidification of a weak binary alloy (specifically, Bi - 1 at% Sn) based on the fixed grid single domain approach is being undertaken. The enthalpy method is used to solve for the temperature field over the computational domain including both the solid and liquid phases; latent heat evolution is treated with the aid of an effective specific heat coefficient. A source term accounting for the release of solute into the liquid during solidification has been incorporated into the solute transport equation. The vorticity-stream function formulation is used to describe thermosolutal convection in the liquid region. In this paper we numerically investigate the effects of g-jitter on directional solidification. A background gravity of 1 micro-g has been assumed, and new results for the effects of periodic disturbances over a range of amplitudes and frequencies on solute field and segregation have been presented.
Benjapiyaporn, C.↗
Effect of g-jitter on Directional Solidification of a Binary Alloy
A study of directional solidification of a weak binary alloy (specifically, Bi - 1 at% Sn) based on the fixed grid single domain approach is being undertaken. The enthalpy method is used to solve for the temperature field over the computational domain including both the solid and liquid phases; latent heat evolution is treated with the aid of an effective specific heat coefficient. A source term accounting for the release of solute into the liquid during solidification has been incorporated into the solute transport equation. The vorticity-stream function formulation is used to describe thermosolutal convection in the liquid region. In this paper we present a numerical simulation of g-jitter: the small, rapid fluctuations in gravitational acceleration which may be experienced in an orbiting space vehicle. A background gravity of 1 micro-g has been assumed, and new results for the effects of orientation angle of the periodic disturbances over a range of amplitudes and frequencies on solute field and segregation have been presented.
Santiviriyapanich, P.↗
Assessment of Strengths & Weaknesses for Major IC Packages
Explore the source record for details and available documents.
BGA CSP Ball Grid Array Chip Scale Package Reliabi↗
Weak influence of anthropogenic emissions on aerosol, cloud, and rain in the wet season of the Amazon rainforest
Anthropogenic emissions have been shown to affect new particle formation, aerosol concentrations, and clouds. Such effects vary with region, environmental conditions, and cloud types. In the wet season of Amazonia, anthropogenic emissions emitted from Manaus, Brazil, can significantly increase the cloud condensation nuclei (CCN) concentrations compared to the background of mainly natural aerosols. However, the regional response of cloud and rain to anthropogenic emissions in Amazonia remains very uncertain. Here, we aim to quantify how aerosol concentration, cloud, and rain respond to changes in anthropogenic emissions through parameterized new particle formation and primary aerosol emission in the Manaus region and to understand the underlying mechanisms. We ran the atmosphere-only configuration of the HadGEM3 climate model with a nested regional domain that covers most of the rainforest region (720 km by 1200 km with 3 km resolution) under scaled regional emissions. The 7 d simulations show that, in the areas that are affected by anthropogenic emissions, when aerosol and precursor gas emissions are doubled from the baseline emission inventories, aerosol number concentrations increase by 13 %. The nucleation rate that involves sulfuric acid and biogenic compounds generally increases with pollution levels. However, nucleation is suppressed very close to the pollution source, resulting in lower nucleation and soluble Aitken mode aerosol number concentrations. We also found that doubling the anthropogenic emission can increase the cloud droplet number concentrations ( N d ) by 9 %, but cloud water and rain mass mixing ratios do not change significantly. Even very strong reductions in aerosol number concentrations by a factor of 4, which is an unrealistic condition, cause only a 4 % increase in rain over the domain. If we assume our simulation has a fine enough grid resolution and an accurate representation of the relevant atmospheric processes, the simulated weak and non-linear response of cloud and rain properties to linearly scaled anthropogenic emissions suggests that the interactions among aerosol, cloud, and precipitation in the Amazonian convective environment are buffered by microphysical processes. It also implies that the convective environment is resilient to the changes in Nd that occur in response to localized anthropogenic aerosol perturbations.
Wang, Xuemei [University of Leeds (United Kingdom)↗
Mechanical Design of Carbon Ion Optics
Carbon Ion Optics are expected to provide much longer thruster life due to their resistance to sputter erosion. There are a number of different forms of carbon that have been used for fabricating ion thruster optics. The mechanical behavior of carbon is much different than that of most metals, and poses unique design challenges. In order to minimize mission risk, the behavior of carbon must be well understood, and components designed within material limitations. Thermal expansion of the thruster structure must be compatible with thermal expansion of the carbon ion optics. Specially designed interfaces may be needed so that grid gap and aperture alignment are not adversely affected by dissimilar material properties within the thruster. The assembled thruster must be robust and tolerant of launch vibration. The following paper lists some of the characteristics of various carbon materials. Several past ion optics designs are discussed, identifying strengths and weaknesses. Electrostatics and material science are not emphasized so much as the mechanical behavior and integration of grid electrodes into an ion thruster.
Haag, Thomas↗
Seismic waves in medium with poroelastic/elastic interfaces: a two-dimensional P-SV finite-difference modelling
SUMMARY We present a new methodology of the finite-difference (FD) modelling of seismic wave propagation in a strongly heterogeneous medium composed of poroelastic (P) and (strictly) elastic (E) parts. The medium can include P/P, P/E and E/E material interfaces of arbitrary shapes. The poroelastic part can be with (i) zero resistive friction, (ii) non-zero constant resistive friction or (iii) JKD model of the frequency-dependent permeability and resistive friction. Our FD scheme is capable of subcell resolution: a material interface can have an arbitrary position in the spatial grid. The scheme keeps computational efficiency of the scheme for a smoothly and weakly heterogeneous medium (medium without material interfaces). Numerical tests against independent analytical, semi-analytical and spectral-element methods prove the efficiency and accuracy of our FD modelling. In numerical examples, we indicate effect of the P/E interfaces for the poroelastic medium with a constant resistive friction and medium with the JKD model of the frequency-dependent permeability and resistive friction. We address the 2-D P-SV problem. The approach can be readily extended to the 3-D problem.
58 GEOSCIENCES↗
An analysis of finite-difference and finite-volume formulations of conservation laws
Finite-difference and finite-volume formulations are analyzed in order to clear up the confusion concerning their application to the numerical solution of conservation laws. A new coordinate-free formulation of systems of conservation laws is developed, which clearly distinguishes the role of physical vectors from that of algebraic vectors which characterize the system. The analysis considers general types of equations--potential, Euler, and Navier-Stokes. Three-dimensional unsteady flows with time-varying grids are described using a single, consistent nomeclature for both formulations. Grid motion due to a non-inertial reference frame as well as flow adaptation is covered. In comparing the two formulations, it is found useful to distinguish between differences in numerical methods and differences in grid definition. The former plays a role for non-Cartesian grids, and results in only cosmetic differences in the manner in which geometric terms are handled. The differences in grid definition for the two formulations is found to be more important, since it affects the manner in which boundary conditions, zonal procedures, and grid singularities are handled at computational boundaries. The proper interpretation of strong and weak conservation-law forms for quasi-one-dimensional and axisymmetric flows is brought out.
Vinokur, Marcel↗