Search NASA⌕ Search

SEARCH · Search NASA

Results for “numerical stabilization”

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 505 records · Page 28

Subleading power corrections for N -jettiness subtractions

Here, the N -jettiness observable T N provides a way of describing the leading singular behavior of the N -jet cross section in the τ = T N = Q → 0 limit, where Q is a hard interaction scale. We consider subleading power corrections in the τ << 1 expansion, and employ soft-collinear effective theory to obtain analytic results for the dominant α s τlnτ and α 2 s τln 3 τ subleading terms for thrust in e + e – collisions and 0-jettiness for q¯q-initiated Drell-Yan–like processes at hadron colliders. These results can be used to significantly improve the numerical accuracy and stability of the N -jettiness subtraction technique for performing fixed order calculations at next-to-leading order and next-to-next-to-leading order. They reduce the size of missing power corrections in the subtractions by an order of magnitude. We also point out that the precise definition of N -jettiness has an important impact on the size of the power corrections and thus the numerical accuracy of the subtractions. The sometimes employed definition of N -jettiness in the hadronic center-ofmass frame suffers from power corrections that grow exponentially with rapidity, causing the power expansion to deteriorate away from central rapidity. This degradation does not occur for the original N -jettiness definition, which explicitly accounts for the boost of the Born process relative to the frame of the hadronic collision, and has a well-behaved power expansion throughout the entire phase space. Integrated over rapidity, using this N -jettiness definition in the subtractions yields another order of magnitude improvement compared to employing the hadronic-frame definition.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Magnetization stability analysis of the Stoner-Wohlfarth model under a spin-polarized current with a tilted polarization

The stationary-state solutions of magnetization dynamics under a spin-polarized current that was polarized in an arbitrary direction were investigated by solving the Landau-Lifshitz-Gilbert-Slonczewski equation for a single-domain magnet. Taking into consideration the uniaxial magnetic anisotropy, the equilibrium directions of the magnetization vectors were analytically obtained by solving an algebraic cubic equation. It was found that one to three pairs of magnetization equilibrium states existed, depending on the current intensity and the direction of the spin polarization. By numerically analyzing the stabilities of these equilibrium states, the threshold switching current for the reversing the magnetic vector was obtained under different current polarization configurations, which may be useful for use in future spintronics devices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Algorithm Stability and the Long-Term Geospace Data Record from TIMED/SABER

The ability of satellite instruments to accurately observe long-term changes in atmospheric temperature depends on many factors including the absolute accuracy of the measurement, the stability of the calibration of the instrument, the stability of the satellite orbit, and the stability of the numerical algorithm that produces the temperature data. We present an example of algorithm instability recently discovered in the temperature dataset from the SABER instrument on the NASA TIMED satellite. The instability resulted in derived temperatures that were substantially colder than anticipated from mid-December 2019 to mid-2022. This algorithm-induced change in temperature over one to two years corresponded to the expected change over several decades from increasing anthropogenic CO2. This paper highlights the importance of algorithm stability in developing Geospace Data Records (GDRs) for Earth’s mesosphere and lower thermosphere. A corrected version (Version 2.08) of the temperatures from SABER is described.

M G Mlynczak↗

Operation of an ungated diamond field-emission array cathode in a L-band radiofrequency electron source

We report on the operation of a field-emitter-array cathode in a conventional L-band radio-frequency electron source. The cathode consisted of an array of ∼10{sup 6} diamond tips on pyramids. Maximum current on the order of 15 mA was reached and the cathode did not show appreciable signs of fatigue after weeks of operation. The measured Fowler-Nordheim characteristics, transverse beam density, and current stability are discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Ideal magnetohydrodynamic stability of configurations without nested flux surfaces

Existing numerical tools for calculating the MHD stability of magnetically confined plasmas generally assume the existence of nested flux surfaces. These tools are therefore not immediately applicable to configurations with magnetic islands or regions with an ergodic magnetic field. However, in practice, these islands or ergodic regions are often small, and their effect on MHD stability can then be evaluated using a perturbation theory developed in the present paper. This procedure allows the effect of the broken magnetic topology on the stability of each eigenmode to be calculated without requiring any knowledge about the perturbed eigenfunctions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Effect of negative triangularity on tearing mode stability in tokamak plasmas

