Results for “Numerical stability” — Page 40

NASA NTRS1998

Model Refinement Using Eigensystem Assignment

This paper presented a novel approach for the refinement of finite-element-based analytical models of flexible structures is presented. The proposed approach models the possible refinements in the mass, damping, and stiffness matrices of the finite element model in the form of a constant gain feedback with acceleration, velocity, and displacement measurements, respectively. Once, the free elements of the structural matrices have been defined, the problem of model refinement reduces to obtaining position, velocity, and acceleration gain matrices, which reassign a desired subset of the eigenvalues of the model, along with partial mode shapes, from their baseline values to those obtained from system identification test data. A sequential procedure is used to assign one self-conjugate pair of closed-loop eigenvalues at each step using symmetric output feedback gain matrices, and the closed-loop eigenvectors are partially assigned, while ensuring that the eigenvalues assigned in the previous steps are not disturbed. The procedure can also impose that gain matrices be dissipative in order to guarantee the stability of the refined model. A numerical example, involving finite element model refinement for a structural testbed at NASA Langley (CSI Evolutionary Model) is presented to demonstrate the feasibility of the proposed approach.

Maghami, Pieman G.↗

NASA NTRS1995

Effect of Far-Field Boundary Conditions on Boundary-Layer Transition

The effect of far-field boundary conditions on the evolution of a finite-amplitude two-dimensional wave in the Blasius boundary layer is assessed. With the use of the parabolized stability equations (PSE) theory for the numerical computations, either asymptotic, Dirichlet, Neumann or mixed boundary conditions are imposed at various distances from the wall. The results indicate that asymptotic and mixed boundary conditions yield the most accurate mean-flow distortion and unsteady instability modes in comparison with the results obtained with either Dirichlet or Neumann conditions.

Bertolotti, Fabio P.↗

NASA NTRS1994

Effect of Far-Field Boundary Conditions on Boundary-Layer Transition

The effect of far-field boundary conditions on the evolution of a finite-amplitude two-dimensional wave in the Blasius boundary layer is assessed. With the use of the parabolized stability equations (PSE) theory for the numerical computations, either asymptotic, Dirichlet, Neumann or mixed boundary conditions are imposed at various distances from the wall. The results indicate that asymptotic and mixed boundary conditions yield the most accurate mean-flow distortion and unsteady instability modes in comparison with the results obtained with either Dirichlet or Neumann conditions.

Bertolotti, Fabio P.↗

NASA NTRS1992

Activation of maternal centrosomes in unfertilized sea urchin eggs

Centrosomes are undetectable in unfertilized sea urchin eggs, and normally the sperm introduces the cell's microtubule-organizing center (MTOC) at fertilization. However, artificial activation or parthenogenesis triggers microtubule assembly in the unfertilized egg, and this study explores the reappearance and behavior of the maternal centrosome. During activation with A23187 or ammonia, microtubules appear first at the cortex; centrosomal antigen is detected diffusely throughout the entire cytoplasm. Later, the centrosome becomes more distinct and organizes a radial microtubule shell, and eventually a compact centrosome at the egg center organizes a monaster. In these activated eggs, centrosomes undergo cycles of compaction and decompaction in synchrony with the chromatin, which also undergoes cycles of condensation and decondensation. Parthenogenetic activation with heavy water (50% D2O) or the microtubule-stabilizing drug taxol (10 microM) induces numerous centrosomal foci in the unfertilized sea urchin egg. Within 15 min after incubation in D2O, numerous fine centrosomal foci are detected, and they organize a connected network of numerous asters which fill the entire egg. Taxol induces over 100 centrosomal foci by 15 min after treatment, which organize a corresponding number of asters. The centrosomal material in either D2O- or taxol-treated eggs aggregates with time to form fewer but denser foci, resulting in fewer and larger asters. Fertilization of eggs pretreated with either D2O or taxol shows that the paternal centrosome is dominant over the maternal centrosome. The centrosomal material gradually becomes associated with the enlarged sperm aster. These experiments demonstrate that maternal centrosomal material is present in the unfertilized egg, likely as dispersed undetectable material, which can be activated without paternal contributions. At fertilization, paternal centrosomes become dominant over the maternal centrosomal material.

