Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical continuation”

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 163 records · Page 9

The Impact of Current and Future Polar Orbiting Satellite Data on Numerical Weather Prediction at NASA/GSFC

The lack of adequate observational data continues to be recognized as a major factor limiting both atmospheric research and numerical prediction on a variety of temporal and spatial scales. Since the advent of meteorological satellites in the 1960's, a considerable research effort has been directed toward the design of space-borne meteorological sensors, the development of optimal methods for the utilization of these data, (and an assessment of the influence of existing satellite data and the potential influence of future satellite observations on numerical weather prediction. This has included both Observing System Experiments (OSEs) and Observing System Simulation Experiments (OSSEs). OSEs are conducted to evaluate the impact of specific observations or classes of observations on analyses and forecasts. While OSEs are performed with existing data, OSSEs are conducted to evaluate the potential for future observing systems to improve-NWP, as well as to evaluate trade-offs in observing system design, and to develop and test improved methods for data assimilation. At the conference, results from OSEs to evaluate satellite data sets that have recently become available to the global observing system, such as AIRS and Seawinds, and results from OSSEs to determine the potential impact of space-based lidar winds will be presented.

Atlas, Robert↗

Orientation reversal and the Chern-Simons natural boundary

We show that the fundamental property of preservation of relations, underlying resurgent analysis, provides a new perspective on crossing a natural boundary, an important general problem in theoretical and mathematical physics. This reveals a deeper rigidity aspect of resurgence in a quantum field theory path integral. The physical context here is the non-perturbative completion of complex Chern-Simons theory that associates to a 3-manifold a collection of q-series invariants labeled by Spinc structures, for which crossing the natural boundary corresponds to orientation reversal of the 3-manifold. Our new resurgent perspective leads to a practical numerical algorithm that generates q-series which are dual to unary q-series composed of false theta functions. Until recently, these duals were only known in a limited number of cases, essentially based on Ramanujan’s mock theta functions, and the common belief was that the duals might not even exist in the general case. Resurgence analysis identifies as primary objects Mordell integrals: up to changes of variables, they are Laplace transforms of resurgent functions. Their unique Borel summed transseries decomposition on either side of the Stokes line is simply the unique decomposition into real and imaginary parts. In turn, the latter are combinations of unary q-series in terms of q and its modular counterpart $\overset{\sim }{q}$ , and are resurgent by construction. The Mordell integral is analytic across the natural boundary of the q and $\overset{\sim }{q}$ series, and uniqueness of a similar decomposition which preserves algebraic relations on the other side of the boundary defines the unique boundary crossing of the q series. We demonstrate that this continuation can be efficiently implemented numerically. In the cases where unique mock modular identities are known, they are found by this numerical procedure, but the procedure can go well beyond the known list of identities. A particularly interesting feature of the resurgent approach is that it reveals new aspects, and is very different from other known approaches based on indefinite theta series, Appell-Lerch sums, and representation theory of logarithmic vertex operator algebras.

Chern-Simons theories↗

A New Compendium of Unsteady Aerodynamic Test Cases for CFD: Summary of AVT WG-003 Activities

