Hybrid Reynolds-Averaged/Large Eddy Simulation of the Flow in a Model SCRamjet Cavity Flameholder

NASA Technical Reports Server (NTRS)

Baurle, R. A.

2016-01-01

Steady-state and scale-resolving simulations have been performed for flow in and around a model scramjet combustor flameholder. Experimental data available for this configuration include velocity statistics obtained from particle image velocimetry. Several turbulence models were used for the steady-state Reynolds-averaged simulations which included both linear and non-linear eddy viscosity models. The scale-resolving simulations used a hybrid Reynolds-averaged/large eddy simulation strategy that is designed to be a large eddy simulation everywhere except in the inner portion (log layer and below) of the boundary layer. Hence, this formulation can be regarded as a wall-modeled large eddy simulation. This e ort was undertaken to not only assess the performance of the hybrid Reynolds-averaged / large eddy simulation modeling approach in a flowfield of interest to the scramjet research community, but to also begin to understand how this capability can best be used to augment standard Reynolds-averaged simulations. The numerical errors were quantified for the steady-state simulations, and at least qualitatively assessed for the scale-resolving simulations prior to making any claims of predictive accuracy relative to the measurements. The steady-state Reynolds-averaged results displayed a high degree of variability when comparing the flameholder fuel distributions obtained from each turbulence model. This prompted the consideration of applying the higher-fidelity scale-resolving simulations as a surrogate "truth" model to calibrate the Reynolds-averaged closures in a non-reacting setting prior to their use for the combusting simulations. In general, the Reynolds-averaged velocity profile predictions at the lowest fueling level matched the particle imaging measurements almost as well as was observed for the non-reacting condition. However, the velocity field predictions proved to be more sensitive to the flameholder fueling rate than was indicated in the measurements.

Optimization-Based Calibration of FAST.Farm Parameters Against SOWFA: Preprint

DOE Office of Scientific and Technical Information (OSTI.GOV)

Moreira, Paula D; Annoni, Jennifer; Jonkman, Jason

2018-01-04

FAST.Farm is a medium-delity wind farm modeling tool that can be used to assess power and loads contributions of wind turbines in a wind farm. The objective of this paper is to undertake a calibration procedure to set the user parameters of FAST.Farm to accurately represent results from large-eddy simulations. The results provide an in- depth analysis of the comparison of FAST.Farm and large-eddy simulations before and after calibration. The comparison of FAST.Farm and large-eddy simulation results are presented with respect to streamwise and radial velocity components as well as wake-meandering statistics (mean and standard deviation) in the lateral andmore » vertical directions under different atmospheric and turbine operating conditions.« less

Detached Eddy Simulations of Hypersonic Transition

NASA Technical Reports Server (NTRS)

Yoon, S.; Barnhardt, M.; Candler, G.

2010-01-01

This slide presentation reviews the use of Detached Eddy Simulation (DES) of hypersonic transistion. The objective of the study was to investigate the feasibility of using CFD in general, DES in particular, for prediction of roughness-induced boundary layer transition to turbulence and the resulting increase in heat transfer.

The formation processes of phytoplankton growth and decline in mesoscale eddies in the western North Pacific Ocean

NASA Astrophysics Data System (ADS)

Chang, Yu-Lin; Miyazawa, Yasumasa; Oey, Lie-Yauw; Kodaira, Tsubasa; Huang, Shihming

2017-05-01

In this study, we investigate the processes of phytoplankton growth and decline in mesoscale eddies in the western North Pacific Ocean based on the in situ chlorophyll data obtained from 52 cruises conducted by the Japan Meteorological Agency together with idealized numerical simulations. Both the observation and model results suggest that chlorophyll/phytoplankton concentrations are higher in cold than in warm eddies in near-surface water (z > -70 m). In the idealized simulation, the isopycnal movements associated with upwelling/downwelling transport phytoplankton and nutrients to different vertical depths during eddy formation (stage A). Phytoplankton and nutrients in cold eddies is transported toward shallower waters while those in warm eddies move toward deeper waters. In the period after the eddy has formed (stage B), sunlight and initially upwelled nutrients together promote the growth of phytoplankton in cold eddies. Phytoplankton in warm eddies decays due to insufficient sunlight in deeper waters. In stage B, upwelling and downwelling coexist in both warm and cold eddies, contributing nearly equally to vertical displacement. The upwelling/downwelling-induced nitrate flux accounts for a small percentage (˜3%) of the total nitrate flux in stage B. The vertical velocity caused by propagating eddies, therefore, is not the primary factor causing differences in phytoplankton concentrations between stage-B warm and cold eddies.

Role of varying interface conditions on the eddy current response from cracks in multilayer structures

NASA Astrophysics Data System (ADS)

Cherry, Aaron; Knopp, Jeremy; Aldrin, John C.; Sabbagh, Harold A.; Boehnlein, Thomas; Mooers, Ryan

2013-01-01

There is a need to improve the understanding of the role of interface conditions on eddy current inspections for cracks in multilayer aircraft structures. This paper presents initial experimental and simulated results studying the influence of gaps and contact conditions between two plates with a notch in the second layer. Simulations show an amplification of the eddy current signal for a subsurface notch adjacent to an air gap as opposed to a submerged notch in a solid plate.

Evaluation of Subgrid-Scale Models for Large Eddy Simulation of Compressible Flows

NASA Technical Reports Server (NTRS)

Blaisdell, Gregory A.

1996-01-01

The objective of this project was to evaluate and develop subgrid-scale (SGS) turbulence models for large eddy simulations (LES) of compressible flows. During the first phase of the project results from LES using the dynamic SGS model were compared to those of direct numerical simulations (DNS) of compressible homogeneous turbulence. The second phase of the project involved implementing the dynamic SGS model in a NASA code for simulating supersonic flow over a flat-plate. The model has been successfully coded and a series of simulations has been completed. One of the major findings of the work is that numerical errors associated with the finite differencing scheme used in the code can overwhelm the SGS model and adversely affect the LES results. Attached to this overview are three submitted papers: 'Evaluation of the Dynamic Model for Simulations of Compressible Decaying Isotropic Turbulence'; 'The effect of the formulation of nonlinear terms on aliasing errors in spectral methods'; and 'Large-Eddy Simulation of a Spatially Evolving Compressible Boundary Layer Flow'.

The Role of Ocean Eddies in the Southern Ocean Response to Observed Greenhouse Gas Forcing

NASA Astrophysics Data System (ADS)

Bilgen, S. I.; Kirtman, B. P.

2017-12-01

The Southern Ocean (SO) is crucial to understanding the possible future response to a changing climate. This is a principal region where energy is conveyed to the ocean by the westerly winds and it is here that mesoscale ocean eddies field dominate meridional heat and momentum transport. Compared to the Arctic, the Antarctic and the surrounding SO have a "delayed warming" anthropogenic greenhouse gas (GHG) response. Understanding the role of the ocean dynamics in modulating the mesoscale atmosphere-ocean interactions in the SO in a fully coupled regime is crucial to efforts aimed at predicting the consequences of the warming and variability to the climate system. The response of model run at multiple resolutions (eddy permitting, eddy resolving) to both GHG forcing and historical forcing are examined in NCAR CCSM4 with four experiments. The first simulation, 0.5° atmosphere coupled to ocean and sea ice components with 1° resolution (LR). The second simulation uses the identical atmospheric model but coupled to 0.1° ocean and sea ice component models (HR). For the third and fourth experiments, the global ocean is simulated for LR an HR models, and a climate change scenario are produced by applying a fixed (present-day) CO2 concentration. The analysis focuses on the last 55 years of two individual 155 year simulations. We discuss results from a set of state-of-art model experiments in comparison with observational estimates and explore mechanisms by examining sea surface temperature, westerly winds, surface heat flux, ocean heat transport. In LR simulations, the patterns and mechanisms of SO changes under GHG forcing are similar to those over the historical period: warming is damped southward of the ACC and enhanced to the north, however major changes between the HR simulations are explored. We find that in recent decades the Southern Annual Mode has shown a distinct upward trend, the result of an anthropogenic global warming. Also, HR simulations show that strengthening of the SAM and associated surface wind stress have been invoked to posit enhancement in the strength of the upwelling of the MOC, and increases eddy activity of the ACC. The results also indicate that eddy-permitting models are not able to capture the eddy-driven SST response since ocean dynamics is playing crucial role in the HR simulation but not in the LR models.

Towards a Fine-Resolution Global Coupled Climate System for Prediction on Decadal/Centennial Scales

DOE Office of Scientific and Technical Information (OSTI.GOV)

McClean, Julie L.

The over-arching goal of this project was to contribute to the realization of a fully coupled fine resolution Earth System Model simulation in which a weather-scale atmosphere is coupled to an ocean in which mesoscale eddies are largely resolved. Both a prototype fine-resolution fully coupled ESM simulation and a first-ever multi-decadal forced fine-resolution global coupled ocean/ice simulation were configured, tested, run, and analyzed as part of this grant. Science questions focused on the gains from the use of high horizontal resolution, particularly in the ocean and sea-ice, with respect to climatically important processes. Both these fine resolution coupled ocean/sea icemore » and fully-coupled simulations and precedent stand-alone eddy-resolving ocean and eddy-permitting coupled ocean/ice simulations were used to explore the high resolution regime. Overall, these studies showed that the presence of mesoscale eddies significantly impacted mixing processes and the global meridional overturning circulation in the ocean simulations. Fourteen refereed publications and a Ph.D. dissertation resulted from this grant.« less

Francis-99 turbine numerical flow simulation of steady state operation using RANS and RANS/LES turbulence model

NASA Astrophysics Data System (ADS)

Minakov, A.; Platonov, D.; Sentyabov, A.; Gavrilov, A.

2017-01-01

We performed numerical simulation of flow in a laboratory model of a Francis hydroturbine at three regimes, using two eddy-viscosity- (EVM) and a Reynolds stress (RSM) RANS models (realizable k-ɛ, k-ω SST, LRR) and detached-eddy-simulations (DES), as well as large-eddy simulations (LES). Comparison of calculation results with the experimental data was carried out. Unlike the linear EVMs, the RSM, DES, and LES reproduced well the mean velocity components, and pressure pulsations in the diffusor draft tube. Despite relatively coarse meshes and insufficient resolution of the near-wall region, LES, DES also reproduced well the intrinsic flow unsteadiness and the dominant flow structures and the associated pressure pulsations in the draft tube.

Eddy current simulation in thick cylinders of finite length induced by coils of arbitrary geometry.

PubMed

Sanchez Lopez, Hector; Poole, Michael; Crozier, Stuart

2010-12-01

Eddy currents are inevitably induced when time-varying magnetic field gradients interact with the metallic structures of a magnetic resonance imaging (MRI) scanner. The secondary magnetic field produced by this induced current degrades the spatial and temporal performance of the primary field generated by the gradient coils. Although this undesired effect can be minimized by using actively and/or passively shielded gradient coils and current pre-emphasis techniques, a residual eddy current still remains in the MRI scanner structure. Accurate simulation of these eddy currents is important in the successful design of gradient coils and magnet cryostat vessels. Efficient methods for simulating eddy currents are currently restricted to cylindrical-symmetry. The approach presented in this paper divides thick conducting cylinders into thin layers (thinner than the skin depth) and expresses the current density on each as a Fourier series. The coupling between each mode of the Fourier series with every other is modeled with an inductive network method. In this way, the eddy currents induced in realistic cryostat surfaces by coils of arbitrary geometry can be simulated. The new method was validated by simulating a canonical problem and comparing the results against a commercially available software package. An accurate skin depth of 2.76 mm was calculated in 6 min with the new method. The currents induced by an actively shielded x-gradient coil were simulated assuming a finite length cylindrical cryostat consisting of three different conducting materials. Details of the temporal-spatial induced current diffusion process were simulated through all cryostat layers, which could not be efficiently simulated with any other method. With this data, all quantities that depend on the current density, such as the secondary magnetic field, are simply evaluated. Copyright © 2010 Elsevier Inc. All rights reserved.

A numerical study of the acoustic radiation due to eddy current-cryostat interactions.

PubMed

Wang, Yaohui; Liu, Feng; Zhou, Xiaorong; Li, Yu; Crozier, Stuart

2017-06-01

To investigate the acoustic radiation due to eddy current-cryostat interactions and perform a qualitative analysis on noise reduction methods. In order to evaluate the sound pressure level (SPL) of the eddy current induced warm bore wall vibration, a Finite Element (FE) model was created to simulate the noises from both the warm bore wall vibration and the gradient coil assembly. For the SPL reduction of the warm bore wall vibration, we first improved the active shielding of the gradient coil, thus reducing the eddy current on the warm bore wall. A damping treatment was then applied to the warm bore wall to control the acoustic radiation. Initial simulations show that the SPL of the warm bore wall is higher than that of the gradient assembly with typical design shielding ratios at many frequencies. Subsequent simulation results of eddy current control and damping treatment application show that the average SPL reduction of the warm bore wall can be as high as 9.6 dB, and even higher in some frequency bands. Combining eddy current control and suggested damping scheme, the noise level in a MRI system can be effectively reduced. © 2017 American Association of Physicists in Medicine.

Large-Eddy Atmosphere-Land-Surface Modelling over Heterogeneous Surfaces: Model Development and Comparison with Measurements

NASA Astrophysics Data System (ADS)

Shao, Yaping; Liu, Shaofeng; Schween, Jan H.; Crewell, Susanne

2013-08-01

A model is developed for the large-eddy simulation (LES) of heterogeneous atmosphere and land-surface processes. This couples a LES model with a land-surface scheme. New developments are made to the land-surface scheme to ensure the adequate representation of atmosphere-land-surface transfers on the large-eddy scale. These include, (1) a multi-layer canopy scheme; (2) a method for flux estimates consistent with the large-eddy subgrid closure; and (3) an appropriate soil-layer configuration. The model is then applied to a heterogeneous region with 60-m horizontal resolution and the results are compared with ground-based and airborne measurements. The simulated sensible and latent heat fluxes are found to agree well with the eddy-correlation measurements. Good agreement is also found in the modelled and observed net radiation, ground heat flux, soil temperature and moisture. Based on the model results, we study the patterns of the sensible and latent heat fluxes, how such patterns come into existence, and how large eddies propagate and destroy land-surface signals in the atmosphere. Near the surface, the flux and land-use patterns are found to be closely correlated. In the lower boundary layer, small eddies bearing land-surface signals organize and develop into larger eddies, which carry the signals to considerably higher levels. As a result, the instantaneous flux patterns appear to be unrelated to the land-use patterns, but on average, the correlation between them is significant and persistent up to about 650 m. For a given land-surface type, the scatter of the fluxes amounts to several hundred W { m }^{-2}, due to (1) large-eddy randomness; (2) rapid large-eddy and surface feedback; and (3) local advection related to surface heterogeneity.

Direct-Numerical and Large-Eddy Simulations of a Non-Equilibrium Turbulent Kolmogorov Flow

NASA Technical Reports Server (NTRS)

Woodruff, S. L.; Shebalin, J. V.; Hussaini, M. Y.

1999-01-01

A non-equilibrium form of turbulent Kolmogorov flow is set up by making an instantaneous change in the amplitude of the spatially-periodic forcing. It is found that the response of the flow to this instantaneous change becomes more dramatic as the wavenumber of the forcing is increased, and, at the same time, that the faithfulness with which the large-eddy-simulation results agree with the direct-numerical results decreases.

Improved turbulence models based on large eddy simulation of homogeneous, incompressible turbulent flows

NASA Technical Reports Server (NTRS)

Bardino, J.; Ferziger, J. H.; Reynolds, W. C.

1983-01-01

The physical bases of large eddy simulation and subgrid modeling are studied. A subgrid scale similarity model is developed that can account for system rotation. Large eddy simulations of homogeneous shear flows with system rotation were carried out. Apparently contradictory experimental results were explained. The main effect of rotation is to increase the transverse length scales in the rotation direction, and thereby decrease the rates of dissipation. Experimental results are shown to be affected by conditions at the turbulence producing grid, which make the initial states a function of the rotation rate. A two equation model is proposed that accounts for effects of rotation and shows good agreement with experimental results. In addition, a Reynolds stress model is developed that represents the turbulence structure of homogeneous shear flows very well and can account also for the effects of system rotation.

Dual optimization method of radiofrequency and quasistatic field simulations for reduction of eddy currents generated on 7T radiofrequency coil shielding.

PubMed

Zhao, Yujuan; Zhao, Tiejun; Raval, Shailesh B; Krishnamurthy, Narayanan; Zheng, Hai; Harris, Chad T; Handler, William B; Chronik, Blaine A; Ibrahim, Tamer S

2015-11-01

To optimize the design of radiofrequency (RF) shielding of transmit coils at 7T and reduce eddy currents generated on the RF shielding when imaging with rapid gradient waveforms. One set of a four-element, 2 × 2 Tic-Tac-Toe head coil structure was selected and constructed to study eddy currents on the RF coil shielding. The generated eddy currents were quantitatively studied in the time and frequency domains. The RF characteristics were studied using the finite difference time domain method. Five different kinds of RF shielding were tested on a 7T MRI scanner with phantoms and in vivo human subjects. The eddy current simulation method was verified by the measurement results. Eddy currents induced by solid/intact and simple-structured slotted RF shielding significantly distorted the gradient fields. Echo-planar images, B1+ maps, and S matrix measurements verified that the proposed slot pattern suppressed the eddy currents while maintaining the RF characteristics of the transmit coil. The presented dual-optimization method could be used to design RF shielding and reduce the gradient field-induced eddy currents while maintaining the RF characteristics of the transmit coil. © 2014 Wiley Periodicals, Inc.

Stationary Waves of the Ice Age Climate.

NASA Astrophysics Data System (ADS)

Cook, Kerry H.; Held, Isaac M.

1988-08-01

