Search NASA⌕ Search

SEARCH · Search NASA

Results for “fully linear model”

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 379 records · Page 21

Self-irradiated cooling condensations - The source of the optical line emission from cooling flows

The optical filaments seen at the centers of clusters of galaxies, which cannot be explained by standard photoionization or shock models, can be straightforwardly interpreted as self-irradiated cooling condensations. Fully developed (about 10,000 K) condensations embedded in thermally unstable 10 million K gas will be bathed in a powerful EUV/soft X-ray flux emanating from the surrounding condensing regions. The present models of the condensation/irradiation process give line ratios and luminosities similar to those observed. Because it is assumed that the condensations arise from linear perturbations, which grow significantly only when the local cooling time is less than the free-fall time, condensing regions of the type modeled here should arise only in cooling flows of central pressure less than about 1 billion ergs/cu cm and should extend only over the inner regions of the cluster.

Voit, G. Mark↗

Response analysis of a laminar premixed M-flame to flow perturbations using a linearized compressible Navier-Stokes solver

The response of a laminar premixed methane-air flame subjected to flow perturbations around a steady state is examined experimentally and using a linearized compressible Navier-Stokes solver with a one-step chemistry mechanism to describe combustion. The unperturbed flame takes an M-shape stabilized both by a central bluff body and by the external rim of a cylindrical nozzle. This base flow is computed by a nonlinear direct simulation of the steady reacting flow, and the flame topology is shown to qualitatively correspond to experiments conducted under comparable conditions. The flame is then subjected to acoustic disturbances produced at different locations in the numerical domain, and its response is examined using the linearized solver. This linear numerical model then allows the componentwise investigation of the effects of flow disturbances on unsteady combustion and the feedback from the flame on the unsteady flow field. It is shown that a wrinkled reaction layer produces hydrodynamic disturbances in the fresh reactant flow field that superimpose on the acoustic field. This phenomenon, observed in several experiments, is fully interpreted here. The additional perturbations convected by the mean flow stem from the feedback of the perturbed flame sheet dynamics onto the flow field by a mechanism similar to that of a perturbed vortex sheet. The different regimes where this mechanism prevails are investigated by examining the phase and group velocities of flow disturbances along an axis oriented along the main direction of the flow in the fresh reactant flow field. It is shown that this mechanism dominates the low-frequency response of the wrinkled shape taken by the flame and, in particular, that it fully determines the dynamics of the flame tip from where the bulk of noise is radiated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

New Results on Communication- and Memory-Aware Load Balancing Model and Algorithms

While load balancing in distributed-memory computing has been well-studied, we present an innovative approach to this problem: a unified, reduced-order model that combines three key components to describe “work” in a distributed system: computation, communication, and memory. Our model enables an optimizer to explore complex tradeoffs in task placement, such as augmented parallelism, at the expense of data replication increasing memory usage. We propose a fully distributed, heuristic-based load balancing optimization algorithm, and demonstrate that it quickly finds close-to-optimal solutions. We formalize the complex optimization problem as a mixed-integer linear program, and compare it to our strategy. Finally, we show that when applied to an electromagnetics code, our approach obtains up to 2.3x speedups for the imbalanced execution.

97 MATHEMATICS AND COMPUTING↗

Non-WKB evolution of solar wind fluctuations - A turbulence modeling approach

Recent observations indicate systematic changes in the interplanetary inertial range velocity-magnetic field correlation with increasing heliocentric distance. Here, a multiple length scale model for the evolution of the small scale fluctuations due to local nonlinear couplings as well as couplings to large scale gradients is outlined, which is similar to that used in turbulence modeling of inhomogeneous shear flows. A simple special case solution is given, indicating that couplings to large scale flow gradients can cause radial evolution of the type seen in observations. The strongest couplings of this type are linear and do not appear in the usual WKB orderings. The relationship of the present approach to WKB theory is not fully understood at present.