Non-NASA Center↗

NASA NTRS2018

Nonlinear Evolution and Breakdown of Azimuthally Compact Crossflow Vortex Pattern over a Yawed Cone

Hypersonic boundary-layer flows over a circular cone at moderate angle of incidence can support strong crossflow instability in between the windward and leeward rays on the plane of symmetry. Due to the more efficient excitation of stationary crossflow vortices by surface roughness, a possible path to transition in such flows corresponds to rapid amplification of the high-frequency secondary instabilities of finite amplitude stationary crossflow vortices. In the present paper, the previous analyses of crossflow instability over a 7- degree half-angle, yawed circular cone in a Mach 6 free stream have been extended to the nonlinear evolution of azimuthally localized crossflow vortex packets and the amplification characteristics and nonlinear breakdown of high-frequency secondary instabilities associated with those packets. A comparison between plane marching PSE and direct Navier-Stokes simulations (DNS) reveals favorable agreement in regard to mode shapes, most amplified disturbance frequencies, and N-factor evolution. In contrast, the quasi-parallel predictions are found to result in severe underprediction of the N-factors. The direct numerical simulations also indicate that the breakdown of secondary instabilities in a 3D hypersonic boundary layer shares certain common features with the previous computations of crossflow transition over subsonic swept wings.

Vortex structure↗

NASA NTRS2020

Streak Instabilities on HIFiRE-5 Elliptic Cone

The Hypersonic International Flight Research and Experimentation (HIFiRE) 5 flight experiment by the Air Force Research Laboratories and Australian Defense Science and Technology Organization yielded inflight, hypersonic transition data for a canonical, 3D configuration in the form of an elliptic cone at various combinations of angle-of-attack and yaw. Previous analysis has shown that boundary-layer stability calculations for the as-flown, zero-yaw flight conditions are able to capture the second mode dominated transition along the attachment line as well as the crossflow induced transition in between the attachment line and the minor axis symmetry plane. However, there is still a need to understand the transition in the vicinity of the minor axis plane of symmetry, where the boundary layer thickness exhibits an azimuthally compact bulge due to the collision of the secondary flow approaching from the attachment lines on both sides. This paper reports a computational analysis of the stability of the boundary layer flow along the bulge regions on both the leeward and windward sides of the vehicle surface. At a representative flight condition corresponding to a time instant of t = 514.813 seconds, computations are performed to identify several different unstable modes within the bulge regions, with the eventual goal of developing a linear stability correlation for the transition onset near the minor axis symmetry plane. The most amplified flow instabilities in the pretransitional region correspond to a “stem” mode that is concentrated within a pair of nearly vertical shear layers below the cap of the mushroom shaped bulge in the mean velocity contours. An excellent agreement is found between the predictions of instability amplification based on the plane marching parabolized stability equations (PSE) and the direct numerical solutions (DNS). Given the large number of unstable modes within the bulge region, the observed agreement demonstrates, for first time, that the PSE can provide accurate predictions for the development of flow instabilities in streak-like flow features over three-dimensional bodies. The peak N-factor of approximately 15 at the transition location estimated from the thermocouple measurements along the leeward centerline is close to the transition N-factor for Mack mode transition during the HIFiRE-1 circular cone flight experiment with nearly zero degrees angle of attack.

Meelan Choudhari↗

NASA NTRS2021

A Computational Analysis of Boundary Layer Instability over the BOLT Configuration

