Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “population flow dynamics”

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 109 records · Page 6

The role of O-type neutral lines in magnetic merging during substorms and solar flares

An analytical model of O-type configurations was derived, for visualizing the geometry and for numerical treatment of plasma flows. It was found that at a certain distance from the neutral line the mean particle motion became decoupled from that of magnetic field lines (which obey the MHD condition). The decoupling distance depended on initial conditions in momentum space, suggesting that the MHD approximation which averages out such conditions may not suffice for describing plasma dynamics near the neutral line. It was also found that after inflowing particles are decoupled from the field line motion, they go over to a mode of runaway acceleration along the neutral line. It is concluded that if merging occurs at an X-O pair, two particle populations may be expected - low energy particles accelerated adiabatically by earthward convection past the X-type line and high-energy particles convected towards the O-type line. The second acceleration process depends critically on the rapidity of merging and is therefore expected to vary considerably from event to event.

Stern, D. P.↗

Magnetospheric imaging of high latitude ion outflows

High latitude ion outflows mostly consist of upward streaming O(+) and He(+) emanating from the ionosphere. At heights above 1000 km, these flows consist of cold and hot components which resonantly scatter solar extreme ultraviolet (EUV) light, however, the ion populations respond differently to Doppler shifting resulting from the large relative velocities between the ions and the Sun. The possibility of optical detection of the Doppler effect on the scattering rate will be discussed for the O(+) (83.4 nm) ions. We have contrasted the EUV solar resonance images of these outflows by simulations of the 30.4 nm He(+) and 83.4 nm O(+) emissions for both quiet and disturbed geomagnetic conditions. Input data for the 1000 km level has been obtained from the EICS instrument aboard the Dynamics Explorer (DE) satellite. Our results show emission rates of 50 and 56 milli-Rayleighs at 30.4 nm for quiet and disturbed conditions and 65 and 75 milli-Rayleighs at 83.4 nm for quiet and disturbed conditions, respectively, obtained for a polar orbiting satellite and viewing radially outward. We also find that an imager at an equatorial distance of 9 R(sub E) or more is in a favorable position for detecting ion outflows, particularly when the plasmapause is depressed in latitude. However, an occultation disk is necessary to obscure the bright plasmaspheric emissions.

Garrido, D. E.↗

Particle injection and the structure of energetic-particle-modified shocks

A macroscopic 'self-consistent' nonlinear two-fluid model is developed for energetic-particle-modeled shocks. The model incorporates particle exchange between a thermal gas and an energetic population via a 'thermal leakage' mechanism. In the model, injection is regarded as a process whereby particles in the cooler thermal core of a total particle distribution are energized sufficiently to cross a 'momentum boundary' into the wings of the distribution. Particles in the wings are identified as 'energetic particles' or 'cosmic rays' and are assumed to propagate according to the diffusive transport equation of cosmic ray theory. Thermal particles are energized either as a result of adiabatic compression in a decelerating flow or directly via anomalous heating within a subshock. It is shown that the nature and model of particle injection directly affects the nature and dynamics of the shock, and hence the efficiency of particle acceleration, and that the shock itself regulates particle injection dynamically.

Zank, G. P.↗

Energetic ion beams at the edge of the plasma sheet - ISEE 1 observations plus a simple explanatory model

It has been found that energetic ions with energies in the range from 10 keV to 2 MeV are an important part of the geomagnetic tail particle population. A description is presented of energetic ion, three-dimensional distribution observations from ISEE 1 during 0400-2000 UT, May 3, 1978, a time when several encounters with the geomagnetic tail plasma sheet occurred. Pronounced energetic ion streaming is observed at each entry to and exit from the plasma sheet. Azimuthal gradients observed during these times indicate that this streaming occurs within two gyroradii of the plasma sheet edge. Streaming characteristics show both narrow and broad angular flows at arbitrary angles to B. These flows are observed simultaneously and individually. The time, energy, and pitch angle evolution of the ion fluxes show that ISEE 1 is encountering a dynamic and not a static equilibrium situation at the edge of the plasma sheet.

Williams, D. J.↗

