SEARCH · Engineering Papers
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.
One of the most striking phenomena which accompany the flow of fluids such as air and water about bodies, is that of changes in character. This phenomena is still very little understood. We will first discuss the nature of these changes and then show that pure theory leads to a multiplicity of characters of flow, among which we will endeavor to indicate those bearing some analogy to experimental results.
This paper is concerned with the issue of potential non-unique solutions in three-phase state estimation. Theory of observability analysis for positive sequence power system state estimation is based on certain assumptions that avoid possibility of multiple solutions. Also, it is shown that observability of a positive sequence network remains independent of the network parameters or the operating state. When extending single-phase observability analysis directly to the three-phase case, this paper considers the possibility of converging to multiple solutions, i.e. solution non-uniqueness, even for cases where state estimator successfully converges. The study illustrates via numerical examples the likelihood of converging to entirely different solutions for certain network parameters. It also examines how the operating state, particularly under unbalanced loading, leads to solution non-uniqueness. The paper then describes an alternative approach to ensure a unique solution in three-phase state estimation. This method aims to accurately and uniquely estimate the state of any unbalanced three-phase system, irrespective of load imbalance, network configuration, existence of synchronous generators or transformers.
The High-Lift Common Research Model (HL-CRM) and the JAXA Standard Model (JSM) were analyzed computationally using both the OVERFLOW and LAVA codes for the third AIAA High-Lift Prediction Workshop. Geometry descriptions and the test cases simulated are described. With the HL-CRM, the effects of surface smoothness during grid projection and the effect of partially sealing a flap gap were studied. Grid refinement studies were performed at two angles of attack using both codes. For the JSM, simulations were performed with and without the nacelle/pylon. Without the nacelle/pylon, evidence of multiple solutions was observed when a quadratic constitutive relation is used in the turbulence modeling; however, using time-accurate simulation seemed to alleviate this issue. With the nacelle/pylon, no evidence of multiple solutions was observed. Laminar-turbulent transition modeling was applied to both JSM configuration, and had an overall favorable impact on the lift predictions.
The modeling of compressible multiphase flows is a decades-old area of study with many applications across various fields. Many of these application areas use stiff pressure relaxation. This process involves the solution of a nonlinear system with N + 1 equations and N + 1 unknowns, where N is the number of phases. The resolution of this system with general equations of state (EOSs) is difficult. Furthermore, nonlinear systems can admit multiple solutions, and current solution methods do not address this possibility. Very recently, a thermodynamic relaxation method was introduced, which effectively maps a relatively simple predictor equation of state onto a more complex target equation of state. In this context, the target EOSs are the chosen EOSs for the thermodynamic model. Furthermore, this thermodynamic relaxation has the benefit of simplifying the stiff pressure relaxation system of equations. In this article, we show this system reduces to a polynomial of degree N, which can be recast as an eigenvalue problem through the use of the associated companion matrix. We show that although this eigenvalue method is generally less efficient than Newton–Raphson iteration, it does not suffer from convergence issues and finds all N roots of the polynomial. Hence, the method provides a fail-safe for root-finding iterative methods and a way to address the issue of multiple solutions to the nonlinear system of equations in stiff pressure relaxation.
A daunting challenge in the design of lithium ion batteries (LiBs) is enabling 10-min extreme fast charging (XFC) while achieving appreciable charge acceptance and cycle life. This desirable outcome requires both a comprehensive understanding of LiB operation and aging behavior at different length scales and careful optimization. Lithium plating has been a critical bottleneck because, at XFC rates, it consumes cyclable lithium causing distinct aging and safety concerns even in moderate-loading LiBs. Here, we propose combining multiple solutions, including materials-to-electrode design-to-charging protocols, that are intended to overcome limitations in lithium-ion transport in the electrolyte phase, thus enabling 10-min XFC in LiBs. Some implemented strategies include cathode chemistry, optimized carbon binder domain in the cathode, dual layer anode design, improved separator and advanced electrolyte. Further, innovative charging protocols in moderately loading (~3 mAh/cm 2 anode/2.7 mAh/cm 2 cathode) single-layer pouch cells are proposed, together with demonstrated 10-min XFC with higher charge acceptance between 87.3 and 92.1% (or 2–2.1 mAh/cm 2 ) for 600 cycles without lithium plating. This methodical study with well-defined cells shows promise in combining multiple solution strategies to enable 10-min XFC, charting a pathway to achieve XFC in higher-loading energy-optimized LiBs.
We develop a new methodology for extracting Compton form factors (CFFs) from deeply virtual exclusive reactions such as the unpolarized DVCS cross section using a specialized inverse problem solver, a variational autoencoder inverse mapper (VAIM). The VAIM-CFF framework not only allows us access to a fitted solution set possibly containing multiple solutions in the extraction of all 8 CFFs from a single cross section measurement, but also accesses the lost information contained in the forward mapping from CFFs to cross section. We investigate various assumptions and their effects on the predicted CFFs such as cross section organization, number of extracted CFFs, use of uncertainty quantification technique, and inclusion of prior physics information. We then use dimensionality reduction techniques such as principal component analysis to visualize the missing physics information tracked in the latent space of the VAIM framework. Through re-framing the extraction of CFFs as an inverse problem, we gain access to fundamental properties of the problem not comprehensible in standard fitting methodologies: exploring the limits of the information encoded in deeply virtual exclusive experiments.
Plates can have more than one buckled solution for a fixed set of boundary conditions. The theory for the identification and the computation of multiple solutions in buckled plates is examined. The theory predicts modal interaction (which is also called change in buckle pattern or secondary buckling) in experiments on certain plates with multiple theoretical solutions. A set of coordinate functions is defined for Galerkin's method so that the von Karman plate equations are reduced to a coupled set of cubic equations in generalized coordinates that are uncoupled in the linear terms. An iterative procedure for solving modal interaction problems is suggested based on this cubic form.
An approach using a two-dimensional thin aerofoil, allied with the theory of viscous bluff-body separation, is used to study the initial cross-over from massive separation to an attached flow in a single-row unstaggered cascade. Analytic solutions are developed for the limit of small cascade-spacing. From the analytic solutions several interesting features of the cascade are examined, including multiple-solution branches and multiple regions of hysteresis. In addition, numerical results are presented for several selected aerofoils. Some of the aerofoils are found to contain markedly enlarged regions of hysteresis for certain critical cascade spacings.
The present invention includes an apparatus, system and method for stochastic sensing of an analyte to a protein pore. The protein pore may be an engineer protein pore, such as an ion channel at temperatures above 55.degree. C. and even as high as near 100.degree. C. The analyte may be any reactive analyte, including chemical weapons, environmental toxins and pharmaceuticals. The analyte covalently bonds to the sensor element to produce a detectable electrical current signal. Possible signals include change in electrical current. Detection of the signal allows identification of the analyte and determination of its concentration in a sample solution. Multiple analytes present in the same solution may also be detected.
Study results are presented showing performance capability of a spaceborne scatterometer to operationally measure ocean surface wind speed and direction. In addition, a research mode is described which will allow development of improved radar signatures for ocean, sea ice, and land targets. The study results show that a scatterometer can meet the operational requirements of + or - 2 m/s wind speed accuracy (or + or - 10%, whichever is greater) and + or - 20 deg wind direction accuracy over most of the expected ocean surface conditions. The six beam scatterometer design evaluated is shown to be skillful (greater than 90% correct) in specifying the correct wind vector solution (with a 180 deg ambiguity) from the multiple solutions derived; further improvement must rely on meteorological and pattern recognition techniques now under study.
Potential Solutions to Wildfires This project is looking at multiple aspects of the use of UAV’s during wildfire season. This project identifies multiple solutions to communication problems, early action response, and innovative tactics to fight wildfires during the dark hours of the second shift. This project will identify the features and innovations needed to create a second shift that can more effectively fight wildland fires at night. In addition, this project will look to find which satellite systems work best for detection of wildfires, with a special focus on remote locations with more rural communities. There will also be a focus on ways to improve current UAVs and sensors for aerial fire surveillance as well as the understanding on how to establish a safe airspace for these UAVs during wildfires. More specifically, unmanned ariel vehicles (UAV) or drones can be used during natural disasters and other situations to help first responders be a more efficient team, ultimately saving hundreds of more lives. Depending on the situation, drones can be equipped with many different accessories, each with their own advantages. When it comes to operating a drone, communication between the first responder, dispatcher, and the person flying the drone is very important when exchanging useful information. A focus on existing systems, and what changes could be made to create more novel systems will also be shown. Moreover, wildfires not only pose a hazard to human safety and a dent on communication methods, but it also has detrimental effects on the environment. Specifically, the use of common fire retardants have dangerous health effects on human health and the environment as a whole. Therefore, this project will also explore the use of alternative eco-friendly fire retardants and recent innovations to the use of sustainable retardants.
In this paper, we study the stability, accuracy and convergence behavior of various numerical schemes for phase-field modeling through a simple ODE model. Both theoretical analysis and numerical experiments are carried out on this ODE model to demonstrate the limitation of most numerical schemes that have been used in practice. One main conclusion is that the first-order fully implicit scheme is the only robust algorithm for phase-field simulations while all other schemes (that have been analyzed) may have convergence issue if the time step size is not exceedingly small. More specifically, by rigorous analysis in most cases, we have the following conclusions: (i) The first-order fully implicit scheme converges to the correct steady state solution for all time step sizes. In the case of multiple solutions, one of the solution branches always converges to the correct steady state solution. (ii) The first-order convex splitting scheme, which is equivalent to the first-order fully implicit scheme with a different time scaling, always converges to the correct steady state solution but may seriously lack numerical accuracy for transient solutions. (iii) For the second-order fully implicit and convex splitting schemes, for any time step size $δt > 0$, there exists an initial condition $u_0$ , with $|u_0| > 1$, such that the numerical solution converges to the wrong steady state solution. (iv) For $|u_0|$ ≤ 1, all second-order schemes studied in this paper converge to the correct steady state solution although severe numerical oscillations occur for most of them if the time step size is not sufficiently small. (v) An unconditionally energy-stable scheme (such as the modified Crank–Nicolson scheme) is not necessarily better than a conditionally energy-stable scheme (such as the Crank–Nicolson scheme). Finally, most, if not all, of the above conclusions are expected to be true for more general Allen–Cahn and other phase-field models.
Multiple solutions superimposed to resolve flows about complex configurations. Use of multiple, overset grids in computational fluid dynamics extends application of finite-difference methods to more complex configurations. Rather than trying to generate single mesh about all components of configuration, multiple, individual meshes used, then overset on major grid. Major grid used to resolve flow field or wrapped around main component.
The dual-frequency ratio of radar reflectivity factors (DFR) has been shown to be a useful quantity as it is independent of the number concentration of the particle size distribution and primarily a function of the mass-weighted particle diameter, Dm. A drawback of DFR-related methods for rain estimation, however, is the non-unique relationship between Dm and DFR. At Ku- and Ka-band frequencies, two solutions for Dm exist when DFR is less than zero. This ambiguity generates multiple solutions for the range profiles of the particle size parameters. We investigate characteristics of these solutions for both the initial-value (forward) and final-value (backward) forms of the equations. To choose one of among many possible range profiles of Dm, number concentration and rain rate, R, independently measured path attenuations are used. For the backward approach the possibility exists of dispensing with externally measured path attenuations by achieving consistency between the input and output path attenuations. The methods are tested by means of a simulation based on disdrometer-measured raindrop size distributions and the results are compared with a simplified version of the operational R-Dm method.
Procedures and solution strategies are described to solve the conventional structural optimization problem using the Lagrange multiplier technique. The multipliers, obtained through solution of an auxiliary nonlinear optimization problem, lead to optimality criteria to determine the design variables. It is shown that this procedure is essentially equivalent to an alternative formulation using a dual method Lagrangian function objective. Although mathematical formulations are straight-forward, successful applications and computational efficiency depend upon execution procedure strategies. Strategies examined, with application examples, include selection of active constraints, move limits, line search procedures, and side constraint boundaries.
Stationary linear perturbation responses to Northern Hemisphere orography are calculated in a quasi-geostrophic barotropic model in solid-body rotation. The stationary mountain torque induced by these perturbations is then used to construct graphical solutions to the steady-state wave, mean-flow interaction problem. It is shown that multiple solutions exist in the system and are near either the forcing equilibrium of the zonal forcing or near the resonance points in the system. Some of these near-resonance solutions have block-like configurations with a confluence zone upstream from a large-amplitude structure consisting of a high at high latitudes and a low at low latitudes. These block-like configurations are shown to be near stable solutions of the system. Time-dependent calculations show that the initial state and the zonal forcing equilibrium are important in determining the long-term time evolution of the system.
Discretization error is a common source of uncertainty in Computational Fluid Dynamics (CFD) analyses. Traditional means of controlling discretization error through fixed-mesh refinement studies has proven to be difficult particularly when modeling complex geometries and flow fields. One reason for this is that mesh generation in today’s production CFD workflow is often a labor intensive process that is heavily dependent on user judgment. Unstructured mesh adaptation is known to be an efficient way to control discretization errors in CFD. Adaptive methods replace user based decision making with automated processes that optimize a mesh to reduce discretization error. This paper compares the application of multiple solution adaptive techniques in combination with multiple flow solvers to solve for the flow field about a 2D airfoil section of the NASA High-Lift Common Research Model (HL-CRM). By driving the adaptive mesh processes to a similar level of mesh convergence, the ability to achieve consistent results between multiple adaptive techniques and flow solvers is demonstrated. Mesh convergence for the various adaptive mesh approaches is compared identifying potential areas for improvement and providing mesh generation guidance for future workshops.
Discretization error is a common source of uncertainty in Computational Fluid Dynamics (CFD) analyses. Traditional means of controlling discretization error through fixed-mesh refinement studies has proven to be difficult particularly when modeling complex geometries and flow fields. One reason for this is that mesh generation in today’s production CFD workflow is often a labor intensive process that is heavily dependent on user judgment. Unstructured mesh adaptation is known to be an efficient way to control discretization errors in CFD. Adaptive methods replace user based decision making with automated processes that optimize a mesh to reduce discretization error. This paper compares the application of multiple solution adaptive techniques in combination with multiple flow solvers to solve for the flow field about a 2D airfoil section of the NASA High-Lift Common Research Model (HL-CRM). By driving the adaptive mesh processes to a similar level of mesh convergence, the ability to achieve consistent results between multiple adaptive techniques and flow solvers is demonstrated. Mesh convergence for the various adaptive mesh approaches is compared identifying potential areas for improvement and providing mesh generation guidance for future workshops.