A linearized, steady state, primitive equation model is used to simulate the climatological zonal asymmetries (stationary eddies) in the wind and temperature fields of the 18 000 YBP climate during winter. We compare these results with the eddies simulated in the ice age experiments of Broccoli and Manabe, who used CLIMAP boundary conditions and reduced atmospheric CO2 in an atmospheric general circulation model (GCM) coupled with a static mixed layer ocean model. The agreement between the models is good, indicating that the linear model can be used to evaluate the relative influences of orography, diabatic heating, and transient eddy heat and momentum transports in generating stationary waves. We find that orographic forcing dominates in the ice age climate. The mechanical influence of the continental ice sheets on the atmosphere is responsible for most of the changes between the present day and ice age stationary eddies. This concept of the ice age climate is complicated by the sensitivity of the stationary eddies to the large increase in the magnitude of the zonal mean meridional temperature gradient simulated in the ice age GCM.

Model Validation for Propulsion - On the TFNS and LES Subgrid Models for a Bluff Body Stabilized Flame

NASA Technical Reports Server (NTRS)

Wey, Thomas

2017-01-01

This paper summarizes the reacting results of simulating a bluff body stabilized flame experiment of Volvo Validation Rig using a releasable edition of the National Combustion Code (NCC). The turbulence models selected to investigate the configuration are the sub-grid scaled kinetic energy coupled large eddy simulation (K-LES) and the time-filtered Navier-Stokes (TFNS) simulation. The turbulence chemistry interaction used is linear eddy mixing (LEM).

Large-Eddy Simulation. Guidelines for Its Application to Planetary Boundary Layer Research

DTIC Science & Technology

1984-08-01

34 engineering application of L98 was Deardorff’s simulation of turbulent channel flow, which was carried out at the National Center for Atmospheric...over the past 20 years, and yet in the perception of some observers * the applications of the resulting basic science to practical problem remain...COVERED -- Large Eddy Simulation: Guidelines for its .0 application to planetary boundary layer research Final Report Oct 83-Aug 84 S. PERFORMING ORG

Eddy Current Assessment of Engineered Components Containing Nanofibers

NASA Astrophysics Data System (ADS)

Ko, Ray T.; Hoppe, Wally; Pierce, Jenny

2009-03-01

The eddy current approach has been used to assess engineered components containing nanofibers. Five specimens with different programmed defects were fabricated. A 4-point collinear probe was used to verify the electrical resistivity of each specimen. The liftoff component of the eddy current signal was used to test two extreme cases with different nano contents. Additional eddy current measurements were also used in detecting a missing nano layer simulating a manufacturing process error. The results of this assessment suggest that eddy current liftoff measurement can be a useful tool in evaluating the electrical properties of materials containing nanofibers.

Nesting large-eddy simulations within mesoscale simulations for wind energy applications

DOE Office of Scientific and Technical Information (OSTI.GOV)

Lundquist, J K; Mirocha, J D; Chow, F K

2008-09-08

With increasing demand for more accurate atmospheric simulations for wind turbine micrositing, for operational wind power forecasting, and for more reliable turbine design, simulations of atmospheric flow with resolution of tens of meters or higher are required. These time-dependent large-eddy simulations (LES), which resolve individual atmospheric eddies on length scales smaller than turbine blades and account for complex terrain, are possible with a range of commercial and open-source software, including the Weather Research and Forecasting (WRF) model. In addition to 'local' sources of turbulence within an LES domain, changing weather conditions outside the domain can also affect flow, suggesting thatmore » a mesoscale model provide boundary conditions to the large-eddy simulations. Nesting a large-eddy simulation within a mesoscale model requires nuanced representations of turbulence. Our group has improved the Weather and Research Forecasting model's (WRF) LES capability by implementing the Nonlinear Backscatter and Anisotropy (NBA) subfilter stress model following Kosovic (1997) and an explicit filtering and reconstruction technique to compute the Resolvable Subfilter-Scale (RSFS) stresses (following Chow et al, 2005). We have also implemented an immersed boundary method (IBM) in WRF to accommodate complex terrain. These new models improve WRF's LES capabilities over complex terrain and in stable atmospheric conditions. We demonstrate approaches to nesting LES within a mesoscale simulation for farms of wind turbines in hilly regions. Results are sensitive to the nesting method, indicating that care must be taken to provide appropriate boundary conditions, and to allow adequate spin-up of turbulence in the LES domain.« less

Detached Eddy Simulation of Flap Side-Edge Flow

NASA Technical Reports Server (NTRS)

Balakrishnan, Shankar K.; Shariff, Karim R.

2016-01-01

Detached Eddy Simulation (DES) of flap side-edge flow was performed with a wing and half-span flap configuration used in previous experimental and numerical studies. The focus of the study is the unsteady flow features responsible for the production of far-field noise. The simulation was performed at a Reynolds number (based on the main wing chord) of 3.7 million. Reynolds Averaged Navier-Stokes (RANS) simulations were performed as a precursor to the DES. The results of these precursor simulations match previous experimental and RANS results closely. Although the present DES simulations have not reached statistical stationary yet, some unsteady features of the developing flap side-edge flowfield are presented. In the final paper it is expected that statistically stationary results will be presented including comparisons of surface pressure spectra with experimental data.

Processes influencing formation of low-salinity high-biomass lenses near the edge of the Ross Ice Shelf

NASA Astrophysics Data System (ADS)

Li, Yizhen; McGillicuddy, Dennis J.; Dinniman, Michael S.; Klinck, John M.

2017-02-01

Both remotely sensed and in situ observations in austral summer of early 2012 in the Ross Sea suggest the presence of cold, low-salinity, and high-biomass eddies along the edge of the Ross Ice Shelf (RIS). Satellite measurements include sea surface temperature and ocean color, and shipboard data sets include hydrographic profiles, towed instrumentation, and underway acoustic Doppler current profilers. Idealized model simulations are utilized to examine the processes responsible for ice shelf eddy formation. 3-D model simulations produce similar cold and fresh eddies, although the simulated vertical lenses are quantitatively thinner than observed. Model sensitivity tests show that both basal melting underneath the ice shelf and irregularity of the ice shelf edge facilitate generation of cold and fresh eddies. 2-D model simulations further suggest that both basal melting and downwelling-favorable winds play crucial roles in forming a thick layer of low-salinity water observed along the edge of the RIS. These properties may have been entrained into the observed eddies, whereas that entrainment process was not captured in the specific eddy formation events studied in our 3-D model-which may explain the discrepancy between the simulated and observed eddies, at least in part. Additional sensitivity experiments imply that uncertainties associated with background stratification and wind stress may also explain why the model underestimates the thickness of the low-salinity lens in the eddy interiors. Our study highlights the importance of incorporating accurate wind forcing, basal melting, and ice shelf irregularity for simulating eddy formation near the RIS edge. The processes responsible for generating the high phytoplankton biomass inside these eddies remain to be elucidated. Appendix B. Details for the basal melting and mechanical forcing by the ice shelf edge.

Eddy Fluxes and Sensitivity of the Water Cycle to Spatial Resolution in Idealized Regional Aquaplanet Model Simulations

DOE Office of Scientific and Technical Information (OSTI.GOV)

Hagos, Samson M.; Leung, Lai-Yung R.; Gustafson, William I.

2014-02-28

A multi-scale moisture budget analysis is used to identify the mechanisms responsible for the sensitivity of the water cycle to spatial resolution using idealized regional aquaplanet simulations. In the higher resolution simulations, moisture transport by eddies fluxes dry the boundary layer enhancing evaporation and precipitation. This effect of eddies, which is underestimated by the physics parameterizations in the low-resolution simulations, is found to be responsible for the sensitivity of the water cycle both directly, and through its upscale effect, on the mean circulation. Correlations among moisture transport by eddies at adjacent ranges of scales provides the potential for reducing thismore » sensitivity by representing the unresolved eddies by their marginally resolved counterparts.« less

Mapping B(1)-induced eddy current effects near metallic structures in MR images: a comparison of simulation and experiment.

PubMed

Vashaee, S; Goora, F; Britton, M M; Newling, B; Balcom, B J

2015-01-01

Magnetic resonance imaging (MRI) in the presence of metallic structures is very common in medical and non-medical fields. Metallic structures cause MRI image distortions by three mechanisms: (1) static field distortion through magnetic susceptibility mismatch, (2) eddy currents induced by switched magnetic field gradients and (3) radio frequency (RF) induced eddy currents. Single point ramped imaging with T1 enhancement (SPRITE) MRI measurements are largely immune to susceptibility and gradient induced eddy current artifacts. As a result, one can isolate the effects of metal objects on the RF field. The RF field affects both the excitation and detection of the magnetic resonance (MR) signal. This is challenging with conventional MRI methods, which cannot readily separate the three effects. RF induced MRI artifacts were investigated experimentally at 2.4 T by analyzing image distortions surrounding two geometrically identical metallic strips of aluminum and lead. The strips were immersed in agar gel doped with contrast agent and imaged employing the conical SPRITE sequence. B1 mapping with pure phase encode SPRITE was employed to measure the B1 field around the strips of metal. The strip geometry was chosen to mimic metal electrodes employed in electrochemistry studies. Simulations are employed to investigate the RF field induced eddy currents in the two metallic strips. The RF simulation results are in good agreement with experimental results. Experimental and simulation results show that the metal has a pronounced effect on the B1 distribution and B1 amplitude in the surrounding space. The electrical conductivity of the metal has a minimal effect. Copyright © 2014 Elsevier Inc. All rights reserved.

Wind turbine wakes in forest and neutral plane wall boundary layer large-eddy simulations

NASA Astrophysics Data System (ADS)

Schröttle, Josef; Piotrowski, Zbigniew; Gerz, Thomas; Englberger, Antonia; Dörnbrack, Andreas

2016-09-01

Wind turbine wake flow characteristics are studied in a strongly sheared and turbulent forest boundary layer and a neutral plane wall boundary layer flow. The reference simulations without wind turbine yield similar results as earlier large-eddy simulations by Shaw and Schumann (1992) and Porte-Agel et al. (2000). To use the fields from the homogeneous turbulent boundary layers on the fly as inflow fields for the wind turbine wake simulations, a new and efficient methodology was developed for the multiscale geophysical flow solver EULAG. With this method fully developed turbulent flow fields can be achieved upstream of the wind turbine which are independent of the wake flow. The large-eddy simulations reproduce known boundary-layer statistics as mean wind profile, momentum flux profile, and eddy dissipation rate of the plane wall and the forest boundary layer. The wake velocity deficit is more asymmetric above the forest and recovers faster downstream compared to the velocity deficit in the plane wall boundary layer. This is due to the inflection point in the mean streamwise velocity profile with corresponding turbulent coherent structures of high turbulence intensity in the strong shear flow above the forest.

Effects of Eddy Viscosity on Time Correlations in Large Eddy Simulation

NASA Technical Reports Server (NTRS)

He, Guowei; Rubinstein, R.; Wang, Lian-Ping; Bushnell, Dennis M. (Technical Monitor)

2001-01-01

Subgrid-scale (SGS) models for large. eddy simulation (LES) have generally been evaluated by their ability to predict single-time statistics of turbulent flows such as kinetic energy and Reynolds stresses. Recent application- of large eddy simulation to the evaluation of sound sources in turbulent flows, a problem in which time, correlations determine the frequency distribution of acoustic radiation, suggest that subgrid models should also be evaluated by their ability to predict time correlations in turbulent flows. This paper compares the two-point, two-time Eulerian velocity correlation evaluated from direct numerical simulation (DNS) with that evaluated from LES, using a spectral eddy viscosity, for isotropic homogeneous turbulence. It is found that the LES fields are too coherent, in the sense that their time correlations decay more slowly than the corresponding time. correlations in the DNS fields. This observation is confirmed by theoretical estimates of time correlations using the Taylor expansion technique. Tile reason for the slower decay is that the eddy viscosity does not include the random backscatter, which decorrelates fluid motion at large scales. An effective eddy viscosity associated with time correlations is formulated, to which the eddy viscosity associated with energy transfer is a leading order approximation.

Numerical Investigations of Subduction of Eighteen Degree Water in the Subtropical Northwest Atlantic Ocean

NASA Astrophysics Data System (ADS)

Zhai, P.; He, R.

2016-02-01

Mode waters are upper-ocean water masses with nearly uniform water properties over a thickness of a few hundred meters. Subduction of mode waters plays an important role in changing atmospheric and oceanic long-term variability because they store "memory" of wintertime air-sea interaction. In this study, we investigated dynamic processes associated with subduction of the Eighteen Degree Water (EDW, the principal mode water) in the subtropical Northwest Atlantic during January to June 2007. Numerical simulations of the temporal and spatial evolutions of EDW were performed using both uncoupled (ocean only) and air-sea coupled configurations and results were contrasted. We find the coupled simulation produced deeper mixed layer depth, stronger eddy kinetic energy, and larger subduction areas than their counterparts in the uncoupled ocean simulation. In both configurations, mesoscale eddies enhance the total subduction and eddy-induced subduction has the same order as the mean component. Resolving strong air-sea coupling and mesoscale eddies is therefore important for understanding EDW dynamics.

Simulated Design Strategies for SPECT Collimators to Reduce the Eddy Currents Induced by MRI Gradient Fields

NASA Astrophysics Data System (ADS)

Samoudi, Amine M.; Van Audenhaege, Karen; Vermeeren, Günter; Verhoyen, Gregory; Martens, Luc; Van Holen, Roel; Joseph, Wout

2015-10-01

Combining single photon emission computed tomography (SPECT) with magnetic resonance imaging (MRI) requires the insertion of highly conductive SPECT collimators inside the MRI scanner, resulting in an induced eddy current disturbing the combined system. We reduced the eddy currents due to the insert of a novel tungsten collimator inside transverse and longitudinal gradient coils. The collimator was produced with metal additive manufacturing, that is part of a microSPECT insert for a preclinical SPECT/MRI scanner. We characterized the induced magnetic field due to the gradient field and adapted the collimators to reduce the induced eddy currents. We modeled the x-, y-, and z-gradient coil and the different collimator designs and simulated them with FEKO, a three-dimensional method of moments / finite element methods (MoM/FEM) full-wave simulation tool. We used a time analysis approach to generate the pulsed magnetic field gradient. Simulation results show that the maximum induced field can be reduced by 50.82% in the final design bringing the maximum induced magnetic field to less than 2% of the applied gradient for all the gradient coils. The numerical model was validated with measurements and was proposed as a tool for studying the effect of a SPECT collimator within the MRI gradient coils.

Large Eddy Simulation of a Film Cooling Flow Injected from an Inclined Discrete Cylindrical Hole into a Crossflow with Zero-Pressure Gradient Turbulent Boundary Layer

NASA Technical Reports Server (NTRS)

Johnson, Perry L.; Shyam, Vikram

2012-01-01

A Large Eddy Simulation (LES) is performed of a high blowing ratio (M = 1.7) film cooling flow with density ratio of unity. Mean results are compared with experimental data to show the degree of fidelity achieved in the simulation. While the trends in the LES prediction are a noticeable improvement over Reynolds-Averaged Navier-Stokes (RANS) predictions, there is still a lack a spreading on the underside of the lifted jet. This is likely due to the inability of the LES to capture the full range of influential eddies on the underside of the jet due to their smaller structure. The unsteady structures in the turbulent coolant jet are also explored and related to turbulent mixing characteristics

Tidal Residual Eddies and their Effect on Water Exchange in Puget Sound

DOE Office of Scientific and Technical Information (OSTI.GOV)

Yang, Zhaoqing; Wang, Taiping

Tidal residual eddies are one of the important hydrodynamic features in tidally dominant estuaries and coastal bays, and they could have significant effects on water exchange in a tidal system. This paper presents a modeling study of tides and tidal residual eddies in Puget Sound, a tidally dominant fjord-like estuary in the Pacific Northwest coast, using a three-dimensional finite-volume coastal ocean model. Mechanisms of vorticity generation and asymmetric distribution patterns around an island/headland were analyzed using the dynamic vorticity transfer approach and numerical experiments. Model results of Puget Sound show that a number of large twin tidal residual eddies existmore » in the Admiralty Inlet because of the presence of major headlands in the inlet. Simulated residual vorticities near the major headlands indicate that the clockwise tidal residual eddy (negative vorticity) is generally stronger than the anticlockwise eddy (positive vorticity) because of the effect of Coriolis force. The effect of tidal residual eddies on water exchange in Puget Sound and its sub-basins were evaluated by simulations of dye transport. It was found that the strong transverse variability of residual currents in the Admiralty Inlet results in a dominant seaward transport along the eastern shore and a dominant landward transport along the western shore of the Inlet. A similar transport pattern in Hood Canal is caused by the presence of tidal residual eddies near the entrance of the canal. Model results show that tidal residual currents in Whidbey Basin are small in comparison to other sub-basins. A large clockwise residual circulation is formed around Vashon Island near entrance of South Sound, which can potentially constrain the water exchange between the Central Basin and South Sound.« less

The nonlinear Galerkin method: A multi-scale method applied to the simulation of homogeneous turbulent flows

NASA Technical Reports Server (NTRS)

Debussche, A.; Dubois, T.; Temam, R.

1993-01-01

Using results of Direct Numerical Simulation (DNS) in the case of two-dimensional homogeneous isotropic flows, the behavior of the small and large scales of Kolmogorov like flows at moderate Reynolds numbers are first analyzed in detail. Several estimates on the time variations of the small eddies and the nonlinear interaction terms were derived; those terms play the role of the Reynolds stress tensor in the case of LES. Since the time step of a numerical scheme is determined as a function of the energy-containing eddies of the flow, the variations of the small scales and of the nonlinear interaction terms over one iteration can become negligible by comparison with the accuracy of the computation. Based on this remark, a multilevel scheme which treats differently the small and the large eddies was proposed. Using mathematical developments, estimates of all the parameters involved in the algorithm, which then becomes a completely self-adaptive procedure were derived. Finally, realistic simulations of (Kolmorov like) flows over several eddy-turnover times were performed. The results are analyzed in detail and a parametric study of the nonlinear Galerkin method is performed.