Variances of the components and magnitude of the polar heliospheric magnetic field

The heliolatitude dependences of the variances in the components and the magnitude of the heliospheric magnetic field have been analysed, using the Ulysses magnetic field observations from close to the ecliptic plane to 80 southern solar latitude. The normalized variances in the components of the field increased significantly (by a factor about 5) as Ulysses entered the purely polar flows from the southern coronal hole. At the same time, there was at most a small increase in the variance of the field magnitude. The analysis of the different components indicates that the power in the fluctuations is not isotropically distributed: most of the power is in the components of the field transverse to the radial direction. Examining the variances calculated over different time scales from minutes to hours shows that the anisotropy of the field variances is different on different scales, indicating the influence of the two distinct populations of fluctuations in the polar solar wind which have been previously identified. We discuss these results in terms of evolutionary, dynamic processes as a function of heliocentric distance and as a function of the large scale geometry of the magnetic field associated with the polar coronal hole.

Balogh, A.↗

statemodify: a Python framework to facilitate accessible exploratory modeling for discovering drought vulnerabilities

The Colorado River Basin (CRB) is experiencing an unprecedented water shortage crisis brought upon by a combination of factors arising from interactions across the region’s coupled human and natural systems. Allocation of water to the seven states that rely on the Colorado River was settled in the Colorado River Compact of 1922 during a period now known to be characterized by atypically high flows (Christensen et al., 2004). Since then, aridification due to anthropogenic-driven warming has steadily reduced the overall water supply available in the basin, with a 10% decrease in the river’s flow occurring over just the past two decades (Bass et al., 2023). The river is further strained by increasing demands associated with a growing population and diverse multi-sectoral demands. Navigating these challenges also requires accounting for the complex prior appropriation water rights system governing water allocation across the region’s diverse users.

54 ENVIRONMENTAL SCIENCES↗

Comprehensive Characterization of Water Group Ion Composition and Distributions in Saturn's Magnetosphere With Cassini Plasma Spectrometer Data

Saturn's magnetosphere is continuously supplied with neutrals from the Enceladus plume and the icy rings, which undergo ionization and charge-exchange to form a complex water-group plasma environment. While the Cassini Plasma Spectrometer (CAPS) instrument has provided extensive compositional information, detailed separation of individual water-group ion species in time-of-flight (TOF) data has not previously been achieved. In this study, we perform forward modeling of CAPS-IMS energy-per-charge (E/Q) and TOF spectra obtained between 2004 and 2012 to resolve O + , OH + , H 2 O + , and H 3 O + and to characterize their plasma properties, including number density, temperature, and thermodynamic κ. Our results demonstrate that O + is the dominant thermal ion species throughout Saturn's magnetosphere, comprising up to ∼70% of the total ion population beyond ∼5 Saturn radii (R S ). In contrast, molecular ions such as OH + , H 2 O + , and H 3 O + dominate closer to Enceladus but rapidly dissociate into atomic ions between ∼5 and 10 R S . This radial region is also characterized by the steepest increase in plasma flow speed, which rises from ∼40% to ∼80% of rigid corotation. Simultaneously, ion velocity distributions approach Maxwell–Boltzmann equilibrium, as indicated by high kappa values. These findings provide new constraints on the ion–neutral chemistry that regulates the balance between molecular and atomic ions in Saturn's magnetosphere. They also emphasize the critical role of the 5–10 R S region as a transition zone for both plasma composition and dynamics. Our results refine previous CAPS-based studies and underscore the need to incorporate seasonal variability and ionospheric coupling into future global models of Saturn's plasma environment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Interrelated structures of the transport shock and collisional relaxation layer in a multitemperature, multilevel ionized gas

The gas dynamic structures of the transport shock and the downstream collisional relaxation layer are evaluated for partially ionized monatomic gases. Elastic and inelastic collisional nonequilibrium effects are taken into consideration. In the microscopic model of the atom, three electronic levels are accounted for. By using an asymptotic technique, the shock morphology is found on a continuum flow basis. This procedure gives two distinct layers in which the nonequilibrium effects to be considered are different. A transport shock appears as the inner solution to an outer collisional relaxation layer. The results show four main interesting points: (1) on structuring the transport shock, ionization and excitation rates must be included in the formulation, since the flow is not frozen with respect to the population of the different electronic levels; (2) an electron temperature precursor appears at the beginning of the transport shock; (3) the collisional layer is rationally reduced to quadrature for special initial conditions, which (4) are obtained from new Rankine-Hugoniot relations for the inner shock.

