Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Conditional distributions”

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 1,207 records · Page 67

The Dependence of Cirrus Gamma Size Distributions Expressed as Volumes in N(sub 0)-Lambda-Mu Phase Space and Bulk Cloud Properties on Environmental Conditions: Results from the Small Ice Particles in Cirrus Experiment (SPARTICUS)

The variability of cirrus ice microphysical properties is investigated using observations obtained during the Small Particles in Cirrus (SPARTICUS) campaign. An existing approach that represents a size distribution (SD) as a single gamma function using an ellipsoid of equally realizable solutions in (N(sub 0), lambda, mu) phase space is modified to automatically identify multiple modes in SDs and characterize each mode by such an ellipsoid. The modified approach is applied to ice crystals with maximum dimension D greater than15 micrometers collected by the 2-D stereo and 2-D precipitation probes on the Stratton Park Engineering Company Learjet. The dependencies of N(sub 0), mu, and lambda from each mode, total number concentration, bulk extinction, ice water content (IWC), and mass median maximum dimension D(sub mm) as a function of temperature T and cirrus type are then analyzed. The changes in the observed codependencies between N(sub 0), mu, and lambda, bulk extinction, IWC, and D(sub mm) with environmental conditions indicate that particles were larger at higher T during SPARTICUS. At most two modes were observed in any SD during SPARTICUS, with the average boundary between them at 115 micrometers, similar to past studies not using probes with shatter mitigating tips and artifact removal algorithms. The bimodality of the SDs increased with T. This and the differences in N(sub 0), mu, and lambda between the modes suggest that particles with smaller D nucleated more recently than particles with larger D, which grew via vapor deposition and aggregation. Because smaller crystals, whose concentrations are uncertain, make marginal contributions to higher order moments, the use of higher moments for evaluating model fields is suggested.

Ellipsoids↗

Matrix Completion Using Alternating Minimization for Distribution System State Estimation

This paper examines the problem of state estimation in power distribution systems under low-observability conditions. The recently proposed constrained matrix completion method which combines the standard matrix completion method and power flow constraints has been shown to be effective in estimating voltage phasors under low-observability conditions using single-snapshot information. However, the method requires solving a semidefinite programming (SDP) problem, which becomes computationally infeasible for large systems and if multiple-snapshot (time-series) information is used. This paper proposes an efficient algorithm to solve the constrained matrix completion problem with time-series data. This algorithm is based on reformulating the matrix completion problem as a bilinear (non-convex) optimization problem, and applying the alternating minimization algorithm to solve this problem. This paper proves the summable convergence of the proposed algorithm, and demonstrates its efficacy and scalability via IEEE 123-bus system and a real utility feeder system. This paper also explores the value of adding more data from the history in terms of computation time and estimation accuracy.

41 EE - Solar Energy Technologies Office (EE-4S)↗

Matrix Completion Using Alternating Minimization for Distribution System State Estimation: Preprint

This paper examines the problem of state estimation in power distribution systems under low-observability conditions. The recently proposed constrained matrix completion method which combines the standard matrix completion method and power flow constraints has been shown to be effective in estimating voltage phasors under low-observability conditions using single-snapshot information. However, the method requires solving a semidefinite programming (SDP) problem, which becomes computationally infeasible for large systems and if multiple-snapshot (time-series) information is used. This paper proposes an efficient algorithm to solve the constrained matrix completion problem with time-series data. This algorithm is based on reformulating the matrix completion problem as a bilinear (non-convex) optimization problem, and applying the alternating minimization algorithm to solve this problem. This paper proves the summable convergence of the proposed algorithm, and demonstrates its efficacy and scalability via IEEE 123-bus system and a real utility feeder system. This paper also explores the value of adding more data from the history in terms of computation time and estimation accuracy.

41 EE - Solar Energy Technologies Office (EE-4S)↗

A statistical study of EMIC waves observed by Cluster: 2. Associated plasma conditions