With the continuous progress in hardware and numerical schemes, Computational Unsteady Aerodynamics (CUA), that is, the application of Computational Fluid Dynamics (CFD) to unsteady flowfields, is slowly finding its way as a useful and reliable tool (turbulence and transition modeling permitting) in the aircraft, helicopter, engine and missile design and development process. Before a specific code may be used with confidence it is essential to validate its capability to describe the physics of the flow correctly, or at least to the level of approximation required, for which purpose a comparison with accurate experimental data is needed. Unsteady wind tunnel testing is difficult and expensive; two factors which dramatically limit the number of organizations with the capability and/or resources to perform it. Thus, unsteady experimental data is scarce, often classified and scattered in diverse documents. Additionally, access to the reports does not necessarily assure access to the data itself. The collaborative effort described in this paper was conceived with the aim of collecting into a single easily accessible document as much quality data as possible. The idea is not new. In the early 80's NATO's AGARD (Advisory Group for Aerospace Research & Development) Structures and Material Panel (SMP) produced AGARD Report No. 702 "Compendium of Unsteady Aerodynamic Measurements", which has found and continues to find extensive use within the CUA Community. In 1995 AGARD's Fluid Dynamics Panel (FDP) decided to update and expand the former database with new geometries and physical phenomena, and launched Working Group WG-22 on "Validation Data for Computational Unsteady Aerodynamic Codes". Shortly afterwards AGARD was reorganized as the RTO (Research and Technology Organization) and the WG was renamed as AVT (Applied Vehicle Technolology) WG-003. Contributions were received from AEDC, BAe, DLR, DERA, Glasgow University, IAR, NAL, NASA, NLR, and ONERA. The final publication with the results of the exercise is expected in the second part of 1999. The aim of the present paper is to announce and present the new database to the Aeroelasticity community. It is also intended to identify, together with one of the groups of end users it targets, deficiencies in the compendium that should be addressed by means of new wind tunnel tests or by obtaining access to additionally existing data.

Ruiz-Calavera, Luis P.↗

MODIS Reflective Solar Bands On-Orbit Calibration and Performance

The design of the MODIS instrument was driven by the scientific community's desire to have near-daily global coverage at moderate resolution (1 km) with a comprehensive spectral coverage from visible to long-wave infrared wavelengths. Since their launches in 1999 and 2002, respectively, the Terra and Aqua MODIS instruments have made continuous global observations and generated numerous data products to help users worldwide with their studies of the Earth's system and its short- and long-term changes. The 20 reflective solar bands (RSB) with wavelengths from 0.41 to 2.2 m collect data at three nadir spatial resolutions: 250 m, 500 m, and 1 km. The solar diffuser (SD) coupled with the solar diffuser stability monitor (SDSM) provide a reflectance-based calibration on-orbit. In addition, lunar observations and response trends from pseudo-invariant desert sites are used to characterize the response versus scan- angle changes on-orbit. This paper provides a brief overview of MODIS RSB calibration algorithms, as implemented in the latest Level 1B version 6.1, operational activities, on-orbit performance, remaining challenges and potential improvements. Results from the SD and SDSM measurements show a wavelength and mirror side dependent degradation in RSB responses, with the largest degradation at the shortest wavelengths, particularly for Terra MODIS. Aqua MODIS has experienced far less degradation of its optics and on-board calibrators compared to Terra MODIS, resulting in an overall better performance. With the exception of Aqua band 6, there have been no new noisy or inoperable detectors in the RSB of either instrument during post-launch operations. As the instruments age and continue to endure the space environment, the detectors and optical systems degrade. The challenges associated with incorporating these on-orbit changes to ensure a production of high-quality calibrated L1B data products are also discussed in this paper.

solar diffuser stability monitor↗

Thermodynamic analysis of journal bearings operating under steady state loading in laminar regime

Thermohydrodynamic effects in journal bearings operating under steady load in laminar regime are investigated. An analytical model for the finite and infinitely long journal bearings is formulated. The model includes correction factors for the cavitation effects in the unloaded region of the bearing and the mixing of the recirculating oil and supply oil at the oil inlet. A finite difference computer program is developed to numerically solve the governing equations of the continuity, Reynolds, energy, Laplace heat conduction, and a viscosity-temperature relation simultaneously. The program includes a numerical technique for obtaining an isothermal shaft temperature. The numerical results of temperature distribution and the heat effects on the bearing load carrying capacity agree closely with those of experimental findings. Several different sets of simpler boundary conditions for the energy equation are studied.

Khonsari, M. M.↗

Reynold stress closure in jet flows using wave models

Research program efforts have continued to concentrate on the development of the numerical methods that will form the computational part of the turbulence closure scheme. Studies have continued on the wave model for the two dimensional shear layer. This configuration is being used as a test case for the closure schemes. Several numerical schemes for the solution of the non-separable Rayleigh equation were developed. This solution is required for the closure scheme in more complex geometries. The most efficient method found is a Hybrid scheme that combines both pseudospectral and finite difference techniques. In addition, conformal transformation techniques were developed to transform the arbitrary geometry of the jet to a simple computational domain. The study of the shock structure in arbitrary geometry jets and multiple jets. These developments are described briefly.