Ye, Zou↗

Vascular mechanics of the coronary artery

This paper describes our research into the vascular mechanics of the coronary artery and plaque. The three sections describe the determination of arterial mechanical properties using intravascular ultrasound (IVUS), a constitutive relation for the arterial wall, and finite element method (FEM) models of the arterial wall and atheroma. METHODS: Inflation testing of porcine left anterior descending coronary arteries was conducted. The changes in the vessel geometry were monitored using IVUS, and intracoronary pressure was recorded using a pressure transducer. The creep and quasistatic stress/strain responses were determined. A Standard Linear Solid (SLS) was modified to reproduce the non-linear elastic behavior of the arterial wall. This Standard Non-linear Solid (SNS) was implemented into an axisymetric thick-walled cylinder numerical model. Finite element analysis models were created for five age groups and four levels of stenosis using the Pathobiological Determinants of Atherosclerosis Youth (PDAY) database. RESULTS: The arteries exhibited non-linear elastic behavior. The total tissue creep strain was epsilon creep = 0.082 +/- 0.018 mm/mm. The numerical model could reproduce both the non-linearity of the porcine data and time dependent behavior of the arterial wall found in the literature with a correlation coefficient of 0.985. Increasing age had a strong positive correlation with the shoulder stress level, (r = 0.95). The 30% stenosis had the highest shoulder stress due to the combination of a fully formed lipid pool and a thin cap. CONCLUSIONS: Studying the solid mechanics of the arterial wall and the atheroma provide important insights into the mechanisms involved in plaque rupture.

NASA Program Biomedical Research and Countermeasur↗

Inflatable Softgoods Design of an Articulating Crew Transfer Tunnel

Future NASA exploration plans call for lunar and Martian surface systems that form a base camp of multiple, pressurized elements. These discrete components, including habitats and pressurized rovers, require interoperability to meet the Artemis Accords and standard docking systems to physically connect elements together. A pressurized, articulating crew transfer tunnel can be used between a rover and habitat to enable shirt-sleeve transfer of crew and cargo, saving valuable crew time and resources. While transfer tunnels have been described in the past, this work details the design of a structural softgoods system that has compliant capability through a proposed docking range of motion. The inflatable softgoods design is based on a zero-hoop stress shape, known as a Taylor surface, that is stacked and truncated to form a unique and flexible configuration. Analytical, non-linear models have been developed to examine and predict the behavior of the structure, and material testing was used to determine the properties that were used in the model. Finally, a sub-scale test article was constructed using the baseline design and pressurized testing is in work. Additional full-scale testing is planned for future years to fully demonstrate the capability of the system.

Tunnel↗

Inflatable Softgoods Design of an Articulating Crew Transfer Tunnel

Future NASA exploration plans call for lunar and Martian surface systems that form a base camp of multiple, pressurized elements. These discrete components, including habitats and pressurized rovers, require interoperability to meet the Artemis Accords and standard docking systems to physically connect elements together. A pressurized, articulating crew transfer tunnel can be used between a rover and habitat to enable shirt-sleeve transfer of crew and cargo, saving valuable crew time and resources. While transfer tunnels have been described in the past, this work details the design of a structural softgoods system that has compliant capability through a proposed docking range of motion. The inflatable softgoods design is based on a zero-hoop stress shape, known as a Taylor surface, that is stacked and truncated to form a unique and flexible configuration. Analytical, non-linear models have been developed to examine and predict the behavior of the structure, and material testing was used to determine the properties that were used in the model. Finally, a sub-scale test article was constructed using the baseline design and pressurized testing is in work. Additional full-scale testing is planned for future years to fully demonstrate the capability of the system.

Tunnel↗

Surface and Atmospheric Parameter Retrieval From AVIRIS Data: The Importance of Non-Linear Effects