Vinolo, A. R.↗

Cold plasma distribution above a few thousand kilometers at high latitudes

Great advances with respect to the study of the earth's magnetosphere could be achieved in connection with the launch of the Dynamics Explorer (DE) 1 and 2 spacecraft. These advances were made possible partly because of the unique combination of coplanar orbits which simultaneously sample the low altitude ionospheric and atmospheric signature (DE-2) and the high altitude phenomena of the inner magnetosphere (DE-1). It was, thus, feasible to obtain new measurements of the coupling of plasmas and fields between these fundamentally important regions. One basic element of the coupling involves the interchange of low energy plasma between the ionosphere and magnetosphere. The flow characteristics of the low energy plasma are indicators of the magnetospheric and ionospheric electric fields which drive the current system. The Retarding Ion Mass Spectrometer (RIMS) instrument has been designed for conducting measurements regarding this plasma population. Attention is given to details regarding RIMS, an example RIMS orbit, plasma trough characteristics, and polar cap phenomena.

Chappell, C. R.↗

Effect of Dietary Countermeasures and Impact of Gravity on Renal Calculi Size Distributions Predicted by PBE-System and PBE-CFD Models

Renal stone disease is not only a concern on earth but can conceivably pose a serious risk to the astronauts health and safety in Space. In this work, two different deterministic models based on a Population Balance Equation (PBE) analysis of renal stone formation are developed to assess the risks of critical renal stone incidence for astronauts during space travel. In the first model, the nephron is treated as a continuous mixed suspension mixed product removal crystallizer and the PBE for the nucleating, growing and agglomerating renal calculi is coupled to speciation calculations performed by JESS. Predictions of stone size distributions in the kidney using this model indicate that the astronaut in microgravity is at noticeably greater but still subcritical risk and recommend administration of citrate and augmented hydration as effective means of minimizing and containing this risk. In the second model, the PBE analysis is coupled to a Computational Fluid Dynamics (CFD) model for flow of urine and transport of Calcium and Oxalate in the nephron to predict the impact of gravity on the stone size distributions. Results presented for realistic 3D tubule and collecting duct geometries, clearly indicate that agglomeration is the primary mode of size enhancement in both 1g and microgravity. 3D numerical simulations seem to further indicate that there will be an increased number of smaller stones developed in microgravity that will likely pass through the nephron in the absence of wall adhesion. However, upon reentry to a 1g (Earth) or 38g (Mars) partial gravitational fields, the renal calculi can lag behind the urinary flow in tubules that are adversely oriented with respect to the gravitational field and grow agglomerate to large sizes that are sedimented near the wall with increased propensity for wall adhesion, plaque formation, and risk to the astronauts.

Renal Stones↗

Valued peaks: Sustainable water allocation for small hydropower plants in an era of explicit ecological needs

Optimizing hydropower operations to balance economic profitability and support functioning ecosystem services is integral to river management policy. In this article, we propose a dynamic, constrained optimization framework for small hydropower plants (SHPs) to evaluate trade-offs between economic profitability and socio-ecological requirements. Specifically, we examine the balance between short-term losses in hydropower generation and the potential for compensatory benefits in the form of revenue from recreational ecosystem services, irrespective of the direct beneficiary. Our framework integrates a fish habitat model, a hydropower optimization model, and a recreational ecosystem service estimate to evaluate different environmental flow scenarios. The optimization process gives three outflow release scenarios, informed by previous streamflow realisations (dam inflow), and designed environmental flow constraints. The framework is applied and tested for the river Kuusinkijoki in North-eastern Finland, which is a habitat for migratory brown trout and grayling populations. We show that the revenue loss due to the environmental flow constraints arises through a reduction in revenue per generated energy unit and through a reduction in turbine efficiency. Additionally, the simulation results reveal that all the designed environmental flow constraints cannot be met simultaneously. Under the environmental flow scenario with both minimum flow and flow ramping rate constraints, the annual hydropower revenue decreases by 16.5 %. An annual increase of 8 % in recreational fishing visits offsets the revenue loss. In conclusion, the developed framework provides knowledge of the costs and benefits of hydropower environmental flow constraints and guides the prioritizing process of environmental measures.

