Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical continuation”

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 649 records · Page 36

Numerical simulations of trailing vortex bursting

Solutions of the steady-state Navier-Stokes equations for the axisymmetric bursting of a laminar trailing vortex are computed with Newton's method and the pseudo-arc length continuation method for wide ranges of vortex strength and Reynolds number. The results indicate that a trailing vortex can undergo a transition from a state in which the core slowly diffuses to a state marked by large amplitude, spatial oscillations of core radius and core axial velocity. At the transition point the core grows rapidly in size. This event is interpreted as vortex bursting. The results also suggest that when the maximum core swirl velocity is sufficiently large the centerline axial flow downstream of transition will be reversed.

Beran, Philip S.↗

Numerical studies of motion of vortex filaments - Implementing the asymptotic analysis

A computational code is developed for the integro-differential equations governing the motion of the centerlines of vortex filaments submerged in a background potential flow. These equations, which are derived from the method of matched asymptotic analysis, include the effect of the decaying large-magnitude circumferential and axial velocity components in the vortical cores. Numerical examples are presented to assess the effect of a large axial velocity and that of nonsimilar initial profiles in the vortical cores. The initial configurations of the filaments are chosen so as to fulfill the basic assumption of the asymptotic analysis, which is that the effective vortical core size is much smaller than all the other length scales in the flowfield, e.g., the radius of curvature and the interfilament distance. The computations are continued until the basic assumption is no longer valid, that is when the merging or intersection of filaments has begun. A classification of the various types of local or global merging or intersection of filaments is made and demonstrated by numerical examples. It is then shown that the asymptotic solution not only provides the initial data but also can be used to formulate the appropriate boundary conditions for the numerical solution of a merged region.

Liu, C. H.↗

A computer code for multiphase all-speed transient flows in complex geometries. MAST version 1.0

The operation of the MAST code, which computes transient solutions to the multiphase flow equations applicable to all-speed flows, is described. Two-phase flows are formulated based on the Eulerian-Lagrange scheme in which the continuous phase is described by the Navier-Stokes equation (or Reynolds equations for turbulent flows). Dispersed phase is formulated by a Lagrangian tracking scheme. The numerical solution algorithms utilized for fluid flows is a newly developed pressure-implicit algorithm based on the operator-splitting technique in generalized nonorthogonal coordinates. This operator split allows separate operation on each of the variable fields to handle pressure-velocity coupling. The obtained pressure correction equation has the hyperbolic nature and is effective for Mach numbers ranging from the incompressible limit to supersonic flow regimes. The present code adopts a nonstaggered grid arrangement; thus, the velocity components and other dependent variables are collocated at the same grid. A sequence of benchmark-quality problems, including incompressible, subsonic, transonic, supersonic, gas-droplet two-phase flows, as well as spray-combustion problems, were performed to demonstrate the robustness and accuracy of the present code.

Chen, C. P.↗

Multigrid calculation of internal flows in complex geometries

The development, validation, and application of a general purpose multigrid solution algorithm and computer program for the computation of elliptic flows in complex geometries is presented. This computer program combines several desirable features including a curvilinear coordinate system, collocated arrangement of the variables, and Full Multi-Grid/Full Approximation Scheme (FMG/FAS). Provisions are made for the inclusion of embedded obstacles and baffles inside the flow domain. The momentum and continuity equations are solved in a decoupled manner and a pressure corrective equation is used to update the pressures such that the fluxes at the cell faces satisfy local mass continuity. Despite the computational overhead required in the restriction and prolongation phases of the multigrid cycling, the superior convergence results in reduced overall CPU time. The numerical scheme and selected results of several validation flows are presented. Finally, the procedure is applied to study the flowfield in a side-inlet dump combustor and twin jet impingement from a simulated aircraft fuselage.

Smith, K. M.↗

Diffusion Limited Supercritical Water Oxidation (SCWO) in Microgravity Environments