Morris, P. J.↗

A flexible and differentiable coil proxy for stellarator equilibrium optimization

Balancing plasma performance and coil cost is a significant challenge when designing a stellarator power plant. Most current stellarator designs are produced through two-stage optimization: stage-1 for the equilibrium and stage-2 for a coil design that reproduces its magnetic configuration. Because few proxies connect both stages, two-stage optimization can produce plasmas that have high-quality physical properties but overly complex coils. In recent years, single-stage optimization has increasingly been used to optimize the plasma and coils simultaneously in order to improve the plasma–coil balance. However, all existing single-stage tools are specialized for filament coils, cannot model coil systems containing permanent magnets or dipole arrays, and continue to be challenged by numerical problems. The quasi-single-stage optimization finds a middle ground by integrating a coil optimization subproblem into stage-1 optimization. We present a flexible, differentiable coil complexity proxy based on the newly developed QUADCOIL coil optimization code. QUADCOIL is fast and can target realistic coil metrics and constraints that are unavailable to codes with comparable speed. We demonstrate the effectiveness and flexibility of the QUADCOIL proxy by presenting two QSS optimization studies. The first study produces a permanent magnet solution for the MUSE stellarator with 34% fewer magnets than previous solutions. The second study produces a coil solution for the ARIES-CS stellarator with 30% reductions in both peak and root-mean-square force.