13 HYDRO ENERGY↗

Determining the impacts of venoarterial extracorporeal membrane oxygenation on cerebral oxygenation using a one-dimensional blood flow simulator

Venoarterial extracorporeal membrane oxygenation (VA-ECMO) is a mechanical system that provides rapid and short-term support for patients with cardiac failure. In many patients, pulmonary function is also impaired, resulting in poorly-oxygenated cardiac outflow competing against well-oxygenated VA-ECMO outflow, a condition known as North-South syndrome. North-South syndrome is a primary concern because of its potential to cause cerebral hypoxia, which has a critical influence on neurological complications often seen in this patient population. In order to reduce ischemic neurological complications, it is important to understand how clinical decisions regarding VA-ECMO parameters influence blood oxygenation. In our work, we studied the impacts of flow rate and cannulation site on oxygenation using a one-dimensional (1D) model to simulate blood flow. Our model was initially tested by comparing blood flow results to those observed from experimental work in VA-ECMO patients. The 1D model was combined with a two-phase flow model to simulate oxygenation. Additionally, the influence of various other clinician-tunable parameters on oxygenation in the common carotid arteries (CCAs) were tested, including, blood viscosity, cannula position within the insertion artery, heart rate, and systemic vascular resistance (SVR), as well as geometrical changes such as arterial radius and length. Our results indicated that blood oxygenation to the brain strongly depended on the cannula insertion site and the VA-ECMO flow rate with a weaker but potentially significant dependence on arterial radius. During femoral cannulation, VA-ECMO flow rates greater than ~4.9 L/min were needed to perfuse the CCAs. However, axillary and central cannulation began to perfuse the CCAs at significantly lower flow (~1 L/min). These results may help explain the incidence of cerebral hypoxia in this patient population and the common need to change cannulation strategies during treatment to address this clinical problem. While this work describes patient-averaged results, determining these relationships between VA-ECMO parameters and cerebral hypoxia is an important step towards future work to develop patient-specific models that clinicians can use to improve outcomes.

1D blood flow↗

Distribution of centrality measures on undirected random networks via the cavity method

The Katz centrality of a node in a complex network is a measure of the node’s importance as far as the flow of information across the network is concerned. For ensembles of locally tree-like undirected random graphs, this observable is a random variable. Its full probability distribution is of interest but difficult to handle analytically because of its “global” character and its definition in terms of a matrix inverse. Leveraging a fast Gaussian Belief Propagation-Cavity algorithm to solve linear systems on tree-like structures, we show that i) the Katz centrality of a single instance can be computed recursively in a very fast way, and ii) the probability P ( K ) that a random node in the ensemble of undirected random graphs has centrality K satisfies a set of recursive distributional equations, which can be analytically characterized and efficiently solved using a population dynamics algorithm. We test our solution on ensembles of Erdős-Rényi and Scale Free networks in the locally tree-like regime, with excellent agreement. The analytical distribution of centrality for the configuration model conditioned on the degree of each node can be employed as a benchmark to identify nodes of empirical networks with over- and underexpressed centrality relative to a null baseline. We also provide an approximate formula based on a rank- 1 projection that works well if the network is not too sparse, and we argue that an extension of our method could be efficiently extended to tackle analytical distributions of other centrality measures such as PageRank for directed networks in a transparent and user-friendly way.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Extreme elevational migration spurred cryptic speciation in giant hummingbirds