This is the second in a pair of papers discussing a statistical study of electromagnetic ion cyclotron (EMIC) waves detected during 10 years (2001–2010) of Cluster observations. In the first paper, an analysis of EMIC wave properties (i.e., wave power, polarization, normal angle, and wave propagation angle) is presented in both the magnetic latitude (MLAT)-distance as well as magnetic local time (MLT)-L frames. In addition, this paper focuses on the distribution of EMIC wave-associated plasma conditions as well as two EMIC wave generation proxies (the electron plasma frequency to gyrofrequency ratio proxy and the linear theory proxy) in these same frames. Based on the distributions of hot H + anisotropy, electron and hot H+ density measurements, hot H + parallel plasma beta, and the calculated wave generation proxies, three source regions of EMIC waves appear to exist: (1) the well-known overlap between cold plasmaspheric or plume populations with hot anisotropic ring current populations in the postnoon to dusk MLT region; (2) regions all along the dayside magnetosphere at high L shells related to dayside magnetospheric compression and drift shell splitting; and (3) off-equator regions possibly associated with the Shabansky orbits in the dayside magnetosphere.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Criteria for stochastic pinning control of networks of chaotic maps

This paper investigates the controllability of discrete-time networks of coupled chaotic maps through stochastic pinning. In this control scheme, the network dynamics are steered towards a desired trajectory through a feedback control input that is applied stochastically to the network nodes. The network controllability is studied by analyzing the local mean square stability of the error dynamics with respect to the desired trajectory. Through the analysis of the spectral properties of salient matrices, a toolbox of conditions for controllability are obtained, in terms of the dynamics of the individual maps, algebraic properties of the network, and the probability distribution of the pinning control. We demonstrate the use of these conditions in the design of a stochastic pinning control strategy for networks of Chirikov standard maps. To elucidate the applicability of the approach, we consider different network topologies and compare five different stochastic pinning strategies through extensive numerical simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Self-field measurements of an HTS twisted stacked-tape cable conductor

For a twisted stacked-tape cable (TSTC) conductor composed of REBCO tapes, self-field measurements were conducted with Hall sensors. In the measurements, a 650 mm diameter single turn coil wound with the TSTC conductor, which was made with 48 REBCO tapes whose width was 6 mm, was utilized as a test sample. Based on the measurement results, the current distribution of the TSTC conductor was investigated with analytical models. The analytical results indicate the current distribution of the TSTC is uniform under the condition that the operating current is 10 kA and the sample temperature is approximately 30 K. On the other hand, the current distribution is not uniform at the excitation and the degauss of the TSTC conductor with the ramp rate of 50 A/s.

Obana, T.↗

Laser aberration signatures in expelled electrons from a tenuous gas

We describe a mechanism by which the aberration content of a focused, multiterawatt laser is imprinted upon the forward angular distribution of electrons ionized within and ponderomotively expelled from the focal volume. In our experiments, the laser aberration type and magnitude are controllably varied, and the measured electron distributions are correspondingly modified in a way consistent with predictions from numerical simulations. This imprint mechanism shows potential for enabling the development of accurate focal-spot characterization of intense lasers when fired at full power and is being developed for extension to petawatt laser systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Multiscale Computational Model for Particle Size Evolution during Coprecipitation of Li-Ion Battery Cathode Precursors

Next generation lithium ion batteries require higher energy and power density, which can be achieved by tailoring the cathode particle morphology, such as particle size, size distribution, and internal porosity. All these morphological features are determined during the cathode synthesis process, which consists of two steps, (i) coprecipitation and (ii) calcination. Transition metal hydroxide precursors are synthesized during the coprecipitation process, whereas their oxidation and lithiation occur during calcination. The size and size distribution of crystalline primary and aggregated secondary particles and their internal porosity are determined during coprecipitation. Operating conditions of the chemical reactor, such as solution pH, ammonia concentration, and stirring speed control the final morphological features. Here, a multiscale computational model has been developed to capture the nucleation and growth of crystalline primary particles and their aggregation into secondary transition metal hydroxide precursor particles. The simulations indicate that increasing solution pH and decreasing ammonia concentration lead to smaller sizes of the secondary particles. Finally, a phase map has been developed that can help identify the synthesis conditions needed for a specified particle size and size distribution.