The complex boundary layer flow over the BOLT flight configuration is known to exhibit multiple and potentially interacting instability mechanisms. This paper represents a continuation of our numerical investigation of the flow instabilities over the main test surface of the BOLT configuration by using state-of-the-art tools in multidimensional instability analysis. Specifically, the paper extends our previous computations by considering the separate effects of a nonzero angle-of-attack and a nonzero yaw on the modal instability characteristics of the boundary layer streaks adjacent to the minor-axis symmetry plane. At a selected condition from the ascent part of the anticipated flight trajectory, either type of departure from the design condition is shown to have a considerable impact on the structure of the laminar basic state and, hence, also on the amplification characteristics of the instability waves. A broader set of computations, including N-factor predictions for the relevant instabilities, will be presented in the final paper.

Boundary layer transition↗

NASA NTRS2021

A Computational Analysis of Boundary Layer Instability over the BOLT Configuration

The complex boundary layer flow over the BOLT flight configuration is known to exhibit multiple and potentially interacting instability mechanisms. This paper represents a continuation of our numerical investigation of the flow instabilities over the main test surface of the BOLT configuration by using state-of-the-art tools in multidimensional instability analysis. Specifically, the paper extends our previous computations by considering the separate effects of a nonzero angle of attack and a nonzero yaw on the modal instability characteristics of the boundary layer streaks adjacent to the minor-axis symmetry plane, specifically near both ends of the azimuthal region of a thick boundary layer in the middle, where this region rapidly changes to a thinner boundary layer on either side. At the t = 28.8767 s condition from the ascent part of the anticipated flight trajectory with a flight Mach number of M ∞ = 5.53 and unit Reynolds number of 4.25 x 10 6 /m, either type of departure from the design condition is shown to have a considerable impact on the structure of the basic state rollup within the region of interest and, hence, also on the amplification characteristics of instability waves within the resulting streaks. Yet, for a yaw angle of β = 4 degrees, the computations indicate only a slight reduction with respect to the peak N-factor of nearly 11 at the design condition of zero degrees yaw and zero degrees angle of attack. In contrast, an angle of attack equal to α = 4 degrees, the peak N-factor decreases to nearly 6 on the leeward side and increases above 16 on the windward side, making the onset of transition highly likely on the windward side. Computations also highlight the role of streak instabilities that originate as Mack mode disturbances and also demonstrate the potential pitfalls in using surface pressure sensors alone to gauge the magnitude of instability amplification.

Boundary layer transition↗

NASA NTRS2022

Linear Disturbance Amplification Over Blunted Flat Plates in High-Speed Flows

Modal and nonmodal instability characteristics of cylindrically blunted flat plates with varying leading edge radii are described for Mach 4 and Mach 6 freestream conditions. The selection of leading edge radii and freestream parameters is informed by experimental conditions. The investigation of this 2D problem provides a slow entropy layer swallowing which allows for an isolated development of perturbations seeded upstream within different wall-normal regions of the flow. At both Mach numbers, a decrease in modal instability amplification was seen as the leading edge radius was increased. Nonmodal analysis reveals amplifying perturbations in the boundary layer as well as in the entropy layer. The medium bluntness regime exhibits the strongest amplification of nonmodal disturbances that is nonmonotonic in character. Small amplitude boundary forcing at the plate surface or volumetric forcing at various wall-normal heights was used to account for receptivity effects. While wall forcing effectively induced modal instabilities, only an actuation above the boundary layer captured disturbances that amplify within the entropy layer. The optimal nonmodal theory’s entropy-layer disturbance evolution exhibited outstanding agreement with controlled forcing, including receptivity effects. The evolution of entropy-layer disturbances from the optimal nonmodal theory showed excellent agreement with the results of receptivity to controlled forcing. Therefore, the nonmodal optimal growth analysis may provide a useful as well as efficient technique to identify the complete disturbance spectrum in blunt hypersonic configurations, where both modal and nonmodal disturbances can amplify in the boundary-layer and entropy-layer regions.

Boundary layer transition↗

DOE OSTI2020

Enforcing constraints for time series prediction in supervised, unsupervised and reinforcement learning

