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.

SourceResearch connectionFrom yearThrough year

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

At least 1,387 records · Page 77

All sources NASA NTRS DOE OSTI

NASA NTRS1990

Numerical study of shock-wave/boundary layer interactions in premixed hydrogen-air hypersonic flows

A computational study of shock wave/boundary layer interactions involving premixed combustible gases, and the resulting combustion processes is presented. The analysis is carried out using a new fully implicit, total variation diminishing (TVD) code developed for solving the fully coupled Reynolds-averaged Navier-Stokes equations and species continuity equations in an efficient manner. To accelerate the convergence of the basic iterative procedure, this code is combined with vector extrapolation methods. The chemical nonequilibrium processes are simulated by means of a finite-rate chemistry model for hydrogen-air combustion. Several validation test cases are presented and the results compared with experimental data or with other computational results. The code is then applied to study shock wave/boundary layer interactions in a ram accelerator configuration. Results indicate a new combustion mechanism in which a shock wave induces combustion in the boundary layer, which then propagates outwards and downstream. At higher Mach numbers, spontaneous ignition in part of the boundary layer is observed, which eventually extends along the entire boundary layer at still higher values of the Mach number.

Yungster, Shaye↗

NASA NTRS1991

Numerical study of shock-wave/boundary layer interactions in premixed hydrogen-air hypersonic flows

Yungster, Shaye↗

NASA NTRS2001

The Continuing Evolution of Land Surface Parameterizations

Land surface models (LSMs) play a critical role in the simulation of climate, for they determine the character of a large fraction of the atmosphere's lower boundary. The LSM partitions the net radiative energy at the land surface into sensible heat, latent heat, and energy storage, and it partitions incident precipitation water into evaporation, runoff, and water storage. Numerous modeling experiments and the existing (though very scant) observational evidence suggest that variations in these partitionings can feed back on the atmospheric processes that induce them. This land-atmosphere feedback can in turn have a significant impact on the generation of continental precipitation. For this and other reasons (including the role of the land surface in converting various atmospheric quantities, such as precipitation, into quantities of perhaps higher societal relevance, such as runoff), many modeling groups are placing a high emphasis on improving the treatment of land surface processes in their models. LSMs have evolved substantially from the original bucket model of Manabe et al. This evolution, which is still ongoing, has been documented considerably. The present paper also takes a look at the evolution of LSMs. The perspective here, though, is different - the evolution is considered strictly in terms of the 'balance' between the formulations of evaporation and runoff processes. The paper will argue that a proper balance is currently missing, largely due to difficulties in treating subgrid variability in soil moisture and its impact on the generation of runoff.

Koster, Randal↗

DOE OSTI2023

Calculation of Dynamical Response Functions Using a Bound-State Method

Abstract We investigate a method to extract response functions (dynamical polarisabilities) directly from a bound-state approach applied to calculations of perturbation-induced reactions. The use of a square-integrable basis leads to a response in the form of a sum of $$\delta $$ δ functions. We integrate this over energy and fit a smooth function to the resulting stepwise-continuous one. Its derivative gives the final approximation to the physical response function. We show that the method reproduces analytical results where known, and analyse the details for a variety of models. We apply it to some simple models, using the stochastic variational method as the numerical method. Albeit we find that this approach, and other numerical techniques, have some difficulties with the threshold behavior in coupled-channel problems with multiple thresholds, its stochastic nature allows us to extract robust results even for such cases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2025

Generative learning for slow manifolds and bifurcation diagrams