Fu, Lanke [Princeton Plasma Physics Laboratory (PP↗

Broken Luttinger theorem in the two-dimensional Fermi-Hubbard model

One of the fundamental questions about high-temperature cuprate superconductors is the size of the Fermi surface underlying the superconducting state. By analyzing the single-particle spectral function for the Fermi-Hubbard model as a function of repulsion U and chemical potential μ, we find that the Fermi surface in the normal state undergoes a transition from a large Fermi surface matching the Luttinger volume as expected in a Fermi liquid, to a Fermi surface that encloses fewer electrons that we dub the “Luttinger breaking” phase, as the Mott insulator is approached. This transition into a non-Fermi-liquid phase that violates the Luttinger count occurs at a critical density in the absence of any other broken symmetry. We obtain the Fermi-surface contour from the spectral weight A k (ω = 0) and from an analysis of the singularities of the Green's function Re G k (E = 0), calculated using determinantal quantum Monte Carlo and analytic continuation methods. We discuss our numerical results in connection with experiments on Hall measurements, scanning tunneling spectroscopy, and angle-resolved photoemission spectroscopy.

36 MATERIALS SCIENCE↗

Degradation of temporal contrast from post-pedestal interference with a chirped pulse in an optical parametric amplifier

Pre-pedestal generation is observed in a 0.35-PW laser front end coming from a post-pedestal via instantaneous gain and pump depletion in an optical parametric amplifier during chirped-pulse amplification. Generalized simulations show how this effect arises from gain nonlinearity and applies to all optical parametric chirped-pulse–amplification systems with a post-pedestal. An experiment minimizing the effect of B -integral is used to isolate and study the newly observed conversion of a continuous post-pedestal into a continuous pre-pedestal. Matching numerical simulations confirm experimental results and additionally reveal how third-order dispersion largely controls the slope of the generated pre-pedestal.

47 OTHER INSTRUMENTATION↗

Air Force Institute of Technology Dust Cloud Model Vertical Dispersion

The Air Force Institute of Technology (AFIT) has developed and refined a nuclear dust cloud model for modeling nuclear fallout, and the dust mass and radiation encountered by aircraft flying through a nuclear dust cloud. Recently, a question was raised requesting an explanation of why the vertical mass and activity distributions had an anomalous shape for yields above 1 Mt. This report explains that the anomalous shape occurs because the dust cloud model is a discrete numeric model, and not a continuous model.

54 ENVIRONMENTAL SCIENCES↗

Microgrid Assisted Design for Remote Areas

In this work, we present a three-stage multiobjective mixed-integer linear programming (MILP) for the optimal expansion planning and operation of isolated multienergy microgrids in remote areas. By selecting the optimal distributed generators (DGs) and energy storage systems (ESSs) mix selection, siting, sizing, and scheduling in the remote microgrid, the proposed model is targeted to minimize the annualized total cost of microgrids while enhancing the performance of the system, i.e., minimizing the voltage deviations and line power loss. To represent the electricity and heat flow between generation resources and various electrical, heating, and cooling loads in the isolated microgrid, linearized power flow, and heat flow constraints are employed in the proposed optimization model. The available capacity of DGs and ESSs are modeled as discrete constants instead of continuous variables for practical purpose. Numerical simulation results on a remote microgrid consisting of DGs, ESSs, and various loads validate the proposed method.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Calculation of turbulent shear stress in supersonic boundary layer flows

Turbulent shear stress distributions for supersonic boundary layer flows have been computed from experimental mean boundary layer data. The computations have been made by numerically integrating the time averaged continuity and streamwise momentum equations. Distributions have been obtained for flows upstream and downstream of shock-wave-boundary layer interactions and for both two-dimensional and axisymmetric flows. The computed results are compared with recently reported shear stress measurements which were obtained by hot wire anemometer and laser velocimeter techniques.

Sun, C. C.↗

Seasonal and magnetic storm related changes in the thermosphere induced by eddy mixing

The possibility of explaining the seasonal and magnetic storm related changes in the thermosphere and the ionosphere through variation in eddy diffusion coefficient in the lower thermosphere is investigated theoretically by obtaining simultaneous numerical solutions of the relevant continuity, momentum, and energy balance equations in a self-consistent manner. It is shown that various important features of thermospheric seasonal behavior, including the winter helium bulge phenomenon, can be explained by assuming effective reduction in the winter-time eddy diffusion coefficient by about one order of magnitude. Storm related changes in the thermal structure and composition are seen to arise by assuming a downward shift of turbulence. Physical implications and relative significance of variation in thermospheric turbulence and large scale inter-hemispheric circulation are briefly discussed.

Sinha, A. K.↗

The NASA Ames Research Center one- and two-dimensional stratospheric models. Part 1: The one-dimensional model

A one-dimensional model of stratospheric trace constituents is described in detail. Specifically, the numerical solution of the species continuity equations, including a technique for treating the stiff differential equations representing the chemical kinetic terms, and an appropriate method for simulating the diurnal variations of the species concentrations, are discussed. A specialized treatment of atmospheric photodissociation rates is outlined in the text. The choice of a vertical eddy diffusivity profile and its success in predicting the vertical tracer distributions (carbon 14, methane, and nitrous oxide) are also discussed.

Turco, R. P.↗

The NASA Ames Research Center one- and two-dimensional stratospheric models. Part 2: The two-dimensional model

The two-dimensional model of stratospheric constituents is presented in detail. The derivation of pertinent transport parameters and the numerical solution of the species continuity equations, including a technique for treating the stiff differential equations that represent the chemical kinetic terms, and appropriate methods for simulating the diurnal variations of the solar zenith angle and species concentrations are discussed. Predicted distributions of tracer constituents (ozone, carbon 14, nitric acid) are compared with observed distributions.

Whitten, R. C.↗

An analytic formulation for heating source memory in the thermospheric composition

Numerical calculations from a spectral circulation model are utilized to construct an analytic Green's function formulation describing the meridional, time-dependent thermospheric composition and temperature response during magnetic storms. The purpose is to develop a formulation that embodies source memory while being sufficiently simple to serve as a heuristic guide for empirical modeling. By passing from the discrete Fourier series representation, utilized for the numerical circulation model, to a continuous Fourier integral representation, explicit waves are obtained for the thermospheric response times. The response times are altitude and species dependent and can exceed two days below 200 km. Thus, for certain storm scenarios, pronounced source memory signatures for the composition and temperature are predicted. Response times obtained from the formulation are shown to give a response consistent with previously published neutral composition data from AE-C for the February 1974 storm when an ap dependent heat source is employed.

Porter, H. S.↗