Toward large eddy simulation of turbulent flow over an airfoil

NASA Technical Reports Server (NTRS)

Choi, Haecheon

1993-01-01

The flow field over an airfoil contains several distinct flow characteristics, e.g. laminar, transitional, turbulent boundary layer flow, flow separation, unstable free shear layers, and a wake. This diversity of flow regimes taxes the presently available Reynolds averaged turbulence models. Such models are generally tuned to predict a particular flow regime, and adjustments are necessary for the prediction of a different flow regime. Similar difficulties are likely to emerge when the large eddy simulation technique is applied with the widely used Smagorinsky model. This model has not been successful in correctly representing different turbulent flow fields with a single universal constant and has an incorrect near-wall behavior. Germano et al. (1991) and Ghosal, Lund & Moin have developed a new subgrid-scale model, the dynamic model, which is very promising in alleviating many of the persistent inadequacies of the Smagorinsky model: the model coefficient is computed dynamically as the calculation progresses rather than input a priori. The model has been remarkably successful in prediction of several turbulent and transitional flows. We plan to simulate turbulent flow over a '2D' airfoil using the large eddy simulation technique. Our primary objective is to assess the performance of the newly developed dynamic subgrid-scale model for computation of complex flows about aircraft components and to compare the results with those obtained using the Reynolds average approach and experiments. The present computation represents the first application of large eddy simulation to a flow of aeronautical interest and a key demonstration of the capabilities of the large eddy simulation technique.

Continuous assimilation of simulated Geosat altimetric sea level into an eddy-resolving numerical ocean model. I - Sea level differences. II - Referenced sea level differences

NASA Technical Reports Server (NTRS)

White, Warren B.; Tai, Chang-Kou; Holland, William R.

1990-01-01

The optimal interpolation method of Lorenc (1981) was used to conduct continuous assimilation of altimetric sea level differences from the simulated Geosat exact repeat mission (ERM) into a three-layer quasi-geostrophic eddy-resolving numerical ocean box model that simulates the statistics of mesoscale eddy activity in the western North Pacific. Assimilation was conducted continuously as the Geosat tracks appeared in simulated real time/space, with each track repeating every 17 days, but occurring at different times and locations within the 17-day period, as would have occurred in a realistic nowcast situation. This interpolation method was also used to conduct the assimilation of referenced altimetric sea level differences into the same model, performing the referencing of altimetric sea sevel differences by using the simulated sea level. The results of this dynamical interpolation procedure are compared with those of a statistical (i.e., optimum) interpolation procedure.

Wind Forced Variability in Eddy Formation, Eddy Shedding, and the Separation of the East Australian Current

NASA Astrophysics Data System (ADS)

Bull, Christopher Y. S.; Kiss, Andrew E.; Jourdain, Nicolas C.; England, Matthew H.; van Sebille, Erik

2017-12-01

The East Australian Current (EAC), like many other subtropical western boundary currents, is believed to be penetrating further poleward in recent decades. Previous observational and model studies have used steady state dynamics to relate changes in the westerly winds to changes in the separation behavior of the EAC. As yet, little work has been undertaken on the impact of forcing variability on the EAC and Tasman Sea circulation. Here using an eddy-permitting regional ocean model, we present a suite of simulations forced by the same time-mean fields, but with different atmospheric and remote ocean variability. These eddy-permitting results demonstrate the nonlinear response of the EAC to variable, nonstationary inhomogeneous forcing. These simulations show an EAC with high intrinsic variability and stochastic eddy shedding. We show that wind stress variability on time scales shorter than 56 days leads to increases in eddy shedding rates and southward eddy propagation, producing an increased transport and southward reach of the mean EAC extension. We adopt an energetics framework that shows the EAC extension changes to be coincident with an increase in offshore, upstream eddy variance (via increased barotropic instability) and increase in subsurface mean kinetic energy along the length of the EAC. The response of EAC separation to regional variable wind stress has important implications for both past and future climate change studies.

The prospect of using large eddy and detached eddy simulations in engineering design, and the research required to get there

PubMed Central

Larsson, Johan; Wang, Qiqi

2014-01-01

In this paper, we try to look into the future to envision how large eddy and detached eddy simulations will be used in the engineering design process about 20–30 years from now. Some key challenges specific to the engineering design process are identified, and some of the critical outstanding problems and promising research directions are discussed. PMID:25024421

Benchmarking the mesoscale variability in global ocean eddy-permitting numerical systems

NASA Astrophysics Data System (ADS)

Cipollone, Andrea; Masina, Simona; Storto, Andrea; Iovino, Doroteaciro

2017-10-01

The role of data assimilation procedures on representing ocean mesoscale variability is assessed by applying eddy statistics to a state-of-the-art global ocean reanalysis (C-GLORS), a free global ocean simulation (performed with the NEMO system) and an observation-based dataset (ARMOR3D) used as an independent benchmark. Numerical results are computed on a 1/4 ∘ horizontal grid (ORCA025) and share the same resolution with ARMOR3D dataset. This "eddy-permitting" resolution is sufficient to allow ocean eddies to form. Further to assessing the eddy statistics from three different datasets, a global three-dimensional eddy detection system is implemented in order to bypass the need of regional-dependent definition of thresholds, typical of commonly adopted eddy detection algorithms. It thus provides full three-dimensional eddy statistics segmenting vertical profiles from local rotational velocities. This criterion is crucial for discerning real eddies from transient surface noise that inevitably affects any two-dimensional algorithm. Data assimilation enhances and corrects mesoscale variability on a wide range of features that cannot be well reproduced otherwise. The free simulation fairly reproduces eddies emerging from western boundary currents and deep baroclinic instabilities, while underestimates shallower vortexes that populate the full basin. The ocean reanalysis recovers most of the missing turbulence, shown by satellite products , that is not generated by the model itself and consistently projects surface variability deep into the water column. The comparison with the statistically reconstructed vertical profiles from ARMOR3D show that ocean data assimilation is able to embed variability into the model dynamics, constraining eddies with in situ and altimetry observation and generating them consistently with local environment.

Large Eddy Simulation of Wake Vortices in the Convective Boundary Layer

NASA Technical Reports Server (NTRS)

Lin, Yuh-Lang; Han, Jongil; Zhang, Jing; Ding, Feng; Arya, S. Pal; Proctor, Fred H.

2000-01-01

The behavior of wake vortices in a convective boundary layer is investigated using a validated large eddy simulation model. Our results show that the vortices are largely deformed due to strong turbulent eddy motion while a sinusoidal Crow instability develops. Vortex rising is found to be caused by the updrafts (thermals) during daytime convective conditions and increases with increasing nondimensional turbulence intensity eta. In the downdraft region of the convective boundary layer, vortex sinking is found to be accelerated proportional to increasing eta, with faster speed than that in an ideal line vortex pair in an inviscid fluid. Wake vortices are also shown to be laterally transported over a significant distance due to large turbulent eddy motion. On the other hand, the decay rate of the, vortices in the convective boundary layer that increases with increasing eta, is larger in the updraft region than in the downdraft region because of stronger turbulence in the updraft region.

Eulerian and Lagrangian Parameterization of the Oceanic Mixed Layer using Large Eddy Simulation and MPAS-Ocean

DOE Office of Scientific and Technical Information (OSTI.GOV)

Van Roekel, Luke

We have conducted a suite of Large Eddy Simulation (LES) to form the basis of a multi-model comparison (left). The results have led to proposed model improvements. We have verified that Eulerian-Lagrangian effective diffusivity estimates of mesoscale mixing are consistent with traditional particle statistics metrics (right). LES and Lagrangian particles will be utilized to better represent the movement of water into and out of the mixed layer.

Hybrid Solution-Adaptive Unstructured Cartesian Method for Large-Eddy Simulation of Detonation in Multi-Phase Turbulent Reactive Mixtures

DTIC Science & Technology

2012-03-27

pulse- detonation engines ( PDE ), stage separation, supersonic cav- ity oscillations, hypersonic aerodynamics, detonation induced structural...ADAPTIVE UNSTRUCTURED CARTESIAN METHOD FOR LARGE-EDDY SIMULATION OF DETONATION IN MULTI-PHASE TURBULENT REACTIVE MIXTURES 5b. GRANT NUMBER FA9550...CCL Report TR-2012-03-03 Hybrid Solution-Adaptive Unstructured Cartesian Method for Large-Eddy Simulation of Detonation in Multi-Phase Turbulent

Large eddy simulation of piloted pulverised coal combustion using extended flamelet/progress variable model

NASA Astrophysics Data System (ADS)

Wen, Xu; Luo, Kun; Jin, Hanhui; Fan, Jianren

2017-09-01

An extended flamelet/progress variable (EFPV) model for simulating pulverised coal combustion (PCC) in the context of large eddy simulation (LES) is proposed, in which devolatilisation, char surface reaction and radiation are all taken into account. The pulverised coal particles are tracked in the Lagrangian framework with various sub-models and the sub-grid scale (SGS) effects of turbulent velocity and scalar fluctuations on the coal particles are modelled by the velocity-scalar joint filtered density function (VSJFDF) model. The presented model is then evaluated by LES of an experimental piloted coal jet flame and comparing the numerical results with the experimental data and the results from the eddy break up (EBU) model. Detailed quantitative comparisons are carried out. It is found that the proposed model performs much better than the EBU model on radial velocity and species concentrations predictions. Comparing against the adiabatic counterpart, we find that the predicted temperature is evidently lowered and agrees well with the experimental data if the conditional sampling method is adopted.

Large eddy simulation for atmospheric boundary layer flow over flat and complex terrains

NASA Astrophysics Data System (ADS)

Han, Yi; Stoellinger, Michael; Naughton, Jonathan

2016-09-01

In this work, we present Large Eddy Simulation (LES) results of atmospheric boundary layer (ABL) flow over complex terrain with neutral stratification using the OpenFOAM-based simulator for on/offshore wind farm applications (SOWFA). The complete work flow to investigate the LES for the ABL over real complex terrain is described including meteorological-tower data analysis, mesh generation and case set-up. New boundary conditions for the lateral and top boundaries are developed and validated to allow inflow and outflow as required in complex terrain simulations. The turbulent inflow data for the terrain simulation is generated using a precursor simulation of a flat and neutral ABL. Conditionally averaged met-tower data is used to specify the conditions for the flat precursor simulation and is also used for comparison with the simulation results of the terrain LES. A qualitative analysis of the simulation results reveals boundary layer separation and recirculation downstream of a prominent ridge that runs across the simulation domain. Comparisons of mean wind speed, standard deviation and direction between the computed results and the conditionally averaged tower data show a reasonable agreement.

Assessment of the Partially Resolved Numerical Simulation (PRNS) Approach in the National Combustion Code (NCC) for Turbulent Nonreacting and Reacting Flows

NASA Technical Reports Server (NTRS)

Shih, Tsan-Hsing; Liu, Nan-Suey

2008-01-01

This paper describes an approach which aims at bridging the gap between the traditional Reynolds-averaged Navier-Stokes (RANS) approach and the traditional large eddy simulation (LES) approach. It has the characteristics of the very large eddy simulation (VLES) and we call this approach the partially-resolved numerical simulation (PRNS). Systematic simulations using the National Combustion Code (NCC) have been carried out for fully developed turbulent pipe flows at different Reynolds numbers to evaluate the PRNS approach. Also presented are the sample results of two demonstration cases: nonreacting flow in a single injector flame tube and reacting flow in a Lean Direct Injection (LDI) hydrogen combustor.

Multilayer integral method for simulation of eddy currents in thin volumes of arbitrary geometry produced by MRI gradient coils.

PubMed

Sanchez Lopez, Hector; Freschi, Fabio; Trakic, Adnan; Smith, Elliot; Herbert, Jeremy; Fuentes, Miguel; Wilson, Stephen; Liu, Limei; Repetto, Maurizio; Crozier, Stuart

2014-05-01

This article aims to present a fast, efficient and accurate multi-layer integral method (MIM) for the evaluation of complex spatiotemporal eddy currents in nonmagnetic and thin volumes of irregular geometries induced by arbitrary arrangements of gradient coils. The volume of interest is divided into a number of layers, wherein the thickness of each layer is assumed to be smaller than the skin depth and where one of the linear dimensions is much smaller than the remaining two dimensions. The diffusion equation of the current density is solved both in time-harmonic and transient domain. The experimentally measured magnetic fields produced by the coil and the induced eddy currents as well as the corresponding time-decay constants were in close agreement with the results produced by the MIM. Relevant parameters such as power loss and force induced by the eddy currents in a split cryostat were simulated using the MIM. The proposed method is capable of accurately simulating the current diffusion process inside thin volumes, such as the magnet cryostat. The method permits the priori-calculation of optimal pre-emphasis parameters. The MIM enables unified designs of gradient coil-magnet structures for an optimal mitigation of deleterious eddy current effects. Copyright © 2013 Wiley Periodicals, Inc.

Large eddy simulation of transitional flow in an idealized stenotic blood vessel: evaluation of subgrid scale models.

PubMed

Pal, Abhro; Anupindi, Kameswararao; Delorme, Yann; Ghaisas, Niranjan; Shetty, Dinesh A; Frankel, Steven H

2014-07-01

In the present study, we performed large eddy simulation (LES) of axisymmetric, and 75% stenosed, eccentric arterial models with steady inflow conditions at a Reynolds number of 1000. The results obtained are compared with the direct numerical simulation (DNS) data (Varghese et al., 2007, "Direct Numerical Simulation of Stenotic Flows. Part 1. Steady Flow," J. Fluid Mech., 582, pp. 253-280). An inhouse code (WenoHemo) employing high-order numerical methods for spatial and temporal terms, along with a 2nd order accurate ghost point immersed boundary method (IBM) (Mark, and Vanwachem, 2008, "Derivation and Validation of a Novel Implicit Second-Order Accurate Immersed Boundary Method," J. Comput. Phys., 227(13), pp. 6660-6680) for enforcing boundary conditions on curved geometries is used for simulations. Three subgrid scale (SGS) models, namely, the classical Smagorinsky model (Smagorinsky, 1963, "General Circulation Experiments With the Primitive Equations," Mon. Weather Rev., 91(10), pp. 99-164), recently developed Vreman model (Vreman, 2004, "An Eddy-Viscosity Subgrid-Scale Model for Turbulent Shear Flow: Algebraic Theory and Applications," Phys. Fluids, 16(10), pp. 3670-3681), and the Sigma model (Nicoud et al., 2011, "Using Singular Values to Build a Subgrid-Scale Model for Large Eddy Simulations," Phys. Fluids, 23(8), 085106) are evaluated in the present study. Evaluation of SGS models suggests that the classical constant coefficient Smagorinsky model gives best agreement with the DNS data, whereas the Vreman and Sigma models predict an early transition to turbulence in the poststenotic region. Supplementary simulations are performed using Open source field operation and manipulation (OpenFOAM) ("OpenFOAM," http://www.openfoam.org/) solver and the results are inline with those obtained with WenoHemo.

Information content of in situ and remotely sensed chlorophyll-a: Learning from size-structured phytoplankton model

NASA Astrophysics Data System (ADS)

Laiolo, Leonardo; Matear, Richard; Baird, Mark E.; Soja-Woźniak, Monika; Doblin, Martina A.

2018-07-01

Chlorophyll-a measurements in the form of in situ observations and satellite ocean colour products are commonly used in data assimilation to calibrate marine biogeochemical models. Here, a two size-class phytoplankton biogeochemical model, with a 0D configuration, was used to simulate the surface chlorophyll-a dynamics (simulated surface Chl-a) for cyclonic and anticyclonic eddies off East Australia. An optical model was then used to calculate the inherent optical properties from the simulation and convert them into remote-sensing reflectance (Rrs). Subsequently, Rrs was used to produce a satellite-like estimate of the simulated surface Chl-a concentrations through the MODIS OC3M algorithm (simulated OC3M Chl-a). Identical parameter optimisation experiments were performed through the assimilation of the two separate datasets (simulated surface Chl-a and simulated OC3M Chl-a), with the purpose of investigating the contrasting information content of simulated surface Chl-a and remotely-sensed data sources. The results we present are based on the analysis of the distribution of a cost function, varying four parameters of the biogeochemical model. In our idealized experiments the simulated OC3M Chl-a product is a poor proxy for the total simulated surface Chl-a concentration. Furthermore, our result show the OC3M algorithm can underestimate the simulated chlorophyll-a concentration in offshore eddies off East Australia (Case I waters), because of the weak relationship between large-sized phytoplankton and remote-sensing reflectance. Although Case I waters are usually characteristic of oligotrophic environments, with a photosynthetic community typically represented by relatively small-sized phytoplankton, mesoscale features such as eddies can generate seasonally favourable conditions for a photosynthetic community with a greater proportion of large phytoplankton cells. Furthermore, our results show that in mesoscale features such as eddies, in situ chlorophyll-a observations and the ocean colour products can carry different information related to phytoplankton sizes. Assimilating both remote-sensing reflectance and measurements of in situ chlorophyll-a concentration reduces the uncertainty of the parameter values more than either data set alone, thus reducing the spread of acceptable solutions, giving an improved simulation of the natural environment.

Detached-Eddy Simulations of Attached and Detached Boundary Layers

NASA Astrophysics Data System (ADS)

Caruelle, B.; Ducros, F.

2003-12-01

This article presents Detached-Eddy Simulations (DESs) of attached and detached turbulent boundary layers. This hybrid Reynolds Averaged Navier-Stokes (RANS) / Large Eddy Simulation (LES) model goes continuously from RANS to LES according to the mesh definition. We propose a parametric study of the model over two "academic" configurations, in order to get information on the influence of the mesh to correctly treat complex flow with attached and detached boundary layers.

Finite element analysis of gradient z-coil induced eddy currents in a permanent MRI magnet.

PubMed

Li, Xia; Xia, Ling; Chen, Wufan; Liu, Feng; Crozier, Stuart; Xie, Dexin

2011-01-01

In permanent magnetic resonance imaging (MRI) systems, pulsed gradient fields induce strong eddy currents in the conducting structures of the magnet body. The gradient field for image encoding is perturbed by these eddy currents leading to MR image distortions. This paper presents a comprehensive finite element (FE) analysis of the eddy current generation in the magnet conductors. In the proposed FE model, the hysteretic characteristics of ferromagnetic materials are considered and a scalar Preisach hysteresis model is employed. The developed FE model was applied to study gradient z-coil induced eddy currents in a 0.5 T permanent MRI device. The simulation results demonstrate that the approach could be effectively used to investigate eddy current problems involving ferromagnetic materials. With the knowledge gained from this eddy current model, our next step is to design a passive magnet structure and active gradient coils to reduce the eddy current effects. Copyright © 2010 Elsevier Inc. All rights reserved.

Model-Simulated Northern Winter Cyclone and Anticyclone Activity under a Greenhouse Warming Scenario.

NASA Astrophysics Data System (ADS)

Zhang, Yi; Wang, Wei-Chyung

1997-07-01

Two 100-yr equilibrium simulations from the NCAR Community Climate Model coupled to a nondynamic slab ocean are used to investigate the activity of northern winter extratropical cyclones and anticyclones under a greenhouse warming scenario. The first simulation uses the 1990 observed CO2, CH4, N2O, CFC-11, and CFC-12 concentrations, and the second adopts the year 2050 concentrations according to the Intergovernmental Panel on Climate Change business-as-usual scenario. Variables that describe the characteristic properties of the cyclone-scale eddies, such as surface cyclone and anticyclone frequency and the bandpassed root-mean-square of 500-hPa geopotential height, along with the Eady growth rate maximum, form a framework for the analysis of the cyclone and anticyclone activity.Objective criteria are developed for identifying cyclone and anticyclone occurrences based on the 1000-hPa geopotential height and vorticity fields and tested using ECMWF analyses. The potential changes of the eddy activity under the greenhouse warming climate are then examined. Results indicate that the activity of cyclone-scale eddies decreases under the greenhouse warming scenario. This is not only reflected in the surface cyclone and anticyclone frequency and in the bandpassed rms of 500-hPa geopotential height, but is also discerned from the Eady growth rate maximum. Based on the analysis, three different physical mechanisms responsible for the decreased eddy activity are discussed: 1) a decrease of the extratropical meridional temperature gradient from the surface to the midtroposphere, 2) a reduction in the land-sea thermal contrast in the east coastal regions of the Asian and North American continents, and 3) an increase in the eddy meridional latent heat fluxes. Uncertainties in the results related to the limitations of the model and the model equilibrium simulations are discussed.