In dynamical systems characterized by separation of time scales, the approximation of so called “slow manifolds”, on which the long term dynamics lie, is a useful step for model reduction. Initializing on such slow manifolds is a useful step in modeling, since it circumvents fast transients, and is crucial in multiscale algorithms (like the equation-free approach) alternating between fine scale (fast) and coarser scale (slow) simulations. In a similar spirit, when one studies the infinite time dynamics of systems depending on parameters, the system attractors (e.g., its steady states) lie on bifurcation diagrams (curves for one-parameter continuation, and more generally, on manifolds in state parameter space. Sampling these manifolds gives us representative attractors (here, steady states of ODEs or PDEs) at different parameter values. Algorithms for the systematic construction of these manifolds (slow manifolds, bifurcation diagrams) are required parts of the “traditional” numerical nonlinear dynamics toolkit. In more recent years, as the field of Machine Learning develops, conditional score-based generative models (cSGMs) have been demonstrated to exhibit remarkable capabilities in generating plausible data from target distributions that are conditioned on some given label. It is tempting to exploit such generative models to produce samples of data distributions (points on a slow manifold, steady states on a bifurcation surface) conditioned on (consistent with) some quantity of interest (QoI, observable). In this work, we present a framework for using cSGMs to quickly (a) initialize on a low-dimensional (reduced-order) slow manifold of a multi-time-scale system consistent with desired value(s) of a QoI (a “label”) on the manifold, and (b) approximate steady states in a bifurcation diagram consistent with a (new, out-of-sample) parameter value. This conditional sampling can help uncover the geometry of the reduced slow-manifold and/or approximately “fill in” missing segments of steady states in a bifurcation diagram. Finally, the quantity of interest, which determines how the sampling is conditioned, is either known a priori or identified using manifold learning-based dimensionality reduction techniques applied to the training data.

Dynamical systems↗

NASA NTRS1980

Modeling of growth, evaporation and sedimentation effects on transmission of visible and IR laser beams in artificial fogs

The dense polydisperse aerosol particles in a quiet chamber may spontaneously go through different microphysical processes including gravitational sedimentation, thermal coagulation, and growth or evaporation. In an earlier paper, we presented the results of a parametric study of the combined and separate effects of thermal coagulation and sedimentation on the time dependence of extinction of four visible and IR laser beams traversing an aerosol medium. As a continuation of this series of studies, the separate and combined effects of growth or evaporation and gravitational sedimentation on the time dependence of extinction of the same four visible and IR laser beams traversing in artificial fogs will be reported in this paper. The method of numerically modeling the change of water droplet size distribution with time due to growth/evaporation and the cutoff of larger aerosols due to gravitational sedimentation is described in detail. Factors governing the relative importance of these two processes are discussed. Results of this study show that the relative humidity or ambient temperature is a crucial parameter in determining the optical depth of the water droplet and aerosol media undergoing microphysical processes.

Yue, G. K.↗

NASA NTRS1983

Helicopter rotor loads using discretized matched asymptotic expansions

The numerical practicality of a matched asymptotic expansion approach for the computation of unsteady three dimensional airloads on a helicopter rotor was improved. This effort utilizes a discretized repesentation of the doublet strength distribution and helical streamlines to decrease the computational requirements of the original analysis. The continuous variation of the doublet strength was approximated by piecewise constant or piecewise quadratic distributions, and the helical trajectory of a fluid particle was approximated by connected straight line segments. As a direct result of these simplified representations the computational time required for the execution of a typical flight condition was reduced by an order of magnitude with respect to the requirements of the original analysis. Airloads which were computed using the discretized method for a two bladed model rotor and a full scale four bladed rotor are in close agreement with measured results and airloads from the original asymptotic analysis. For conditions characterized by significant rotor/wake interaction the piecewise constant representation requires a reduced azimuth spacing to maintain acceptable accuracy.

Pierce, G. A.↗

NASA NTRS1984

Numerical simulation of reconnection in an emerging magnetic flux region

Numerical solutions in two dimensions are presented for the resistive MHD equations of an initial boundary value problem, simulating reconnection between an emerging magnetic flux region and an overlying coronal magnetic field. This numerical model displays four main phases, which are interpreted in terms of (1) a slowly evolving quasi-steady phase during which most of the magnetic flux emerges, with approximate equilibrium between magnetic and pressure forces; (2) an impulsive phase in which either the reconnection or continuing emergence of the first phase comes to disrupt the equilibrium, and extensive acceleration occurs as the high pressure region within the emerging region drives the fluid upwards and outwards; (3) a second quasi-steady phase; and (4) a potential-static phase in which continuing reconnection, ohmic dissipation, and fluid transport through boundaries depletes the system of all currents and flows.

Forbes, T. G.↗

NASA NTRS1986

Transport of plasma across the magnetopause

It is here demonstrated numerically, by tracing test particle motions, that solar wind plasma particles can cross the magnetopause and get into the boundary layer if the tangential electric field in the magnetopause rest frame changes in concert with changes in the magnetic field as E = V x B for continuous earthward flow, whether or not the field lines are closed. Such a simultaneous change in the electric and magnetic fields is proposed as a key feature of flux transfer events (FTEs). It is argued that the FTE process could more aptly be described as a transient or impulsive transfer of plasma particles on a single flux tube over a localized region of the magnetopause, generally ending up on different flux tubes after the crossing. Thus, it is proposed that FTE should denote flux 'trade' event, rather than 'transfer'.

Heikkila, Walter J.↗

NASA NTRS1989

Reflection and trapping of transient Alfven waves propagating in an isothermal atmosphere with constant gravity and uniform magnetic field

A time-dependent linear magnetohydrodynamic numerical model was used to investigate the propagation of Alfven waves in an isothermal and stratified atmosphere with constant gravity and uniform vertical magnetic field. Results show that the Alfven wave transit time for the wave source to infinity is finite and that the wave exhibits continuous partial reflection which becomes total reflection as the front approaches infinity. The total reflection causes the waves to be trapped in the cavity that extends from the wave source to infinity and in which the wave energy is stored. The results suggest that the reflection of Alfven waves (of sufficiently long period) from the outer corona is an intrinsic phenomenon for any stellar atmosphere stratified by gravity and an open magnetic field, and that, therefore, such waves may be trapped in the stellar atmosphere.

An, C.-H.↗

NASA NTRS1991

Predicting chaos for infinite dimensional dynamical systems: The Kuramoto-Sivashinsky equation, a case study

The results of extensive computations are presented in order to accurately characterize transitions to chaos for the Kuramoto-Sivashinsky equation. In particular, the oscillatory dynamics in a window that supports a complete sequence of period doubling bifurcations preceding chaos is followed. As many as thirteen period doublings are followed and used to compute the Feigenbaum number for the cascade and so enable, for the first time, an accurate numerical evaluation of the theory of universal behavior of nonlinear systems, for an infinite dimensional dynamical system. Furthermore, the dynamics at the threshold of chaos exhibit a fractal behavior which is demonstrated and used to compute a universal scaling factor that enables the self-similar continuation of the solution into a chaotic regime.

Smyrlis, Yiorgos S.↗

NASA NTRS1991

Convection in vertical Bridgman configurations

We are interested in tracing the convective profiles of vertical Bridgman growth in two configurations, the pure Rayleigh convective mode and the combined Rayleigh-Marangoni mode. In order to do so, we conducted a numerical investigation that involved a finite volume calculation. The governing equations were integrated about a cell volume, using the Gauss Theorem and the volume variables like temperature and velocity were related to the surface variables. In order to solve for the pressure field, we employed the continuity equation and the residuals resulted in a Poisson equation. Results and comments for the Rayleigh and Marangoni problems in a vertical cylinder or Bridgman configuration are given.

Narayanan, Ranga↗

NASA NTRS1991

A far wing line shape theory and its application to the water continuum absorption in the infrared region. I

The present theory for the continuous absorption that is due to the far-wing contribution of allowed lines is based on the quasistatic approximation for the far wing limit and the binary collision approximation of one absorber molecule and one bath molecule. The validity of the theory is discussed, and numerical results of the water-continuum absorption in the IR region are presented for comparison with experimental data. Good agreement is obtained for both the magnitude and temperature dependence of the absorption coefficients.

Ma, Q.↗

NASA NTRS2001

The Application of Time-Frequency Methods to HUMS

This paper reports the study of four time-frequency transforms applied to vibration signals and presents a new metric for comparing them for fault detection. The four methods to be described and compared are the Short Time Frequency Transform (STFT), the Choi-Williams Distribution (WV-CW), the Continuous Wavelet Transform (CWT) and the Discrete Wavelet Transform (DWT). Vibration data of bevel gear tooth fatigue cracks, under a variety of operating load levels, are analyzed using these methods. The new metric for automatic fault detection is developed and can be produced from any systematic numerical representation of the vibration signals. This new metric reveals indications of gear damage with all of the methods on this data set. Analysis with the CWT detects mechanical problems with the test rig not found with the other transforms. The WV-CW and CWT use considerably more resources than the STFT and the DWT. More testing of the new metric is needed to determine its value for automatic fault detection and to develop methods of setting the threshold for the metric.

Pryor, Anna H.↗

NASA NTRS2010

An Overview of the Micro Pulse Lidar Network (MPLNET)

The NASA Micro Pulse Lidar Network (MPLNET) is a federated network of Micro Pulse Lidar (MPL) systems designed to measure aerosol and cloud vertical structure continuously, day and night, over long time periods required to contribute to climate change studies and provide ground validation for models and satellite sensors in the NASA Earth Observing System (FOS). At present, there are eighteen active sites worldwide, and several more in the planning stage. Numerous temporary sites are deployed in support of various field campaigns. Most MPLNET sites are co-located with sites in the NASA Aerosol Robotic Network (AERONET) to provide both column and vertically resolved aerosol and cloud data. MPLNET data and more information on the project are available at http://mpinet.gsfc.nasa.gov . Here we present a summary of the first ten years of MPLNET, along with an overview of our current status, specifically our version two data products and applications. Future network plans will be presented, with a focus on our activities in South East Asia.

Welton, Ellsworth↗

NASA NTRSDate not supplied

Demonstration of Linear Covariance Analysis Techniques to Evaluate Entry Descent and Landing Guidance Algorithms, Vehicle Configurations, Analysis Techniques, and Trajectory Profiles

Linear covariance analysis techniques have been previously developed to analyze closed-loop entry, descent, and landing (EDL) scenarios and the initial validation efforts are under-way confirming the generated GN&C system performance results. Given both the theoretical foundation and previous conceptual demonstration, this work begins to flex the potential of linear covariance analysis for atmospheric flight and highlight its versatility and reliability by evaluating multiple entry guidance algorithms, vehicle configurations, trajectory profiles, environment conditions, and analysis techniques for a variety of trade studies. To demonstrate the benefit linear covariance analysis can provide in producing rapid yet accurate performance data, two entry profiles are adopted including the NASA Mars Science Laboratory (MSL) and Exploration Flight Test-1 (EFT-1) while utilizing two different guidance algorithms, the Apollo Final Phase (AFP) and the Fully Numeric Predictor-Corrector Entry Guidance (FNPEG) with different navigation sensor suites in a 6 degree-of-freedom (6-DOF) simulation environment. Results are shown using both linear covariance and Monte Carlo analysis techniques to high-light the consistency between the two methodologies and continue the validation maturation of linear covariance analysis for entry, descent, and landing.

EDL↗

NASA NTRS1973

Numerical simulation of noise propagation in jet engine ducts

A finite difference formulation is presented which is useful in the study of acoustically treated inlet and exhaust ducts used in turbofan engines. The difference formulation can readily handle acoustic flow field complications, such as axial variations in wall impedance and cross-sectional area, that would occur in a sonic inlet. In formulating the difference solutions, the continuous acoustic field is lumped into a series of grid points spread uniformly throughout the field. At each point, the pressure is separated into its real and imaginary terms. Example solutions are presented for sound propagation in a one-dimensional straight hard-wall duct and in a two-dimensional straight soft-wall duct without steady flow.

Baumeister, K. J.↗

DOE OSTI2021

Adjoint DSMC for nonlinear Boltzmann equation constrained optimization

Applications for kinetic equations such as optimal design and inverse problems often involve finding unknown parameters through gradient-based optimization algorithms. Based on the adjoint-state method, we derive two different frameworks for approximating the gradient of an objective functional constrained by the nonlinear Boltzmann equation. While the forward problem can be solved by the DSMC method, it is difficult to efficiently solve the high-dimensional continuous adjoint equation obtained by the “optimize-then-discretize” approach. This challenge motivates us to propose an adjoint DSMC method following the “discretize-then-optimize” approach for Boltzmann-constrained optimization. We also analyze the properties of the two frameworks and their connections. Here, several numerical examples are presented to demonstrate their accuracy and efficiency.

97 MATHEMATICS AND COMPUTING↗