Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Gradient Limiters”

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

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

At least 1,729 records · Page 96

MAS: A versatile Landau-fluid eigenvalue code for plasma stability analysis in general geometry

We have developed a new global eigenvalue code, Multiscale Analysis for plasma Stabilities (MAS), for studying plasma problems with wave toroidal mode number (n) and frequency (ω) in a broad range of interest in general tokamak geometry, based on a five-field Landau-fluid description of thermal plasmas. Beyond keeping the necessary plasma fluid response, we further retain the important kinetic effects including diamagnetic drift, ion finite Larmor radius, finite parallel electric field (E||), ion and electron Landau resonances in a self-consistent and non-perturbative manner without sacrificing the attractive efficiency in computation. The physical capabilities of the code are evaluated and examined in the aspects of both theory and simulation. In theory, the comprehensive Landau-fluid model implemented in MAS can be reduced to the well-known ideal MHD model, electrostatic ion-fluid model, and drift-kinetic model in various limits, which clearly delineates the physics validity regime. In simulation, MAS has been well benchmarked with theory and other gyrokinetic and kinetic-MHD hybrid codes in a manner of adopting the unified physical and numerical framework, which covers the kinetic Alfv\'en wave (KAW), ion sound wave (ISW), low-n kink, high-n ion temperature gradient mode (ITG) and kinetic ballooning mode (KBM). Moreover, MAS is successfully applied to model the Alfv\'en eigenmode (AE) activities in DIII-D discharge #159243, which faithfully captures the frequency sweeping of reversed shear Alfv\'en eigenmode (RSAE), the tunneling damping of toroidal Alfv\'en eigenmode (TAE), as well as the polarization characteristics of kinetic beta-induced Alfv\'en eigenmode (KBAE) and beta-induced Alfv\'en-acoustic eigenmode (BAAE) being consistent with former gyrokinetic theory and simulation. With respect to the key progress contributed to the community, MAS has the advantage of combining rich physics ingredients, realistic global geometry and high computation efficiency together for plasma stability analysis in linear regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Toward the Dynamical Convergence on the Jet Stream in Aquaplanet AGCMs

Systematic sensitivity of the jet position and intensity to horizontal model resolution is identified in several aquaplanet AGCMs, with the coarser resolution producing a more equatorward eddy-driven jet and a stronger upper tropospheric jet intensity. As the resolution of the models increases to 50km or finer, the jet position and intensity show signs of convergence within each model group. The mechanism for this convergence behavior is investigated using a hybrid Eulerian-Lagrangian finite-amplitude wave activity budget developed for the upper tropospheric absolute vorticity. The result suggests that the poleward shift of the 33 eddy-driven jet with higher resolution can be attributed to the smaller effective diffusivity of the model in the midlatitudes that allows more wave activity to survive the dissipation and to reach the subtropical critical latitude for wave breaking. The enhanced subtropical wave breaking and associated irreversible vorticity mixing act to maintain a more poleward peak of the vorticity gradient, and thus a more poleward jet. Being over-dissipative, the coarse-resolution AGCMs misrepresent the nuanced nonlinear aspect of the midlatitude eddy-mean flow interaction, giving rise to the equatorward bias of the eddy-driven jet. In accordance with the asymptotic behavior of effective diffusivity of Batchelor turbulence in the large Peclet number limit, the upper tropospheric effective diffusivity of the aquaplanet AGCMs displays signs of convergence in the midlatitude towards a value of ~107 m2 s-1 for the ?2 diffusion. This provides a dynamical underpinning for the convergence of the jet stream observed in these AGCMs at high resolution.

jet position, horizontal model resolution, AGCM,↗

Best estimate of the planetary boundary layer height from multiple remote sensing measurements