Abstract The influence of negative triangularity (NT) of the plasma shape on the n = 1 (n is the toroidal mode number) tearing mode (TM) stability has been numerically investigated, with results compared to that of the positive triangularity (PT) counterpart. By matching the safety factor profile for a series of toroidal equilibria, several important plasma parameters, including the triangularity, the plasma equilibrium pressure, the plasma resistivity as well as the toroidal rotation, have been varied. The TM localized near the plasma edge is found to be more unstable in the NT plasmas as compared to the PT counterpart. The fundamental reason for this difference is the lack of favorable average curvature stabilization in NT configurations. Direct comparison of the Mercier index corroborates this conclusion. For the core-localized mode, where the difference in the local triangularity between NT and PT becomes small and the curvature stabilization is significantly reduced, larger Shafranov shift in the plasma core associated with the NT configuration results in more stable TM. The plasma toroidal flow generally stabilizes the TM in plasmas with both NTs and PTs. The flow stabilization is however weaker in the case of negative triangularity with finite plasma pressure.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Input-output-controlled nonlinear equation solvers

To upgrade the efficiency and stability of the successive substitution (SS) and Newton-Raphson (NR) schemes, the concept of input-output-controlled solvers (IOCS) is introduced. By employing the formal properties of the constrained version of the SS and NR schemes, the IOCS algorithm can handle indefiniteness of the system Jacobian, can maintain iterate monotonicity, and provide for separate control of load incrementation and iterate excursions, as well as having other features. To illustrate the algorithmic properties, the results for several benchmark examples are presented. These define the associated numerical efficiency and stability of the IOCS.

Padovan, Joseph↗

Stable calcium-free myocilin olfactomedin domain variants reveal challenges in differentiating between benign and glaucoma-causing mutations

Nonsynonymous gene mutations can be beneficial, neutral, or detrimental to the stability, structure, and biological function of the encoded protein, but the effects of these mutations are often not readily predictable. For example, the β-propeller olfactomedin domain of myocilin (mOLF) exhibits a complex interrelationship among structure(s), stability, and aggregation. Numerous mutations within mOLF are linked to glaucoma; the resulting variants are less stable, aggregation-prone, and sequestered intracellularly, causing cytotoxicity. In this study, we report the first stable mOLF variants carrying substitutions in the calcium-binding site that exhibit solution characteristics indistinguishable from those of glaucoma variants. Crystal structures of these stable variants at 1.8–2.0-Å resolution revealed features that we could not predict by molecular dynamics simulations, including loss of loop structure, helix unwinding, and a blade shift. Double mutants that combined a stabilizing substitution and a selected glaucoma-causing single-point mutant rescued in vitro folding and stability defects. In the context of full-length myocilin, secretion of stable single variants was indistinguishable from that of the WT protein, and the double mutants were secreted to varying extents. In summary, our finding that mOLF can tolerate particular substitutions that render the protein stable despite a conformational switch emphasizes the complexities in differentiating between benign and glaucoma-causing variants and provides new insight into the possible biological function of myocilin.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

An analysis of the stability of the compressible Ekman boundary layer

The linear stability problem for the compressible Ekman boundary layer common to rotating fluids is formulated and the stability properties determined numerically. Three classes of unstable waves are identified (called class A, B, and C), their properties are described. The class C waves have only recently been reported in the literature and are present only in compressible Ekman boundary layers. Most of the calculations presented here are for uranium hexafluoride gas; however, critical Reynolds numbers are also computed for air and ammonia gas. Compressibility is generally found to decrease the critical Reynolds number for each class of wave. A comparison of results for the three different gases shows the stability to be largely unaffected by changes in the gas properties. Maximum growth rate calculations for each wave show the class A and B waves to be the dominant instabilities.

Spall, J. R.↗

Comparison of momentum transport models for numerical relativity

The main problems of nonvacuum numerical relativity, compact binary mergers and stellar collapse, involve hydromagnetic instabilities and turbulent flows, so that kinetic energy at small scales leads to mean effects at large scale that drive the secular evolution. Notable among these effects is momentum transport. We investigate two models of this transport effect, a relativistic Navier-Stokes system and a turbulent mean stress model, that are similar to all of the prescriptions that have been attempted to date for treating subgrid effects on binary neutron star mergers and their aftermath. Our investigation involves both stability analysis and numerical experimentation on star and disk systems. We also begin the investigation of the effects of particle and heat transport on postmerger simulations. Here, we find that correct handling of turbulent heating is crucial for avoiding unphysical instabilities. Given such appropriate handling, the evolution of a differentially rotating star and the accretion rate of a disk are reassuringly insensitive to the choice of prescription. However, disk outflows can be sensitive to the choice of method, even for the same effective viscous strength. We also consider the effects of eddy diffusion in the evolution of an accretion disk and show that it can interestingly affect the composition of outflows.

79 ASTRONOMY AND ASTROPHYSICS↗

