Search NASA⌕ Search

SEARCH · Search NASA

Results for “Domain wall motion”

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.

116 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Incorporating Full Elastodynamic Effects and Dipping Fault Geometries in Community Code Verification Exercises for Simulations of Earthquake Sequences and Aseismic Slip (SEAS)

ABSTRACT Numerical modeling of earthquake dynamics and derived insight for seismic hazard relies on credible, reproducible model results. The sequences of earthquakes and aseismic slip (SEAS) initiative has set out to facilitate community code comparisons, and verify and advance the next generation of physics-based earthquake models that reproduce all phases of the seismic cycle. With the goal of advancing SEAS models to robustly incorporate physical and geometrical complexities, here we present code comparison results from two new benchmark problems: BP1-FD considers full elastodynamic effects, and BP3-QD considers dipping fault geometries. Seven and eight modeling groups participated in BP1-FD and BP3-QD, respectively, allowing us to explore these physical ingredients across multiple codes and better understand associated numerical considerations. With new comparison metrics, we find that numerical resolution and computational domain size are critical parameters to obtain matching results. Codes for BP1-FD implement different criteria for switching between quasi-static and dynamic solvers, which require tuning to obtain matching results. In BP3-QD, proper remote boundary conditions consistent with specified rigid body translation are required to obtain matching surface displacements. With these numerical and mathematical issues resolved, we obtain excellent quantitative agreements among codes in earthquake interevent times, event moments, and coseismic slip, with reasonable agreements made in peak slip rates and rupture arrival time. We find that including full inertial effects generates events with larger slip rates and rupture speeds compared to the quasi-dynamic counterpart. For BP3-QD, both dip angle and sense of motion (thrust versus normal faulting) alter ground motion on the hanging and foot walls, and influence event patterns, with some sequences exhibiting similar-size characteristic earthquakes, and others exhibiting different-size events. These findings underscore the importance of considering full elastodynamics and nonvertical dip angles in SEAS models, as both influence short- and long-term earthquake behavior and are relevant to seismic hazard.

Geochemistry & Geophysics↗
DOE OSTI2024

Stochastic modal velocity field in rough-wall turbulence

Stochastically generated instantaneous velocity profiles are used to reproduce the outer region of rough-wall turbulent boundary layers in a range of Reynolds numbers extending from the wind tunnel to field conditions. Each profile consists in a sequence of steps, defined by the modal velocities and representing uniform momentum zones (UMZs), separated by velocity jumps representing the internal shear layers. Height-dependent UMZ is described by a minimal set of attributes: thickness, mid-height elevation, and streamwise (modal) and vertical velocities. These are informed by experimental observations and reproducing the statistical behaviour of rough-wall turbulence and attached eddy scaling, consistent with the corresponding experimental datasets. Sets of independently generated profiles are reorganized in the streamwise direction to form a spatially consistent modal velocity field, starting from any randomly selected profile. The operation allows one to stretch or compress the velocity field in space, increases the size of the domain and adjusts the size of the largest emerging structures to the Reynolds number of the simulated flow. By imposing the autocorrelation function of the modal velocity field to be anchored on the experimental measurements, we obtain a physically based spatial resolution, which is employed in the computation of the velocity spectrum, and second-order structure functions. The results reproduce the Kolmogorov inertial range extending from the UMZ and their attached-eddy vertical organization to the very-large-scale motions (VLSMs) introduced with the reordering process. The dynamic role of VLSM is confirmed in the –u'w' co-spectra and in their vertical derivative, representing a scale-dependent pressure gradient contribution.

42 ENGINEERING↗
DOE OSTI2020

Evaluating the Influence of Hemorheological Parameters on Circulating Tumor Cell Trajectory and Simulation Time