We assume that we are given a time series of data from a dynamical system and our task is to learn the flow map of the dynamical system. We present a collection of results on how to enforce constraints coming from the dynamical system in order to accelerate the training of deep neural networks to represent the flow map of the system as well as increase their predictive ability. In particular, we provide ways to enforce constraints during training for all three major modes of learning, namely supervised, unsupervised and reinforcement learning. In general, the dynamic constraints need to include terms which are analogous to memory terms in model reduction formalisms. Such memory terms act as a restoring force which corrects the errors committed by the learned flow map during prediction. For supervised learning, the constraints are added to the objective function. For the case of unsupervised learning, in particular generative adversarial networks, the constraints are introduced by augmenting the input of the discriminator. Finally, for the case of reinforcement learning and in particular actor-critic methods, the constraints are added to the reward function. In addition, for the reinforcement learning case, we present a novel approach based on homotopy of the action-value function in order to stabilize and accelerate training. We use numerical results for the Lorenz system to illustrate the various constructions.

Stinis, Panagiotis↗

DOE OSTI2021

The Concept and Applications of a Dual Energy Storage Ring

A dual energy electron storage ring configuration is initially proposed as an electron cooler to cool the ion beam in a collider. It consists of two energy loops, the electron beam in the high energy loop undergoes the synchrotron radiation damping to obtain the desired beam property and the beam in the low energy loop is for cooling of the ion beam. The two different energy loops are connected by an energy recovery linac. A lattice design of such a dual energy storage ring has been completed and beam stability conditions are established. We performed numerical simulations to demonstrate the beam qualities and evaluated the cooling performance. In this paper, we present the study results and discuss possible applications of such a concept in many physics research and medical fields.

Dhital, B.↗

NASA NTRS1986

Reynolds stress closure in jet flows using wave models

A wave model was developed for the two dimensional shear layer. This configuration is being used as a test case for the closure schemes. Numerical methods are under development to solve the nonseparable Rayleigh equation. A model problem is being used to assist in the algorithm development. An analytic solution of the Rayleigh equation for a basic elliptic flow was obtained. This will be used to verify the stability codes developed for arbitrary geometries. Other numerical methods for solving the Rayleigh equation based on the boundary element technique are being examined.

Morris, P. J.↗

NASA NTRS1997

A Thermal and Electrical Analysis of Power Semiconductor Devices

The state-of-art power semiconductor devices require a thorough understanding of the thermal behavior for these devices. Traditional thermal analysis have (1) failed to account for the thermo-electrical interaction which is significant for power semiconductor devices operating at high temperature, and (2) failed to account for the thermal interactions among all the levels involved in, from the entire device to the gate micro-structure. Furthermore there is a lack of quantitative studies of the thermal breakdown phenomenon which is one of the major failure mechanisms for power electronics. This research work is directed towards addressing. Using a coupled thermal and electrical simulation, in which the drift-diffusion equations for the semiconductor and the energy equation for temperature are solved simultaneously, the thermo-electrical interactions at the micron scale of various junction structures are thoroughly investigated. The optimization of gate structure designs and doping designs is then addressed. An iterative numerical procedure which incorporates the thermal analysis at the device, chip and junction levels of the power device is proposed for the first time and utilized in a BJT power semiconductor device. In this procedure, interactions of different levels are fully considered. The thermal stability issue is studied both analytically and numerically in this research work in order to understand the mechanism for thermal breakdown.

Vafai, Kambiz↗

NASA NTRS1984

Stability of finite difference models containing two boundaries or interfaces

The stability of finite difference models of hyperbolic initial boundary value problems is connected with the propagation and reflection of parasitic waves. Wave propagation ideas are applied to models containing two boundaires or interfaces, where repeated reflection of trapped wave packets is a potential new source of instability. Various known instability phenomena are accounted for in a unified way. Results show: (1) dissipativity does not ensure stability when three or more formulas are concatenated at a boundary or internal interface; (2) algebraic GKS instabilities can be converted by a second boundary to exponential instabilities only when an infinite numerical reflection coefficient is present; and (3) GKS-stability and P-stability can be established in certain problems by showing that all numerical reflection coefficients have modulus less than 1.