South Atlantic Ocean circulation: Simulation experiments with a quasi-geostrophic model and assimilation of TOPEX/POSEIDON and ERS 1 altimeter data

NASA Astrophysics Data System (ADS)

Florenchie, P.; Verron, J.

1998-10-01

Simulation experiments of South Atlantic Ocean circulations are conducted with a 1/6°, four-layered, quasi-geostrophic model. By means of a simple nudging data assimilation procedure along satellite tracks, TOPEX/POSEIDON and ERS 1 altimeter measurements are introduced into the model to control the simulation of the basin-scale circulation for the period from October 1992 to September 1994. The model circulation appears to be strongly influenced by the introduction of altimeter data, offering a consistent picture of South Atlantic Ocean circulations. Comparisons with observations show that the assimilating model successfully simulates the kinematic behavior of a large number of surface circulation components. The assimilation procedure enables us to produce schematic diagrams of South Atlantic circulation in which patterns ranging from basin-scale currents to mesoscale eddies are portrayed in a realistic way, with respect to their complexity. The major features of the South Atlantic circulation are described and analyzed, with special emphasis on the Brazil-Malvinas Confluence region, the Subtropical Gyre with the formation of frontal structures, and the Agulhas Retroflection. The Agulhas eddy-shedding process has been studied extensively. Fourteen eddies appear to be shed during the 2-year experiment. Because of their strong surface topographic signature, Agulhas eddies have been tracked continuously during the assimilation experiment as they cross the South Atlantic basin westward. Other effects of the assimilation procedure are shown, such as the intensification of the Subtropical Gyre, the appearance of a strong seasonal cycle in the Brazil Current transport, and the increase of the mean Brazil Current transport. This last result, combined with the westward oriention of the Agulhas eddies' trajectories, leads to a southward transport of mean eddy kinetic energy across 30°S.

Large-eddy simulation of flow in a plane, asymmetric diffuser

NASA Technical Reports Server (NTRS)

Kaltenbach, Hans-Jakob

1993-01-01

Recent improvements in subgrid-scale modeling as well as increases in computer power make it feasible to investigate flows using large-eddy simulation (LES) which have been traditionally studied with techniques based on Reynolds averaging. However, LES has not yet been applied to many flows of immediate technical interest. Preliminary results from LES of a plane diffuser flow are described. The long term goal of this work is to investigate flow separation as well as separation control in ducts and ramp-like geometries.

Large-Eddy Simulation of Internal Flow through Human Vocal Folds

NASA Astrophysics Data System (ADS)

Lasota, Martin; Šidlof, Petr

2018-06-01

The phonatory process occurs when air is expelled from the lungs through the glottis and the pressure drop causes flow-induced oscillations of the vocal folds. The flow fields created in phonation are highly unsteady and the coherent vortex structures are also generated. For accuracy it is essential to compute on humanlike computational domain and appropriate mathematical model. The work deals with numerical simulation of air flow within the space between plicae vocales and plicae vestibulares. In addition to the dynamic width of the rima glottidis, where the sound is generated, there are lateral ventriculus laryngis and sacculus laryngis included in the computational domain as well. The paper presents the results from OpenFOAM which are obtained with a large-eddy simulation using second-order finite volume discretization of incompressible Navier-Stokes equations. Large-eddy simulations with different subgrid scale models are executed on structured mesh. In these cases are used only the subgrid scale models which model turbulence via turbulent viscosity and Boussinesq approximation in subglottal and supraglottal area in larynx.

Improved Climate Simulations through a Stochastic Parameterization of Ocean Eddies

NASA Astrophysics Data System (ADS)

Williams, Paul; Howe, Nicola; Gregory, Jonathan; Smith, Robin; Joshi, Manoj

2016-04-01

In climate simulations, the impacts of the sub-grid scales on the resolved scales are conventionally represented using deterministic closure schemes, which assume that the impacts are uniquely determined by the resolved scales. Stochastic parameterization relaxes this assumption, by sampling the sub-grid variability in a computationally inexpensive manner. This presentation shows that the simulated climatological state of the ocean is improved in many respects by implementing a simple stochastic parameterization of ocean eddies into a coupled atmosphere-ocean general circulation model. Simulations from a high-resolution, eddy-permitting ocean model are used to calculate the eddy statistics needed to inject realistic stochastic noise into a low-resolution, non-eddy-permitting version of the same model. A suite of four stochastic experiments is then run to test the sensitivity of the simulated climate to the noise definition, by varying the noise amplitude and decorrelation time within reasonable limits. The addition of zero-mean noise to the ocean temperature tendency is found to have a non-zero effect on the mean climate. Specifically, in terms of the ocean temperature and salinity fields both at the surface and at depth, the noise reduces many of the biases in the low-resolution model and causes it to more closely resemble the high-resolution model. The variability of the strength of the global ocean thermohaline circulation is also improved. It is concluded that stochastic ocean perturbations can yield reductions in climate model error that are comparable to those obtained by refining the resolution, but without the increased computational cost. Therefore, stochastic parameterizations of ocean eddies have the potential to significantly improve climate simulations. Reference PD Williams, NJ Howe, JM Gregory, RS Smith, and MM Joshi (2016) Improved Climate Simulations through a Stochastic Parameterization of Ocean Eddies. Journal of Climate, under revision.

3D analysis of eddy current loss in the permanent magnet coupling.

PubMed

Zhu, Zina; Meng, Zhuo

2016-07-01

This paper first presents a 3D analytical model for analyzing the radial air-gap magnetic field between the inner and outer magnetic rotors of the permanent magnet couplings by using the Amperian current model. Based on the air-gap field analysis, the eddy current loss in the isolation cover is predicted according to the Maxwell's equations. A 3D finite element analysis model is constructed to analyze the magnetic field spatial distributions and vector eddy currents, and then the simulation results obtained are analyzed and compared with the analytical method. Finally, the current losses of two types of practical magnet couplings are measured in the experiment to compare with the theoretical results. It is concluded that the 3D analytical method of eddy current loss in the magnet coupling is viable and could be used for the eddy current loss prediction of magnet couplings.

Simulation of Atmospheric-Entry Capsules in the Subsonic Regime

NASA Technical Reports Server (NTRS)

Murman, Scott M.; Childs, Robert E.; Garcia, Joseph A.

2015-01-01

The accuracy of Computational Fluid Dynamics predictions of subsonic capsule aerodynamics is examined by comparison against recent NASA wind-tunnel data at high-Reynolds-number flight conditions. Several aspects of numerical and physical modeling are considered, including inviscid numerical scheme, mesh adaptation, rough-wall modeling, rotation and curvature corrections for eddy-viscosity models, and Detached-Eddy Simulations of the unsteady wake. All of these are considered in isolation against relevant data where possible. The results indicate that an improved predictive capability is developed by considering physics-based approaches and validating the results against flight-relevant experimental data.

Impact of Preferred Eddy Tracks on Transport and Mixing in the Eastern South Pacific

NASA Astrophysics Data System (ADS)

Belmadani, A.; Donoso, D.; Auger, P. A.; Chaigneau, A.

2017-12-01

Mesoscale eddies, which play a fundamental role in the transport of mass, heat, nutrients, and biota across the oceans, have been suggested to propagate preferently along specific tracks. These preferred pathways, also called eddy trains, are near-zonal due to westward drift of individual vortices, and tend to be polarized (ie alternatively dominated by anticyclonic/cyclonic eddies), coinciding with the recently discovered latent striations (quasi-zonal mesoscale jet-like features). While significant effort has been made to understand the dynamics of striations and their interplay with mesoscale eddies, the impact of repeated eddy tracks on physical (temperature, salinity), biogeochemical (oxygen, carbon, nutrients) and other tracers (e.g. chlorophyll, marine debris) has received little attention. Here we report on the results of numerical modeling experiments that simulate the impact of preferred eddy tracks on the transport and mixing of water particles in the Eastern South Pacific off Chile. A 30-year interannual simulation of the oceanic circulation in this region has been performed over 1984-2013 with the ROMS (Regional Oceanic Modeling System) at an eddy-resolving resolution (10 km). Objective tracking of mesoscale coherent vortices is obtained using automated methods, allowing to compute the contribution of eddies to the ocean circulation. Preferred eddy tracks are further isolated from the more random eddies, by comparing the distances between individual tracks and the striated pattern in long-term mean eddy polarity with a least-squares approach. The remaining non-eddying flow may also be decomposed into time-mean and anomalous circulation, and/or small- and large-scale circulation. Neutrally-buoyant Lagrangian floats are then released uniformly into the various flow components as well as the total flow, and tracked forward in time with the ARIANE software. The dispersion patterns of water particles are used to estimate the respective contributions of organized and random eddies, mean flow, large-scale perturbations etc. to mixing properties and transport pathways. Float release into the full flow inside selected vortices is also used to document the impact of eddy trains on the transformation of water masses inferred from changes in temperature/salinity along float trajectories.

Application of renormalization group theory to the large-eddy simulation of transitional boundary layers

NASA Technical Reports Server (NTRS)

Piomelli, Ugo; Zang, Thomas A.; Speziale, Charles G.; Lund, Thomas S.

1990-01-01

An eddy viscosity model based on the renormalization group theory of Yakhot and Orszag (1986) is applied to the large-eddy simulation of transition in a flat-plate boundary layer. The simulation predicts with satisfactory accuracy the mean velocity and Reynolds stress profiles, as well as the development of the important scales of motion. The evolution of the structures characteristic of the nonlinear stages of transition is also predicted reasonably well.

Immersed transient eddy current flow metering: a calibration-free velocity measurement technique for liquid metals

NASA Astrophysics Data System (ADS)

Krauter, N.; Stefani, F.

2017-10-01

Eddy current flow meters are widely used for measuring the flow velocity of electrically conducting fluids. Since the flow induced perturbations of a magnetic field depend both on the geometry and the conductivity of the fluid, extensive calibration is needed to get accurate results. Transient eddy current flow metering has been developed to overcome this problem. It relies on tracking the position of an impressed eddy current system that is moving with the same velocity as the conductive fluid. We present an immersed version of this measurement technique and demonstrate its viability by numerical simulations and a first experimental validation.

Evaluation of scale-aware subgrid mesoscale eddy models in a global eddy-rich model

NASA Astrophysics Data System (ADS)

Pearson, Brodie; Fox-Kemper, Baylor; Bachman, Scott; Bryan, Frank

2017-07-01

Two parameterizations for horizontal mixing of momentum and tracers by subgrid mesoscale eddies are implemented in a high-resolution global ocean model. These parameterizations follow on the techniques of large eddy simulation (LES). The theory underlying one parameterization (2D Leith due to Leith, 1996) is that of enstrophy cascades in two-dimensional turbulence, while the other (QG Leith) is designed for potential enstrophy cascades in quasi-geostrophic turbulence. Simulations using each of these parameterizations are compared with a control simulation using standard biharmonic horizontal mixing.Simulations using the 2D Leith and QG Leith parameterizations are more realistic than those using biharmonic mixing. In particular, the 2D Leith and QG Leith simulations have more energy in resolved mesoscale eddies, have a spectral slope more consistent with turbulence theory (an inertial enstrophy or potential enstrophy cascade), have bottom drag and vertical viscosity as the primary sinks of energy instead of lateral friction, and have isoneutral parameterized mesoscale tracer transport. The parameterization choice also affects mass transports, but the impact varies regionally in magnitude and sign.

Hybrid Reynolds-Averaged/Large-Eddy Simulations of a Coaxial Supersonic Free-Jet Experiment

NASA Technical Reports Server (NTRS)

Baurle, Robert A.; Edwards, Jack R.

2010-01-01

Reynolds-averaged and hybrid Reynolds-averaged/large-eddy simulations have been applied to a supersonic coaxial jet flow experiment. The experiment was designed to study compressible mixing flow phenomenon under conditions that are representative of those encountered in scramjet combustors. The experiment utilized either helium or argon as the inner jet nozzle fluid, and the outer jet nozzle fluid consisted of laboratory air. The inner and outer nozzles were designed and operated to produce nearly pressure-matched Mach 1.8 flow conditions at the jet exit. The purpose of the computational effort was to assess the state-of-the-art for each modeling approach, and to use the hybrid Reynolds-averaged/large-eddy simulations to gather insight into the deficiencies of the Reynolds-averaged closure models. The Reynolds-averaged simulations displayed a strong sensitivity to choice of turbulent Schmidt number. The initial value chosen for this parameter resulted in an over-prediction of the mixing layer spreading rate for the helium case, but the opposite trend was observed when argon was used as the injectant. A larger turbulent Schmidt number greatly improved the comparison of the results with measurements for the helium simulations, but variations in the Schmidt number did not improve the argon comparisons. The hybrid Reynolds-averaged/large-eddy simulations also over-predicted the mixing layer spreading rate for the helium case, while under-predicting the rate of mixing when argon was used as the injectant. The primary reason conjectured for the discrepancy between the hybrid simulation results and the measurements centered around issues related to the transition from a Reynolds-averaged state to one with resolved turbulent content. Improvements to the inflow conditions were suggested as a remedy to this dilemma. Second-order turbulence statistics were also compared to their modeled Reynolds-averaged counterparts to evaluate the effectiveness of common turbulence closure assumptions.

Large-eddy simulation of laminar-turbulent breakdown at high speeds with dynamic subgrid-scale modeling

NASA Technical Reports Server (NTRS)

El-Hady, Nabil M.

1993-01-01

The laminar-turbulent breakdown of a boundary-layer flow along a hollow cylinder at Mach 4.5 is investigated with large-eddy simulation. The subgrid scales are modeled dynamically, where the model coefficients are determined from the local resolved field. The behavior of the dynamic-model coefficients is investigated through both an a priori test with direct numerical simulation data for the same case and a complete large-eddy simulation. Both formulations proposed by Germano et al. and Lilly are used for the determination of unique coefficients for the dynamic model and their results are compared and assessed. The behavior and the energy cascade of the subgrid-scale field structure are investigated at various stages of the transition process. The investigations are able to duplicate a high-speed transition phenomenon observed in experiments and explained only recently by the direct numerical simulations of Pruett and Zang, which is the appearance of 'rope-like' waves. The nonlinear evolution and breakdown of the laminar boundary layer and the structure of the flow field during the transition process were also investigated.

Development of eddy current probe for fiber orientation assessment in carbon fiber composites

NASA Astrophysics Data System (ADS)

Wincheski, Russell A.; Zhao, Selina

2018-04-01

Measurement of the fiber orientation in a carbon fiber composite material is crucial in understanding the load carrying capability of the structure. As manufacturing conditions including resin flow and molding pressures can alter fiber orientation, verification of the as-designed fiber layup is necessary to ensure optimal performance of the structure. In this work, the development of an eddy current probe and data processing technique for analysis of fiber orientation in carbon fiber composites is presented. A proposed directional eddy current probe is modeled and its response to an anisotropic multi-layer conductor simulated. The modeling results are then used to finalize specifications of the eddy current probe. Experimental testing of the fabricated probe is presented for several samples including a truncated pyramid part with complex fiber orientation draped to the geometry for resin transfer molding. The inductively coupled single sided measurement enables fiber orientation characterization through the thickness of the part. The fast and cost-effective technique can be applied as a spot check or as a surface map of the fiber orientations across the structure. This paper will detail the results of the probe design, computer simulations, and experimental results.