The ecoevolutionary drivers of species niche expansion or contraction are critical for biodiversity but challenging to infer. Niche expansion may be promoted by local adaptation or constrained by physiological performance trade-offs. For birds, evolutionary shifts in migratory behavior permit the broadening of the climatic niche by expansion into varied, seasonal environments. Broader niches can be short-lived if diversifying selection and geography promote speciation and niche subdivision across climatic gradients. To illuminate niche breadth dynamics, we can ask how “outlier” species defy constraints. Of the 363 hummingbird species, the giant hummingbird (Patagona gigas) has the broadest climatic niche by a large margin. To test the roles of migratory behavior, performance trade-offs, and genetic structure in maintaining its exceptional niche breadth, we studied its movements, respiratory traits, and population genomics. Satellite and light-level geolocator tracks revealed an >8,300-km loop migration over the Central Andean Plateau. This migration included a 3-wk, ~4,100-m ascent punctuated by upward bursts and pauses, resembling the acclimatization routines of human mountain climbers, and accompanied by surging blood-hemoglobin concentrations. Extreme migration was accompanied by deep genomic divergence from high-elevation resident populations, with decisive postzygotic barriers to gene flow. The two forms occur side-by-side but differ almost imperceptibly in size, plumage, and respiratory traits. The high-elevation resident taxon is the world’s largest hummingbird, a previously undiscovered species that we describe and name here. The giant hummingbirds demonstrate evolutionary limits on niche breadth: when the ancestral niche expanded due to evolution (or loss) of an extreme migratory behavior, speciation followed.

Science & Technology - Other Topics↗

The Characterization of Secondary Interstellar Neutral Oxygen Beyond the Heliopause: A Detailed Analysis of the IBEX-Lo Oxygen Observations

In this study, we analyze the directional distribution of the secondary interstellar neutral (ISN) O population observed by the IBEX-Lo neutral atom camera on the Interstellar Boundary EXplorer (IBEX) via the comparison with simulated ISN O intensity maps produced by an analytical model. In the analytical model, we assume that there are primary and secondary ISN populations at the heliopause. We further assume that each population is represented by a Maxwellian velocity distribution function with its own flow parameters. For the viewing directions of IBEX-Lo, we compute the incoming atom speeds at the heliopause with a Keplerian equation of motion in the solar gravity field. Then, we calculate analytically the distribution function to obtain the ISN intensities at Earth’s orbit. We compare the simulated O intensity maps with the IBEX-Lo O sky map to determine the most likely flow parameters of the secondary ISN O population. Using this method, we find the most likely flow parameters of the secondary ISN O population: V(sub SecISNO) = 11 ± 2.2 km s(exp -1), λ(sub SecISNO) = 67° ± 1°.5, β(sub SecISNO) = -12° ± 1°.6, and T(sub SecISNO) = 10,000 ± 1500 K. The results indicate that the secondary ISN O flow direction is deflected toward lower ecliptic longitude and higher negative ecliptic latitude from the ISN gas flow direction at the heliopause. The secondary ISN O flow direction is more deflected from the ISN gas flow direction than the secondary ISN He flow direction.

Interstellar matter clouds↗

Radiative Feedback from Primordial Protostars and Final Mass of the First Stars

In this contribution, we review our efforts toward understanding the typical mass-scale of primordial stars. Our direct numerical simulations show that, in both of Population III.1 and III.2 cases, strong UV stellar radiative feedback terminatesmass accretion onto a protostar.AnHII region formed around the protostar very dynamically expands throughout the gas accreting envelope, which cuts off the gas supply to a circumstellar disk. The disk is exposed to the stellar UV radiation and loses its mass by photoevaporation. The derived final masses are 43 Stellar Mass and 17 Stellar Mass in our fiducial Population III.1 and III.2 cases. Much more massive stars should form in other exceptional conditions. In atomic-cooling halos where H2 molecules are dissociated, for instance, a protostar grows via very rapid mass accretion with the rates M* approx. 0.1 - 1 Stellar Mass/yr. Our newstellar evolution calculations show that the protostar significantly inflates and never contracts to reach the ZAMS stage in this case. Such the "supergiant protostars" have very low UV luminosity, which results in weak radiative feedback against the accretion flow. In the early universe, supermassive stars formed through this process might provide massive seeds of supermassive black holes.