Tests designed to quantify the gravitational effects on thermal mixing and reactant injection in a Supercritical Water Oxidation (SCWO) reactor have recently been performed in the Zero Gravity Facility (ZGF) at NASA s Glenn Research Center. An artificial waste stream, comprising aqueous mixtures of methanol, was pressurized to approximately 250 atm and then heated to 450 C. After uniform temperatures in the reactor were verified, a controlled injection of air was initiated through a specially designed injector to simulate diffusion limited reactions typical in most continuous flow reactors. Results from a thermal mapping of the reaction zone in both 1-g and 0-g environments are compared. Additionally, results of a numerical model of the test configuration are presented to illustrate first order effects on reactant mixing and thermal transport in the absence of gravity.

Hicks, M. C.↗

Advanced Intensity-Modulation Continuous-Wave Lidar Techniques for ASCENDS O2 Column Measurements

Global atmospheric carbon dioxide (CO2) measurements for the NASA Active Sensing of CO2 Emissions over Nights, Days, and Seasons (ASCENDS) space mission are critical for improving our understanding of global CO2 sources and sinks. Advanced Intensity- Modulated Continuous-Wave (IM-CW) lidar techniques are investigated as a means of facilitating CO2 measurements from space to meet the ASCENDS measurement requirements. In recent numerical, laboratory and flight experiments we have successfully used the Binary Phase Shift Keying (BPSK) modulation technique to uniquely discriminate surface lidar returns from intermediate aerosol and cloud contamination. We demonstrate the utility of BPSK to eliminate sidelobes in the range profile as a means of making Integrated Path Differential Absorption (IPDA) column CO2 measurements in the presence of optically thin clouds, thereby eliminating the need to correct for sidelobe bias errors caused by the clouds. Furthermore, high accuracy and precision ranging to the surface as well as to the top of intermediate cloud layers, which is a requirement for the inversion of column CO2 number density measurements to column CO2 mixing ratios, has been demonstrated using new hyperfine interpolation techniques that takes advantage of the periodicity of the modulation waveforms. This approach works well for both BPSK and linear swept-frequency modulation techniques. The BPSK technique under investigation has excellent auto-correlation properties while possessing a finite bandwidth. A comparison of BPSK and linear swept-frequency is also discussed in this paper. These results are extended to include Richardson-Lucy deconvolution techniques to extend the resolution of the lidar beyond that implied by limit of the bandwidth of the modulation, where it is shown useful for making tree canopy measurements.

Campbell, Joel F.↗

initMIP-Antarctica: an Ice Sheet Model Initialization Experiment of ISMIP6

Ice sheet numerical modeling is an important tool to estimate the dynamic contribution of the Antarctic ice sheet to sea level rise over the coming centuries. The influence of initial conditions on ice sheet model simulations, however, is still unclear. To better understand this influence, an initial state intercomparison exercise (initMIP) has been developed to compare, evaluate, and improve initialization procedures and estimate their impact on century-scale simulations. initMIP is the first set of experiments of the Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6), which is the primary Coupled Model Intercomparison Project Phase 6 (CMIP6) activity focusing on the Greenland and Antarctic ice sheets. Following initMIP-Greenland, initMIP-Antarctica has been designed to explore uncertainties associated with model initialization and spin-up and to evaluate the impact of changes in external forcings. Starting from the state of the Antarctic ice sheet at the end of the initialization procedure, three forward experiments are each run for 100 years: a control run, a run with a surface mass balance anomaly, and a run with a basal melting anomaly beneath floating ice. This study presents the results of initMIP-Antarctica from 25 simulations performed by 16 international modeling groups. The submitted results use different initial conditions and initialization methods, as well as ice flow model parameters and reference external forcings. We find a good agreement among model responses to the surface mass balance anomaly but large variations in responses to the basal melting anomaly. These variations can be attributed to differences in the extent of ice shelves and their upstream tributaries, the numerical treatment of grounding line, and the initial ocean conditions applied, suggesting that ongoing efforts to better represent ice shelves in continental-scale models should continue.