Numerical simulation of flow around the NREL S826 airfoil at moderate Reynolds number using delayed detached Eddy simulation (DDES)

NASA Astrophysics Data System (ADS)

Prytz, Erik R.; Huuse, Øyvind; Müller, Bernhard; Bartl, Jan; Sætran, Lars Roar

2017-07-01

Turbulent flow at Reynolds numbers 5 . 104 to 106 around the NREL S826 airfoil used for wind turbine blades is simulated using delayed detached eddy simulation (DDES). The 3D domain is built as a replica of the low speed wind tunnel at the Norwegian University of Science and Technology (NTNU) with the wind tunnel walls considered as slip walls. The subgrid turbulent kinetic energy is used to model the sub-grid scale in the large eddy simulation (LES) part of DDES. Different Reynoldsaveraged Navier-Stokes (RANS) models are tested in ANSYS Fluent. The realizable k - ∈ model as the RANS model in DDES is found to yield the best agreement of simulated pressure distributions with the experimental data both from NTNU and the Technical University of Denmark (DTU), the latter for a shorter spanwise domain. The present DDES results are in excellent agreement with LES results from DTU. Since DDES requires much fewer cells in the RANS region near the wing surface than LES, DDES is computationally much more efficient than LES. Whereas DDES is able to predict lift and drag in close agreement with experiment up to stall, pure 2D RANS simulations fail near stall. After testing different numerical settings, time step sizes and grids for DDES, a Reynolds number study is conducted. Near stall, separated flow structures, so-called stall cells, are observed in the DDES results.

Calculations of High-Temperature Jet Flow Using Hybrid Reynolds-Average Navier-Stokes Formulations

NASA Technical Reports Server (NTRS)

Abdol-Hamid, Khaled S.; Elmiligui, Alaa; Giriamaji, Sharath S.

2008-01-01

Two multiscale-type turbulence models are implemented in the PAB3D solver. The models are based on modifying the Reynolds-averaged Navier Stokes equations. The first scheme is a hybrid Reynolds-averaged- Navier Stokes/large-eddy-simulation model using the two-equation k(epsilon) model with a Reynolds-averaged-Navier Stokes/large-eddy-simulation transition function dependent on grid spacing and the computed turbulence length scale. The second scheme is a modified version of the partially averaged Navier Stokes model in which the unresolved kinetic energy parameter f(sub k) is allowed to vary as a function of grid spacing and the turbulence length scale. This parameter is estimated based on a novel two-stage procedure to efficiently estimate the level of scale resolution possible for a given flow on a given grid for partially averaged Navier Stokes. It has been found that the prescribed scale resolution can play a major role in obtaining accurate flow solutions. The parameter f(sub k) varies between zero and one and is equal to one in the viscous sublayer and when the Reynolds-averaged Navier Stokes turbulent viscosity becomes smaller than the large-eddy-simulation viscosity. The formulation, usage methodology, and validation examples are presented to demonstrate the enhancement of PAB3D's time-accurate turbulence modeling capabilities. The accurate simulations of flow and turbulent quantities will provide a valuable tool for accurate jet noise predictions. Solutions from these models are compared with Reynolds-averaged Navier Stokes results and experimental data for high-temperature jet flows. The current results show promise for the capability of hybrid Reynolds-averaged Navier Stokes and large eddy simulation and partially averaged Navier Stokes in simulating such flow phenomena.

Evaluation of machine learning algorithms for prediction of regions of high Reynolds averaged Navier Stokes uncertainty

DOE PAGES

Ling, Julia; Templeton, Jeremy Alan

2015-08-04

Reynolds Averaged Navier Stokes (RANS) models are widely used in industry to predict fluid flows, despite their acknowledged deficiencies. Not only do RANS models often produce inaccurate flow predictions, but there are very limited diagnostics available to assess RANS accuracy for a given flow configuration. If experimental or higher fidelity simulation results are not available for RANS validation, there is no reliable method to evaluate RANS accuracy. This paper explores the potential of utilizing machine learning algorithms to identify regions of high RANS uncertainty. Three different machine learning algorithms were evaluated: support vector machines, Adaboost decision trees, and random forests.more » The algorithms were trained on a database of canonical flow configurations for which validated direct numerical simulation or large eddy simulation results were available, and were used to classify RANS results on a point-by-point basis as having either high or low uncertainty, based on the breakdown of specific RANS modeling assumptions. Classifiers were developed for three different basic RANS eddy viscosity model assumptions: the isotropy of the eddy viscosity, the linearity of the Boussinesq hypothesis, and the non-negativity of the eddy viscosity. It is shown that these classifiers are able to generalize to flows substantially different from those on which they were trained. As a result, feature selection techniques, model evaluation, and extrapolation detection are discussed in the context of turbulence modeling applications.« less

Time-Domain Filtering for Spatial Large-Eddy Simulation

NASA Technical Reports Server (NTRS)

Pruett, C. David

1997-01-01

An approach to large-eddy simulation (LES) is developed whose subgrid-scale model incorporates filtering in the time domain, in contrast to conventional approaches, which exploit spatial filtering. The method is demonstrated in the simulation of a heated, compressible, axisymmetric jet, and results are compared with those obtained from fully resolved direct numerical simulation. The present approach was, in fact, motivated by the jet-flow problem and the desire to manipulate the flow by localized (point) sources for the purposes of noise suppression. Time-domain filtering appears to be more consistent with the modeling of point sources; moreover, time-domain filtering may resolve some fundamental inconsistencies associated with conventional space-filtered LES approaches.

Large Eddy Simulations using oodlesDST

DTIC Science & Technology

2016-01-01

Research Agency DST-Group-TR-3205 ABSTRACT The oodlesDST code is based on OpenFOAM software and performs Large Eddy Simulations of......maritime platforms using a variety of simulation techniques. He is currently using OpenFOAM software to perform both Reynolds Averaged Navier-Stokes

An efficient 3-D eddy-current solver using an independent impedance method for transcranial magnetic stimulation.

PubMed

De Geeter, Nele; Crevecoeur, Guillaume; Dupre, Luc

2011-02-01

In many important bioelectromagnetic problem settings, eddy-current simulations are required. Examples are the reduction of eddy-current artifacts in magnetic resonance imaging and techniques, whereby the eddy currents interact with the biological system, like the alteration of the neurophysiology due to transcranial magnetic stimulation (TMS). TMS has become an important tool for the diagnosis and treatment of neurological diseases and psychiatric disorders. A widely applied method for simulating the eddy currents is the impedance method (IM). However, this method has to contend with an ill conditioned problem and consequently a long convergence time. When dealing with optimal design problems and sensitivity control, the convergence rate becomes even more crucial since the eddy-current solver needs to be evaluated in an iterative loop. Therefore, we introduce an independent IM (IIM), which improves the conditionality and speeds up the numerical convergence. This paper shows how IIM is based on IM and what are the advantages. Moreover, the method is applied to the efficient simulation of TMS. The proposed IIM achieves superior convergence properties with high time efficiency, compared to the traditional IM and is therefore a useful tool for accurate and fast TMS simulations.

How does the Red Sea outflow water interact with Gulf of Aden Eddies?

NASA Astrophysics Data System (ADS)

Ilıcak, Mehmet; Özgökmen, Tamay M.; Johns, William E.

As the Red Sea overflow water (RSOW) enters the Gulf of Aden (GOA), it interacts with a sequence of nearly barotropic, mesoscale eddies originating in the Indian Ocean. To investigate how these eddies impact the dispersal and eastward transport of the RSOW toward the Indian Ocean, a high resolution 3D regional model is employed to explore systematically the interaction between the RSOW and mesoscale eddies. Two types of experiments are conducted. In the first set, we simulate the behavior of RSOW in the presence of an idealized cyclone and an idealized anticyclone. The second type of simulation involves nesting of the regional model (ROMS) within a data-assimilating global model (HYCOM), in which a sequence of mesoscale eddies entering the Gulf of Aden is realistically captured. This simulation is integrated for one year, and includes a simple representation of the seasonality of the RSOW. Bower et al. (2002) suggest that the Red Sea overflow might be a western boundary undercurrent. Consistent with these expectations, the idealized simulations show that the preferred pathway of the RSOW in the absence of eddies is along the coast of Somalia (southern continental shelf) as a western boundary undercurrent. Simultaneously, a cyclonic circulation is generated in the far western GOA due to vortex stretching by the descending outflow. The presence of a cyclone in the western GOA increases the peak RSOW transport, but the cyclone itself rapidly loses its coherence after interacting with the rough topography in the western GOA. The presence of an anticyclone tends to block the preferred boundary pathway and inhibits the eastward transport of the RSOW. The eddies also result in substantially increased mixing of the RSOW in the western GOA. On the basis of the more realistic ROMS experiment, it is found that the modeled RSOW leaves the western part of the Gulf of Aden in short episodic bursts with transports that are an order of magnitude greater than that associated with the quasi-steady RSOW inflow into GOA. Such enhancement in RSOW transport is shown to be induced by cyclonic eddies that cause a rapid discharge of RSOW from the western part of the GOA. We conclude that mesoscale eddies play a key role in the transport and mixing of the RSOW within GOA.

Hybrid Reynolds-Averaged/Large-Eddy Simulations of a Co-Axial Supersonic Free-Jet Experiment

NASA Technical Reports Server (NTRS)

Baurle, R. A.; Edwards, J. R.

2009-01-01

Reynolds-averaged and hybrid Reynolds-averaged/large-eddy simulations have been applied to a supersonic coaxial jet flow experiment. The experiment utilized either helium or argon as the inner jet nozzle fluid, and the outer jet nozzle fluid consisted of laboratory air. The inner and outer nozzles were designed and operated to produce nearly pressure-matched Mach 1.8 flow conditions at the jet exit. The purpose of the computational effort was to assess the state-of-the-art for each modeling approach, and to use the hybrid Reynolds-averaged/large-eddy simulations to gather insight into the deficiencies of the Reynolds-averaged closure models. The Reynolds-averaged simulations displayed a strong sensitivity to choice of turbulent Schmidt number. The baseline value chosen for this parameter resulted in an over-prediction of the mixing layer spreading rate for the helium case, but the opposite trend was noted when argon was used as the injectant. A larger turbulent Schmidt number greatly improved the comparison of the results with measurements for the helium simulations, but variations in the Schmidt number did not improve the argon comparisons. The hybrid simulation results showed the same trends as the baseline Reynolds-averaged predictions. The primary reason conjectured for the discrepancy between the hybrid simulation results and the measurements centered around issues related to the transition from a Reynolds-averaged state to one with resolved turbulent content. Improvements to the inflow conditions are suggested as a remedy to this dilemma. Comparisons between resolved second-order turbulence statistics and their modeled Reynolds-averaged counterparts were also performed.

Detached eddy simulation for turbulent fluid-structure interaction of moving bodies using the constraint-based immersed boundary method

NASA Astrophysics Data System (ADS)

Nangia, Nishant; Bhalla, Amneet P. S.; Griffith, Boyce E.; Patankar, Neelesh A.

2016-11-01

Flows over bodies of industrial importance often contain both an attached boundary layer region near the structure and a region of massively separated flow near its trailing edge. When simulating these flows with turbulence modeling, the Reynolds-averaged Navier-Stokes (RANS) approach is more efficient in the former, whereas large-eddy simulation (LES) is more accurate in the latter. Detached-eddy simulation (DES), based on the Spalart-Allmaras model, is a hybrid method that switches from RANS mode of solution in attached boundary layers to LES in detached flow regions. Simulations of turbulent flows over moving structures on a body-fitted mesh incur an enormous remeshing cost every time step. The constraint-based immersed boundary (cIB) method eliminates this operation by placing the structure on a Cartesian mesh and enforcing a rigidity constraint as an additional forcing in the Navier-Stokes momentum equation. We outline the formulation and development of a parallel DES-cIB method using adaptive mesh refinement. We show preliminary validation results for flows past stationary bodies with both attached and separated boundary layers along with results for turbulent flows past moving bodies. This work is supported by the National Science Foundation Graduate Research Fellowship under Grant No. DGE-1324585.

Large Eddy Simulations (LES) and Direct Numerical Simulations (DNS) for the computational analyses of high speed reacting flows

NASA Technical Reports Server (NTRS)

Givi, Peyman; Madnia, Cyrus K.; Steinberger, C. J.; Frankel, S. H.

1992-01-01

The principal objective is to extend the boundaries within which large eddy simulations (LES) and direct numerical simulations (DNS) can be applied in computational analyses of high speed reacting flows. A summary of work accomplished during the last six months is presented.

Lateral eddy diffusivity estimates from simulated and observed drifter trajectories: a case study for the Agulhas Current system

NASA Astrophysics Data System (ADS)

Rühs, Siren; Zhurbas, Victor; Durgadoo, Jonathan V.; Biastoch, Arne

2017-04-01

The Lagrangian description of fluid motion by sets of individual particle trajectories is extensively used to characterize connectivity between distinct oceanic locations. One important factor influencing the connectivity is the average rate of particle dispersal, generally quantified as Lagrangian diffusivity. In addition to Lagrangian observing programs, Lagrangian analyses are performed by advecting particles with the simulated flow field of ocean general circulation models (OGCMs). However, depending on the spatio-temporal model resolution, not all scale-dependent processes are explicitly resolved in the simulated velocity fields. Consequently, the dispersal of advective Lagrangian trajectories has been assumed not to be sufficiently diffusive compared to observed particle spreading. In this study we present a detailed analysis of the spatially variable lateral eddy diffusivity characteristics of advective drifter trajectories simulated with realistically forced OGCMs and compare them with estimates based on observed drifter trajectories. The extended Agulhas Current system around South Africa, known for its intricate mesoscale dynamics, serves as a test case. We show that a state-of-the-art eddy-resolving OGCM indeed features theoretically derived dispersion characteristics for diffusive regimes and realistically represents Lagrangian eddy diffusivity characteristics obtained from observed surface drifter trajectories. The estimates for the maximum and asymptotic lateral single-particle eddy diffusivities obtained from the observed and simulated drifter trajectories show a good agreement in their spatial pattern and magnitude. We further assess the sensitivity of the simulated lateral eddy diffusivity estimates to the temporal and lateral OGCM output resolution and examine the impact of the different eddy diffusivity characteristics on the Lagrangian connectivity between the Indian Ocean and the South Atlantic.

Large-Eddy Simulation of Propeller Crashback

NASA Astrophysics Data System (ADS)

Kumar, Praveen; Mahesh, Krishnan

2013-11-01

Crashback is an operating condition to quickly stop a propelled vehicle, where the propeller is rotated in the reverse direction to yield negative thrust. The crashback condition is dominated by the interaction of free stream flow with strong reverse flow. Crashback causes highly unsteady loads and flow separation on blade surface. This study uses Large-Eddy Simulation to predict the highly unsteady flow field in propeller crashback. Results are shown for a stand-alone open propeller, hull-attached open propeller and a ducted propeller. The simulations are compared to experiment, and used to discuss the essential physics behind the unsteady loads. This work is supported by the Office of Naval Research.

Influence of Kuroshio Oceanic Eddies on North Pacific Weather Patterns

NASA Astrophysics Data System (ADS)

Ma, X.; Chang, P.; Saravanan, R.; Montuoro, R.; Hsieh, J. S.; Wu, D.; Lin, X.; Wu, L.; Jing, Z.

2016-02-01

High-resolution satellite observations reveal energetic meso-scale ocean eddy activity and positive correlation between meso-scale sea surface temperature (SST) and surface wind along oceanic frontal zones, such as the Kuroshio and Gulf Stream, suggesting a potential role of meso-scale oceanic eddies in forcing the atmosphere. Using a 27 km horizontal resolution Weather Research Forecasting (WRF) model forced with observed daily SST at 0.09° spatial resolution during boreal winter season, two ensembles of 10 WRF simulations, in one of which meso-scale SST variability induced by ocean eddies was suppressed, were conducted in the North Pacific to study the local and remote influence of meso-scale oceanic eddies in the Kuroshio Extention Region (KER) on the atmosphere. Suppression of meso-scale oceanic eddies results in a deep tropospheric response along and downstream of the KER, including a significant decrease (increase) in winter season mean rainfall along the KER (west coast of US), a reduction of storm genesis in the KER, and a southward shift of the jet stream and North Pacific storm track in the eastern North Pacific. The simulated local and remote rainfall response to meso-scale oceanic eddies in the KER is also supported by observational analysis. A mechanism invoking moist baroclinic instability is proposed as a plausible explanation for the linkage between meso-scale oceanic eddies in the KER and large-scale atmospheric response in the North Pacific. It is argued that meso-scale oceanic eddies can have a rectified effect on planetary boundary layer moisture, the stability of the lower atmosphere and latent heat release, which in turn affect cyclogenesis. The accumulated effect of the altered storm development downstream further contributes to the equivalent barotropic mean flow change in the eastern North Pacific basin.

Impacts of the Mesoscale Ocean-Atmosphere Coupling on the Peru-Chile Ocean Dynamics: The Current-Induced Wind Stress Modulation

NASA Astrophysics Data System (ADS)

Oerder, V.; Colas, F.; Echevin, V.; Masson, S.; Lemarié, F.

2018-02-01

