Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fixed grid”

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

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

At least 343 records · Page 19

A Study of Turbulence-Driven Inward Momentum Flux & Construction of a Difference Scheme for Fokker-Planck-Landau Operator

Momentum transport driven by microinstability is studied via analytic theories, gyrokinetic simulations, and experimental data analysis. A scenario for rotation profile evolution is constructed based on experimental observation that parallel flow initiates at edge region and penetrates into core region of tokamak [1]. In particular, turbulent diffusion and inward convective pinch are studied according to this scenario. For study of convective momentum pinch, toroidal ion temperature gradient mode near marginal stability is examined in kinetic limit. It is found from quasilinear estimation that the inward turbulent equipartition momentum pinch remains as the most robust part of pinch. In addition, it is found that ion temperature gradient driven momentum pinch is inward for typical parameters, while ion density gradient driven momentum pinch is outward. It is shown by theoretical formulas and experimental data analysis that momentum pinch velocity is sensitive to the electron to ion temperature ratio. For study of diffusive momentum transport, parallel shear flow instability is studied in fluid limit. Diffusion coefficient using mixing length estimate from nonlocal theory reveals that it has a strong dependence on the parallel flow shear and is mainly carried by fluctuations in short wavelength regime. Gyrokinetic simulation of collisionless parallel shear flow instability shows that the analytic prediction of diffusion coefficient is comparable to numerical values in the linear phase of the simulations. In the nonlinear phase, fluctuation energy gets transferred to the long wavelength regime. A two-dimensional, uniform velocity-grid based solver for the nonlinear Fokker-Planck-Landau (FPL) collision operator has been developed in the limit of strong magnetic field. The conservative and entropy increasing semidiscretized numerical scheme for the isotropic FPL equation by C. Buet et al. (2002) has been extended to 2D velocity space. A Picard fixed point method is used for semi-implicit time advance. One time evaluation of interpolation coefficients for probability distribution function is introduced so that numerical equilibrium is achieved with Maxwellian. However, positivity property is sacrificed at the saving of computational costs. In addition, a highly practical method is suggested to apply the velocity-grid based collision solver to particle-in-cell code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Numerical simulation of a powered-lift landing, tracking flow features using overset grids, and simulation of high lift devices on a fighter-lift-and-control wing

Attached as appendices to this report are documents describing work performed on the simulation of a landing powered-lift delta wing, the tracking of flow features using overset grids, and the simulation of flaps on the Wright Patterson Lab's fighter-lift-and-control (FLAC) wing. Numerical simulation of a powered-lift landing includes the computation of flow about a delta wing at four fixed heights as well as a simulated landing, in which the delta wing descends toward the ground. Comparison of computed and experimental lift coefficients indicates that the simulations capture the qualitative trends in lift-loss encountered by thrust-vectoring aircraft operating in ground effect. Power spectra of temporal variations of pressure indicate computed vortex shedding frequencies close to the jet exit are in the experimentally observed frequency range; the power spectra of pressure also provide insights into the mechanisms of lift oscillations. Also, a method for using overset grids to track dynamic flow features is described and the method is validated by tracking a moving shock and vortices shed behind a circular cylinder. Finally, Chimera gridding strategies were used to develop pressure coefficient contours for the FLAC wing for a Mach no. of 0.18 and Reynolds no. of 2.5 million.

Chawla, Kalpana↗

Effective HTCondor-based monitoring system for CMS

The CMS experiment at the LHC relies on HTCondor and glideinWMS as its primary batch and pilot-based Grid provisioning systems, respectively. Given the scale of the global queue in CMS, the operators found it increasingly difficult to monitor the pool to find problems and fix them. The operators had to rely on several different web pages, with several different levels of information, and sift tirelessly through log files in order to monitor the pool completely. Therefore, coming up with a suitable monitoring system was one of the crucial items before the beginning of the LHC Run 2 in order to ensure early detection of issues and to give a good overview of the whole pool. Our new monitoring page utilizes the HTCondor ClassAd information to provide a complete picture of the whole submission infrastructure in CMS. The monitoring page includes useful information from HTCondor schedulers, central managers, the glideinWMS frontend, and factories. It also incorporates information about users and tasks making it easy for operators to provide support and debug issues.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Towards Accurate Prediction of Turbulent, Three-Dimensional, Recirculating Flows with the NCC