A Discrete Constraint for Entropy Conservation and Sound Waves in Cloud-Resolving Modeling

Ideal cloud-resolving models contain little-accumulative errors. When their domain is so large that synoptic large-scale circulations are accommodated, they can be used for the simulation of the interaction between convective clouds and the large-scale circulations. This paper sets up a framework for the models, using moist entropy as a prognostic variable and employing conservative numerical schemes. The models possess no accumulative errors of thermodynamic variables when they comply with a discrete constraint on entropy conservation and sound waves. Alternatively speaking, the discrete constraint is related to the correct representation of the large-scale convergence and advection of moist entropy. Since air density is involved in entropy conservation and sound waves, the challenge is how to compute sound waves efficiently under the constraint. To address the challenge, a compensation method is introduced on the basis of a reference isothermal atmosphere whose governing equations are solved analytically. Stability analysis and numerical experiments show that the method allows the models to integrate efficiently with a large time step.

Zeng, Xi-Ping↗

Numerical Simulation of Dynamic Stall Using Near-Body Adaptive Mesh Refinement

Time-dependent Reynolds-averaged Navier-Stokes simulations have been carried out for a NACA 0012 airfoil and a flexible UH-60A rotor undergoing dynamic stall. The OVERFLOW computational fluid dynamics (CFD) code is used to explore the use of near-body adaptive mesh refinement (NB-AMR) for the first time on a flexible helicopter rotor in forward flight. Emphasis is placed on understanding the temporal and spatial convergence of the solutions, their numerical efficiency and stability, and establishing a grid-converged solution. Some of the remaining differences between CFD and flight-test measurements are discussed. Time-dependent flow visualization is used to provide an improved understanding of the physical mechanisms involved with two-dimensional and three-dimensional dynamic stall, and the NB-AMR process.

Near-Body↗

Influence of weakening minerals on ensemble strength and slip stability of faults: Influence of Weakening Minerals on Fault

We explore the impact of phyllosilicate (weak but velocity strengthening) in a majority tectosilicate (strong but velocity weakening) matrix on bulk shear strength and slip stability of faults. Numerical simple shear experiments using a distinct element model (DEM) are conducted on both uniform mixtures of quartz and talc analogs and on textured mixtures consisting of a talc layer embedded in a quartz matrix. The mechanical response of particles is represented by a linear elastic contact model with a slip-weakening constitutive relation representing the essence of rate-state friction. The weight percentage of the talc in the uniform mixtures and the relative thickness of the talc layer in the textured mixtures are varied to investigate the transitional behavior of shear strength and slip stability. Specifically, for uniform mixtures, ~50% reduction on bulk shear strength is observed with 25% talc present, and a dominant influence of talc occurs at 50%; for textured mixtures, a noticeable weakening effect is shown at a relative layer thickness of 1 particle, ~50% shear strength reduction is observed with 3-particles, and a dominant influence occurs at 5 particles. In terms of slip stability, a transition from velocity weakening to velocity strengthening is observed with 10% to 25% talc present in the uniform mixtures or with 3 particles to 5 particles in the textured mixtures. In addition, further analysis suggests that quartz has a high tendency toward dilation, potentially promoting permeability; while talc dilates with increased slip rate but compacts rapidly when slip rate is reduced, potentially destroying permeability. The simulation results match well with previous laboratory observations.

58 GEOSCIENCES↗

Wiggle instability of galactic spiral shocks: Effects of magnetic fields

It has been suggested that the wiggle instability (WI) of spiral shocks in a galactic disk is responsible for the formation of gaseous feathers observed in grand-design spiral galaxies. We perform both a linear stability analysis and numerical simulations to investigate the effect of magnetic fields on the WI. The disk is assumed to be infinitesimally thin, isothermal, and non-self-gravitating. We control the strengths of magnetic fields and spiral-arm forcing using the dimensionless parameters β and F, respectively. By solving the perturbation equations as a boundary-eigenvalue problem, we obtain dispersion relations of the WI for various values of β=1−∞ and F=5% and 10%. We find that the WI arising from the accumulation of potential vorticity at disturbed shocks is suppressed, albeit not completely, by magnetic fields. The stabilizing effect of magnetic fields is not from the perturbed fields but from the unperturbed fields that reduce the density compression factor in the background shocks. When F=5% and β≲10 or F=10% and β ∼ 5–10, the most unstable mode has a wavelength of ∼0.1–0.2 times the arm-to-arm separation, which appears consistent with a mean spacing of observed feathers.

79 ASTRONOMY AND ASTROPHYSICS↗

Nonlinear evolution of the electromagnetic electron-cyclotron instability in bi-Kappa distributed plasma