The ocean dynamical responses to the surface current-wind stress interaction at the oceanic mesoscale are investigated in the South-East Pacific using a high-resolution regional ocean-atmosphere coupled model. Two simulations are compared: one includes the surface current in the wind stress computation while the other does not. In the coastal region, absolute wind velocities are different between the two simulations but the wind stress remains very similar. As a consequence, the mean regional oceanic circulation is almost unchanged. On the contrary, the mesoscale activity is strongly reduced when taking into account the effect of the surface current on the wind stress. This is caused by a weakening of the eddy kinetic energy generation near the coast by the wind work and to intensified offshore eddy damping. We show that, above coherent eddies, the current-stress interaction generates eddy damping through Ekman pumping and eddy kinetic energy dissipation through wind work. This alters significantly the coherent eddy vertical structures compared with the control simulation, weakening the temperature and vorticity anomalies and increasing strongly the vertical velocity anomalies associated to eddies.

Turbulence and pollutant transport in urban street canyons under stable stratification: a large-eddy simulation

NASA Astrophysics Data System (ADS)

Li, X.

2014-12-01

Thermal stratification of the atmospheric surface layer has strong impact on the land-atmosphere exchange of turbulent, heat, and pollutant fluxes. Few studies have been carried out for the interaction of the weakly to moderately stable stratified atmosphere and the urban canopy. This study performs a large-eddy simulation of a modeled street canyon within a weakly to moderately stable atmosphere boundary layer. To better resolve the smaller eddy size resulted from the stable stratification, a higher spatial and temporal resolution is used. The detailed flow structure and turbulence inside the street canyon are analyzed. The relationship of pollutant dispersion and Richardson number of the atmosphere is investigated. Differences between these characteristics and those under neutral and unstable atmosphere boundary layer are emphasized.

Large-eddy simulation, fuel rod vibration and grid-to-rod fretting in pressurized water reactors

DOE PAGES

Christon, Mark A.; Lu, Roger; Bakosi, Jozsef; ...

2016-10-01

Grid-to-rod fretting (GTRF) in pressurized water reactors is a flow-induced vibration phenomenon that results in wear and fretting of the cladding material on fuel rods. GTRF is responsible for over 70% of the fuel failures in pressurized water reactors in the United States. Predicting the GTRF wear and concomitant interval between failures is important because of the large costs associated with reactor shutdown and replacement of fuel rod assemblies. The GTRF-induced wear process involves turbulent flow, mechanical vibration, tribology, and time-varying irradiated material properties in complex fuel assembly geometries. This paper presents a new approach for predicting GTRF induced fuelmore » rod wear that uses high-resolution implicit large-eddy simulation to drive nonlinear transient dynamics computations. The GTRF fluid–structure problem is separated into the simulation of the turbulent flow field in the complex-geometry fuel-rod bundles using implicit large-eddy simulation, the calculation of statistics of the resulting fluctuating structural forces, and the nonlinear transient dynamics analysis of the fuel rod. Ultimately, the methods developed here, can be used, in conjunction with operational management, to improve reactor core designs in which fuel rod failures are minimized or potentially eliminated. Furthermore, robustness of the behavior of both the structural forces computed from the turbulent flow simulations and the results from the transient dynamics analyses highlight the progress made towards achieving a predictive simulation capability for the GTRF problem.« less

Large-eddy simulation, fuel rod vibration and grid-to-rod fretting in pressurized water reactors

DOE Office of Scientific and Technical Information (OSTI.GOV)

Christon, Mark A.; Lu, Roger; Bakosi, Jozsef

Grid-to-rod fretting (GTRF) in pressurized water reactors is a flow-induced vibration phenomenon that results in wear and fretting of the cladding material on fuel rods. GTRF is responsible for over 70% of the fuel failures in pressurized water reactors in the United States. Predicting the GTRF wear and concomitant interval between failures is important because of the large costs associated with reactor shutdown and replacement of fuel rod assemblies. The GTRF-induced wear process involves turbulent flow, mechanical vibration, tribology, and time-varying irradiated material properties in complex fuel assembly geometries. This paper presents a new approach for predicting GTRF induced fuelmore » rod wear that uses high-resolution implicit large-eddy simulation to drive nonlinear transient dynamics computations. The GTRF fluid–structure problem is separated into the simulation of the turbulent flow field in the complex-geometry fuel-rod bundles using implicit large-eddy simulation, the calculation of statistics of the resulting fluctuating structural forces, and the nonlinear transient dynamics analysis of the fuel rod. Ultimately, the methods developed here, can be used, in conjunction with operational management, to improve reactor core designs in which fuel rod failures are minimized or potentially eliminated. Furthermore, robustness of the behavior of both the structural forces computed from the turbulent flow simulations and the results from the transient dynamics analyses highlight the progress made towards achieving a predictive simulation capability for the GTRF problem.« less

Southern Ocean eddy compensation in a forced eddy-resolving GCM

NASA Astrophysics Data System (ADS)

Bruun Poulsen, Mads; Jochum, Markus; Eden, Carsten; Nuterman, Roman

2017-04-01

Contemporary eddy-resolving model studies have demonstrated that the common parameterisation of isopycnal mixing in the ocean is subject to limitations in the Southern Ocean where the mesoscale eddies are of leading order importance to the dynamics. We here present forced simulations from the Community Earth System Model on a global {1/10}° and 1° horizontal grid, the latter employing an eddy parameterisation, where the strength of the zonal wind stress south of 25°S has been varied. With a 50% zonally symmetric increase of the wind stress, we show that the two models arrive at two radically different solutions in terms of the large-scale circulation, with an increase of the deep inflow of water to the Southern Ocean at 40°S by 50% in the high resolution model against 20% at coarse resolution. Together with a weaker vertical displacement of the pycnocline in the 1° model, these results suggest that the parameterised eddies have an overly strong compensating effect on the water mass transformation compared to the explicit eddies. Implications for eddy mixing parameterisations will be discussed.

Evaluation of eddy-current proximity devices for measuring thin potassium film thicknesses

NASA Technical Reports Server (NTRS)

Asadourian, A. S.

1972-01-01

Two eddy current proximity probe systems were tested over a range of 0 to 508 micrometers (0 to 20 mils) of simulated potassium film thicknesses for simulated temperatures of 66 C (150 F), 232 C (450 F), and 666 C (1230 F). The results of short time calibration tests are presented. Instrument drift was a problem throughout the testing and, without correction, may limit the use of such systems to short periods of time. Additional development will be required prior to their being usable as practical instrumentation systems.

Recent Advances in Simulation of Eddy Current Testing of Tubes and Experimental Validations

DOE Office of Scientific and Technical Information (OSTI.GOV)

Reboud, C.; Premel, D.; Lesselier, D.

2007-03-21

Eddy current testing (ECT) is widely used in iron and steel industry for the inspection of tubes during manufacturing. A collaboration between CEA and the Vallourec Research Center led to the development of new numerical functionalities dedicated to the simulation of ECT of non-magnetic tubes by external probes. The achievement of experimental validations led us to the integration of these models into the CIVA platform. Modeling approach and validation results are discussed here. A new numerical scheme is also proposed in order to improve the accuracy of the model.

Recent Advances in Simulation of Eddy Current Testing of Tubes and Experimental Validations

NASA Astrophysics Data System (ADS)

Reboud, C.; Prémel, D.; Lesselier, D.; Bisiaux, B.

2007-03-01

Eddy current testing (ECT) is widely used in iron and steel industry for the inspection of tubes during manufacturing. A collaboration between CEA and the Vallourec Research Center led to the development of new numerical functionalities dedicated to the simulation of ECT of non-magnetic tubes by external probes. The achievement of experimental validations led us to the integration of these models into the CIVA platform. Modeling approach and validation results are discussed here. A new numerical scheme is also proposed in order to improve the accuracy of the model.

Filter size definition in anisotropic subgrid models for large eddy simulation on irregular grids

NASA Astrophysics Data System (ADS)

Abbà, Antonella; Campaniello, Dario; Nini, Michele

2017-06-01

The definition of the characteristic filter size to be used for subgrid scales models in large eddy simulation using irregular grids is still an unclosed problem. We investigate some different approaches to the definition of the filter length for anisotropic subgrid scale models and we propose a tensorial formulation based on the inertial ellipsoid of the grid element. The results demonstrate an improvement in the prediction of several key features of the flow when the anisotropicity of the grid is explicitly taken into account with the tensorial filter size.

Seasonal Variability in Global Eddy Diffusion and the Effect on Thermospheric Neutral Density

NASA Astrophysics Data System (ADS)

Pilinski, M.; Crowley, G.

2014-12-01

We describe a method for making single-satellite estimates of the seasonal variability in global-average eddy diffusion coefficients. Eddy diffusion values as a function of time between January 2004 and January 2008 were estimated from residuals of neutral density measurements made by the CHallenging Minisatellite Payload (CHAMP) and simulations made using the Thermosphere Ionosphere Mesosphere Electrodynamics - Global Circulation Model (TIME-GCM). The eddy diffusion coefficient results are quantitatively consistent with previous estimates based on satellite drag observations and are qualitatively consistent with other measurement methods such as sodium lidar observations and eddy-diffusivity models. The eddy diffusion coefficient values estimated between January 2004 and January 2008 were then used to generate new TIME-GCM results. Based on these results, the RMS difference between the TIME-GCM model and density data from a variety of satellites is reduced by an average of 5%. This result, indicates that global thermospheric density modeling can be improved by using data from a single satellite like CHAMP. This approach also demonstrates how eddy diffusion could be estimated in near real-time from satellite observations and used to drive a global circulation model like TIME-GCM. Although the use of global values improves modeled neutral densities, there are some limitations of this method, which are discussed, including that the latitude-dependence of the seasonal neutral-density signal is not completely captured by a global variation of eddy diffusion coefficients. This demonstrates the need for a latitude-dependent specification of eddy diffusion consistent with diffusion observations made by other techniques.

Mitigation of eddy current heating during magnetic nanoparticle hyperthermia therapy.

PubMed

Stigliano, Robert V; Shubitidze, Fridon; Petryk, James D; Shoshiashvili, Levan; Petryk, Alicia A; Hoopes, P Jack

2016-11-01

Magnetic nanoparticle hyperthermia therapy is a promising technology for cancer treatment, involving delivering magnetic nanoparticles (MNPs) into tumours then activating them using an alternating magnetic field (AMF). The system produces not only a magnetic field, but also an electric field which penetrates normal tissue and induces eddy currents, resulting in unwanted heating of normal tissues. Magnitude of the eddy current depends, in part, on the AMF source and the size of the tissue exposed to the field. The majority of in vivo MNP hyperthermia therapy studies have been performed in small animals, which, due to the spatial distribution of the AMF relative to the size of the animals, do not reveal the potential toxicity of eddy current heating in larger tissues. This has posed a non-trivial challenge for researchers attempting to scale up to clinically relevant volumes of tissue. There is a relative dearth of studies focused on decreasing the maximum temperature resulting from eddy current heating to increase therapeutic ratio. This paper presents two simple, clinically applicable techniques for decreasing maximum temperature induced by eddy currents. Computational and experimental results are presented to understand the underlying physics of eddy currents induced in conducting, biological tissues and leverage these insights to mitigate eddy current heating during MNP hyperthermia therapy. Phantom studies show that the displacement and motion techniques reduce maximum temperature due to eddy currents by 74% and 19% in simulation, and by 77% and 33% experimentally. Further study is required to optimise these methods for particular scenarios; however, these results suggest larger volumes of tissue could be treated, and/or higher field strengths and frequencies could be used to attain increased MNP heating when these eddy current mitigation techniques are employed.

Geometrical optics analysis of atmospheric turbulence

NASA Astrophysics Data System (ADS)

Wu, Chensheng; Davis, Christopher C.

2013-09-01

2D phase screen methods have been frequently applied to estimate atmospheric turbulence in free space optic communication and imaging systems. In situations where turbulence is "strong" enough to cause severe discontinuity of the wavefront (small Fried coherence length), the transmitted optic signal behaves more like "rays" rather than "waves". However, to achieve accurate simulation results through ray modeling requires both a high density of rays and a large number of eddies. Moreover, their complicated interactions require significant computational resources. Thus, we introduce a 3D ray model based on simple characteristics of turbulent eddies regardless of their particular geometry. The observed breakup of a beam wave into patches at a receiver and the theoretical description indicates that rays passing through the same sequence of turbulent eddies show "group" behavior whose wavefront can still be regarded as continuous. Thus, in our approach, we have divided the curved trajectory of rays into finite line segments and intuitively related their redirections to the refractive property of large turbulent eddies. As a result, our proposed treatment gives a quick and effective high-density ray simulation of a turbulent channel which only requires knowledge of the magnitude of the refractive index deviations. And our method points out a potential correction in reducing equivalent Cn2 by applying adaptive optics. This treatment also shows the possibility of extending 2D phase screen simulations into more general 3D treatments.

Simulations of eddy kinetic energy transport in barotropic turbulence

NASA Astrophysics Data System (ADS)

Grooms, Ian

2017-11-01

Eddy energy transport in rotating two-dimensional turbulence is investigated using numerical simulation. Stochastic forcing is used to generate an inhomogeneous field of turbulence and the time-mean energy profile is diagnosed. An advective-diffusive model for the transport is fit to the simulation data by requiring the model to accurately predict the observed time-mean energy distribution. Isotropic harmonic diffusion of energy is found to be an accurate model in the case of uniform, solid-body background rotation (the f plane), with a diffusivity that scales reasonably well with a mixing-length law κ ∝V ℓ , where V and ℓ are characteristic eddy velocity and length scales. Passive tracer dynamics are added and it is found that the energy diffusivity is 75 % of the tracer diffusivity. The addition of a differential background rotation with constant vorticity gradient β leads to significant changes to the energy transport. The eddies generate and interact with a mean flow that advects the eddy energy. Mean advection plus anisotropic diffusion (with reduced diffusivity in the direction of the background vorticity gradient) is moderately accurate for flows with scale separation between the eddies and mean flow, but anisotropic diffusion becomes a much less accurate model of the transport when scale separation breaks down. Finally, it is observed that the time-mean eddy energy does not look like the actual eddy energy distribution at any instant of time. In the future, stochastic models of the eddy energy transport may prove more useful than models of the mean transport for predicting realistic eddy energy distributions.

A High-Sensitivity Flexible Eddy Current Array Sensor for Crack Monitoring of Welded Structures under Varying Environment.

PubMed

Chen, Tao; He, Yuting; Du, Jinqiang

2018-06-01

This paper develops a high-sensitivity flexible eddy current array (HS-FECA) sensor for crack monitoring of welded structures under varying environment. Firstly, effects of stress, temperature and crack on output signals of the traditional flexible eddy current array (FECA) sensor were investigated by experiments that show both stress and temperature have great influences on the crack monitoring performance of the sensor. A 3-D finite element model was established using Comsol AC/DC module to analyze the perturbation effects of crack on eddy currents and output signals of the sensor, which showed perturbation effect of cracks on eddy currents is reduced by the current loop when crack propagates. Then, the HS-FECA sensor was proposed to boost the sensitivity to cracks. Simulation results show that perturbation effect of cracks on eddy currents excited by the HS-FECA sensor gradually grows stronger when the crack propagates, resulting in much higher sensitivity to cracks. Experimental result further shows that the sensitivity of the new sensor is at least 19 times that of the original one. In addition, both stress and temperature variations have little effect on signals of the new sensor.

Scale-Similar Models for Large-Eddy Simulations

NASA Technical Reports Server (NTRS)

Sarghini, F.

1999-01-01

Scale-similar models employ multiple filtering operations to identify the smallest resolved scales, which have been shown to be the most active in the interaction with the unresolved subgrid scales. They do not assume that the principal axes of the strain-rate tensor are aligned with those of the subgrid-scale stress (SGS) tensor, and allow the explicit calculation of the SGS energy. They can provide backscatter in a numerically stable and physically realistic manner, and predict SGS stresses in regions that are well correlated with the locations where large Reynolds stress occurs. In this paper, eddy viscosity and mixed models, which include an eddy-viscosity part as well as a scale-similar contribution, are applied to the simulation of two flows, a high Reynolds number plane channel flow, and a three-dimensional, nonequilibrium flow. The results show that simulations without models or with the Smagorinsky model are unable to predict nonequilibrium effects. Dynamic models provide an improvement of the results: the adjustment of the coefficient results in more accurate prediction of the perturbation from equilibrium. The Lagrangian-ensemble approach [Meneveau et al., J. Fluid Mech. 319, 353 (1996)] is found to be very beneficial. Models that included a scale-similar term and a dissipative one, as well as the Lagrangian ensemble averaging, gave results in the best agreement with the direct simulation and experimental data.

Large-eddy and unsteady RANS simulations of a shock-accelerated heavy gas cylinder

DOE PAGES

Morgan, B. E.; Greenough, J. A.

2015-04-08

Two-dimensional numerical simulations of the Richtmyer–Meshkov unstable “shock-jet” problem are conducted using both large-eddy simulation (LES) and unsteady Reynolds-averaged Navier–Stokes (URANS) approaches in an arbitrary Lagrangian–Eulerian hydrodynamics code. Turbulence statistics are extracted from LES by running an ensemble of simulations with multimode perturbations to the initial conditions. Detailed grid convergence studies are conducted, and LES results are found to agree well with both experiment and high-order simulations conducted by Shankar et al. (Phys Fluids 23, 024102, 2011). URANS results using a k–L approach are found to be highly sensitive to initialization of the turbulence lengthscale L and to the timemore » at which L becomes resolved on the computational mesh. As a result, it is observed that a gradient diffusion closure for turbulent species flux is a poor approximation at early times, and a new closure based on the mass-flux velocity is proposed for low-Reynolds-number mixing.« less

Nonlinear Eddy-Eddy Interactions in Dry Atmospheres Macroturbulence

NASA Astrophysics Data System (ADS)

Ait Chaalal, F.; Schneider, T.

2012-12-01