Extravasation of circulating tumor cells (CTCs) occurs primarily in the microvasculature, where flow and cell interactions significantly affect the blood rheology. Capturing cell trajectory at this scale requires the coupling of several interaction models, leading to increased computational cost that scales as more cells are added or the domain size is increased. In this work, we focus on micro-scale vessels and study the influence of certain hemorheological factors, including the presence of red blood cell aggregation, hematocrit level, microvessel size, and shear rate, on the trajectory of a circulating tumor cell. We determine which of the aforementioned factors significantly affect CTC motion and identify those which can potentially be disregarded, thus reducing simulation time. We measure the effect of these elements by studying the radial CTC movement and runtime at various combinations of these hemorheological parameters. To accurately capture blood flow dynamics and single cell movement, we perform high-fidelity hemodynamic simulations at a sub-micron resolution using our in-house fluid dynamics solver, HARVEY. We find that increasing hematocrit increases the likelihood of tumor cell margination, which is exacerbated by the presence of red blood cell aggregation. As microvessel diameter increases, there is no major CTC movement towards the wall; however, including aggregation causes the CTC to marginate quicker as the vessel size increases. Finally, as the shear rate is increased, the presence of aggregation has a diminished effect on tumor cell margination.

Roychowdhury, Sayan↗
NASA NTRS1992

Direct numerical simulation of instabilities in parallel flow with spherical roughness elements

Results from a direct numerical simulation of laminar flow over a flat surface with spherical roughness elements using a spectral-element method are given. The numerical simulation approximates roughness as a cellular pattern of identical spheres protruding from a smooth wall. Periodic boundary conditions on the domain's horizontal faces simulate an infinite array of roughness elements extending in the streamwise and spanwise directions, which implies the parallel-flow assumption, and results in a closed domain. A body force, designed to yield the horizontal Blasius velocity in the absence of roughness, sustains the flow. Instabilities above a critical Reynolds number reveal negligible oscillations in the recirculation regions behind each sphere and in the free stream, high-amplitude oscillations in the layer directly above the spheres, and a mean profile with an inflection point near the sphere's crest. The inflection point yields an unstable layer above the roughness (where U''(y) is less than 0) and a stable region within the roughness (where U''(y) is greater than 0). Evidently, the instability begins when the low-momentum or wake region behind an element, being the region most affected by disturbances (purely numerical in this case), goes unstable and moves. In compressible flow with periodic boundaries, this motion sends disturbances to all regions of the domain. In the unstable layer just above the inflection point, the disturbances grow while being carried downstream with a propagation speed equal to the local mean velocity; they do not grow amid the low energy region near the roughness patch. The most amplified disturbance eventually arrives at the next roughness element downstream, perturbing its wake and inducing a global response at a frequency governed by the streamwise spacing between spheres and the mean velocity of the most amplified layer.

Deanna, R. G.↗
DOE OSTI2022

Understanding the Impacts of Support–Polymer Interactions on the Dynamics of Poly(ethyleneimine) Confined in Mesoporous SBA-15

Supported amines are a promising class of CO 2 sorbents offering large uptake capacities and fast uptake rates. Additionally, among supported amines, poly(ethyleneimine) (PEI) physically impregnated in the mesopores of SBA-15 silica is widely used. Within these composite materials, the chain dynamics and morphologies of PEI strongly influence the CO 2 capture performance, yet little is known about chain and macromolecule mobility in confined pores. Here, we probe the impact of the support–PEI interactions on the dynamics and structures of PEI at the support interface and the corresponding impact on CO 2 uptake performance, which yields critical structure–property relationships. The pore walls of the support are grafted with organosilanes with different chemical end groups to differentiate interaction modes (spanning from strong attraction to repulsion) between the pore surface and PEI. Combinations of techniques, such as quasi-elastic neutron scattering (QENS), 1 H T 1 –T 2 relaxation correlation solid-state NMR, and molecular dynamics (MD) simulations, are used to comprehensively assess the physical properties of confined PEI. We hypothesized that PEI would have faster dynamics when subjected to less attractive or repulsive interactions. However, we discover that complex interfacial interactions resulted in complex structure–property relationships. Indeed, both the chain conformation of the surface-grafted chains and of the PEI around the surface influenced the chain mobility and CO 2 uptake performance. By coupling knowledge of the dynamics and distributions of PEI with CO 2 sorption performance and other characteristics, we determine that the macroscopic structures of the hybrid materials dictate the first rapid CO 2 uptake, and the rate of CO 2 sorption during the subsequent gradual uptake stage is determined by PEI chain motions that promote diffusive jumps of CO 2 through PEI-packed domains.

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