AVIRIS data represent a new and important approach for the retrieval of atmospheric and surface parameters from optical remote sensing data. Not only as a test for future space systems, but also as an operational airborne remote sensing system, the development of algorithms to retrieve information from AVIRIS data is an important step to these new approaches and capabilities. Many things have been learned since AVIRIS became operational, and the successive technical improvements in the hardware and the more sophisticated calibration techniques employed have increased the quality of the data to the point of almost meeting optimum user requirements. However, the potential capabilities of imaging spectrometry over the standard multispectral techniques have still not been fully demonstrated. Reasons for this are the technical difficulties in handling the data, the critical aspect of calibration for advanced retrieval methods, and the lack of proper models with which to invert the measured AVIRIS radiances in all the spectral channels. To achieve the potential of imaging spectrometry, these issues must be addressed. In this paper, an algorithm to retrieve information about both atmospheric and surface parameters from AVIRIS data, by using model inversion techniques, is described. Emphasis is put on the derivation of the model itself as well as proper inversion techniques, robust to noise in the data and an inadequate ability of the model to describe natural variability in the data. The problem of non-linear effects is addressed, as it has been demonstrated to be a major source of error in the numerical values retrieved by more simple, linear-based approaches. Non-linear effects are especially critical for the retrieval of surface parameters where both scattering and absorption effects are coupled, as well as in the cases of significant multiple-scattering contributions. However, sophisticated modeling approaches can handle such non-linear effects, which are especially important over vegetated surfaces. All the data used in this study were acquired during the 1991 Multisensor Airborne Campaign (MAC-Europe), as part of the European Field Experiment on a Desertification-threatened Area (EFEDA), carried out in Spain in June-July 1991.

Green Robert O.↗

A non‐intrusive domain‐decomposition model reduction method for linear steady‐state partial differential equations with random coefficients

Abstract Domain decomposition methods have been proved to be an effective strategy to reduce the dimension of parametric partial differential equations (PDEs). However, existing domain decomposition methods for parametric PDEs are usually intrusive, which means domain decomposition based solvers need to be implemented from scratch for each target parametric PDE. To address this issue, we develop a new non‐intrusive domain‐decomposition model reduction method for linear steady‐state PDEs with random‐field coefficients. As a variant of our previous work by Mu and Zhang, the new method only needs access to the final linear system, that is, the global stiffness matrix and the right hand side, of a deterministic PDE solver, in order to build a domain‐decomposition‐based reduced model without intrusive implementation from scratch. The key idea is to remove the interface condition between sub‐domains and rely on the correlation between columns of the linear system to couple the sub‐domains. The non‐intrusive feature enables the applicability of the proposed method to a broader class of uncertainty quantification problems, where many legacy codes/solvers can be fully reused by our method. Two numerical examples including diffusion equations with random diffusivity and convection‐dominated transport with random velocity, are provided to demonstrate the effectiveness and efficiency of our method.

Zhang, Guannan↗

Aerodynamics and CFD analysis of equal size dual-rotor wind turbine

In this study, the linear momentum theory is used to analyze the power-extraction capability of dual-rotor wind turbines with equal-size rotors. The rotors of a dual-rotor wind turbine are modeled as two separate actuator disks. The stream tube encompassing the front rotor is modeled as two (inner and outer) stream tubes, with the rear rotor being fully enclosed within the front rotor inner stream tube. No assumption is made on airflow pressure in between the rotors. The effect of the front and rear rotor interaction on the airflow within the inner stream tube is included in the analysis. With the results obtained, axial thrusts on front and rear rotors are determined and later used as input for computational fluid dynamics simulation to determine flow characteristics across the rotors. Based on the flow pattern between the rotors, the total power coefficient of a dual-rotor wind turbine is related to the rotor separation distance. A general solution for the dual-rotor wind turbine is developed, which shows that the largest total power coefficient that can be obtained is 0.814, occurring at a rotor separation distance of 2.8 times the rotor diameter. Finally, the results of this study are compared with those obtained by similar studies reported in the literature. Discrepancies in the largest possible power coefficient of dual-rotor wind turbines reported in various investigations are examined and discussed.