ice sheet model↗

Ten Years of VIIRS On-Orbit Geolocation Calibration and Performance

The first innovative Visible Infrared Imaging Radiometer Suite (VIIRS) sensor aboard the Suomi National Polar-orbiting Partnership (SNPP) satellite has been in operation for 10 years since its launch on 28 October 2011. The second VIIRS sensor aboard the first Join Polar Satellite System (JPSS-1) satellite has been in operation for 4 years since its launch on 18 November 2017, which became NOAA-20. Well-geolocated and radiometrically calibrated Level-1 sensor data records (SDRs) from VIIRS are crucial to numerical weather prediction (NWP) and Level-2+ environmental data record (EDR) algorithms and products. The high quality of Level-2+ EDRs is a requirement for the continuity of NASA Earth science data records (ESDRs) and climate data records (CDRs), one of the two objectives of the SNPP mission and one of the three elements in the JPSS mission objective. The other objective of the SNPP mission is risk reduction for the follow-on JPSS missions. This paper summarizes the on-orbit geolocation calibration and validation (Cal/Val) activities for both VIIRS sensors onboard SNPP and NOAA-20 in the past 10 years. These activities include nominal geolocation Cal/Val activities, risk reduction activities, and improvements for the on-orbit VIIRS sensor operations. After these activities, sub-pixel geolocation accuracy is achieved. Nadir equivalent geolocation uncertainty is generally within 75 m (1-σ), or 20% imagery band pixels, in either the along-scan or along-track direction for both SNPP and NOAA-20 VIIRS sensors. The worst 16-day measured geolocation errors (radial, 3-σ) are 280 m and 267 m, respectively, in the latest SNPP and NOAA-20 VIIRS data collections, which are better than the required accuracy of 375 m (radial, 3-σ). The risk reduction activities also improved VIIRS builds for JPSS-3 and JPSS-4 satellites, and provide lessons learned for other VIIRS-like sensor builds.

VIIRS↗

Operational maize yield forecasts for Sub-Saharan Africa using Earth observation data and machine learning

Food insecurity continues to grow in Sub-Saharan Africa (SSA). In 2019, chronically malnourished people numbered nearly 240 million, or 20% of the population in SSA. Globally, numerous efforts have been made to anticipate potential droughts, crop conditions, and food shortages to foster improved food insecurity early warning and risk management. To support this goal, we develop an Earth Observation (EO) and machine-learning-based operational, subnational maize yield forecast system and evaluate its out-of-sample forecast skills during the growing seasons for Kenya, Somalia, Malawi, and Burkina Faso. In general, forecast skills improve substantially during the vegetative growth period (VP) and gradually during the reproductive development period (RP). Thus, mid-season assessment can provide effective early warning months before harvest. Skillful forecasts (Nash Sutcliffe Efficiency (NSE) > 0.6 and Mean Absolute Percentage Error (MAPE) < 20%) appear approximately two dekads after the VP; for example, they appear in May in Kenya and Somalia, January in Malawi, and July in Burkina Faso.

Donghoon Lee↗

A numerical method for the solution of three dimensional, incompressible, viscous flows over slender bodies

A marching iterative method for the solution of the three dimensional, incompressibhle, steady and parabolized Navier-Stokes equations is described. The equations are written in primitive variables and discretized in general axisymmetric orthogonal coordinate systems. The coupled set of finite-difference equations are solved without any splitting or factorization errors. Moreover, the continuity equation and the two crossflow momentum equations are exactly satisfied at every step of the iterative process. The solution scheme is equivalent to the solution of one Poisson equation by the Successive Plane Over Relaxation method and has good convergence properties. Other existing solution methods resemble a Jacobi-type iterative scheme and therefore are less efficient. Numerical experiments include the laminar, incompressible flow over prolate spheroids at incidence.