The National Combustion Code (NCC) was used to calculate the steady state, nonreacting flow field of a prototype Lean Direct Injection (LDI) swirler. This configuration used nine groups of eight holes drilled at a thirty-five degree angle to induce swirl. These nine groups created swirl in the same direction, or a corotating pattern. The static pressure drop across the holes was fixed at approximately four percent. Computations were performed on one quarter of the geometry, because the geometry is considered rotationally periodic every ninety degrees. The final computational grid used was approximately 2.26 million tetrahedral cells, and a cubic nonlinear k - epsilon model was used to model turbulence. The NCC results were then compared to time averaged Laser Doppler Velocimetry (LDV) data. The LDV measurements were performed on the full geometry, but four ninths of the geometry was measured. One-, two-, and three-dimensional representations of both flow fields are presented. The NCC computations compare both qualitatively and quantitatively well to the LDV data, but differences exist downstream. The comparison is encouraging, and shows that NCC can be used for future injector design studies. To improve the flow prediction accuracy of turbulent, three-dimensional, recirculating flow fields with the NCC, recommendations are given.

Iannetti, A.↗

Comparing Anisotropic Output-Based Grid Adaptation Methods by Decomposition

Anisotropic grid adaptation is examined by decomposing the steps of flow solution, ad- joint solution, error estimation, metric construction, and simplex grid adaptation. Multiple implementations of each of these steps are evaluated by comparison to each other and expected analytic results when available. For example, grids are adapted to analytic metric fields and grid measures are computed to illustrate the properties of multiple independent implementations of grid adaptation mechanics. Different implementations of each step in the adaptation process can be evaluated in a system where the other components of the adaptive cycle are fixed. Detailed examination of these properties allows comparison of different methods to identify the current state of the art and where further development should be targeted.

Park, Michael A.↗

Turbulence and deterministic chaos

Several turbulent and nonturbulent solutions of the Navier-Stokes equations are obtained. The unaveraged equations are used numerically in conjunction with tools and concepts from nonlinear dynamics, including time series, phase portraits, Poincare sections, largest Liapunov exponents, power spectra, and strange attractors. Initially neighboring solutions for a low Reynolds number fully developed turbulence are compared. Several flows are noted: fully chaotic, complex periodic, weakly chaotic, simple periodic, and fixed-point. Of these, only fully chaotic is classified as turbulent. Besides the sustained flows, a flow which decays as it becomes turbulent is examined. For the finest grid, 128(exp 3) points, the spatial resolution appears to be quite good. As a final note, the variation of the velocity derivatives skewness of a Navier-Stokes flow as the Reynolds number goes to zero is calculated numerically. The value of the skewness is shown to become small at low Reynolds numbers, in agreement with intuitive arguments that nonlinear terms should be negligible.

Deissler, Robert G.↗

Investigating Multi-Microgrid Black Start Methods Using Grid-Forming Inverters: Preprint

This paper examines state-of-the-art microgrid (MG) black-start technologies with high IBR penetrations and proposes black-start and interconnection methods for 100% grid-forming (GFM) inverter-based MGs. A multiple-microgrid approach is proposed and compared to the existing methods in a 4-bus, 12-GFM inverter simulation test setup. This investigation involves intelligent synchronization units called "synchrobreakers," which enable the synchronization of multiple MGs based on their terminal measurements. MGs in this setup are held at different loading levels and comprise averaged models of GFM inverters with ideal voltage sources. The results of the black-start techniques are compared, and conclusions are drawn to better prepare MG planners and DSOs for the next-generation GFM inverter-based multi-MG black-start procedures.

black start↗

Integrated solution techniques for security constrained unit commitment problem

Apparatus and methods are disclosed for solving Mixed Integer Programming (MIP) problems, such as Security Constrained Unit Commitment (SCUC) problems used by power grid authorities to perform day-ahead market clearing. In certain examples, a plurality of threads of a software tool implementing a concurrent optimizer can be executed concurrently and sequentially to generate new solutions to a SCUC problem for an upcoming planning horizon. Data can be shared among the concurrently executing threads, such as intermediate/incumbent solutions and hints regarding the fixing of variables and constraints to reduce the size of the SCUC problem. In some examples, the threads are seeded with historical solutions from prior planning horizons. The software tool can select a best solution from the solutions generated by the threads, and determine dispatch instructions for a device coupled to the power grid for the upcoming planning horizon based at least in part on the selected solution.

Pan, Feng↗

Continuous Integration Service at Fermilab

This paper describes the current architecture of Continuous Integration (CI) service developed at Fermilab, encountered successes and difficulties, as well as future development plans. The current experiment code has hundreds of contributors that provide new features, bug fixes, and other improvements. Version control systems help developers to collaborate in contributing software for their experiments, while the CI system helps developers keep their code healthy. The Fermilab CI service allows experiments and projects to test and validate their offline production and analysis code on the supported platforms. It is built on top of Jenkins, designed to be set up from a configuration file that provides implementation for each phase of the CI workflow, and able to validate experiment code through grid jobs. This CI service provides a dashboard for easy access to logs and statistical graphs. Since the CI service has been adopted by Fermilab experiments/projects, it proved to be very useful to intercept issues in their code early on and get them fixed before running it in production. Currently the CI service is in use by the ArgoNeuT, DUNE, g-2, LArIAT, MINERvA, mu2e, NOvA, SBND and uBooNE experiments and by the following projects: ART and LArSoft software suites, GENIE, and Glidein-WMS. The CI service is under active development and planning to support code profiling.