17 WIND ENERGY↗

Planning Optimization for Inductively Charged On-demand Automated Electric Shuttles Project at Greenville, South Carolina

Wireless charging technology presents an ideal fit for autonomous electric vehicles for realizing a fully automated system (vehicle and charger). This paper presents a planning optimization analysis for a pilot project of in-route wireless charging infrastructure serving fixed-route on-demand shared automated electric shuttles (SAESs) at Greenville, South Carolina, USA. A single-objective non-linear mixed integer system planning optimization problem is formulated. A comprehensive cost model representing the overall inductively charged SAESs system is developed, considering road construction, power electronics and materials, traction battery, and installation costs. The optimization problem is solved to determine the best combination of the system key design parameters (number and allocations of wireless chargers, charging power level, track length and on-board battery capacity) that show the most cost-effective solution and allow the SAESs to achieve charge-sustaining operation. The planning platform incorporates representative simulated traffic data (driving speed and routes) for four SAESs in the Greenville using the Simulation of Urban Mobility tool. These data are fed to a vehicle powertrain model and a wireless charger power model to predict the battery power, energy and state-of-charge profiles, which are provided to the search algorithm to assess the design objectives under specific constraints. Finally, the results indicate that implementing high-power (100 kW) wireless charger at a few designated stops for fixed-route SAESs with the proper track length allows the vehicles to realize charge-sustaining operation, infinite range, and zero recharge downtime, with a significant reduction in the on-board battery (36%) and road coverage (69%), at minimum cost.

33 ADVANCED PROPULSION SYSTEMS↗

High-resolution microwave images of Saturn

An analysis of high-resolution microwave images of Saturn and Saturn's individual rings is presented. Radio interferometric observations of Saturn taken at the Very Large Array in New Mexico at wavelengths of 2 and 6 centimeters reveal interesting new features in both the atmosphere and rings. The resulting maps show an increase in brightness temperature of about 3 K from equator to pole at both wavelengths, while the 6-centimeter map shows a bright band at northern midlatitudes. The data are consistent with a radiative transfer model of the atmosphere that constrains the well-mixed, fully saturated, NH3 mixing ratio to be 0.00012 in a region just below the NH3 clouds, while the observed bright band indicates a 25 percent relative decrease of NH3 in northern midlatitudes. Brightness temperatures for the classical rings are presented. Ring brightness shows a variation with azimuth and is linearly polarized at an average value of about 5 percent. The variations in ring polarization suggest that at least 20 percent of the ring brightness is the result of a single scattering process.

Grossman, A. W.↗

Rheology of Foam Near the Order-Disorder Transition

The first part of our research results are summarized in the recent journal publication: J. Gregory McDaniel and R. Glynn Holt, 'Measurement of aqueous foam rheology by acoustic levitation', Phys. Rev. E 61, 2204 (2000). This aspect of the work was a combination of experiment and analysis. We built a levitation system capable of acoustically levitating small samples of aqueous foam of arbitrary gas and liquid volume fractions. We then modulated the acoustic field to induce normal mode oscillations of the foam samples. The observables from the experiment were frequency and mode number. For dry (roughly > 70% gas by volume) foams and small deformations, we developed an effective medium, normal-modes analysis which took the frequency and mode number from experiment, and gave us the shear elastic modulus of the foam as a function of Poisson's ratio. The second part of our results may be found in a soon-to-be submitted manuscript 'Dynamics of aqueous foam drops', I.Sh. Akhatov, J.G. McDaniel and R.G. Holt, describing our modeling in the wet foam limit by considering the acoustic problem. This aspect of the research is purely theoretical. Beginning from a mass-conserving mixture law, the fully nonlinear equations of motion for a wet (roughly < 10% gas by volume) foam drop of initially spherical shape were derived. The frequencies for normal mode oscillations were derived in the linear inviscid limit. The nonlinear equations were numerically solved to elicit the motion of a foam drop under acoustic excitation. The role of the time-varying void fraction in breathing-mode oscillations is of particular interest. As of the end of the current (NAG#3-2121) grant, this work was not yet concluded. We continue to work on this aspect in order to extend the analysis to cover the transition regime of gas volume fractions, as well as to compare to experiments in the wet regime.