The statistical moment equations derived from the atmospheric equation of motions are not closed. However neglecting the large-scale eddy-eddy nonlinear interactions in an idealized dry general circulation model (GCM), which is equivalent to truncating the moment equations at the second order, can reproduce some of the features of the general circulation ([1]), highlighting the significance of eddy-mean flow interactions and the weakness of eddy-eddy interactions in atmospheric macroturbulence ([2]). The goal of the present study is to provide new insight into the rôle of these eddy-eddy interactions and discuss the relevance of a simple stochastic parametrization to represent them. We investigate in detail the general circulation in an idealized dry GCM, comparing full simulations with simulations where the eddy-eddy interactions are removed. The radiative processes are parametrized through Newtonian relaxation toward a radiative-equilibrium state with a prescribed equator to pole temperature contrast. A convection scheme relaxing toward a prescribed convective vertical lapse rate mimics some aspects of moist convection. The study is performed over a wide range of parameters covering the planetary rotation rate, the equator to pole temperature contrast and the vertical lapse rate. Particular attention is given to the wave-mean flow interactions and to the spectral budget. It is found that the no eddy-eddy simulations perform well when the baroclinic activity is weaker, for example for lower equator to pole temperature contrasts or higher rotation rates: the mean meridional circulation is well reproduced, with realistic eddy-driven jets and energy-containing eddy length scales of the order of the Rossby deformation radius. For a stronger baroclinic activity the no eddy-eddy model does not achieve a realistic isotropization of the eddies, the meridional circulation is compressed in the meridional direction and secondary eddy-driven jets emerge. In addition, the baroclinic wave activity does not reach the upper troposphere in association with a very weak or absent Rossby wave absorption in the upper subtropical troposphere. Understanding these deficiencies and the rôle of the eddy-eddy nonlinear interactions in determining the mean meridional circulation paves the way to the development of stochastic third order moments parametrizations, to eventually build GCMs that directly solve for the flow statistics and that could provide a deeper understanding of anthropogenic and natural climate changes. [1] O'Gorman, P. A., & Schneider, T. 2007, Geophysical Research Letters, 34, 22801 [2] Schneider, T., and C. C. Walker, 2006, Journal of the Atmospheric Sciences, 63, 1569-1586.

Simulating carbon and water fluxes at Arctic and boreal ecosystems in Alaska by optimizing the modified BIOME-BGC with eddy covariance data

NASA Astrophysics Data System (ADS)

Ueyama, M.; Kondo, M.; Ichii, K.; Iwata, H.; Euskirchen, E. S.; Zona, D.; Rocha, A. V.; Harazono, Y.; Nakai, T.; Oechel, W. C.

2013-12-01

To better predict carbon and water cycles in Arctic ecosystems, we modified a process-based ecosystem model, BIOME-BGC, by introducing new processes: change in active layer depth on permafrost and phenology of tundra vegetation. The modified BIOME-BGC was optimized using an optimization method. The model was constrained using gross primary productivity (GPP) and net ecosystem exchange (NEE) at 23 eddy covariance sites in Alaska, and vegetation/soil carbon from a literature survey. The model was used to simulate regional carbon and water fluxes of Alaska from 1900 to 2011. Simulated regional fluxes were validated with upscaled GPP, ecosystem respiration (RE), and NEE based on two methods: (1) a machine learning technique and (2) a top-down model. Our initial simulation suggests that the original BIOME-BGC with default ecophysiological parameters substantially underestimated GPP and RE for tundra and overestimated those fluxes for boreal forests. We will discuss how optimization using the eddy covariance data impacts the historical simulation by comparing the new version of the model with simulated results from the original BIOME-BGC with default ecophysiological parameters. This suggests that the incorporation of the active layer depth and plant phenology processes is important to include when simulating carbon and water fluxes in Arctic ecosystems.

Improved Climate Simulations through a Stochastic Parameterization of Ocean Eddies

NASA Astrophysics Data System (ADS)

Williams, Paul; Howe, Nicola; Gregory, Jonathan; Smith, Robin; Joshi, Manoj

2017-04-01

In climate simulations, the impacts of the subgrid scales on the resolved scales are conventionally represented using deterministic closure schemes, which assume that the impacts are uniquely determined by the resolved scales. Stochastic parameterization relaxes this assumption, by sampling the subgrid variability in a computationally inexpensive manner. This study shows that the simulated climatological state of the ocean is improved in many respects by implementing a simple stochastic parameterization of ocean eddies into a coupled atmosphere-ocean general circulation model. Simulations from a high-resolution, eddy-permitting ocean model are used to calculate the eddy statistics needed to inject realistic stochastic noise into a low-resolution, non-eddy-permitting version of the same model. A suite of four stochastic experiments is then run to test the sensitivity of the simulated climate to the noise definition by varying the noise amplitude and decorrelation time within reasonable limits. The addition of zero-mean noise to the ocean temperature tendency is found to have a nonzero effect on the mean climate. Specifically, in terms of the ocean temperature and salinity fields both at the surface and at depth, the noise reduces many of the biases in the low-resolution model and causes it to more closely resemble the high-resolution model. The variability of the strength of the global ocean thermohaline circulation is also improved. It is concluded that stochastic ocean perturbations can yield reductions in climate model error that are comparable to those obtained by refining the resolution, but without the increased computational cost. Therefore, stochastic parameterizations of ocean eddies have the potential to significantly improve climate simulations. Reference Williams PD, Howe NJ, Gregory JM, Smith RS, and Joshi MM (2016) Improved Climate Simulations through a Stochastic Parameterization of Ocean Eddies. Journal of Climate, 29, 8763-8781. http://dx.doi.org/10.1175/JCLI-D-15-0746.1

An extended algebraic variational multiscale-multigrid-multifractal method (XAVM4) for large-eddy simulation of turbulent two-phase flow

NASA Astrophysics Data System (ADS)

Rasthofer, U.; Wall, W. A.; Gravemeier, V.

2018-04-01

A novel and comprehensive computational method, referred to as the eXtended Algebraic Variational Multiscale-Multigrid-Multifractal Method (XAVM4), is proposed for large-eddy simulation of the particularly challenging problem of turbulent two-phase flow. The XAVM4 involves multifractal subgrid-scale modeling as well as a Nitsche-type extended finite element method as an approach for two-phase flow. The application of an advanced structural subgrid-scale modeling approach in conjunction with a sharp representation of the discontinuities at the interface between two bulk fluids promise high-fidelity large-eddy simulation of turbulent two-phase flow. The high potential of the XAVM4 is demonstrated for large-eddy simulation of turbulent two-phase bubbly channel flow, that is, turbulent channel flow carrying a single large bubble of the size of the channel half-width in this particular application.

Large Eddy Simulation of Transitional Flow in an Idealized Stenotic Blood Vessel: Evaluation of Subgrid Scale Models

PubMed Central

Pal, Abhro; Anupindi, Kameswararao; Delorme, Yann; Ghaisas, Niranjan; Shetty, Dinesh A.; Frankel, Steven H.

2014-01-01