Remote sensing measurements have been widely used to estimate the planetary boundary layer height (PBLHT). Each remote sensing approach offers unique strengths and faces different limitations. In this study, we use machine learning (ML) methods to produce a best-estimate PBLHT (PBLHT-BE-ML) by integrating four PBLHT estimates derived from remote sensing measurements at the Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) observatory. Three ML models – random forest (RF) classifier, RF regressor, and light gradient-boosting machine (LightGBM) – were trained on a dataset from 2017 to 2023 that included radiosonde, various remote sensing PBLHT estimates, and atmospheric meteorological conditions. Evaluations indicated that PBLHT-BE-ML from all three models improved alignment with the PBLHT derived from radiosonde data (PBLHT-SONDE), with LightGBM demonstrating the highest accuracy under both stable and unstable boundary layer conditions. Feature analysis revealed that the most influential input features at the SGP site were the PBLHT estimates derived from (a) potential temperature profiles retrieved using Raman lidar (RL) and atmospheric emitted radiance interferometer (AERI) measurements (PBLHT-THERMO), (b) vertical velocity variance profiles from Doppler lidar (PBLHT-DL), and (c) aerosol backscatter profiles from micropulse lidar (PBLHT-MPL). The trained models were then used to predict PBLHT-BE-ML at a temporal resolution of 10 min, effectively capturing the diurnal evolution of PBLHT and its significant seasonal variations, with the largest diurnal variation observed over summer at the SGP site. We applied these trained models to data from the ARM Eastern Pacific Cloud Aerosol Precipitation Experiment (EPCAPE) field campaign (EPC), where the PBLHT-BE-ML, particularly with the LightGBM model, demonstrated improved accuracy against PBLHT-SONDE. Analyses of model performance at both the SGP and EPC sites suggest that expanding the training dataset to include various surface types, such as ocean and ice-covered areas, could further enhance ML model performance for PBLHT estimation across varied geographic regions.