Holt, R. Glynn↗

Transition Analysis for the CRM-NLF Wind Tunnel Configuration

This paper reports the results of a comprehensive linear stability analysis of the boundary layer flow over the common research model with natural laminar flow (CRM-NLF) aircraft configuration. The flow conditions match selected test conditions from a recent wind tunnel experiment in the National Transonic Facility at the NASA Langley Research Center. Previous work has shown that the measured onset of laminar-turbulent transition during the experiments can be correlated with the linear amplification of Tollmien-Schlichting (TS) and stationary crossflow (CF) instabilities in the swept wing boundary layer. However, a significant scatter ( N ∈ (4,9)) was observed in the values of the logarithmic amplification factor along the measured transition front. This previous analysis was based on an approximate basic state (based on a boundary layer code with conical flow approximation) and parallel stability computations without surface curvature effects. Here, we examine the effects of these various approximations with the goal of quantifying the resulting changes in the N-factor correlations. Specifically, both linear stability theory (LST) and the parabolized stability equations (PSE)are used in conjunction with an accurate definition of the laminar boundary layer flow as computed with a Navier-Stokes solver with a Reynolds-Averaged-Navier-Stokes (RANS) based turbulence model within the turbulent parts of the flow. Furthermore, the effects of instability wave propagation within a fully three-dimensional boundary layer are also evaluated by integrating the disturbance growth rates along suitably chosen, curvilinear (i.e., nonplanar) propagation trajectories. The results of this analysis are also used in an accompanying paper by Venkatachari et al. to develop improved, physics based transition predictions for the same CRM-NLF configuration.

Boundary layer transition↗

A positive and entropy-satisfying finite volume scheme for the Baer–Nunziato model

We present a relaxation scheme for approximating the entropy dissipating weak solutions of the Baer–Nunziato two-phase flow model. This relaxation scheme is straightforwardly obtained as an extension of the relaxation scheme designed in for the isentropic Baer–Nunziato model and consequently inherits its main properties. To our knowledge, this is the only existing scheme for which the approximated phase fractions, phase densities and phase internal energies are proven to remain positive without any restrictive condition other than a classical fully computable CFL condition. For ideal gas and stiffened gas equations of state, real values of the phasic speeds of sound are also proven to be maintained by the numerical scheme. It is also the only scheme for which a discrete entropy inequality is proven, under a CFL condition derived from the natural sub-characteristic condition associated with the relaxation approximation. This last property, which ensures the non-linear stability of the numerical method, is satisfied for any admissible equation of state. We provide a numerical study for the convergence of the approximate solutions towards some exact Riemann solutions. The numerical simulations show that the relaxation scheme compares well with two of the most popular existing schemes available for the Baer–Nunziato model, namely Schwendeman–Wahle–Kapila's Godunov-type scheme and Tokareva–Toro's HLLC scheme . The relaxation scheme also shows a higher precision and a lower computational cost (for comparable accuracy) than a standard numerical scheme used in the nuclear industry, namely Rusanov's scheme. Finally, we assess the good behavior of the scheme when approximating vanishing phase solutions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A fully implicit, scalable, conservative nonlinear relativistic Fokker–Planck 0D-2P solver for runaway electrons