25 ENERGY STORAGE↗

A Laboratory Model of a Hydrogen/Oxygen Engine for Combustion and Nozzle Studies

A small laboratory diagnostic thruster was developed to augment present low thrust chemical rocket optical and heat flux diagnostics at the NASA Lewis Research Center. The objective of this work was to evaluate approaches for the use of temperature and pressure sensors for the investigation of low thrust rocket flow fields. The nominal engine thrust was 110 N. Tests were performed at chamber pressures of about 255 kPa, 370 kPa, and 500 kPa with oxidizer to fuel mixture ratios between 4.0 and 8.0. Two gaseous hydrogen/gaseous oxygen injector designs were tested with 60 percent and 75 percent fuel film cooling. The thruster and instrumentation designs were proven to be effective via hot fire testing. The thruster diagnostics provided inner wall temperature and static pressure measurements which were compared to the thruster global performance data. For several operating conditions, the performance data exhibited unexpected trends which were correlated with changes in the axial wall temperature distribution. Azimuthal temperature distributions were found to be a function of operating conditions and hardware configuration. The static pressure profiles showed that no severe pressure gradients were present in the rocket. The results indicated that small differences in injector design can result in dramatically different thruster performance and wall temperature behavior, but that these injector effects may be overshadowed by operating at a high fuel film cooling rate.

Morren, Sybil Huang↗

Study of a Single-Power Two-Circuit ESR Process with Current-Carrying Mold: Mathematical Simulation of the Process and Experimental Verification

A novel single-power two-circuit ESR process (ESR-STCCM) with current-carrying mold has been investigated via numerical simulation and experimental research in this paper. A 2D quasi-steady-state mathematical model is developed to describe ESR-STCCM. The electromagnetic field, flow field, slag pool temperature distribution, and the shape of a molten steel pool in ESR-STCCM have been investigated by FLUENT software as well as user-defined functions (UDF). The results indicate that ESR-STCCM is different from the conventional ESR process. The maximum electromagnetic force, current density, Joule heat, and slag pool flow velocity are located in the lower part of the conductor in the ESR-STCCM process. The direction of the maximum electromagnetic force inclines upward. There are two distinct vortices in the slag pool. The larger swirl rotates counterclockwise near the conductor, with a value of 0.0263 m s{sup −1} due to the interaction of the electromagnetic force and gravity. The maximum temperature of the slag pool is 2070 K (1797 °C) and is located in the center of the swirl with a filling ratio of 0.6 and a 20 mm electrode immersion depth. The depth of a molten steel pool is shallower, which is conducive to improving solidification quality. In addition, the filling ratio of 0.6 is conducive to controlling steel solidification quality. Some experiments have been done, and the numerical model is confirmed by experimental results.

36 MATERIALS SCIENCE↗

Tilting Saturn without tilting Jupiter: Constraints on giant planet migration

The migration and encounter histories of the giant planets in our solar system can be constrained by the obliquities of Jupiter and Saturn. We have performed secular simulations with imposed migration and N-body simulations with planetesimals to study the expected obliquity distribution of migrating planets with initial conditions resembling those of the smooth migration model, the resonant Nice model and two models with five giant planets initially in resonance (one compact and one loose configuration). For smooth migration, the secular spin–orbit resonance mechanism can tilt Saturn’s spin axis to the current obliquity if the product of the migration timescale and the orbital inclinations is sufficiently large (exceeding 30 Myr deg). For the resonant Nice model with imposed migration, it is difficult to reproduce today’s obliquity values, because the compactness of the initial system raises the frequency that tilts Saturn above the spin precession frequency of Jupiter, causing a Jupiter spin–orbit resonance crossing. Migration timescales sufficiently long to tilt Saturn generally suffice to tilt Jupiter more than is observed. The full N-body simulations tell a somewhat different story, with Jupiter generally being tilted as often as Saturn, but on average having a higher obliquity. The main obstacle is the final orbital spacing of the giant planets, coupled with the tail of Neptune’s migration. The resonant Nice case is barely able to simultaneously reproduce the orbital and spin properties of the giant planets, with a probability ∼0.15%. The loose five planet model is unable to match all our constraints (probability <0.08%). The compact five planet model has the highest chance of matching the orbital and obliquity constraints simultaneously (probability ∼0.3%).