Zhang, Damao [Pacific Northwest National Laborator↗

Optimal detumbling of a large manned spacecraft using an internal moving mass

This investigation deals with the use of a movable mass control system to stabilize a tumbling asymmetric spacecraft about the maximum inertia axis. A first-order gradient optimization technique is used to minimize angular velocity components along the intermediate and minimum inertia axes, thus, permitting a wide range of initial guesses for mass position history. Motion of the control mass is along a linear track fixed in the vehicle. The control variable is taken as mass acceleration with respect to body coordinates. Motion is limited to defined quantities and a penalty function is used to insure a given range of positions. Numerical solutions of the optimization equations verify that minimum time detumbling is achieved with the largest permissible movable mass, length of linear track, and positions of the mass on the two coordinates perpendicular to the linear motion. The optimal method permits detumbling in about one-fourth the time when compared to a force control law formulation available in the literature.

Kunciw, B. G.↗

Advances in electrophoretic separations

Free fluid electrophoresis is described using laboratory and space experiments combined with extensive mathematical modeling. Buoyancy driven convective flows due to thermal and concentration gradients are absent in the reduced gravity environment of space. The elimination of convection in weightlessness offers possible improvements in electrophoresis and other separation methods which occur in fluid media. The mathematical modeling suggests new ways of doing electrophoresis in space and explains various phenomena observed during past experiments. The extent to which ground based separation techniques are limited by gravity induced convection is investigated and space experiments are designed to evaluate specific characteristics of the fluid/particle environment. A series of experiments are proposed that require weightlessness and apparatus is developed that can be used to carry out these experiments in the near future.

Snyder, R. S.↗

Thermocapillary Bubble Migration: Thermal Boundary Layers for Large Marangoni Numbers

The migration of an isolated gas bubble in an immiscible liquid possessing a temperature gradient is analyzed in the absence of gravity. The driving force for the bubble motion is the shear stress at the interface which is a consequence of the temperature dependence of the surface tension. The analysis is performed under conditions for which the Marangoni number is large, i.e. energy is transferred predominantly by convection. Velocity fields in the limit of both small and large Reynolds numbers are used. The thermal problem is treated by standard boundary layer theory. The outer temperature field is obtained in the vicinity of the bubble. A similarity solution is obtained for the inner temperature field. For both small and large Reynolds numbers, the asymptotic values of the scaled migration velocity of the bubble in the limit of large Marangoni numbers are calculated. The results show that the migration velocity has the same scaling for both low and large Reynolds numbers, but with a different coefficient. Higher order thermal boundary layers are analyzed for the large Reynolds number flow field and the higher order corrections to the migration velocity are obtained. Results are also presented for the momentum boundary layer and the thermal wake behind the bubble, for large Reynolds number conditions.

Balasubramaniam, R.↗

Thin Interface Asymptotics for an Energy/Entropy Approach to Phase-Field Models with Unequal Conductivities

Karma and Rapped recently developed a new sharp interface asymptotic analysis of the phase-field equations that is especially appropriate for modeling dendritic growth at low undercoolings. Their approach relieves a stringent restriction on the interface thickness that applies in the conventional asymptotic analysis, and has the added advantage that interfacial kinetic effects can also be eliminated. However, their analysis focussed on the case of equal thermal conductivities in the solid and liquid phases; when applied to a standard phase-field model with unequal conductivities, anomalous terms arise in the limiting forms of the boundary conditions for the interfacial temperature that are not present in conventional sharp-interface solidification models, as discussed further by Almgren. In this paper we apply their asymptotic methodology to a generalized phase-field model which is derived using a thermodynamically consistent approach that is based on independent entropy and internal energy gradient functionals that include double wells in both the entropy and internal energy densities. The additional degrees of freedom associated with the generalized phased-field equations can be chosen to eliminate the anomalous terms that arise for unequal conductivities.

McFadden, G. B.↗

Microelectrical Mechanical Systems Flow Control Used to Manage Engine Face Distortion in Compact Inlet Systems

Turbofan engine-face flow distortion is one of the most troublesome and least understood problems for designers of modern engine inlet systems. One concern is that there are numerous sources of flow-field distortion that are ingested by the inlet or generated within the inlet duct itself. Among these are: (1) flow separation at the cowl lip during in-flight maneuvering, (2) flow separation on the compression surfaces due to shock-wave/boundary layer interactions, (3) spillage of the fuselage boundary layer into the inlet duct, (4) ingestion of aircraft vortices and wakes emanating from upstream disturbances, and (5) strong secondary flow gradients and flow separation induced by wall curvature within the inlet duct itself. Most developing aircraft (including the B70, F-111, F-14, Mig-25, Tornado, and Airbus A300) have experienced one or more of these types of problems, particularly at high Mach numbers and/or extreme maneuver conditions when flow distortion at the engine face exceeded the allowable limits of the engine.

Anderson, Bernhard H.↗

Wildlife Densities and Habitat Use Across Temporal and Spatial Scales on the Mid-Atlantic Outer Continental Shelf. Final Report to the Department of Energy EERE Wind & Water Power Technologies Office

The Mid-Atlantic Baseline Studies Project helped address environmental barriers to offshore wind energy development in the mid-Atlantic region by providing regulators, developers, and other stakeholders with comprehensive baseline ecological data and analyses. Project funders and collaborators from a range of academic institutions, non-governmental organizations, federal agencies, foundations, and private companies came together to study bird, sea turtle, and marine mammal distributions, densities, and movements on the mid-Atlantic Outer Continental Shelf between 2012 and 2014. Specific project activities and goals included the following: (1) Conduct standardized surveys to quantify bird, sea turtle, and marine mammal densities seasonally and annually throughout the study region and identify important habitat use or aggregation areas. (2) Develop statistical models to help understand the drivers of wildlife distribution and abundance patterns. (3) Use individual tracking data for several focal bird species to provide information on population connectivity and individual movements that is complementary to survey data. (4) Identify species that are likely to be exposed to offshore wind energy development activities in the mid-Atlantic study area. (5) Develop U.S.-based technological resources and assessment methods for future monitoring efforts, including a comparison of high resolution digital video aerial surveys to boat-based surveys. (6) Help meet data needs associated with National Environmental Policy Act (NEPA), Marine Mammal Protection Act, and Endangered Species Act requirements, by contributing several years of data and analysis towards future Environmental Impact Statements. This report consists of six parts: Project overview (executive summary and Chapters 1-2); Examining wildlife distributions and relative abundance from a digital video aerial survey platform (Chapters 3-6); Examining wildlife distributions and abundance using boat-based surveys (Chapters 7-12); Integrating data across survey platforms (Chapters 13-19); Individual movements and habitat use for focal bird species (Chapters 20-25); and Nocturnal avian migration monitoring (Chapters 26-27). Boat-based and digital video aerial surveys each had specific advantages and disadvantages, but were largely complementary (Chapters 1, 5, 13-14). Digital aerial surveys were particularly useful for covering offshore areas at broad scales, where general distributions of taxonomic groups were a priority; boat surveys could provide more detailed data on species identities and behaviors, but were more limited in geographic scope due to their slower survey pace. The mid-Atlantic study area was important for wintering and breeding taxa, and its location also made it a key migratory corridor. There was considerable variation in species composition and spatial patterns by season, largely driven by dynamic environmental conditions (Chapters 12, 15, and 20-22). Habitat gradients in nearshore waters, however, were reliable influences on productivity and patterns of species distributions and abundance. Areas within about 30-40 km of the coast offshore of the mouths of Chesapeake and Delaware Bays, as well as to the south of Delaware Bay along the coast, were consistent hotspots of abundance and species diversity, regardless of survey methodology or analytical approach (Chapters 2, 12, 17). Inter-annual variation was substantial, and the importance of certain environmental variables in predicting animal distributions indicates that these species may well respond to future environmental shifts brought about by anthropogenic effects and climatic change. This study is an important first step, however, towards understanding how bird, marine mammal, and sea turtle populations in the mid-Atlantic may be exposed to offshore wind energy development and other anthropogenic activities. The results of this study provide insight to help address environmental permitting requirements for current and future offshore development projects, and serve as a starting point for more site-specific studies, risk analyses, and evaluation of potential measures to avoid and minimize those risks.

17 WIND ENERGY↗

Test Results of the CERN HL-LHC low- Quadrupole Short Models MQXFS3c and MQXFS4

For the High Luminosity Upgrade of the CERN Large Hadron Collider, lower * quadrupole magnets based on advanced Nb3Sn conductors will be installed on each side of the ATLAS and CMS experiment insertion zones. As part of the technological developments needed to achieve the required field gradient of 132.6 T/m within a 150-mm aperture, short length model magnets, named MQXFS, are tested both at the CERN SM18 and FERMI Lab test facilities. The model magnets rely on two types of Nb3Sn conductors (RRP and PIT) and on an innovative bladders and keys design to provide mechanical support against the Lorentz forces. In 2016 and 2017, the powering tests of the first two models MQXFS3 (RRP) and MQXFS5 (PIT) proved that nominal performance (16.5 kA) could be reached with excellent memory of the quench current after thermal cycle. However both magnets showed a slow training behavior with clear observations of voltage disturbances before the quench. Besides MQXFS5 only could reach ultimate current (17.9 kA) whereas erratic behavior was observed on MQXFS3 due to conductor local degradation at the head of one of the coil. In 2018, this limiting coil was changed and the applied azimuthal pre-stress increased. If ultimate current could then be reached, no stable current could be maintained due to identified defect on the outer layer of the new coil. Finally the outcome of the test of the new model MQXFS4, featuring the final RRP conductors that will be used for the series production and variation on the inner layer quench heater designs are here reported in details.

43 PARTICLE ACCELERATORS↗

AN EXTERNAL ORIGIN FOR CARBON MONOXIDE ON URANUS FROM HERSCHEL/SPIRE?

Previous studies have demonstrated an external source of CO on Jupiter, Saturn, and Neptune. However, it has not been possible to demonstrate this on Uranus because of its low CO abundance, low upper-tropospheric temperatures, and low stratospheric thermal gradient, which make detection very challenging. Here we use 17 Herschel/SPIRE observation sequences spanning 3 yr (2009-2012), which cover 14.6-51.8 cm{sup –1} with a combined integration time of 5 hr. These spectra were originally taken for routine calibration purposes, so were corrected for continuum offsets prior to analysis. The final stacked spectra had an extremely low noise level of 10-50 pW cm{sup –2} sr{sup –1}/cm{sup –1}. Despite this, CO was not observed, but we were able to obtain stringent 3σ upper limits at the 0.1-0.2 bar level of 2.1 ppb for a uniform profile, and 9.4 ppb for a stratosphere-only profile—an order of magnitude improvement over previous studies. Comparison with observed CO fluorescence by Encrenaz et al. suggests the majority of Uranus' stratospheric CO has an external origin. It thus appears that external supply of oxygen species—via comets, micrometeorites, or dust—is an important process on all giant planets in our solar system.

79 ASTRONOMY AND ASTROPHYSICS↗

Multipoint Aerostructural Optimization of Wind Turbine Rotors Using a Coupled Blade‐Resolved Aerostructural Solver

Physics‐based design optimization workflows thread the needle between computational cost limitations and simulation complexity, often compromising between modeling detail and the range of operating design conditions. Multipoint aerostructural optimization of wind turbine rotors has so far been confined to low‐fidelity analyses or to high‐fidelity studies with simplified structural models, leaving the most complex design trade‐offs unexplored. We close this gap by performing the first tightly coupled gradient‐based multipoint aerostructural rotor optimization using 3D aerodynamic and structural solvers with discrete coupled adjoints. The optimizer simultaneously varies blade planform, airfoil shapes, and structural thickness through more than 270 design variables, minimizing a weighted combination of rotor mass and power across multiple wind speeds. Applied to a modified DTU 10‐MW benchmark under conservative structural and aerodynamic constraints, our multipoint optimization reduces rotor mass by up to 36% and increases power by 12%–15% across the main operating conditions; biasing the objective toward power yields power gains up to 18% and a 17% mass reduction. For a nominal wind distribution, 3‐point rotor designs accounting for low RPM and high thrust conditions capture dominant trade‐offs and outperform single‐point designs. Adding two off‐design points changes individual‐condition power by less than 3% but leaves the weighted average within 0.5%, and the mass‐power bias has a stronger effect on the final design than the operating‐point weighting itself. Our framework extends naturally to richer load cases and site‐specific wind distributions, providing a basis for high‐fidelity multipoint design earlier in industrial workflows.

17 WIND ENERGY↗

Cooperative pulses for pseudo-pure state preparation

Using an extended version of the optimal-control-based gradient ascent pulse engineering algorithm, cooperative (COOP) pulses are designed for multi-scan experiments to prepare pseudo-pure states in quantum computation. COOP pulses can cancel undesired signal contributions, complementing and generalizing phase cycles. They also provide more flexibility and, in particular, eliminate the need to select specific individual target states and achieve the fidelity of theoretical limit by flexibly choosing appropriate number of scans and duration of pulses. The COOP approach is experimentally demonstrated for three-qubit and four-qubit systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

X-ray transparent microfluidic chips for high-throughput screening and optimization of in meso membrane protein crystallization

Elucidating and clarifying the function of membrane proteins ultimately requires atomic resolution structures as determined most commonly by X-ray crystallography. Many high impact membrane protein structures have resulted from advanced techniques such as in meso crystallization that present technical difficulties for the set-up and scale-out of high-throughput crystallization experiments. In prior work, we designed a novel, low-throughput X-ray transparent microfluidic device that automated the mixing of protein and lipid by diffusion for in meso crystallization trials. Here, in this work, we report X-ray transparent microfluidic devices for high-throughput crystallization screening and optimization that overcome the limitations of scale and demonstrate their application to the crystallization of several membrane proteins. Two complementary chips are presented: (1) a high-throughput screening chip to test 192 crystallization conditions in parallel using as little as 8 nl of membrane protein per well and (2) a crystallization optimization chip to rapidly optimize preliminary crystallization hits through fine-gradient re-screening. We screened three membrane proteins for new in meso crystallization conditions, identifying several preliminary hits that we tested for X-ray diffraction quality. Further, we identified and optimized the crystallization condition for a photosynthetic reaction center mutant and solved its structure to a resolution of 3.5 Å.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Interparameter trade-off quantification and reduction in isotropic-elastic full-waveform inversion: synthetic experiments and Hussar land data set application

The problem of inverting for multiple physical parameters in the subsurface using seismic full-waveform inversion (FWI) is complicated by interparameter trade-off arising from inherent ambiguities between different physical parameters. Parameter resolution is often characterized using scattering radiation patterns, but these neglect some important aspects of interparameter trade-off. More general analysis and mitigation of interparameter trade-off in isotropic-elastic FWI is possible through judiciously chosen multiparameter Hessian matrix-vector products. We show that products of multiparameter Hessian off-diagonal blocks with model perturbation vectors, referred to as interparameter contamination kernels , are central to the approach. We then apply the multiparameter Hessian to various vectors designed to provide information regarding the strengths and characteristics of interparameter contamination, both locally and within the whole volume. With numerical experiments, we observe that S-wave velocity perturbations introduce strong contaminations into density and phase-reversed contaminations into P-wave velocity, but themselves experience only limited contaminations from other parameters. Based on these findings, we introduce a novel strategy to mitigate the influence of interparameter trade-off with approximate contamination kernels. Furthermore, we recommend that the local spatial and interparameter trade-off of the inverted models be quantified using extended multiparameter point spread functions (EMPSFs) obtained with pre-conditioned conjugate-gradient algorithm. Compared to traditional point spread functions, the EMPSFs appear to provide more accurate measurements for resolution analysis, by de-blurring the estimations, scaling magnitudes and mitigating interparameter contamination. Approximate eigenvalue volumes constructed with stochastic probing approach are proposed to evaluate the resolution of the inverted models within the whole model. With a synthetic Marmousi model example and a land seismic field data set from Hussar, Alberta, Canada, we confirm that the new inversion strategy suppresses the interparameter contamination effectively and provides more reliable density estimations in isotropic-elastic FWI as compared to standard simultaneous inversion approach.

58 GEOSCIENCES↗

Application of the weighted-density approximation to the accurate description of electron-positron correlation effects in materials

Here, we discuss positron-annihilation lifetimes for a set of illustrative bulk materials within the framework of the weighted-density approximation (WDA). The WDA can correctly describe electron-positron correlations in strongly inhomogeneous systems, such as surfaces, where the applicability of (semi-)local approximations is limited. We analyze the WDA in detail and show that the electrons which cannot screen external charges efficiently, such as the core electrons, cannot be treated accurately via the pair correlation of the homogeneous electron gas. We discuss how this problem can be addressed by reducing the screening in the homogeneous electron gas by adding terms depending on the gradient of the electron density. Further improvements are obtained when core electrons are treated within the LDA and the valence electron using the WDA. Finally, we discuss a semiempirical WDA-based approach in which a sum rule is imposed to reproduce the experimental lifetimes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Quasiparticle self-consistent GW band structures and phase transitions of LiAlO 2 in tetrahedrally and octahedrally coordinated structures

In this work, a first-principles computational study is presented of various phases of LiAlO 2 . The relative total energies and equations of state of the α, β, and γ phases are determined after structural relaxation of each phase. The β and γ tetrahedral phases are found to be very close in energy and lattice volume with the γ phase having the lowest energy. The octahedral α phase is a high-pressure phase and the transition pressure from the γ and β phases to α is determined to be about 1 GPa. The electronic band structures of each phase at their own equilibrium volume are determined using the quasiparticle self-consistent (QS) GW method as well as using the 0.8 Σ approach in which the QS GW self-energy is reduced by a factor of 0.8 to correct for the underscreening of W in QS GW. The effective masses of the band edges and the nature of the band gaps are presented. The lowest energy γ phase is found to have a pseudodirect gap of 7.69 eV. The gap is direct at Γ but corresponds to a dipole forbidden transition. The imaginary part of the dielectric function and the absorption coefficient are calculated in the long-wavelength limit and the random phase approximation, without local field or electron-hole interaction effects for each phase, and their anisotropies are discussed. Si doping on the Al site is investigated as a possible n-type dopant in γ-LiAlO 2 using a 128-atom supercell corresponding to 3.125% Si on the Al sublattice in the generalized gradient approximation and a smaller 16-atom cell with 25% Si in the QS GW approximation. The Si is found to significantly perturb the conduction band and lower the gap but a clearly separated deep donor defect level is not found. However, the donor binding energy is still expected to be relatively deep, on the order of a few tenths eV in the hydrogenic effective mass approximation.

36 MATERIALS SCIENCE↗