Di benedetto, Vito↗

Investigating Multi-Microgrid Black Start Methods Using Grid-Forming Inverters

This paper examines state-of-the-art microgrid (MG) black-start technologies with grid-forming (GFM) inverter-based resources (IBRs) and proposes black start and interconnection methods for 100% inverter-based MGs. A multiple-MG approach is proposed and compared to the existing methods in a 4-bus, 12-GFM inverter simulation test setup. This investigation involves intelligent synchronization units that enable the autonomous synchronization of multiple MGs based on their terminal measurements. MGs in this setup are held at different loading levels and comprise averaged models of GFM inverters. The results of the black-start techniques are compared, and conclusions are drawn to better prepare MG planners and distribution system operators for next-generation, multi-MG, GFM inverter-based, black-start procedures.

black start↗

Reconductoring Economic and Financial Analysis Tool (REFA) v1.0

REFA is designed to help transmission planners better understand the financial, environmental and economic benefits of reconductoring upgrades, using traditional or advanced conductors. More specifically, the tool enables grid planners and utilities to demonstrate the entire lifetime value of reconductoring projects, providing potential justifications for projects with higher Capex - the purchase of long-term physical or fixed assets used in a business's operations. Existing tools are designed to compare technical aspects of different conductor applications. They typically do not allow a cost-benefit analysis over the lifetime of the conductor's operations.

Heleno, Miguel↗

Forebody and Inlet Design for the HIFiRE 2 Flight Test

A forebody and inlet have been designed for the HIFiRE 2 scramjet flight test. The test will explore the operating, performance, and stability characteristics of a simple hydrocarbon-fueled scramjet combustor as it transitions from dual-mode to scramjet-mode operation and during supersonic combustion at Mach 8+ flight conditions. Requirements for the compression system were derived from inlet starting and combustor inflow requirements as well as physical size constraints. The design process is described. A planar, fixed geometry, mixed compression concept was used to produce laterally uniform flow at the inlet entrance and a conservative amount of internal contraction with respect to inlet starting. A grid sensitivity study was performed so that important flow physics caused by three-dimensional shock boundary layer interactions could be captured with confidence. Results from low Mach number operability studies, nominal trajectory cases, and high dynamic pressure heat load cases are discussed. The forebody and inlet solutions provide information for on-going combustor calculations, mass capture across the trajectory for fuel system design, and surface heating rates for thermal/structural analysis. The design has a one freestream Mach number margin for inlet starting, exceeds the high Mach number combustor entrance pressure requirement, produces high quality flow at the inlet exit for all Mach numbers and vehicle attitudes in the design space, and fits inside the booster shroud.

Ferlemann, Paul G.↗

Crack Offset Measurement With the Projected Laser Target Device

The device and associated analysis methodology summarized in this report were developed for the purpose of estimating the size of discontinuities in the surface of the foam that covers the Space Shuttle External Tank. These surface offsets are thought to be due to subsurface cracks in the foam insulation. The mathematical analysis and procedure described here provide a method to quantity the dimensions of the crack offset in a direction perpendicular to the surface, making use of the projected laser target device (PLTD) tool and a laser line projector. The keys to the construction and use of the PLTD are the following geometrical design requirements: Laser dots are on a square grid: length on a side. Laser beams are perpendicular to projected surface. Beams are parallel out to the distance being projected. The PLTD can be used to (1) calibrate fixed cameras of unknown magnification and orientation (far-field solution); (2) provide equivalent calibration to multiple cameras, previously achieved only by the use of known target points, for example, in 3.D foreign-object debris tracking on a fixed launch platform; (3) compute scaling for conventional 2.D images, and depth of field for 3.D images (near-field solution); and (4) in conjunction with a laser line projector, achieve accurate measurements of surface discontinuity (cracks) in a direction perpendicular to the surface.

Source record↗

Methods and Results of a Search for Gravitational Waves Associated with Gamma-Ray Bursts Using the GEO 600, LIGO, and Virgo Detectors