79 ASTRONOMY AND ASTROPHYSICS↗

A DLM/FD/IB method for simulating compound vesicle motion under creeping flow condition

In this article we present first a new distributed Lagrange multiplier/fictitious domain (DLM/FD) method for simulating fluid–particle interaction in Stokes flow. A conjugate gradient method driven by both pressure and distributed Lagrange multiplier, called one-shot method, has been developed to solve the discrete Stokes problem while enforcing the rigid body motion within the region occupied by the particle. The methodology is validated by comparing the numerical results of a neutrally buoyant particle of either a circular or elliptic shape with the associated Jeffery's solutions. We have successively combined the above methodology with an immersed boundary (IB) method and an elastic membrane modeled by a spring network to simulate the dynamics of a compound vesicle. In simple shear flow under creeping flow condition, the results are consistent with those obtained in literature. In Poiseuille flow, the compound vesicle motion is dominated by the motion of the vesicle membrane as expected and stays in the central region of the channel.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

In situ measurements of plasma properties during gas-condensation of Cu nanoparticles

Since the mean, standard deviation, and modality of nanoparticle size distributions can vary greatly between similar input conditions (e.g., power and gas flow rate), plasma diagnostics were carried out in situ using a double-sided, planar Langmuir probe to determine the effect the plasma has on the heating of clusters and their final size distributions. The formation of Cu nanoparticles was analyzed using cluster-plasma physics, which relates the processes of condensation and evaporation to internal plasma properties (e.g., electron temperature and density). Monitoring these plasma properties while depositing Cu nanoparticles with different size distributions revealed a negative correlation between average particle size and electron temperature. Furthermore, the modality of the size distributions also correlated with the modality of the electron energy distributions. It was found that the maximum cluster temperature reached during plasma heating and the material's evaporation point regulates the growth process inside the plasma. In the case of Cu, size distributions with average sizes of 8.2, 17.3, and 24.9 nm in diameter were monitored with the Langmuir probe, and from the measurements made, the cluster temperatures for each deposition were calculated to be 1028, 1009, and 863 K. These values are then compared with the onset evaporation temperature of particles of this size, which was estimated to be 1059, 1068, and 1071 K. Thus, when the cluster temperature is too close to the evaporation temperature, less particle growth occurs, resulting in the formation of smaller particles.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Microwave radiometer measurement of tidally induced salinity changes off the Georgia coast

A quasi-synoptic survey of tidally induced salinity changes off the Georgia coast was performed by using a L band microwave radiomater onboard a NASA aircraft. Salinity maps were obtained for ebb and flood conditions in order to define the salinity distributions near rivers and sounds and major changes that occur from ebb flow to flood flow. The Savannah River plume dominated the salinity regime and extended out from the Savannah River mouth about 12 km during ebb tidal conditions. The plume merged into a band of low salinity water extending along the Georgia-South Carolina coast which was produced by the many river sources of freshwater entering the coastal waters. The changes in salinity observed offshore of the river plume area were consistent with estimates of the changes that would occur over a typical tidal excursion perpendicular to the observed gradient.

Kendall, B. M.↗

Sequentially coupled flow and geomechanical simulation with a discrete fracture model for analyzing fracturing fluid recovery and distribution in fractured ultra-low permeability gas reservoirs