This paper presents a numerical study of the linear and nonlinear evolution of the electromagnetic electron-cyclotron (EMEC) instability in a bi-Kappa distributed plasma. Distributions with high energy tails described by the Kappa power-laws are often observed in collision-less plasmas (e.g., solar wind and accelerators), where wave-particle interactions control the plasma thermodynamics and keep the particle distributions out of Maxwellian equilibrium. Under certain conditions, the anisotropic bi-Kappa distribution gives rise to plasma instabilities creating low-frequency EMEC waves in the whistler branch. The instability saturates nonlinearly by reducing the temperature anisotropy until marginal stability is reached. Numerical simulations of the Vlasov-Maxwell system of equations show excellent agreement with the growth-rate and real frequency of the unstable modes predicted by linear theory. The wave-amplitude of the EMEC waves at nonlinear saturation is consistent with magnetic trapping of the electrons.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Electroconvection near an ion-selective surface with Butler–Volmer kinetics

We study the effects of interfacial kinetics on the electro-hydrodynamics of ion transport near an ion-selective surface using a combination of linear stability analysis and numerical simulation. The finite kinetics of the electrolyte–electrode interface affects the ion transfer and electroconvection in many ways. On a surface of fixed topography, such as a metal surface of slow and stable ion deposition or covered by a polymer membrane, the finite kinetics reduces the current in one-dimensional ion diffusion/migration, increases the critical voltage for the onset of the electroconvective instability, changes the dynamics of the electroconvection and the overlimiting current, and enhances the lateral ion diffusion within the interfacial layer. The first three effects are indirectly caused by the reaction kinetics and can be characterized by an effective voltage difference across the liquid electrolyte. In comparison, the last effect is controlled by a direct interplay between kinetics and nonlinear electroconvection. Scaling laws for ion transport and features of electroconvection are proposed. We also analyse the linear stability of a surface which evolves under ion deposition and find that the finite kinetics decreases the growth rate of both electroconvective and morphological instabilities and therefore modifies the wavenumber of the most unstable mode.

Mechanics↗

Fundamental Investigations of Airframe Noise

An extensive numerical and experimental study of airframe noise mechanisms associated with a subsonic high-lift system has been performed at NASA Langley Research Center (LaRC). Investigations involving both steady and unsteady computations and experiments on a small-scale, part-span flap model are presented. Both surface (steady and unsteady pressure measurements, hot films, oil flows, pressure sensitive paint) and off surface (5 hole-probe, particle-imaged velocimetry, laser velocimetry, laser light sheet measurements) were taken in the LaRC Quiet Flow Facility (QFF) and several hard-wall tunnels up to flight Reynolds number. Successful microphone array measurements were also taken providing both acoustic source maps on the model, and quantitative spectra. Critical directivity measurements were obtained in the QFF. NASA Langley unstructured and structured Reynolds- Averaged Navier-Stokes codes modeled the flap geometries excellent comparisons with surface and offsurface experimental data were obtained. Subsequently, these meanflow calculations were utilized in both linear stability and direct numerical simulations of the flap-edge flow field to calculate unsteady surface pressures and farfield acoustic spectra. Accurate calculations were critical in obtaining not only noise source characteristics, but shear layer correction data as well. Techniques utilized in these investigations as well as brief overviews of results will be given.

Macaraeg, M. G.↗

Shear Instabilities as a Probe of Jupiter's Atmosphere

Linear wave patterns in Jupiter clouds with wavelengths strongly clustered around 300 km are commonly observed in the planet's equatorial atmosphere. We propose that the preferred wavelength is related to the thickness of an unstable shear layer within the clouds. We numerically analyze the linear stability of wavelike disturbances that have nonzero horizontal phase speeds in Jupiter's atmosphere and find that. if the static stability in the strongly clustered around 300 km are commonly observed in the planet's equatorial atmosphere. We propose that the preferred wavelength is related to the thickness of an unstable shear layer within the clouds. We numerically analyze the linear stability of wavelike disturbances that have nonzero horizontal phase speeds in Jupiter's atmosphere and find that. if the static stability in the shear layer is very low (but still nonnegative), a deep vertical shear layer like the one measured by the Galileo probe can generate the instabilities. The fastest growing waves grow exponentially within an hour, and their wavelengths match the observations. Close to zero values of static stability that permit the growth of instabilities are within the range of values measured by the Galileo probe in a hot spot. Our model probes Jupiter's equatorial atmosphere below the cloud deck and suggests that thick regions of wind shear and low static stability exist outside hot spots.

Bosak, Tanja↗