Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “local minimum”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Nonlinear control of the minimum safety factor in tokamaks by optimal allocation of spatially moving electron cyclotron current drive

The minimum value of the safety factor profile is related to the magnetohydrodynamic (MHD) stability of the plasma confined in a tokamak. Therefore, active control of the minimum safety factor may mitigate MHD instabilities that can degrade or even terminate plasma confinement. Typically, in most tokamak scenarios, the minimum safety factor evolves spatially with time, i.e., the location at which the safety factor achieves the minimum value changes with time. In addition to the inherent nonlinearities in the minimum safety factor evolution, its spatial variation makes the control design challenging. In particular, complexity in control design may arise from the need for time-dependent nonlinear models that account for spatial variation of the minimum safety factor. Furthermore, the minimum safety factor may drift to locations where the actuator authority is low. The problem of minimum safety factor control with target location tracking and moving electron cyclotron current drive (ECCD) is addressed in this work. A nonlinear time-dependent model that incorporates the spatial variation of the minimum safety factor is presented. A nonlinear controller based on optimal feedback linearization is developed to track a target minimum safety factor. The proposed controller treats the ECCD position as a controllable variable. In other words, the controller prescribes the ECCD position (in addition to the non-inductive powers) in real time based on an optimal criterion that is defined a priori. This work also presents the steps necessary to integrate the minimum safety factor controller with a total energy controller to achieve multiple control objectives simultaneously. In conclusion, the proposed integrated control algorithm is tested using nonlinear simulations in the Control Oriented Transport SIMulator (COTSIM) for a DIII-D tokamak scenario.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1973

Noniterative estimation of a nonlinear parameter

An algorithm is described which solves the parameters X = (x1,x2,...,xm) and p in an approximation problem Ax nearly equal to y(p), where the parameter p occurs nonlinearly in y. Instead of linearization methods, which require an approximate value of p to be supplied as a priori information, and which may lead to the finding of local minima, the proposed algorithm finds the global minimum by permitting the use of series expansions of arbitrary order, exploiting an a priori knowledge that the addition of a particular function, corresponding to a new column in A, will not improve the goodness of the approximation.

Bergstroem, A.↗
NASA NTRS1983

A simple Lagrangian forecast system with aviation forecast potential

A trajectory forecast procedure is developed which uses geopotential tendency fields obtained from a simple, multiple layer, potential vorticity conservative isentropic model. This model can objectively account for short-term advective changes in the mass field when combined with fine-scale initial analyses. This procedure for producing short-term, upper-tropospheric trajectory forecasts employs a combination of a detailed objective analysis technique, an efficient mass advection model, and a diagnostically proven trajectory algorithm, none of which require extensive computer resources. Results of initial tests are presented, which indicate an exceptionally good agreement for trajectory paths entering the jet stream and passing through an intensifying trough. It is concluded that this technique not only has potential for aiding in route determination, fuel use estimation, and clear air turbulence detection, but also provides an example of the types of short range forecasting procedures which can be applied at local forecast centers using simple algorithms and a minimum of computer resources.

Petersen, R. A.↗
NASA NTRS1989

The ion energy spectra in the ring current during the geomagnetic storm of February 1986

The largest geomagnetic storm in several decades occurred in February 1986, with a minimum Dst of -312 nT. In this experiment, the CHEM instrument on the AMPTE/CCE spacecraft has been used to study the development of this storm in more detail. By comparing measurements of phase space density (f) vs magnetic moment (mu) for four ion species H(+), O(+), He (+), and He(++) during consecutive passes through the ring current region, a distinction can be made between flux increases that require the injection of a new population and those that are consistent with the adiabatic acceleration of a pre-existing population. It is shown that the smaller flux observed during the outbound pass compared to the preceding inbound pass at the time of minimum Dst is due to the difference in local times of the two measurements. This local time difference is consistent with a greatly enhanced convection electric field, which brings a new population from the magnetotail to the post-noon, but not to pre-noon local time sector.

Kistler, L. M.↗
NASA NTRS1995

Dissecting the EUV spectrum of Capella