In this paper we report on a search for short-duration gravitational wave bursts in the frequency range 64 Hz-1792 Hz associated with gamma-ray bursts (GRBs), using data from GEO600 and one of the LIGO or Virgo detectors. We introduce the method of a linear search grid to analyze GRB events with large sky localization uncertainties such as the localizations provided by the Fermi Gamma-ray Burst Monitor (GBM). Coherent searches for gravitational waves (GWs) can be computationally intensive when the GRB sky position is not well-localized, due to the corrections required for the difference in arrival time between detectors. Using a linear search grid we are able to reduce the computational cost of the analysis by a factor of O(10) for GBM events. Furthermore, we demonstrate that our analysis pipeline can improve upon the sky localization of GRBs detected by the GBM, if a high-frequency GW signal is observed in coincidence. We use the linear search grid method in a search for GWs associated with 129 GRBs observed satellite-based gamma-ray experiments between 2006 and 2011. The GRBs in our sample had not been previously analyzed for GW counterparts. A fraction of our GRB events are analyzed using data from GEO600 while the detector was using squeezed-light states to improve its sensitivity; this is the first search for GWs using data from a squeezed-light interferometric observatory. We find no evidence for GW signals, either with any individual GRB in this sample or with the population as a whole. For each GRB we place lower bounds on the distance to the progenitor, assuming a fixed GW emission energy of 10(exp -2)Stellar Mass sq c, with a median exclusion distance of 0.8 Mpc for emission at 500 Hz and 0.3 Mpc at 1 kHz. The reduced computational cost associated with a linear search grid will enable rapid searches for GWs associated with Fermi GBM events in the Advanced detector era.

methods↗

Power System Recovery Coordinated with (Non-)Black-Start Generators

Power restoration is an urgent task after a black-out, and recovery efficiency is critical when quantifying system resilience. Multiple elements should be considered to restore the power system quickly and safely. This paper proposes a recovery model to solve a direct-current optimal power flow (DCOPF) based on mixed-integer linear programming (MILP). Since most of the generators cannot start independently, the interaction between black-start (BS) and non-black-start (NBS) generators must be modeled appropriately. The energization status of the NBS is coordinated with the recovery status of transmission lines, and both of them are modeled as binary variables. Also, only after an NBS unit receives the cranking power through connected transmission lines, will it be allowed to participate in the following system dispatch. The amount of cranking power is estimated as a fixed proportion to the maximum generation capacity. The proposed model is validated on several test systems, as well as a 1393-bus representation system of the Puerto Rican electric power grid. Test results demonstrate how the recovery of NBS units and damaged transmission lines can be optimized, resulting in an efficient and well-coordinated recovery procedure.

Zhao, Meng↗

Numerical simulation of the effect of spatial disturbances on vortex asymmetry

The steady asymmetric vortex pattern observed on slender bodies of revolution at large angle of attack was investigated using fine-grid Navier-Stokes computations. The computed results demonstrate the marked asymmetry which has been observed in experiments. To obtain asymmetry, it was found essential to introduce a space-fixed time-invariant perturbation into the computation. The computational results suggest that vortex asymmetry is forced by amplification of small disturbances, such as those due to surface roughness, occurring within the body viscous boundary layer.

Degani, David↗

Calculation of unsteady linearized Euler flows in cascades using harmonically deforming grids

A method for calculating unsteady, inviscid, compressible flows in cascades is presented. Using the linearized Euler technique, the flow is decomposed into a steady or mean flow plus a harmonically varying small disturbance flow. The equations that describe the small disturbance flow are linear variable coefficient equations, and are solved using a pseudo-time time marching Lax-Wendroff technique. Unlike previous linearized methods, however, the solution is computed on a harmonically deforming computational grid that conforms to the motion of the vibrating airfoils. The mean flow and perturbation flow solutions are defined in the deforming coordinate system rather than in a coordinate system fixed in space. Hence, no extrapolation terms are required to implement the upwash boundary conditions at the airfoil surfaces, significantly improving the accuracy of the method. For transonic flow calculations, unsteady shock motions are modelled using shock capturing. The unsteady loads due to the shock motion are then seen as pressure impulses. Representative computational results are presented for transonic channel flows and subsonic and transonic cascade flows.

Hall, Kenneth C.↗

Supersonic vortex breakdown over a delta wing in transonic flow

The effects of freestream Mach number and angle of attack on the leading-edge vortex breakdown due to the terminating shock on a 65-degree, sharp-edged, cropped delta wing are investigated computationally, using the time-accurate solution of the laminar unsteady compressible full Navier-Stokes equations with the implicit upwind flux-difference splitting, finite-volume scheme. A fine O-H grid consisting of 125 x 85 x 84 points in the wrap-around, normal, and axial directions, respectively, is used for all the flow cases. Keeping the Reynolds number fixed at 3.23 x 10 exp 6, the Mach number is varied from 0.85 to 0.9 and the angle of attack is varied from 20 to 24 deg. The results show that, at 20-deg angle of attack, the increase of the Mach number from 0.85 to 0.9 results in moving the location of the terminating shock downstream. The results also show that, at 0.85 Mach number, the increase of the angle of attack from 20 to 24 deg results in moving the location of the terminating shock upstream. The results are in good agreement with the experimental data.

Kandil, Hamdy A.↗