Upon application of a sufficiently strong electric field, electrons break away from thermal equilibrium and approach relativistic speeds. These highly energetic ‘runaway’ electrons (~ MeV) play a significant role in tokamak disruption physics, and therefore their accurate understanding is essential to develop reliable mitigation strategies. As such, we have developed a fully implicit solver for the 0D-2P (i.e., including two momenta coordinates) relativistic nonlinear Fokker–Planck equation (rFP). As in earlier implicit rFP studies (NORSE, CQL3D), electron–ion interactions are modeled using the Lorentz operator, and synchrotron damping using the Abraham–Lorentz–Dirac reaction term. However, our implementation improves on these earlier studies by (1) ensuring exact conservation properties for electron collisions, (2) strictly preserving positivity, and (3) being scalable algorithmically and in parallel. Key to our proposed approach is an efficient multigrid preconditioner for the linearized rFP equation, a multigrid elliptic solver for the Braams–Karney potentials, and a novel adaptive technique to determine the associated boundary values. We verify the accuracy and efficiency of the proposed scheme with numerical results ranging from small electric-field electrical conductivity measurements to the accurate reproduction of runaway tail dynamics when strong electric fields are applied.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Simulation and Regression Modeling of X-59 Low-Boom Carpets Across America

The NASA X-59 aircraft is predicted to produce a significantly quieter cruise sonic boom than traditional N-wave-producing aircraft. A propagation simulation study was undertaken to quantify loudness levels, exposure size, and variability of the X-59 low-boom carpet using realistic atmospheric profiles across the contiguous United States of America (CONUS). Near-field pressure data of the X-59 in supersonic cruise from NASA’s fully unstructured Navier–Stokes three-dimensional (known as FUN3D) computational fluid dynamics code were propagated using NASA’s PCBoom code, which solves an enhanced Burgers equation along acoustic rays. Atmospheric profiles from the National Oceanic and Atmospheric Administration’s Climate Forecast System Version 2 database were used for propagation at 138 locations across the CONUS. Carpets at each location were generated for aircraft headings in the four cardinal directions. Over one million X-59 carpets were generated in total. The effects of the heading, season, geography, and climate zone on boom levels and exposure size are presented. Multiple linear regression models were developed to estimate carpet width and loudness metrics across the CONUS. These results inform regulators and mission planners on expected variations in boom levels and carpet extent from atmospheric variations. Understanding potential carpet variability is important when planning community noise surveys using the X-59.

X-59↗

Mass and energy transfer across the Earth's magnetopause caused by vortex-induced reconnection: Mass and energy transfer by K-H vortex

When the interplanetary magnetic field (IMF) is strongly northward, a boundary layer that contains a considerable amount of plasma of magnetosheath origin is often observed along and earthward of the low-latitude magnetopause. Such a pre-existing boundary layer, with a higher density than observed in the adjacent magnetosphere, reduces the local Alfvén speed and allows the Kelvin-Helmholtz instability (KHI) to grow more strongly. We employ a three-dimensional fully kinetic simulation to model an event observed by the Magnetospheric Multiscale (MMS) mission in which the spacecraft detected substantial KH waves between a pre-existing boundary layer and the magnetosheath during strong northward IMF. Initial results of this simulation [Nakamura et al., 2017] have successfully demonstrated ion-scale signatures of magnetic reconnection induced by the non-linearly developed KH vortex, which are quantitatively consistent with MMS observations. Furthermore, we quantify the simulated mass and energy transfer processes driven by this vortex-induced reconnection (VIR) and show that during this particular MMS event (i) mass enters a new mixing layer formed by the VIR more efficiently from the pre-existing boundary layer side than from the magnetosheath side, (ii) mixed plasmas within the new mixing layer convect tailward along the magnetopause at more than half the magnetosheath flow speed, and (iii) energy dissipation in localized VIR dissipation regions results in a strong parallel electron heating within the mixing layer. Finally, the quantitative agreements between the simulation and MMS observations allow new predictions that elucidate how the mass and energy transfer processes occur near the magnetopause during strong northward IMF.

58 GEOSCIENCES↗