Trefethen, L. N.↗

NASA NTRS1985

Stability of finite-difference models containing two boundaries or interfaces

The stability of finite difference models of hyperbolic initial boundary value problems is connected with the propagation and reflection of parasitic waves. Wave propagation ideas are applied to models containing two boundaries or interfaces, where repeated reflection of trapped wave packets is a potential new source of instability. Various known instability phenomena are accounted for in a unified way. Results show: (1) dissipativity does not ensure stability when three or more formulas are concatenated at a boundary or internal interface; (2) algebraic GKS instabilities can be converted by a second boundary to exponential instabilities only when an infinite numerical reflection coefficient is present; and (3) GKS-stability and P-stability can be established in certain problems by showing that all numerical reflection coefficients have modulus less than 1.

Trefethen, L. N.↗

NASA NTRS2022

Stability Analysis of Streaks Induced by Optimized Vortex Generators

Numerical computations are performed to investigate the potential for transition control in an axisymmetric boundary layer via fully realizable, streamwise stationary streaks induced by an azimuthally periodic array of surface mounted vortex generators (VGs). Previous work has shown that suitable streaks of this type can significantly reduce the growth of Mack’s second mode instabilities, but large streak amplitudes can make the flow susceptible to previously absent streak instabilities that can become the leading cause of transition. Here, we use the adjoint capabilities of the SU2 flow solver to optimize the VG shape to maximize the reduction in the growth of second-mode disturbances while also preventing the streak amplitudes from reaching large enough values to precipitate an earlier onset of transition via streak instabilities. The geometry and the freestream flow conditions are selected to match a relevant trajectory lo-cation from the HIFiRE-1 flight experiment. Results show that the optimized VGs can increase the mean streak amplitude by 117% with respect to a manually developed baseline design. The stability of this optimized basic state is analyzed via the plane-marching parabolized stability equations, predicting a fully laminar flow over the entire cone, or equivalently, yielding transition delay of 130% versus the 17% for the baseline VGs.

Boundary layer transition↗

NASA NTRS2022

Stability Analysis of Streaks Induced By Optimized Vortex Generators

Numerical computations are performed to investigate the potential for transition control in an axisymmetric boundary layer via fully realizable, streamwise stationary streaks induced by an azimuthally periodic array of surface mounted vortex generators (VGs). Previous work has shown that suitable streaks of this type can significantly reduce the growth of Mack’s second mode instabilities, but large streak amplitudes can make the flow susceptible to previously absent streak instabilities that can become the leading cause of transition. Here, we use the adjoint capabilities of the SU2 flow solver to optimize the VG shape to maximize the reduction in the growth of second-mode disturbances while also preventing the streak amplitudes from reaching large enough values to precipitate an earlier onset of transition via streak instabilities. The geometry and the freestream flow conditions are selected to match a relevant trajectory lo-cation from the HIFiRE-1 flight experiment. Results show that the optimized VGs can increase the mean streak amplitude by 117% with respect to a manually developed baseline design. The stability of this optimized basic state is analyzed via the plane-marching parabolized stability equations, predicting a fully laminar flow over the entire cone, or equivalently, yielding transition delay of 130% versus the 17% for the baseline VGs.

Boundary layer transition↗

NASA NTRS1979

Stability of sequences generated by nonlinear differential systems

A local stability analysis is presented for both the analytic and numerical solutions of the initial value problem for a system of ordinary differential equations. It is shown that, using a proper choice of Liapunov function, a connected region of stable initial values of both the analytic solution and the one-leg k-step numerical solution can be approximated. Attention is given to the example of the two-dimensional problem involving the stability of the longitudinal equations of motion of a gliding jet aircraft.

Brown, R. L.↗