Extreme ultraviolet spectra of Capella, obtained at various orbital phases over the past two years by the EUVE satellite, show strong emission lines from a continuous distribution of temperatures (approximately 10(exp 5 -10(exp 7.3) K). In addition to the strong He II lambda 303.8, the spectra are dominated by emission lines of highly ionized iron. Strong lines of Fe IX, XV, XVI, and XVIII-XXIV are used to construct emission measure distributions for the individual pointings, which show several striking features, including a minimum near 10(exp 6) K and a local maximum at 10(exp 6.8) K. Furthermore, intensities of the highest temperature lines (Te is greater than 10(exp 7) K) show variations (factors of 2-3) at different orbital phases, while the lower temperature Fe lines show variations of about 30% or less. The low variability of most of the strong low temperature features motivates a detailed analysis of the summed spectrum. With approximately 280 ksec of total exposure time, we have measured over 200 emission features with S/N greater than or equal to 3.0 in the summed spectrum. We report here initial results from the analysis of this spectrum. We can now identify lines of Fe VIII and X-XIV, as well as a number of electron density and abundance diagnostic lines. We also report here the first direct measurement of the continuum flux around approximately 100 A in a cool star atmosphere with EUVE. The continuum flux can be predicted from the emission measure model based on Fe line emission, and demonstrates that the Fe/H abundance ratio is dose to the solar photospheric value.

Brickhouse, Nancy S.↗
NASA NTRS2001

A Shell/3D Modeling Technique for the Analyses of Delaminated Composite Laminates

A shell/3D modeling technique was developed for which a local three-dimensional solid finite element model is used only in the immediate vicinity of the delamination front. The goal was to combine the accuracy of the full three-dimensional solution with the computational efficiency of a plate or shell finite element model. Multi-point constraints provided a kinematically compatible interface between the local three-dimensional model and the global structural model which has been meshed with plate or shell finite elements. Double Cantilever Beam (DCB), End Notched Flexure (ENF), and Single Leg Bending (SLB) specimens were modeled using the shell/3D technique to study the feasibility for pure mode I (DCB), mode II (ENF) and mixed mode I/II (SLB) cases. Mixed mode strain energy release rate distributions were computed across the width of the specimens using the virtual crack closure technique. Specimens with a unidirectional layup and with a multidirectional layup where the delamination is located between two non-zero degree plies were simulated. For a local three-dimensional model, extending to a minimum of about three specimen thicknesses on either side of the delamination front, the results were in good agreement with mixed mode strain energy release rates obtained from computations where the entire specimen had been modeled with solid elements. For large built-up composite structures modeled with plate elements, the shell/3D modeling technique offers a great potential for reducing the model size, since only a relatively small section in the vicinity of the delamination front needs to be modeled with solid elements.

Krueger, Ronald↗
NASA NTRS2015

Benchmark Data Set for Wheat Growth Models: Field Experiments and AgMIP Multi-Model Simulations.

The data set includes a current representative management treatment from detailed, quality-tested sentinel field experiments with wheat from four contrasting environments including Australia, The Netherlands, India and Argentina. Measurements include local daily climate data (solar radiation, maximum and minimum temperature, precipitation, surface wind, dew point temperature, relative humidity, and vapor pressure), soil characteristics, frequent growth, nitrogen in crop and soil, crop and soil water and yield components. Simulations include results from 27 wheat models and a sensitivity analysis with 26 models and 30 years (1981-2010) for each location, for elevated atmospheric CO2 and temperature changes, a heat stress sensitivity analysis at anthesis, and a sensitivity analysis with soil and crop management variations and a Global Climate Model end-century scenario.

climate change impact↗
NASA NTRS2017

The Mechanism for the Energy Buildup Driving Solar Eruptive Events

The underlying origin of solar eruptive events (SEEs), ranging from giant coronal mass ejections to small coronalhole jets, is that the lowest-lying magnetic flux in the Sun's corona undergoes continual buildup of stress and free energy. This magnetic stress has long been observed as the phenomenon of "filament channels:" strongly sheared magnetic field localized around photospheric polarity inversion lines. However, the mechanism for the stress buildup-formation of filament channels-is still debated. We present magnetohydrodynamic simulations of a coronal volume that is driven by transient, cellular boundary flows designed to model the processes by which the photosphere drives the corona. The key feature of our simulations is that they accurately preserve magnetic helicity, the topological quantity that is conserved even in the presence of ubiquitous magnetic reconnection. Although small-scale random stress is injected everywhere at the photosphere, driving stochastic reconnection throughout the corona, the net result of the magnetic evolution is a coherent shearing of the lowest-lying field lines. This highly counterintuitive result-magnetic stress builds up locally rather than spreading out to attain a minimum energy state-explains the formation of filament channels and is the fundamental mechanism underlying SEEs. Furthermore, this process is likely to be relevant to other astrophysical and laboratory plasmas.

Knizhnik, K. J.↗
NASA NTRS1994

Nonlinear equilibration of localized instabilities on a baroclinic jet

Dynamical mechanisms underlying the equilibration of absolute instability are examined in a nonlinear, quasigeostrophic, two-layer model. The key to understanding the nonlinear equilibration is in recognizing that linear absolute instabilities can be stabilized both by a reduction of the vertical shear and by enhancement of the mean barotropic velocity. In a localized domain, the equilibration process proceeds with the creation of locally convectively unstable regions downstream, which encroach onto the locally absolutely unstable region until the local instability is suppressed. That local instabilities exist only if absolutely unstable regions span a minimum size is verified by eigenvalue calculations of three-dimensional flows. Numerical examples suggest that this critical size is at least 9000 km for a wide range of parameter values chosen to investigate the midlatitude storm tracks. Fluctuations arising from local absolute instability obtain maximum amplitude in the downstream convectively unstable regions rather than in the absolutely unstable regions themselves. Together, these results suggest that if an equilibrated absolute instability were to occur in midlatitudes, a zonal band of surface easterlies exceeding 900 km would be required and the associated enhanced variances would not be found coincident with the regions of absolute instability.

Delsole, Timothy M.↗
DOE OSTI2020

Temperature estimation in a spatially inhomogeneous flame by diode laser absorption spectroscopy

A two-temperature (2T) model is proposed for estimating the temperature of a spatially inhomogeneous hot zone using single-beam diode laser absorption spectrometry. The proposed algorithm is based on fitting experimental absorption spectra by a linear combination of two database-simulated single-temperature spectra with different temperatures. The model efficiency is experimentally demonstrated when determining the temperature in different flame zones of a Wolfgard – Parker slot burner. The proposed 2T algorithm is used to find the maximum and minimum temperatures in different flame sections, which are compared with local temperatures in the same flame zones measured by the method of coherent anti-Stokes Raman scattering (CARS). It is shown that the maximum temperatures determined by the proposed 2T algorithm are in good agreement with CARS data in flame zones where high-temperature regions prevail. A good agreement between the minimum temperatures obtained by these two methods is also observed for predominantly cold zones. (paper)

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS1982

Surface waves on solar wind tangential discontinuities

It is demonstrated that (tangential) discontinuities in the magnetic field direction can support MHD surface waves. The surface waves are similar to the usual Alfven wave, but there are seven important differences. The first is that the surface waves exhibit a low-frequency cutoff; the second is that the velocity and magnetic field fluctuations are elliptically, and sometimes circularly, polarized. It is noted that they may account for the solar wind helicity spectrum. The third difference is that the surface waves are compressive, although there are special cases where they are noncompressive. The fourth is that the wave vector k, the local normals to the surfaces of constant phase, and the magnetic minimum variance direction do not all coincide. The fifth is that there is a tendency for the minimum variance direction to align itself with the mean magnetic field direction. The sixth difference is that the waves can be intrinsically nonplanar, and the seventh is that equipartition between magnetic and kinetic energies is not obeyed locally. These properties of the surface waves are interpreted to mean that surface waves may be common in the solar wind.

Hollweg, J. V.↗
NASA NTRS1988

Computation of steady and unsteady vortex-dominated flows with shock waves

The unsteady Euler equations have been derived in the conservation form for the flow relative motion with respect to a rotating frame of reference. The resulting equations are solved by using a central-difference finite-volume scheme with four-state Runge-Kutta time stepping. For steady flow applications local time stepping is used, and for unsteady applications the minimum global time stepping is used. A three-dimensional fully vectorized computer program has been developed and applied to steady and unsteady maneuvering delta wings. The capability of the three-dimensional program has been demonstrated for a rigid sharp-edged delta wing undergoing uniform rolling in a conical flow and rolling oscillations in a locally conical flow.

Kandil, Osama A.↗
NASA NTRS1990

Global models of Ne and Te at solar maximum based on DE-2 measurements

Newly developed global models of the electron density of (Ne) and (Te) in the F-region at solar maximum, based on Langmuir probe measurements from the Dynamics Explorer-2 satellite, are compared with solar minimum models that were developed earlier from Atmosphere Explorer data. Spherical harmonics are used in both models to describe the variations with geomagnetic latitude and local time, but the solar maximum model also includes longitudinal variations. The solar minimum models were for the fixed altitudes of 300 km and 400 km, while the solar maximum model covers all altitudes between 300 and 1000 km. In this paper, the global patterns of Ne and Te at 400 km at solar maximum and minimum are compared with the IRI model for the corresponding parts of the solar cycle. In most respects, the IRI model describes the empirical models quite well at both solar maximum and solar minimum.

Brace, L. H.↗
NASA NTRS1990

Global optimization methods for engineering design

The problem is to find a global minimum for the Problem P. Necessary and sufficient conditions are available for local optimality. However, global solution can be assured only under the assumption of convexity of the problem. If the constraint set S is compact and the cost function is continuous on it, existence of a global minimum is guaranteed. However, in view of the fact that no global optimality conditions are available, a global solution can be found only by an exhaustive search to satisfy Inequality. The exhaustive search can be organized in such a way that the entire design space need not be searched for the solution. This way the computational burden is reduced somewhat. It is concluded that zooming algorithm for global optimizations appears to be a good alternative to stochastic methods. More testing is needed; a general, robust, and efficient local minimizer is required. IDESIGN was used in all numerical calculations which is based on a sequential quadratic programming algorithm, and since feasible set keeps on shrinking, a good algorithm to find an initial feasible point is required. Such algorithms need to be developed and evaluated.

Arora, Jasbir S.↗
DOE OSTI2022

Generating Bright Emissive States by Modulating the Bandgap of Monolayer Tungsten Diselenide

Transition metal dichalcogenides (TMDs) are essential due to their fascinating electronic and optical properties, strong exciton binding energy, and layer-dependent bandgap. They can be tuned to function as a single-photon emitter, but the quantum yield of photoluminescence of single-layer WSe 2 is low. There is some evidence that a bright emissive state can be produced by introducing local defects through functionalization. In this paper we use spin-polarized periodic density functional theory (DFT) to study the effect of functionalization for the specific case of cyclic carbenes. We find that the simplest of these molecules, cyclopentadiene carbene (Cyc), binds to the surface by a covalent bond leading to a substantial change in the bandgap (1.24 eV compared to 1.64 eV for the pristine surface). Moreover, there are semiflat bands below the Fermi level that originate from σ bonding and the interaction between the lone-pair p orbitals on the carbon of Cyc that bonds with Se on the TMD. Cyc is found to form a Type IIa heterojunction before and after contact with the surface, where the “a” refers to ordering of the energy levels where the TMD levels have the larger splitting and “b” to the reverse. The effect of various electron-withdrawing and electron-donating groups on Cyc is investigated, and it is found that the direct bandgap and heterointerfaces can be chemically tuned with covalently bound functional groups. With an electron-withdrawing group (EWG) such as –CHO, –COCl, and –CN attached to Cyc, a Type III//Type IIa interface is formed, whereas all other EWGs used in this study form a Type IIa junction, before and after contact to the surface. The electron-donating groups (EDG) form Type Ib//Type Ia junctions, and in particular, we find that the system WSe 2 + Cyc-Me band structure consists of semiflat bands at valence band maxima localized on Cyc-Me and the conduction band minimum is coupled between the surface and functional group with a direct bandgap of 0.88 eV. Hence, we predict that the Cyc-Me-functionalized monolayer WSe 2 will produce a red-shifted bright emissive state. Furthermore, we also find that carbenes with the Cyc all-carbon ring have (1) a triplet ground state and (2) form covalent bonds to WSe 2 , while heterocyclic carbenes (1) have a singlet ground state and (2) do not form covalent bonds to WSe 2 .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Hybrid Density Functional Theory Comparison of Oxygen Release and Solvent Decomposition Kinetics on Li x NiO 2 Surfaces

High-nickel-content layered oxides are among the most promising electric vehicle battery cathode materials. However, their interfacial reactivity with electrolytes and tendency toward oxygen release (possibly yielding reactive 1 O 2 ) remain degradation concerns. Elucidating the most relevant (i.e., fastest) interfacial degradation mechanism will facilitate future mitigation strategies. We apply screened hybrid density functional (HSE06) calculations to compare the reaction kinetics of Li x NiO 2 surfaces with ethylene carbonate (EC) with those of O 2 release. Further, on both the (001) and (104) facets, EC oxidative decomposition exhibits lower activation energies than O 2 release. Our calculations, coupled with previously computed liquid-phase reaction rates of 1 O 2 with EC, strongly question the role of “reactive 1 O 2 ” species in electrolyte oxidative degradation. The possible role of other oxygen species is discussed. To deal with the challenges of modeling Li x NiO 2 surface reactivity, we emphasize a “local structure” approach instead of pursuing the global energy minimum.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Multi‐Scale Seismic Imaging of the Ridgecrest, CA, Region With Waveform Inversion of Regional and Dense Array Data

Abstract We develop an inversion procedure for deriving multi‐scale velocity models with waveform inversions of earthquake and ambient noise data at multi‐frequency bands recorded by regional and dense sensor configurations. The method is applied for the area around the 2019 Ridgecrest earthquake rupture zones, utilizing data recorded by regional stations and dense 2D and 1D arrays with station spacings of ∼5 km and ∼100 m, respectively. Starting with regional Vp, Vs models and locations of Ridgecrest aftershocks, the velocity models and event locations are improved iteratively by inversions of waveforms recorded by regional stations and the 2D array, using a minimum spectral element size of ∼600 m. Waveforms from local events recorded by dense 1D arrays across the M7.1 rupture zone with frequencies of up to 10 Hz are used to resolve small‐scale features of the rupture zone and shallow crust with a local spectral element size of 80 m. The refined models provide self‐consistent descriptions of the rupture zone and the shallow crust embedded in the regional structures. The results reveal pronounced low Vs and high Vp/Vs in the M6.4 and M7.1 rupture zones coinciding with concentrations of seismicity, and also around the Garlock fault and in several local basins. We also observe clear velocity contrasts across the Garlock fault with polarity reversals along strike and with depth. The obtained multi‐scale velocity models can be used to improve derivations of earthquake source properties, simulations of dynamic ruptures and ground motions, and the understanding of fault and tectonic processes in the region.

Geochemistry & Geophysics↗
DOE OSTI2024

Diffusion quantum Monte Carlo approach to the polaritonic ground state

Making and using polaritonic states (i.e., hybrid electron-photon states) for chemical applications has recently become one of the most prominent and active fields that connects the communities of chemistry and quantum optics. Modeling of such polaritonic phenomena using ab initio approaches calls for new methodologies, leading to the reinvention of many commonly used electronic structure methods, such as Hartree-Fock, density functional, and coupled cluster theories. Here in this work, we explore the formally exact diffusion quantum Monte Carlo approach to obtain numerical solutions to the polaritonic ground state during the dissociation of the H 2 molecular system. Here, we examine various electron-nuclear-photon properties throughout the dissociation, such as changes to the minimum of the cavity Born-Oppenheimer surface, the localization of the electronic wave function, and the average mode occupation. Finally, we directly compare our results to that obtained with state-of-the-art, yet approximate, polaritonic coupled cluster approaches.

74 ATOMIC AND MOLECULAR PHYSICS↗