star formation↗

Semi-Supervised Learning of Lift Optimization of Multi-Element Three-Segment Variable Camber Airfoil

This chapter describes a new intelligent platform for learning optimal designs of morphing wings based on Variable Camber Continuous Trailing Edge Flaps (VCCTEF) in conjunction with a leading edge flap called the Variable Camber Krueger (VCK). The new platform consists of a Computational Fluid Dynamics (CFD) methodology coupled with a semi-supervised learning methodology. The CFD component of the intelligent platform comprises of a full Navier-Stokes solution capability (NASA OVERFLOW solver with Spalart-Allmaras turbulence model) that computes flow over a tri-element inboard NASA Generic Transport Model (GTM) wing section. Various VCCTEF/VCK settings and configurations were considered to explore optimal design for high-lift flight during take-off and landing. To determine globally optimal design of such a system, an extremely large set of CFD simulations is needed. This is not feasible to achieve in practice. To alleviate this problem, a recourse was taken to a semi-supervised learning (SSL) methodology, which is based on manifold regularization techniques. A reasonable space of CFD solutions was populated and then the SSL methodology was used to fit this manifold in its entirety, including the gaps in the manifold where there were no CFD solutions available. The SSL methodology in conjunction with an elastodynamic solver (FiDDLE) was demonstrated in an earlier study involving structural health monitoring. These CFD-SSL methodologies define the new intelligent platform that forms the basis for our search for optimal design of wings. Although the present platform can be used in various other design and operational problems in engineering, this chapter focuses on the high-lift study of the VCK-VCCTEF system. Top few candidate design configurations were identified by solving the CFD problem in a small subset of the design space. The SSL component was trained on the design space, and was then used in a predictive mode to populate a selected set of test points outside of the given design space. The new design test space thus populated was evaluated by using the CFD component by determining the error between the SSL predictions and the true (CFD) solutions, which was found to be small. This demonstrates the proposed CFD-SSL methodologies for isolating the best design of the VCK-VCCTEF system, and it holds promise for quantitatively identifying best designs of flight systems, in general.

Intelligent Systems↗

Numerical Model of the Plasma Sheath Generated by the Plasma Source Instrument Aboard the Polar Satellite

The plasma sheath generated by the operation of the Plasma Source Instrument (PSI) aboard the Polar satellite is studied by using a three-dimensional particle-in-cell (PIC) code. When the satellite passes through the region of low-density plasma, the satellite charges to positive potentials as high as 40-50 V, owing to the photoelectron emission. In such a case, ambient core ions cannot accurately be measured or detected. The goal of the onboard PSI is to reduce the floating potential of the satellite to a sufficiently low value so that the ions in the polar wind become detectable. When the PSI is operated, ion-rich xenon plasma is ejected from the satellite, such that the floating potential of the satellite is reduced and is maintained at approximately 2 V. Accordingly, in our three-dimensional PIC simulation we considered that the potential of the satellite is 2 V as a fixed bias. Considering the relatively high density of the xenon plasma in the sheath (10-10(exp 3)/cc), the ambient plasma of low density (<1/cc) is neglected. In the simulations the electric fields and plasma dynamics are calculated self-consistently. We found that an 'apple'-shape positive potential sheath forms surrounding the satellite. In the region near the PSI emission a high positive potential hill develops. Near the Thermal Ion Dynamics Experiment detector away from the PSI, the potentials are sufficiently low for the ambient polar wind ions to reach it. In the simulations it takes only about a couple of tens of electron gyroperiods for the sheath to reach a quasi steady state. This time is approximately the time taken by the heavy Xe(+) ions to expand up to about one average Larmor radius of electrons from the satellite surface. After this time the expansion of the sheath in directions transverse to the ambient magnetic field slows down because the electrons are magnetized. Using the quasi steady sheath, we performed trajectory calculations to characterize the detector response to a highly supersonic polar wind flow. The detected ions' velocity distribution shows significant deviations from a shifted Maxwellian in the ambient polar wind population. The deviations are caused by the effects of electric fields on the ions' motion as they traverse the sheath.

Singh, N.↗