ESR Dipole Power Supply Current Ripple and Noise Specifications

This note presents key findings for the ESR main magnet dipole power supplies (PS), where we find the current ripple specification to be close to or beyond the state-of-the-art. These specifications originate from beam-beam considerations, with the requirement to limit the ripple-induced hadron emittance growth to below 10%/hour. Beam dynamics that drive this PS ripple specification arise from the beam motions at the Interaction Point (IP). The frequency of the motions can be separated into "low", compared to the betatron frequency, and "high", i.e. around the betatron frequency and harmonics. In terms of the driving frequency, "low" implies ƒ<<ƒ 0 v x,y and "fast" means ƒ≈{ƒ 0 v x,y , ƒ 0 (1- v x,y , etc.}, where ƒ 0 =1/ T 0 =78.2 kHz is the revolution frequency, and v x,y are the fractional parts of the betatron tunes. Frequencies higher than ƒ 0 /2 will be folded back due to the particles sampling the field once per turn. To provide flexibility for future lattice adjustments and working point variations, we do not consider the tunes as fixed. Instead, we assume a certain margin and allow them to potentially fall within the range of 0.1< v x,y <0.5. In other words, the high-frequency region spans approximately from 8kHz to 40 kHz. Consequently,, we definte the dipole PS ripple in two distinct frequency ranges: the low-frequency range of [1-8000] Hz and the high-frequency range of [8-40] kHz. For the physics effects we analyzed in this note, there is no distinction between the ripple (which can be approximately reproduced in the frequency domain) and the random noise if both have some power within the frequency bandwidth of interest. Therefore, while we will use the term "ripple" for short, it should always be understood that we are referring to "ripple plus noise." Except for the lower end of the low-frequency range, the impact of the rippling PS current on the beam will be considerably reduced due to the eddy currents induced in the walls of the vacuum chamber. We will account for this effect in the PS ripple specifications to follow. The remaining sections of this note are structure as follows: Section II outlines the beam-beam physics requirements for the positional stability of the beam at the IP. Section III describes the anticipated shielding effect of the eddy currents induced in the vacuum chamber. Section IV derives the ripple requirement for the low-frequency range by propagating the closed orbit ripple resulting from the rippling dipoles to the IP (relevant lattice simulation results are summarized in the Appendix). In Section V, we present the analytical criterion for the ripple in the high-frequency range by considering resonant oscillations of the electron beam around a stable closed orbit near the betatron frequency. Finally, Section VI provides a summary of our findings and discusses related work.

43 PARTICLE ACCELERATORS↗
NASA NTRS2013

Constraints on Decreases in Eta Carinae's Mass-loss from 3D Hydrodynamic Simulations of Its Binary Colliding Winds

Recent work suggests that the mass-loss rate of the primary star Eta-A in the massive colliding wind binary Eta Carinae dropped by a factor of 2-3 between 1999 and 2010. We present result from large- (+/- 1545 au) and small- (+/- 155 au) domain, 3D smoothed particle hydrodynamics (SPH) simulations of Eta Car's colliding winds for three Eta-A mass-loss rates ( (dot-M(sub Eta-A) = 2.4, 4.8 and 8.5 × 10(exp −4) M(solar)/ yr), investigating the effects on the dynamics of the binary wind-wind collision (WWC). These simulations include orbital motion, optically thin radiative cooling and radiative forces. We find that dot-M Eta-A greatly affects the time-dependent hydrodynamics at all spatial scales investigated. The simulations also show that the post-shock wind of the companion star Eta-B switches from the adiabatic to the radiative-cooling regime during periastron passage (Phi approx.= 0.985-1.02). This switchover starts later and ends earlier the lower the value of dot-M Eta-A and is caused by the encroachment of the wind of Eta-A into the acceleration zone of Eta-B's wind, plus radiative inhibition of Eta-B's wind by Eta-A. The SPH simulations together with 1D radiative transfer models of Eta-A's spectra reveal that a factor of 2 or more drop in dot-M EtaA should lead to substantial changes in numerous multiwavelength observables. Recent observations are not fully consistent with the model predictions, indicating that any drop in dot- M Eta-A was likely by a factor of approx. < 2 and occurred after 2004. We speculate that most of the recent observed changes in Eta Car are due to a small increase in the WWC opening angle that produces significant effects because our line of sight to the system lies close to the dense walls of the WWC zone. A modest decrease in dot-M Eta-A may be responsible, but changes in the wind/stellar parameter of Eta-B, while less likely, cannot yet be fully ruled out. We suggest observations during Eta-Car's next periastron in 2014 to further test for decreases in dot-M Eta-A. If dot-M Eta-A is declining and continues to do so, the 2014 X-ray minimum should be even shorter than that of 2009.

Eta Carinae↗
NASA NTRS2016

Applied Time Domain Stability Margin Assessment for Nonlinear Time-Varying Systems

The baseline stability margins for NASA's Space Launch System (SLS) launch vehicle were generated via the classical approach of linearizing the system equations of motion and determining the gain and phase margins from the resulting frequency domain model. To improve the fidelity of the classical methods, the linear frequency domain approach can be extended by replacing static, memoryless nonlinearities with describing functions. This technique, however, does not address the time varying nature of the dynamics of a launch vehicle in flight. An alternative technique for the evaluation of the stability of the nonlinear launch vehicle dynamics along its trajectory is to incrementally adjust the gain and/or time delay in the time domain simulation until the system exhibits unstable behavior. This technique has the added benefit of providing a direct comparison between the time domain and frequency domain tools in support of simulation validation. This technique was implemented by using the Stability Aerospace Vehicle Analysis Tool (SAVANT) computer simulation to evaluate the stability of the SLS system with the Adaptive Augmenting Control (AAC) active and inactive along its ascent trajectory. The gains for which the vehicle maintains apparent time-domain stability defines the gain margins, and the time delay similarly defines the phase margin. This method of extracting the control stability margins from the time-domain simulation is relatively straightforward and the resultant margins can be compared to the linearized system results. The sections herein describe the techniques employed to extract the time-domain margins, compare the results between these nonlinear and the linear methods, and provide explanations for observed discrepancies. The SLS ascent trajectory was simulated with SAVANT and the classical linear stability margins were evaluated at one second intervals. The linear analysis was performed with the AAC algorithm disabled to attain baseline stability margins. At each time point, the system was linearized about the current operating point using Simulink's built-in solver. Each linearized system in time was evaluated for its rigid-body gain margin (high frequency gain margin), rigid-body phase margin, and aero gain margin (low frequency gain margin) for each control axis. Using the stability margins derived from the baseline linearization approach, the time domain derived stability margins were determined by executing time domain simulations in which axis-specific incremental gain and phase adjustments were made to the nominal system about the expected neutral stability point at specific flight times. The baseline stability margin time histories were used to shift the system gain to various values around the zero margin point such that a precise amount of expected gain margin was maintained throughout flight. When assessing the gain margins, the gain was applied starting at the time point under consideration, thereafter following the variation in the margin found in the linear analysis. When assessing the rigid-body phase margin, a constant time delay was applied to the system starting at the time point under consideration. If the baseline stability margins were correctly determined via the linear analysis, the time domain simulation results should contain unstable behavior at certain gain and phase values. Examples will be shown from repeated simulations with variable added gain and phase lag. Faithfulness of margins calculated from the linear analysis to the nonlinear system will be demonstrated.

Kiefer, J. M.↗