Rosenfeld, Moshe↗

Climate Change, Extreme Weather Events, and Fungal Disease Emergence and Spread

Empirical evidence from multiple sources show the Earth has been warming since the late 19th century. More recently, evidence for this warming trend is strongly supported by satellite data since the late 1970s from the cryosphere, atmosphere, oceans, and land that confirms increasing temperature trends and their consequences (e.g., reduced Arctic sea ice, rising sea level, ice sheet mass loss, etc.). At the same time, satellite observations of the Sun show remarkably stable solar cycles since the late 1970s, when direct observations of the Sun's total solar irradiance began. Numerical simulation models, driven in part by assimilated satellite data, suggest that future-warming trends will lead to not only a warmer planet, but also a wetter and drier climate depending upon location in a fashion consistent with large-scale atmospheric processes. Continued global warming poses new opportunities for the emergence and spread of fungal disease, as climate systems change at regional and global scales, and as animal and plant species move into new niches. Our contribution to this proceedings is organized thus: First, we review empirical evidence for a warming Earth. Second, we show the Sun is not responsible for the observed warming. Third, we review numerical simulation modeling results that project these trends into the future, describing the projected abiotic environment of our planet in the next 40 to 50 years. Fourth, we illustrate how Rift Valley fever outbreaks have been linked to climate, enabling a better understanding of the dynamics of these diseases, and how this has led to the development of an operational predictive outbreak model for this disease in Africa. Fifth, We project how this experience may be applicable to predicting outbreaks of fungal pathogens in a warming world. Lastly, we describe an example of changing species ranges due to climate change, resulting from recent warming in the Andes and associated glacier melt that has enabled amphibians to colonize higher elevation lakes, only to be followed shortly by the emergence of fungal disease in the new habitats.

Tucker, Compton J.↗

Electromagnetic Field Penetration Studies

A numerical method is presented to determine electromagnetic shielding effectiveness of rectangular enclosure with apertures on its wall used for input and output connections, control panels, visual-access windows, ventilation panels, etc. Expressing EM fields in terms of cavity Green's function inside the enclosure and the free space Green's function outside the enclosure, integral equations with aperture tangential electric fields as unknown variables are obtained by enforcing the continuity of tangential electric and magnetic fields across the apertures. Using the Method of Moments, the integral equations are solved for unknown aperture fields. From these aperture fields, the EM field inside a rectangular enclosure due to external electromagnetic sources are determined. Numerical results on electric field shielding of a rectangular cavity with a thin rectangular slot obtained using the present method are compared with the results obtained using simple transmission line technique for code validation. The present technique is applied to determine field penetration inside a Boeing-757 by approximating its passenger cabin as a rectangular cavity filled with a homogeneous medium and its passenger windows by rectangular apertures. Preliminary results for, two windows, one on each side of fuselage were considered. Numerical results for Boeing-757 at frequencies 26 MHz, 171-175 MHz, and 428-432 MHz are presented.

Deshpande, M.D.↗

Numerical Evaluation of the "Dual-Kernel Counter-flow" Matric Convolution Integral that Arises in Discrete/Continuous (D/C) Control Theory

Discrete/Continuous (D/C) control theory is a new generalized theory of discrete-time control that expands the concept of conventional (exact) discrete-time control to create a framework for design and implementation of discretetime control systems that include a continuous-time command function generator so that actuator commands need not be constant between control decisions, but can be more generally defined and implemented as functions that vary with time across sample period. Because the plant/control system construct contains two linear subsystems arranged in tandem, a novel dual-kernel counter-flow convolution integral appears in the formulation. As part of the D/C system design and implementation process, numerical evaluation of that integral over the sample period is required. Three fundamentally different evaluation methods and associated algorithms are derived for the constant-coefficient case. Numerical results are matched against three available examples that have closed-form solutions.

Nixon, Douglas D.↗

