Fully 3D Attenuation and Scatter Correction in SPECT Using Monte Carlo Simulated System Matrices that Model for Scatter from Outside the Field-of-View
Explore the source record for details and available documents.
97 MATHEMATICS AND COMPUTING↗
SEARCH · Search NASA
Results for “3D Simulation”
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.
Velocity gradients as a tracer for magnetic fields
Strong Alfvénic turbulence develops eddy-like motions perpendicular to the local direction of magnetic fields. This local alignment induces velocity gradients perpendicular to the local direction of the magnetic field. We use this fact to propose a new technique of studying the direction of magnetic fields from observations, which we call the velocity gradient technique. We test our idea by employing the synthetic observations obtained via 3D magnetohydrodynamical (MHD) numerical simulations for different sonic and Alfvén Mach numbers. We calculate the velocity gradient, Ω, using the velocity centroids. We find that Ω traces the projected magnetic field best for the synthetic maps obtained with sub-Alfvénic simulations and provides good point-wise correspondence between the magnetic field direction and the direction of Ω. The reported alignment is much better than the alignment between the density gradients and the magnetic field, and we demonstrate that it can be used to find the magnetic field strength with an analog of the Chandrasekhar–Fermi method. This new technique does not require dust polarimetry, and our study opens up a new way of studying magnetic fields using spectroscopic data.
79 ASTRONOMY AND ASTROPHYSICS↗
Origin and Evolution of Magnetic Field in PMS Stars: Influence of Rotation and Structural Changes
During stellar evolution, especially in the pre-main-sequence phase, stellar structure and rotation evolve significantly, causing major changes in the dynamics and global flows of the star. We wish to assess the consequences of these changes on stellar dynamo, internal magnetic field topology, and activity level. To do so, we have performed a series of 3D HD and MHD simulations with the ASH code. We choose five different models characterized by the radius of their radiative zone following an evolutionary track computed by a 1D stellar evolution code. These models characterized stellar evolution from 1 to 50 Myr. By introducing a seed magnetic field in the fully convective model and spreading its evolved state through all four remaining cases, we observe systematic variations in the dynamical properties and magnetic field amplitude and topology of the models. The five MHD simulations develop a strong dynamo field that can reach an equipartition state between the kinetic and magnetic energies and even superequipartition levels in the faster-rotating cases. We find that the magnetic field amplitude increases as it evolves toward the zero-age main sequence. Moreover, the magnetic field topology becomes more complex, with a decreasing axisymmetric component and a nonaxisymmetric one becoming predominant. The dipolar components decrease as the rotation rate and the size of the radiative core increase. The magnetic fields possess a mixed poloidal-toroidal topology with no obvious dominant component. Moreover, the relaxation of the vestige dynamo magnetic field within the radiative core is found to satisfy MHD stability criteria. Hence, it does not experience a global reconfiguration but slowly relaxes by retaining its mixed stable poloidal-toroidal topology.
79 ASTRONOMY AND ASTROPHYSICS↗
Comparison of Solar Fine Structure Observed Simultaneously in Ly α and Mg ii h
The Chromospheric Lyman Alpha Spectropolarimeter (CLASP) observed the Sun in H i Ly α during a suborbital rocket flight on 2015 September 3. The Interface Region Imaging Telescope ( IRIS ) coordinated with the CLASP observations and recorded nearly simultaneous and co-spatial observations in the Mg ii h and k lines. The Mg ii h and Ly α lines are important transitions, energetically and diagnostically, in the chromosphere. The canonical solar atmosphere model predicts that these lines form in close proximity to each other and so we expect that the line profiles will exhibit similar variability. In this analysis, we present these coordinated observations and discuss how the two profiles compare over a region of quiet Sun at viewing angles that approach the limb. In addition to the observations, we synthesize both line profiles using a 3D radiation-MHD simulation. In the observations, we find that the peak width and the peak intensities are well correlated between the lines. For the simulation, we do not find the same relationship. We have attempted to mitigate the instrumental differences between IRIS and CLASP and to reproduce the instrumental factors in the synthetic profiles. The model indicates that formation heights of the lines differ in a somewhat regular fashion related to magnetic geometry. This variation explains to some degree the lack of correlation, observed and synthesized, between Mg ii and Ly α . Our analysis will aid in the definition of future observatories that aim to link dynamics in the chromosphere and transition region.
79 ASTRONOMY AND ASTROPHYSICS↗
Improved Measurements of Galaxy Star Formation Stochasticity from the Intrinsic Scatter of Burst Indicators
Abstract Measurements of short-timescale star formation variations (i.e., “burstiness” or star formation stochasticity) are integral to our understanding of star formation feedback mechanisms and the assembly of stellar populations in galaxies. We expand upon the work of Broussard et al. by introducing a new analysis of galaxy star formation burstiness that accounts for variations in the Q sg = E B − V stars / E B − V gas distribution, a major confounding factor. We use Balmer decrements from the MOSFIRE Deep Evolution Field (MOSDEF) survey to measure Q sg , which we use to construct mock catalogs from the Santa Cruz Semi-Analytic Models and Mufasa cosmological hydrodynamical simulation based on 3D-HST, Fiber Multi-Object Spectrograph (FMOS)-COSMOS, and MOSDEF galaxies with H α detections. The results of the mock catalogs are compared against observations using the burst indicator η = log 10 ( SFR H α / SFR NUV ) , with the standard deviation of the η distribution indicating burstiness. We find decent agreement between mock and observed η distribution shapes; however, the FMOS-COSMOS and MOSDEF mocks show a systematically low median and scatter in η in comparison to the observations. This work also presents the novel approach of analytically deriving the relationship between the intrinsic scatter in η , scatter added by measurement uncertainties, and observed scatter, resulting in an intrinsic burstiness measurement of 0.06–0.16 dex.
79 ASTRONOMY AND ASTROPHYSICS↗
Characterizing the Directionality of Gravitational Wave Emission from Matter Motions within Core-collapse Supernovae
Abstract We analyze the directional dependence of the gravitational wave (GW) emission from 15 3D neutrino radiation hydrodynamic simulations of core-collapse supernovae (CCSNe). Using spin weighted spherical harmonics, we develop a new analytic technique to quantify the evolution of the distribution of GW emission over all angles. We construct a physics-informed toy model that can be used to approximate GW distributions for general ellipsoid-like systems, and use it to provide closed form expressions for the distribution of GWs for different CCSN phases. Using these toy models, we approximate the protoneutron star (PNS) dynamics during multiple CCSN stages and obtain similar GW distributions to simulation outputs. When considering all viewing angles, we apply this new technique to quantify the evolution of preferred directions of GW emission. For nonrotating cases, this dominant viewing angle drifts isotropically throughout the supernova, set by the dynamical timescale of the PNS. For rotating cases, during core bounce and the following tens of milliseconds, the strongest GW signal is observed along the equator. During the accretion phase, comparable—if not stronger—GW amplitudes are generated along the axis of rotation, which can be enhanced by the low T /∣ W ∣ instability. We show two dominant factors influencing the directionality of GW emission are the degree of initial rotation and explosion morphology. Lastly, looking forward, we note the sensitive interplay between GW detector site and supernova orientation, along with its effect on detecting individual polarization modes.
Astronomy & Astrophysics↗
The Threatening Magnetic and Plasma Environment of the TRAPPIST-1 Planets
Recently, four additional Earth-mass planets were discovered orbiting the nearby ultracool M8 dwarf, TRAPPIST-1, making a remarkable total of seven planets with equilibrium temperatures compatible with the presence of liquid water on their surface. Temperate terrestrial planets around an M-dwarf orbit close to their parent star, rendering their atmospheres vulnerable to erosion by the stellar wind and energetic electromagnetic and particle radiation. Here, we use state-of-the-art 3D magnetohydrodynamic models to simulate the wind around TRAPPIST-1 and study the conditions at each planetary orbit. All planets experience a stellar wind pressure between 10{sup 3} and 10{sup 5} times the solar wind pressure on Earth. All orbits pass through wind pressure changes of an order of magnitude and most planets spend a large fraction of their orbital period in the sub-Alfvénic regime. For plausible planetary magnetic field strengths, all magnetospheres are greatly compressed and undergo much more dynamic change than that of the Earth. The planetary magnetic fields connect with the stellar radial field over much of the planetary surface, allowing the direct flow of stellar wind particles onto the planetary atmosphere. These conditions could result in strong atmospheric stripping and evaporation and should be taken into account for any realistic assessment of the evolution and habitability of the TRAPPIST-1 planets.
79 ASTRONOMY AND ASTROPHYSICS↗
Reverse Current Model for Coronal Mass Ejection Cavity Formation
In this work, we present a new model for explaining the three-part structure of coronal mass ejections (CMEs). The model proposes that the cavity in a CME forms because a rising electric current in the core prominence induces an oppositely directed electric current in the background plasma; this eddy current is required to satisfy the frozen-in magnetic flux condition in the background plasma. The magnetic force between the inner-core electric current and the oppositely directed induced eddy current propels the background plasma away from the core, creating a cavity and a density pileup at the cavity edge. The cavity radius saturates when an inward restoring force from magnetic and hydrodynamic pressure in the region outside the cavity edge balances the outward magnetic force. The model is supported by (i) laboratory experiments showing the development of a cavity as a result of the repulsion of an induced reverse current by a rising inner-core flux-rope current, (ii) 3D numerical magnetohydrodynamic (MHD) simulations that reproduce the laboratory experiments in quantitative detail, and (iii) an analytic model that describes cavity formation as a result of the plasma containing the induced reverse current being repelled from the inner core. This analytic model offers broad applicability because the predicted cavity widths are relatively independent of both the current injection mechanism and the injection timescale.
79 ASTRONOMY AND ASTROPHYSICS↗
Current Sheet Alignment in Oblique Black Hole Magnetospheres: A Black Hole Pulsar?
Abstract We study the magnetospheric evolution of a nonaccreting spinning black hole (BH) with an initially inclined split monopole magnetic field by means of 3D general relativistic magnetohydrodynamic simulations. This serves as a model for a neutron star (NS) collapse or a BH–NS merger remnant after the inherited magnetosphere has settled into a split monopole field creating a striped wind. We show that the initially inclined split monopolar current sheet aligns over time with the BH equatorial plane. The inclination angle evolves exponentially toward alignment, with an alignment timescale that is inversely proportional to the square of the BH angular velocity, where higher spin results in faster alignment. Furthermore, magnetic reconnection in the current sheet leads to exponential decay of event-horizon-penetrating magnetic flux with nearly the same timescale for all considered BH spins. In addition, we present relations for the BH mass and spin in terms of the period and alignment timescale of the striped wind. The explored scenario of a rotating, aligning, and reconnecting current sheet can potentially lead to multimessenger electromagnetic counterparts to a gravitational-wave event due to the acceleration of particles powering high-energy radiation, plasmoid mergers resulting in coherent radio signals, and pulsating emission due to the initial misalignment of the BH magnetosphere.
Astronomy & Astrophysics↗
Investigation on Underhood Thermal Analysis of Truck Platooning
This paper presents a combined aero-thermal computational fluid dynamic (CFD) evaluation of platooning medium duty commercial vehicles in two highway configurations. Thermal analysis comparison is made between an approach that includes vehicle drag reduction on engine heat rejection and one that does not by assuming a constant heat rejection based on open road conditions. The paper concludes that accounting for aerodynamic drag reduction on engine heat load provides a more real world evaluation than assuming a constant heat load based on open road conditions. A 3D CFD underhood thermal simulations are performed in two different vehicle platooning configurations; (i) single-lane and (ii) two-lane traffic conditions. The vehicle platooning consists of two identical vehicles, i.e. leading and trailing vehicle. In this work, heat exchangers are modeled by two different heat rejection rate models. In the first model, a constant heat rejection rates are considered as similar to no-traffic vehicle condition. In the other model, a varied heat rejection rates are implemented by considering an aerodynamic influence on fuel consumption rates. In a constant heat rejection rate model, the trailing vehicle thermal performance is significantly dropped in single-lane traffic due to reduced oncoming cold mass air flow velocities from leading vehicle. Also, the similar observations are found in two-lane traffic but at higher vehicle separation distances. In a varied heat rejection rate model, significant drop in temperature raise is observed in both leading and trailing vehicles when compared to a constant heat rejection rate model in the single-lane and the two-lane traffic conditions. In leading and trailing vehicles, the varied heat rejection rates are obtained from fuel consumption rates.
Platooning, underhood thermal analysis↗
Multimuons events and primary composition
Nucleon decay detectors at large depths offers now a total area larger than 1000 sq m to registrate muons of energy exceeding 1 TeV. Near complete high energy muon families are detected in those arrays. An extensive 3D Monte-Carlo simulation was conducted in view to understand the spatial distribution of those events and the possible link with elementary act or primary composition. As pion or kaon parents have a very small decay probability at so high energy, multimuon phenomena occurs at high altitude where the atmospheric density is small after the most energetic collisions.
Acharya, B. S.↗
Conservative boundary conditions for 3D gas dynamics problems
A method is described for 3D-gas dynamics computer simulation in regions of complicated shape by means of nonadjusted rectangular grids providing unified treatment of various problems. Some test problem computation results are given.
Gerasimov, B. P.↗
Io's volcanic and sublimation atmospheres
Fully 3D axisymmetric gasdynamic equations simulating SO2 and H2S frost sublimation and SO2 dayside and nightside volcanic atmospheres on Io are numerically solved, using a time-explicit finite-volume formulation. Both the sublimation and volcanic atmospheres generate horizontal supersonic winds away from the subsolar point or the volcanic vent. While the sublimation atmosphere is primarily driven by horizontal pressure gradients determined by surface temperatures, the volcanic atmosphere is driven by pressure gradients that are determined by the source rate. Sublimation and condensation produce patterns of surface deposits which are characteristic of the two types of atmospheres. The volcanic model is quantitatively consistent with Voyager observations of ring deposits.
Moreno, Miguel A.↗
Multidisciplinary tailoring of hot composite structures
A computational simulation procedure is described for multidisciplinary analysis and tailoring of multilayered multimaterial hot composite engine structural components subjected to simultaneous multiple discipline-specific thermal, structural, vibration, and acoustic loadings including the effect of aggressive environments. The simulation is based on a 3D finite element analysis technique in conjunction with structural mechanics codes, thermal/acoustic analysis methods, and tailoring procedures. The integrated multidisciplinary simulation procedure is general-purpose including the coupled effects of nonlinearities in structure geometry, material, loading, and environmental complexities. The composite material behavior is assessed at all composite scales, i.e., the laminate/ply/constituents (fiber/matrix), via a nonlinear material characterization hygro-thermomechanical model. Sample tailoring cases exhibiting nonlinear material/loading/environmental behavior of aircraft engine fan blades, are presented. The various multidisciplinary loadings lead to different tailored designs, even those opposite of each other, as in the case of minimum material cost versus minimum structure weight and in the case of minimum vibration frequency versus minimum acoustic noise.
Singhal, Surendra N.↗
Three-dimensional calculations of supersonic reacting flows using an LU scheme
An implicit finite volume lower-upper time-marching method which efficiently solves the complete Navier-Stokes and specied equations in a fully coupled fashion is the basis of the present 3D numerical program for simulating the supersonic reacting flows of H2 in air. The chemistry model incorporated has nine species and 18 reaction steps. Calculations are presented for flowfields of underexpanded hydrogen jets that are transversely injected into the supersonic airstream within scramjet combustors; the shock structure, separated flow regions around the injector, and combustion-product distributions are clearly represented.
Yu, Sheng-Tao↗
Variations of a global constraint factor in cracked bodies under tension and bending loads
Elastic-plastic finite-element analyses were used to calculate stresses and displacements around a crack in finite-thickness plates for an elastic-perfectly plastic material. Middle- and edge-crack specimens were analyzed under tension and bending loads. Specimens were 1.25 to 20 mm thick with various widths and crack lengths. A global constraint factor alpha(sub g), an averaged normal-stress to flow-stress ratio over the plastic region, was defined to simulate three-dimensional (3D) effects in two-dimensional (2D) models. For crack lengths and uncracked ligament lengths greater than four times the thickness, the global constraint factor was found to be nearly a unique function of a normalized stress-intensity factor (related to plastic-zone size to thickness ratio) from small- to large-scale yielding conditions for various specimen types and thickness. For crack length-to-thickness ratios less than four, the global constraint factor was specimen type, crack length and thickness dependent. Using a 2D strip-yield model and the global constraint factors, plastic-zone sizes and crack-tip displacements agreed reasonably well with the 3D analyses. For a thin sheet aluminum alloy, the critical crack-tip-opening angle during stable tearing was found to be independent of specimen type and crack length for crack length-to-thickness ratios greater than 4.
Newman, J. C., Jr.↗
Experimental investigation of the relationship between impact crater morphology and impacting particle velocity and direction
Interpretation of the wealth of impact data available from the Long Duration Exposure Facility, in terms of the absolute and relative populations of space debris and natural micrometeoroids, requires three dimensional models of the distribution of impact directions, velocities and masses of such particles, as well as understanding of the impact processes. Although the stabilized orbit of LDEF provides limited directional information, it is possible to determine more accurate impact directions from detailed crater morphology. The applicability of this technique has already been demonstrated but the relationship between crater shape and impactor direction and velocity has not been derived in detail. We present the results of impact experiments and simulations: (1) impacts at micron dimensions using the Unit's 2MV Van de Graaff accelerator; (2) impacts at mm dimensions using a Light Gas Gun; and (3) computer simulations using AUTODYN-3D from which an empirical relationship between crater shape and impactor velocity, direction and particle properties we aim to derive. Such a relationship can be applied to any surface exposed to space debris or micrometeoroid particles for which a detailed pointing history is available.
Mackay, N. G.↗
Study of a Wake Recovery Mechanism in a High-Speed Axial Compressor Stage
This work addresses the significant differences in compressor rotor wake mixing loss which exist in a stage environment relative to a rotor in isolation. The wake decay for a rotor in isolation is due solely to viscous dissipation which is an irreversible process and thus leads to a loss in both total pressure and efficiency. Rotor wake decay in the stage environment is due to both viscous mixing and the inviscid strain imposed on the wake fluid particles by the stator velocity field. This straining process, referred to by Smith (1993) as recovery, is reversible and for a 2D rotor wake leads to an inviscid reduction of the velocity deficit of the wake. A model for the rotor wake decay process is developed and used to quantify the viscous dissipation effects relative to those of inviscid wake stretching. The model is verified using laser anemometer measurements acquired in the wake of a transonic rotor operated in isolation and in a stage configuration at near peak efficiency and near stall operating conditions. Additional insight is provided by a time-accurate 3D Navier-Stokes simulation of the compressor stator flow field at the corresponding stage loading levels. Results from the wake decay model exhibit good agreement with the experimental data. Data from the model, laser anemometer measurements, and numerical simulations indicate that for the rotor/stator spacing used in this work, which is typical of core compressors, rotor wake straining (stretching) is the primary decay process in the stator passage with viscous mixing playing only a minor role. The implications of these results on compressor stage design are discussed.
VanZante, Dale E.↗