More accurate characterization and prediction of the in-situ distribution of fracturing fluid in fractured reservoirs are needed for enhancing well productivity. In this study, an implicit-sequentially coupled flow/geomechanics simulator incorporating an efficient discrete fracture model is developed to model fluid distribution and recovery performance of ultra-low permeability gas reservoirs. The finite-volume and finite-element methods are used for space discretization of the flow and geomechanics equations, respectively, while the backward Euler method is employed for time discretization. The flow and geomechanics equations are solved sequentially based on fixed-stress splitting. An efficient discrete-fracture model is used to explicitly model the fractured system. Flexible unstructured gridding is employed to model arbitrarily-oriented fractures. The interrelations among pore volume, permeability and geomechanical conditions are considered dynamically using two-way coupled flow and geomechanics computations. The geometry of fracture (networks) due to hydraulic fracturing has significant impacts on the fracturing fluid recovery efficiency and ensuing fluid distribution. Under the same injection volume, the fracturing fluid recovery is higher when the fracture geometry is planar. Fluid recovery is relatively lower whenever natural fractures are activated during fracturing treatments; flowback time is also shortened when complex fracture network with enlarged fracture interface is present. Fracturing fluid in hydraulic fractures may leak off into the natural fractures and subsequently imbibes into the surrounding matrix due to capillarity effects. The fracturing fluid recovery and in-situ fluid distribution are sensitive to the shut-in duration and fracture closure behavior. This study analyzes the coupled flow-geomechanical responses of fractured gas reservoirs during the post-fracturing periods. Understanding the fate of the fracturing fluid can provide insights on, to some extent, the stimulated fracture volume, size of the water invasion zone, and efficiency of the fracturing design. The simulation predictions can also provide more accurate initial reservoir conditions (e.g. distributions of different phases and pressure) for long-term well performance estimation.

42 ENGINEERING↗

Determination of Particle Size Distribution of Zirconium Dioxide Powders by Static Laser Scattering and Optical Microscopy

Dispersing conditions of powdered zirconium dioxide with addition of 3 mol % yttrium oxide are determined. The particle size distribution of the powders is determined by static laser scattering and optical microscopy. The zeta potential of suspension of zirconium dioxide particles as a function of pH of medium is measured by electrophoresis with titration. The influence of acidity on suspension stability during measurements is studied. It is demonstrated that a stable suspension can be obtained by shift of pH to both the acid and alkaline regions as well as by addition of surfactants, such as Dolapix CE 64.

36 MATERIALS SCIENCE↗

Velocity and extraneous boundary conditions of viscous flow problems

A novel method for treating certain troublesome boundary conditions in the numerical solution of time-dependent incompressible viscous flow problems is presented. This new method is developed on the basis of an integral representation for the velocity vector which contains the entire kinematics of the problem, including the boundary conditions of concern. It is shown that for the exterior flow problem the free-stream condition is satisfied at infinity exactly and the need to treat a far-field condition is removed by the use of the integral representation. The distribution of a nonvelocity variable on the solid boundary - i.e., the 'extraneous' boundary condition needed for both the exterior and the interior flows - is shown to be governed by the kinematics of the problem.

Wu, J. C.↗

Effect of non-equilibrium flow chemistry and surface catalysis on surface heating to AFE

The effect of nonequilibrium flow chemistry on the surface temperature distribution over the forebody heat shield on the Aeroassisted Flight Experiment (AFE) vehicle was investigated using a reacting boundary-layer code. Computations were performed by using boundary-layer-edge properties determined from global iterations between the boundary-layer code and flow field solutions from a viscous shock layer (VSL) and a full Navier-Stokes solution. Surface temperature distribution over the AFE heat shield was calculated for two flight conditions during a nominal AFE trajectory. This study indicates that the surface temperature distribution is sensitive to the nonequilibrium chemistry in the shock layer. Heating distributions over the AFE forebody calculated using nonequilibrium edge properties were similar to values calculated using the VSL program.

Stewart, David A.↗