A block iterative finite element algorithm for numerical solution of the steady-state, compressible Navier-Stokes equations

An iterative method for numerically solving the time independent Navier-Stokes equations for viscous compressible flows is presented. The method is based upon partial application of the Gauss-Seidel principle in block form to the systems of the nonlinear algebraic equations which arise in construction of finite element (Galerkin) models approximating solutions of fluid dynamic problems. The continuous cubic element on triangles is employed for function approximation. Computational results for a free shear flow at Re = 1000 indicate significant achievement of economy in iterative convergence rate over finite element and finite difference models which employ the customary time dependent equations and symptotic time marching procedure to steady solution. Numerical results are in excellent agreement with those obtained for the same test problem employing time marching finite element and finite difference solution techniques.

Cooke, C. H.↗

Comparison of four stable numerical methods for Abel's integral equation

The 3-D image reconstruction from cone-beam projections in computerized tomography leads naturally, in the case of radial symmetry, to the study of Abel-type integral equations. If the experimental information is obtained from measured data, on a discrete set of points, special methods are needed in order to restore continuity with respect to the data. A new combined Regularized-Adjoint-Conjugate Gradient algorithm, together with two different implementations of the Mollification Method (one based on a data filtering technique and the other on the mollification of the kernal function) and a regularization by truncation method (initially proposed for 2-D ray sample schemes and more recently extended to 3-D cone-beam image reconstruction) are extensively tested and compared for accuracy and numerical stability as functions of the level of noise in the data.

Murio, Diego A.↗

A method for calculating unsteady turbulent boundary layers in two- and three-dimensional flows.

The governing equations for an unsteady turbulent boundary layer on a swept infinite cylinder, composed of a continuity equation, a pair of momentum equations and a pair of rate equations for the shear stress, based on the turbulent kinetic-energy equation, are solved numerically. These rate equations for the shear-stress express a balance between the convection, production, dissipation and diffusion of the turbulent shear stress components. Due to the physical model of turbulence employed, this system of equations is hyperbolic. An explicit, second-order accurate, conditionally stable difference scheme for this system, together with appropriate boundary and initial conditions, is formulated and developed into a computer program. Calculations for oscillating free stream flows with no pressure gradient show significant unsteady effects on the turbulent boundary layer.

Singleton, R. E.↗

Status and projections of the NAS program

NASA's Numerical Aerodynamic Simulation (NAS) Program has completed development of the initial operating configuration of the NAS Processing System Network (NPSN). This is the first milestone in the continuing and pathfinding effort to provide state-of-the-art supercomputing for aeronautics research and development. The NPSN, available to a nation-wide community of remote users, provides a uniform UNIX environment over a network of host computers ranging from the Cray-2 supercomputer to advanced scientific workstations. This system, coupled with a vendor-independent base of common user interface and network software, presents a new paradigm for supercomputing environments. Background leading to the NAS program, its programmatic goals and strategies, technical goals and objectives, and the development activities leading to the current NPSN configuration are presented. Program status, near-term plans, and plans for the next major milestone, the extended operating configuration, are also discussed.

Bailey, Frank R.↗

Status and projections of the NAS Program

NASA's Numerical Aerodynamic Simulation (NAS) Program has completed development of the initial operating configuration of the NAS Processing System Network (NPSN). This is the first milestone in the continuing and pathfinding effort to provide state-of-the-art supercomputing for aeronautics research and development. The NPSN, available to a nation-wide community of remote users, provides a uniform UNIX environment over a network of host computers ranging from the new Cray-2 supercomputer to advanced scientific workstations. This system, coupled with a vendor-independent base of common user interface and network software, presents a new paradigm for supercomputing environments. Presented here is the background leading to the NAS Program, its programmatic goals and strategies, technical goals and objectives, and the development activities leading to the current NPSN configuration. Program status, near-term plans and plans for the next major milestone, the extended operating configuration, are also discussed.

Bailey, F. R.↗