In the present study, we performed large eddy simulation (LES) of axisymmetric, and 75% stenosed, eccentric arterial models with steady inflow conditions at a Reynolds number of 1000. The results obtained are compared with the direct numerical simulation (DNS) data (Varghese et al., 2007, “Direct Numerical Simulation of Stenotic Flows. Part 1. Steady Flow,” J. Fluid Mech., 582, pp. 253–280). An inhouse code (WenoHemo) employing high-order numerical methods for spatial and temporal terms, along with a 2nd order accurate ghost point immersed boundary method (IBM) (Mark, and Vanwachem, 2008, “Derivation and Validation of a Novel Implicit Second-Order Accurate Immersed Boundary Method,” J. Comput. Phys., 227(13), pp. 6660–6680) for enforcing boundary conditions on curved geometries is used for simulations. Three subgrid scale (SGS) models, namely, the classical Smagorinsky model (Smagorinsky, 1963, “General Circulation Experiments With the Primitive Equations,” Mon. Weather Rev., 91(10), pp. 99–164), recently developed Vreman model (Vreman, 2004, “An Eddy-Viscosity Subgrid-Scale Model for Turbulent Shear Flow: Algebraic Theory and Applications,” Phys. Fluids, 16(10), pp. 3670–3681), and the Sigma model (Nicoud et al., 2011, “Using Singular Values to Build a Subgrid-Scale Model for Large Eddy Simulations,” Phys. Fluids, 23(8), 085106) are evaluated in the present study. Evaluation of SGS models suggests that the classical constant coefficient Smagorinsky model gives best agreement with the DNS data, whereas the Vreman and Sigma models predict an early transition to turbulence in the poststenotic region. Supplementary simulations are performed using Open source field operation and manipulation (OpenFOAM) (“OpenFOAM,” http://www.openfoam.org/) solver and the results are inline with those obtained with WenoHemo. PMID:24801556

Numerical modeling of turbulent swirling flow in a multi-inlet vortex nanoprecipitation reactor using dynamic DDES

NASA Astrophysics Data System (ADS)

Hill, James C.; Liu, Zhenping; Fox, Rodney O.; Passalacqua, Alberto; Olsen, Michael G.

2015-11-01

The multi-inlet vortex reactor (MIVR) has been developed to provide a platform for rapid mixing in the application of flash nanoprecipitation (FNP) for manufacturing functional nanoparticles. Unfortunately, commonly used RANS methods are unable to accurately model this complex swirling flow. Large eddy simulations have also been problematic, as expensive fine grids to accurately model the flow are required. These dilemmas led to the strategy of applying a Delayed Detached Eddy Simulation (DDES) method to the vortex reactor. In the current work, the turbulent swirling flow inside a scaled-up MIVR has been investigated by using a dynamic DDES model. In the DDES model, the eddy viscosity has a form similar to the Smagorinsky sub-grid viscosity in LES and allows the implementation of a dynamic procedure to determine its coefficient. The complex recirculating back flow near the reactor center has been successfully captured by using this dynamic DDES model. Moreover, the simulation results are found to agree with experimental data for mean velocity and Reynolds stresses.

Detached Eddy Simulation for the F-16XL Aircraft Configuration

NASA Technical Reports Server (NTRS)

Elmiligui, Alaa; Abdol-Hamid, Khaled; Parlette, Edward B.

2015-01-01

Numerical simulations for the flow around the F-16XL configuration as a contribution to the Cranked Arrow Wing Aerodynamic Project International 2 (CAWAPI-2) have been performed. The NASA Langley Tetrahedral Unstructured Software System (TetrUSS) with its USM3D solver was used to perform the unsteady flow field simulations for the subsonic high angle-of-attack case corresponding to flight condition (FC) 25. Two approaches were utilized to capture the unsteady vortex flow over the wing of the F-16XL. The first approach was to use Unsteady Reynolds-Averaged Navier-Stokes (URANS) coupled with standard turbulence closure models. The second approach was to use Detached Eddy Simulation (DES), which creates a hybrid model that attempts to combine the most favorable elements of URANS models and Large Eddy Simulation (LES). Computed surface static pressure profiles are presented and compared with flight data. Time-averaged and instantaneous results obtained on coarse, medium and fine grids are compared with the flight data. The intent of this study is to demonstrate that the DES module within the USM3D solver can be used to provide valuable data in predicting vortex-flow physics on a complex configuration.

Evaluation of subgrid-scale models in large-eddy simulations of turbulent flow in a centrifugal pump impeller

NASA Astrophysics Data System (ADS)

Yang, Zhengjun; Wang, Fujun; Zhou, Peijian

2012-09-01

The current research of large eddy simulation (LES) of turbulent flow in pumps mainly concentrates in applying conventional subgrid-scale (SGS) model to simulate turbulent flow, which aims at obtaining the flow field in pump. The selection of SGS model is usually not considered seriously, so the accuracy and efficiency of the simulation cannot be ensured. Three SGS models including Smagorinsky-Lilly model, dynamic Smagorinsky model and dynamic mixed model are comparably studied by using the commercial CFD code Fluent combined with its user define function. The simulations are performed for the turbulent flow in a centrifugal pump impeller. The simulation results indicate that the mean flows predicted by the three SGS models agree well with the experimental data obtained from the test that detailed measurements of the flow inside the rotating passages of a six-bladed shrouded centrifugal pump impeller performed using particle image velocimetry (PIV) and laser Doppler velocimetry (LDV). The comparable results show that dynamic mixed model gives the most accurate results for mean flow in the centrifugal pump impeller. The SGS stress of dynamic mixed model is decompose into the scale similar part and the eddy viscous part. The scale similar part of SGS stress plays a significant role in high curvature regions, such as the leading edge and training edge of pump blade. It is also found that the dynamic mixed model is more adaptive to compute turbulence in the pump impeller. The research results presented is useful to improve the computational accuracy and efficiency of LES for centrifugal pumps, and provide important reference for carrying out simulation in similar fluid machineries.

Multi-model analysis of terrestrial carbon cycles in Japan: limitations and implications of model calibration using eddy flux observations

NASA Astrophysics Data System (ADS)

Ichii, K.; Suzuki, T.; Kato, T.; Ito, A.; Hajima, T.; Ueyama, M.; Sasai, T.; Hirata, R.; Saigusa, N.; Ohtani, Y.; Takagi, K.

2010-07-01

Terrestrial biosphere models show large differences when simulating carbon and water cycles, and reducing these differences is a priority for developing more accurate estimates of the condition of terrestrial ecosystems and future climate change. To reduce uncertainties and improve the understanding of their carbon budgets, we investigated the utility of the eddy flux datasets to improve model simulations and reduce variabilities among multi-model outputs of terrestrial biosphere models in Japan. Using 9 terrestrial biosphere models (Support Vector Machine - based regressions, TOPS, CASA, VISIT, Biome-BGC, DAYCENT, SEIB, LPJ, and TRIFFID), we conducted two simulations: (1) point simulations at four eddy flux sites in Japan and (2) spatial simulations for Japan with a default model (based on original settings) and a modified model (based on model parameter tuning using eddy flux data). Generally, models using default model settings showed large deviations in model outputs from observation with large model-by-model variability. However, after we calibrated the model parameters using eddy flux data (GPP, RE and NEP), most models successfully simulated seasonal variations in the carbon cycle, with less variability among models. We also found that interannual variations in the carbon cycle are mostly consistent among models and observations. Spatial analysis also showed a large reduction in the variability among model outputs. This study demonstrated that careful validation and calibration of models with available eddy flux data reduced model-by-model differences. Yet, site history, analysis of model structure changes, and more objective procedure of model calibration should be included in the further analysis.

Local dynamic subgrid-scale models in channel flow

NASA Technical Reports Server (NTRS)

Cabot, William H.

1994-01-01

The dynamic subgrid-scale (SGS) model has given good results in the large-eddy simulation (LES) of homogeneous isotropic or shear flow, and in the LES of channel flow, using averaging in two or three homogeneous directions (the DA model). In order to simulate flows in general, complex geometries (with few or no homogeneous directions), the dynamic SGS model needs to be applied at a local level in a numerically stable way. Channel flow, which is inhomogeneous and wall-bounded flow in only one direction, provides a good initial test for local SGS models. Tests of the dynamic localization model were performed previously in channel flow using a pseudospectral code and good results were obtained. Numerical instability due to persistently negative eddy viscosity was avoided by either constraining the eddy viscosity to be positive or by limiting the time that eddy viscosities could remain negative by co-evolving the SGS kinetic energy (the DLk model). The DLk model, however, was too expensive to run in the pseudospectral code due to a large near-wall term in the auxiliary SGS kinetic energy (k) equation. One objective was then to implement the DLk model in a second-order central finite difference channel code, in which the auxiliary k equation could be integrated implicitly in time at great reduction in cost, and to assess its performance in comparison with the plane-averaged dynamic model or with no model at all, and with direct numerical simulation (DNS) and/or experimental data. Other local dynamic SGS models have been proposed recently, e.g., constrained dynamic models with random backscatter, and with eddy viscosity terms that are averaged in time over material path lines rather than in space. Another objective was to incorporate and test these models in channel flow.

On the large eddy simulation of turbulent flows in complex geometry

NASA Technical Reports Server (NTRS)

Ghosal, Sandip

1993-01-01

Application of the method of Large Eddy Simulation (LES) to a turbulent flow consists of three separate steps. First, a filtering operation is performed on the Navier-Stokes equations to remove the small spatial scales. The resulting equations that describe the space time evolution of the 'large eddies' contain the subgrid-scale (sgs) stress tensor that describes the effect of the unresolved small scales on the resolved scales. The second step is the replacement of the sgs stress tensor by some expression involving the large scales - this is the problem of 'subgrid-scale modeling'. The final step is the numerical simulation of the resulting 'closed' equations for the large scale fields on a grid small enough to resolve the smallest of the large eddies, but still much larger than the fine scale structures at the Kolmogorov length. In dividing a turbulent flow field into 'large' and 'small' eddies, one presumes that a cut-off length delta can be sensibly chosen such that all fluctuations on a scale larger than delta are 'large eddies' and the remainder constitute the 'small scale' fluctuations. Typically, delta would be a length scale characterizing the smallest structures of interest in the flow. In an inhomogeneous flow, the 'sensible choice' for delta may vary significantly over the flow domain. For example, in a wall bounded turbulent flow, most statistical averages of interest vary much more rapidly with position near the wall than far away from it. Further, there are dynamically important organized structures near the wall on a scale much smaller than the boundary layer thickness. Therefore, the minimum size of eddies that need to be resolved is smaller near the wall. In general, for the LES of inhomogeneous flows, the width of the filtering kernel delta must be considered to be a function of position. If a filtering operation with a nonuniform filter width is performed on the Navier-Stokes equations, one does not in general get the standard large eddy equations. The complication is caused by the fact that a filtering operation with a nonuniform filter width in general does not commute with the operation of differentiation. This is one of the issues that we have looked at in detail as it is basic to any attempt at applying LES to complex geometry flows. Our principal findings are summarized.

Observed and Simulated Eddy Diffusivity Upstream of the Drake Passage

NASA Astrophysics Data System (ADS)

Tulloch, R.; Ferrari, R. M.; Marshall, J.

2012-12-01

Estimates of eddy diffusivity in the Southern Ocean are poorly constrained due to lack of observations. We compare the first direct estimate of isopycnal eddy diffusivity upstream of the Drake Passage (from Ledwell et al. 2011) with a numerical simulation. The estimate is computed from a point tracer release as part of the Diapycnal and Isopycnal Mixing Experiment in the Southern Ocean (DIMES). We find that the observational diffusivity estimate of about 500m^2/s at 1500m depth is close to that computed in a data-constrained, 1/20th of a degree simulation of the Drake Passage region. This tracer estimate also agrees with Lagrangian float calculations in the model. The role of mean flow suppression of eddy diffusivity at shallower depths will also be discussed.

Seasonal variability in global eddy diffusion and the effect on neutral density

NASA Astrophysics Data System (ADS)

Pilinski, M. D.; Crowley, G.

2015-04-01

We describe a method for making single-satellite estimates of the seasonal variability in global-average eddy diffusion coefficients. Eddy diffusion values as a function of time were estimated from residuals of neutral density measurements made by the Challenging Minisatellite Payload (CHAMP) and simulations made using the thermosphere-ionosphere-mesosphere electrodynamics global circulation model (TIME-GCM). The eddy diffusion coefficient results are quantitatively consistent with previous estimates based on satellite drag observations and are qualitatively consistent with other measurement methods such as sodium lidar observations and eddy diffusivity models. Eddy diffusion coefficient values estimated between January 2004 and January 2008 were then used to generate new TIME-GCM results. Based on these results, the root-mean-square sum for the TIME-GCM model is reduced by an average of 5% when compared to density data from a variety of satellites, indicating that the fidelity of global density modeling can be improved by using data from a single satellite like CHAMP. This approach also demonstrates that eddy diffusion could be estimated in near real-time from satellite observations and used to drive a global circulation model like TIME-GCM. Although the use of global values improves modeled neutral densities, there are limitations to this method, which are discussed, including that the latitude dependence of the seasonal neutral-density signal is not completely captured by a global variation of eddy diffusion coefficients. This demonstrates the need for a latitude-dependent specification of eddy diffusion which is also consistent with diffusion observations made by other techniques.

Comparison of three large-eddy simulations of shock-induced turbulent separation bubbles

NASA Astrophysics Data System (ADS)

Touber, Emile; Sandham, Neil D.

2009-12-01

Three different large-eddy simulation investigations of the interaction between an impinging oblique shock and a supersonic turbulent boundary layer are presented. All simulations made use of the same inflow technique, specifically aimed at avoiding possible low-frequency interferences with the shock/boundary-layer interaction system. All simulations were run on relatively wide computational domains and integrated over times greater than twenty five times the period of the most commonly reported low-frequency shock-oscillation, making comparisons at both time-averaged and low-frequency-dynamic levels possible. The results confirm previous experimental results which suggested a simple linear relation between the interaction length and the oblique-shock strength if scaled using the boundary-layer thickness and wall-shear stress. All the tested cases show evidences of significant low-frequency shock motions. At the wall, energetic low-frequency pressure fluctuations are observed, mainly in the initial part of interaction.

Detached Eddy Simulation Results for a Space Launch System Configuration at Liftoff Conditions and Comparison with Experiment

NASA Technical Reports Server (NTRS)

Krist, Steven E.; Ghaffari, Farhad

2015-01-01

Computational simulations for a Space Launch System configuration at liftoff conditions for incidence angles from 0 to 90 degrees were conducted in order to generate integrated force and moment data and longitudinal lineloads. While the integrated force and moment coefficients can be obtained from wind tunnel testing, computational analyses are indispensable in obtaining the extensive amount of surface information required to generate proper lineloads. However, beyond an incidence angle of about 15 degrees, the effects of massive flow separation on the leeward pressure field is not well captured with state of the art Reynolds Averaged Navier-Stokes methods, necessitating the employment of a Detached Eddy Simulation method. Results from these simulations are compared to the liftoff force and moment database and surface pressure data derived from a test in the NASA Langley 14- by 22-Foot Subsonic Wind Tunnel.

Analysis of the Influence of Construction Insulation Systems on Public Safety in China

PubMed Central

Zhang, Guowei; Zhu, Guoqing; Zhao, Guoxiang

2016-01-01

With the Government of China’s proposed Energy Efficiency Regulations (GB40411-2007), the implementation of external insulation systems will be mandatory in China. The frequent external insulation system fires cause huge numbers of casualties and extensive property damage and have rapidly become a new hot issue in construction evacuation safety in China. This study attempts to reconstruct an actual fire scene and propose a quantitative risk assessment method for upward insulation system fires using thermal analysis tests and large eddy simulations (using the Fire Dynamics Simulator (FDS) software). Firstly, the pyrolysis and combustion characteristics of Extruded polystyrene board (XPS panel), such as ignition temperature, combustion heat, limiting oxygen index, thermogravimetric analysis and thermal radiation analysis were studied experimentally. Based on these experimental data, large eddy simulation was then applied to reconstruct insulation system fires. The results show that upward insulation system fires could be accurately reconstructed by using thermal analysis test and large eddy simulation. The spread of insulation material system fires in the vertical direction is faster than that in the horizontal direction. Moreover, we also find that there is a possibility of flashover in enclosures caused by insulation system fires as the smoke temperature exceeds 600 °C. The simulation methods and experimental results obtained in this paper could provide valuable references for fire evacuation, hazard assessment and fire resistant construction design studies. PMID:27589774

Analysis of the Influence of Construction Insulation Systems on Public Safety in China.

PubMed

Zhang, Guowei; Zhu, Guoqing; Zhao, Guoxiang

2016-08-30

With the Government of China's proposed Energy Efficiency Regulations (GB40411-2007), the implementation of external insulation systems will be mandatory in China. The frequent external insulation system fires cause huge numbers of casualties and extensive property damage and have rapidly become a new hot issue in construction evacuation safety in China. This study attempts to reconstruct an actual fire scene and propose a quantitative risk assessment method for upward insulation system fires using thermal analysis tests and large eddy simulations (using the Fire Dynamics Simulator (FDS) software). Firstly, the pyrolysis and combustion characteristics of Extruded polystyrene board (XPS panel), such as ignition temperature, combustion heat, limiting oxygen index, thermogravimetric analysis and thermal radiation analysis were studied experimentally. Based on these experimental data, large eddy simulation was then applied to reconstruct insulation system fires. The results show that upward insulation system fires could be accurately reconstructed by using thermal analysis test and large eddy simulation. The spread of insulation material system fires in the vertical direction is faster than that in the horizontal direction. Moreover, we also find that there is a possibility of flashover in enclosures caused by insulation system fires as the smoke temperature exceeds 600 °C. The simulation methods and experimental results obtained in this paper could provide valuable references for fire evacuation, hazard assessment and fire resistant construction design studies.

Dynamical analysis of a satellite-observed anticyclonic eddy in the northern Bering Sea

NASA Astrophysics Data System (ADS)

Li, Yineng; Li, Xiaofeng; Wang, Jia; Peng, Shiqiu

2016-05-01

The characteristics and evolution of a satellite-observed anticyclonic eddy in the northern Bering Sea during March and April 1999 are investigated using a three-dimensional Princeton Ocean Model (POM). The anticyclonic-like current pattern and asymmetric feature of the eddy were clearly seen in the synthetic aperture radar (SAR), sea surface temperature, and ocean color images in April 1999. The results from model simulation reveal the three-dimensional structure of the anticyclonic eddy, its movement, and dissipation. Energy analysis indicates that the barotropic instability (BTI) is the main energy source for the growth of the anticyclonic eddy. The momentum analysis further reveals that the larger magnitude of the barotropic pressure gradient in the meridional direction causes the asymmetry of the anticyclonic eddy in the zonal and meridional directions, while the different magnitudes of the meridional baroclinic pressure gradient are responsible for the different intensity of currents between the northern and southern parts of the anticyclonic eddy. This article was corrected on 23 JUL 2016. See the end of the full text for details.

Large Eddy Simulation of a Film Cooling Technique with a Plenum

NASA Astrophysics Data System (ADS)

Dharmarathne, Suranga; Sridhar, Narendran; Araya, Guillermo; Castillo, Luciano; Parameswaran, Sivapathasund

2012-11-01

Factors that affect the film cooling performance have been categorized into three main groups: (i) coolant & mainstream conditions, (ii) hole geometry & configuration, and (iii) airfoil geometry Bogard et al. (2006). The present study focuses on the second group of factors, namely, the modeling of coolant hole and the plenum. It is required to simulate correct physics of the problem to achieve more realistic numerical results. In this regard, modeling of cooling jet hole and the plenum chamber is highly important Iourokina et al. (2006). Substitution of artificial boundary conditions instead of correct plenum design would yield unrealistic results Iourokina et al. (2006). This study attempts to model film cooling technique with a plenum using a Large Eddy Simulation.Incompressible coolant jet ejects to the surface of the plate at an angle of 30° where it meets compressible turbulent boundary layer which simulates the turbine inflow conditions. Dynamic multi-scale approach Araya (2011) is introduced to prescribe turbulent inflow conditions. Simulations are carried out for two different blowing ratios and film cooling effectiveness is calculated for both cases. Results obtained from LES will be compared with experimental results.

Large-eddy simulations of compressible convection on massively parallel computers. [stellar physics

NASA Technical Reports Server (NTRS)

Xie, Xin; Toomre, Juri

1993-01-01

We report preliminary implementation of the large-eddy simulation (LES) technique in 2D simulations of compressible convection carried out on the CM-2 massively parallel computer. The convective flow fields in our simulations possess structures similar to those found in a number of direct simulations, with roll-like flows coherent across the entire depth of the layer that spans several density scale heights. Our detailed assessment of the effects of various subgrid scale (SGS) terms reveals that they may affect the gross character of convection. Yet, somewhat surprisingly, we find that our LES solutions, and another in which the SGS terms are turned off, only show modest differences. The resulting 2D flows realized here are rather laminar in character, and achieving substantial turbulence may require stronger forcing and less dissipation.

Asymmetric lake distribution on Titan mediated by methane transport due to atmospheric eddies

NASA Astrophysics Data System (ADS)

Lora, Juan M.; Mitchell, Jonathan L.

2015-11-01

The observed north-south asymmetry in the distribution of Titan's seas and lakes has been proposed to be a consequence of orbital forcing affecting Titan's hydrologic cycle, as in the present the northern summer is longer but milder than its southern counterpart. Though recent general circulation models have simulated asymmetrical surface liquid distributions, the mechanism that generates this asymmetry has not been explained. In this work, we compare axisymmetric and three-dimensional simulations of Titan's atmospheric circulation with the Titan Atmospheric Model (TAM) [Lora et al. 2015, Icarus 250] to investigate the transport of moisture by the atmosphere. A significant hemispheric asymmetry only develops in the latter case, and we demonstrate that equatorward transport by high-latitude, baroclinic eddies is responsible. Eddies transport moisture from the high latitudes into the low and midlatitude cross-equatorial mean meridional circulation, producing an atmospheric "bucket brigade." The moisture transport by eddies is more intense in the south than in the north as a consequence of the orbital forcing, and therefore the result is net northward transport of methane, explaining the surface buildup in the north.

Parameterization of eddy sensible heat transports in a zonally averaged dynamic model of the atmosphere

NASA Technical Reports Server (NTRS)

Genthon, Christophe; Le Treut, Herve; Sadourny, Robert; Jouzel, Jean

1990-01-01

A Charney-Branscome based parameterization has been tested as a way of representing the eddy sensible heat transports missing in a zonally averaged dynamic model (ZADM) of the atmosphere. The ZADM used is a zonally averaged version of a general circulation model (GCM). The parameterized transports in the ZADM are gaged against the corresponding fluxes explicitly simulated in the GCM, using the same zonally averaged boundary conditions in both models. The Charney-Branscome approach neglects stationary eddies and transient barotropic disturbances and relies on a set of simplifying assumptions, including the linear appoximation, to describe growing transient baroclinic eddies. Nevertheless, fairly satisfactory results are obtained when the parameterization is performed interactively with the model. Compared with noninteractive tests, a very efficient restoring feedback effect between the modeled zonal-mean climate and the parameterized meridional eddy transport is identified.

Large-eddy simulation of the urban boundary layer in the MEGAPOLI Paris Plume experiment

NASA Astrophysics Data System (ADS)

Esau, Igor

2010-05-01

This study presents results from the specific large-eddy simulation study of the urban boundary layer in the MEGAPOLI Paris Plume field campaign. We used LESNIC and PALM codes, MEGAPOLI city morphology database, nudging to the observed meteorological conditions during the Paris Plume campaign and some concentration measurements from that campaign to simulate and better understand the nature of the urban boundary layer on scales larger then the street canyon scales. The primary attention was paid to turbulence self-organization and structure-to-surface interaction. The study has been aimed to demonstrate feasibility and estimate required resources for such research. Therefore, at this stage we do not compare the simulation with other relevant studies as well as we do not formulate the theoretical conclusions.

Detailed characteristics of drop-laden mixing layers: LES predictions compared to DNS

NASA Technical Reports Server (NTRS)

Okong'o, N.; Leboissetier, A.; Bellan, J.

2004-01-01

Results have been compared from Direct Numerical Simulation (DNS) and Large Eddy Simulation (LES) of a temporal mixing layer laden with evaporating drops, to assess the ability of LES to reproduce detailed characteristics of DNS.

Error-growth dynamics and predictability of surface thermally induced atmospheric flow

DOE Office of Scientific and Technical Information (OSTI.GOV)

Zeng, X.; Pielke, R.A.

1993-09-01

Using the CSU Regional Atmospheric Modeling System (RAMS) in its nonhydrostatic and compressible configuration, over 200 two-dimensional simulations with [Delta]x = 2 km and [Delta]x = 100 m are performed to study in detail the initial adjustment process and the error-growth dynamics of surface thermally induced circulation including the sensitivity to initial conditions, boundary conditions, and model parameters, and to study the predictability as a function of the size of surface heat patches under a calm mean wind. It is found that the error growth is not sensitive to the characterisitics of the initial perturbations. The numerical smoothing has amore » strong impact on the initial adjustment process and on the error-growth dynamics. The predictability and flow structures, it is found that the vertical velocity field is strongly affected by the mean wind, and the flow structures are quite sensitive to the initial soil water content. The transition from organized flow to the situation in which fluxes are dominated by noncoherent turbulent eddies under a calm mean wind is quantitatively evaluated and this transition is different for different variables. The relationship between the predictability of a realization and of an ensemble average is discussed. The predictability and the coherent circulations modulated by the surface inhomogeneities are also studied by computing the autocorrelations and the power spectra. The three-dimensional mesoscale and large-eddy simulations are performed to verify the above results. It is found that the two-dimensional mesoscale (or fine resolution) simulation yields very close or similar results regarding the predictability as those from the three-dimensional mesoscale (or large eddy) simulation. The horizontally averaged quantities based on two-dimensional fine-resolution simulations are insensitive to initial perturbations and agree with those based on three-dimensional large-eddy simulations. 87 refs., 25 figs.« less

Subgrid-scale models for large-eddy simulation of rotating turbulent channel flows

NASA Astrophysics Data System (ADS)

Silvis, Maurits H.; Bae, Hyunji Jane; Trias, F. Xavier; Abkar, Mahdi; Moin, Parviz; Verstappen, Roel

2017-11-01

We aim to design subgrid-scale models for large-eddy simulation of rotating turbulent flows. Rotating turbulent flows form a challenging test case for large-eddy simulation due to the presence of the Coriolis force. The Coriolis force conserves the total kinetic energy while transporting it from small to large scales of motion, leading to the formation of large-scale anisotropic flow structures. The Coriolis force may also cause partial flow laminarization and the occurrence of turbulent bursts. Many subgrid-scale models for large-eddy simulation are, however, primarily designed to parametrize the dissipative nature of turbulent flows, ignoring the specific characteristics of transport processes. We, therefore, propose a new subgrid-scale model that, in addition to the usual dissipative eddy viscosity term, contains a nondissipative nonlinear model term designed to capture transport processes, such as those due to rotation. We show that the addition of this nonlinear model term leads to improved predictions of the energy spectra of rotating homogeneous isotropic turbulence as well as of the Reynolds stress anisotropy in spanwise-rotating plane-channel flows. This work is financed by the Netherlands Organisation for Scientific Research (NWO) under Project Number 613.001.212.

Eddy-Kuroshio Interactions: Local and Remote Effects

NASA Astrophysics Data System (ADS)

Jan, Sen; Mensah, Vigan; Andres, Magdalena; Chang, Ming-Huei; Yang, Yiing Jang

2017-12-01

Quasi-geostrophic mesoscale eddies regularly impinge on the Kuroshio in the western North Pacific, but the processes underlying the evolution of these eddy-Kuroshio interactions have not yet been thoroughly investigated in the literature. Here this interaction is examined with results from a semi-idealized three-dimensional numerical model and observations from four pressure-sensor equipped inverted echo sounders (PIESs) in a zonal section east of Taiwan and satellite altimeters. Both the observations and numerical simulations suggest that, during the interaction of a cyclonic eddy with the Kuroshio, the circular eddy is deformed into an elliptic shape with the major axis in the northwest-southeast direction, before being dissipated; the poleward velocity and associated Kuroshio transport decrease and the sea level and pycnocline slopes across the Kuroshio weaken. In contrast, for an anticyclonic eddy during the eddy-Kuroshio interaction, variations in the velocity, sea level, and isopycnal depth are reversed; the circular eddy is also deformed to an ellipse but with the major axis parallel to the Kuroshio. The model results also demonstrate that the velocity field is modified first and consequently the SSH and isopycnal depth evolve during the interaction. Furthermore, due to the combined effect of impingement latitude and realistic topography, some eddy-Kuroshio interactions east of Taiwan are found to have remote effects, both in the Luzon Strait and on the East China Sea shelf northeast of Taiwan.