Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Convergence analysis”

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,261 records · Page 70

Development of a method of analysis and computer program for calculating the inviscid flow about the windward surfaces of space shuttle configurations at large angles of attack

A general method developed for the analysis of inviscid hypersonic shock layers is discussed for application to the case of the shuttle vehicle at high (65 deg) angle of attack. The associated extensive subsonic flow region caused convergence difficulties whose resolution is discussed. It is required that the solution be smoother than anticipated.

Maslen, S. H.↗

Interaction of a two-dimensional strip boundary layer with a three-dimensional transonic swept-wing code

A 3D inviscid transonic analysis code was combined with a 2D strip integral boundary layer technique to form an approximate interaction procedure for analyzing the flow over a high aspect ratio wing near cruise conditions. Converged results obtained using the procedure for an aspect ratio 10.3 supercritical wing are discussed. Angle of attack adjustments were made during the iterative procedure in order to compensate for the viscous lift loss. A comparison of the calculations with experimental data is presented.

Newman, P. A.↗

Exploiting parallel computing with limited program changes using a network of microcomputers

Network computing and multiprocessor computers are two discernible trends in parallel processing. The computational behavior of an iterative distributed process in which some subtasks are completed later than others because of an imbalance in computational requirements is of significant interest. The effects of asynchronus processing was studied. A small existing program was converted to perform finite element analysis by distributing substructure analysis over a network of four Apple IIe microcomputers connected to a shared disk, simulating a parallel computer. The substructure analysis uses an iterative, fully stressed, structural resizing procedure. A framework of beams divided into three substructures is used as the finite element model. The effects of asynchronous processing on the convergence of the design variables are determined by not resizing particular substructures on various iterations.

Rogers, J. L., Jr.↗

North American plate dynamics

Deformation within the North American plate in response to various tectonic processes is modeled using an elastic finite element analysis. The tectonic processes considered in the modeling include ridge forces associated with the normal thermal evolution of oceanic lithosphere, shear and normal stresses transmitted across transforms, normal stresses transmitted across convergent boundaries, stresses due to horizontal density contrasts within the continent, and shear tractions applied along the base of the plate. Model stresses are calculated with respect to a lithostatic reference stress state. Shear stresses transmitted across transform boundaries along the San Andreas and Caribbean are small, of the order of 5-10 MPa. Also, compressive stresses of the order of 5-10 MPa transmitted across the major transforms improve the fit to the data. Compressive stresses across convergent margins along the Aleutians and the Middle America trench are important.

Richardson, Randall M.↗

Consistency of Estimated Global Water Cycle Variations Over the Satellite Era

Motivated by the question of whether recent indications of decadal climate variability and a possible "climate shift" may have affected the global water balance, we examine evaporation minus precipitation (E-P) variability integrated over the global oceans and global land from three points of view-remotely sensed retrievals / objective analyses over the oceans, reanalysis vertically-integrated moisture convergence (MFC) over land, and land surface models forced with observations-based precipitation, radiation and near-surface meteorology. Because monthly variations in area-averaged atmospheric moisture storage are small and the global integral of moisture convergence must approach zero, area-integrated E-P over ocean should essentially equal precipitation minus evapotranspiration (P-ET) over land (after adjusting for ocean and land areas). Our analysis reveals considerable uncertainty in the decadal variations of ocean evaporation when integrated to global scales. This is due to differences among datasets in 10m wind speed and near-surface atmospheric specific humidity (2m qa) used in bulk aerodynamic retrievals. Precipitation variations, all relying substantially on passive microwave retrievals over ocean, still have uncertainties in decadal variability, but not to the degree present with ocean evaporation estimates. Reanalysis MFC and P-ET over land from several observationally forced diagnostic and land surface models agree best on interannual variations. However, upward MFC (i.e. P-ET) reanalysis trends are likely related in part to observing system changes affecting atmospheric assimilation models. While some evidence for a low-frequency E-P maximum near 2000 is found, consistent with a recent apparent pause in sea-surface temperature (SST) rise, uncertainties in the datasets used here remain significant. Prospects for further reducing uncertainties are discussed. The results are interpreted in the context of recent climate variability (Pacific Decadal Oscillation, Atlantic Meridional Overturning), and efforts to distinguish these modes from longer-term trends.

Robertson, F. R.↗

Consistency of Estimated Global Water Cycle Variations Over the Satellite Era

Motivated by the question of whether recent indications of decadal climate variability and a possible "climate shift" may have affected the global water balance, we examine evaporation minus precipitation (E-P) variability integrated over the global oceans and global land from three points of view-remotely sensed retrievals / objective analyses over the oceans, reanalysis vertically-integrated moisture convergence (MFC) over land, and land surface models forced with observations-based precipitation, radiation and near-surface meteorology. Because monthly variations in area-averaged atmospheric moisture storage are small and the global integral of moisture convergence must approach zero, area-integrated E-P over ocean should essentially equal precipitation minus evapotranspiration (P-ET) over land (after adjusting for ocean and land areas). Our analysis reveals considerable uncertainty in the decadal variations of ocean evaporation when integrated to global scales. This is due to differences among datasets in 10m wind speed and near-surface atmospheric specific humidity (2m qa) used in bulk aerodynamic retrievals. Precipitation variations, all relying substantially on passive microwave retrievals over ocean, still have uncertainties in decadal variability, but not to the degree present with ocean evaporation estimates. Reanalysis MFC and P-ET over land from several observationally forced diagnostic and land surface models agree best on interannual variations. However, upward MFC (i.e. P-ET) reanalysis trends are likely related in part to observing system changes affecting atmospheric assimilation models. While some evidence for a low-frequency E-P maximum near 2000 is found, consistent with a recent apparent pause in sea-surface temperature (SST) rise, uncertainties in the datasets used here remain significant. Prospects for further reducing uncertainties are discussed. The results are interpreted in the context of recent climate variability (Pacific Decadal Oscillation, Atlantic Meridional Overturning), and efforts to distinguish these modes from longer-term trends.

Robertson, F. R.↗

Colour-producing [beta]-keratin nanofibres in blue penguin (Eudyptula minor) feathers

The colours of living organisms are produced by the differential absorption of light by pigments (e.g. carotenoids, melanins) and/or by the physical interactions of light with biological nanostructures, referred to as structural colours. Only two fundamental morphologies of non-iridescent nanostructures are known in feathers, and recent work has proposed that they self-assemble by intracellular phase separation processes. Here, we report a new biophotonic nanostructure in the non-iridescent blue feather barbs of blue penguins (Eudyptula minor) composed of parallel {beta}-keratin nanofibres organized into densely packed bundles. Synchrotron small angle X-ray scattering and two-dimensional Fourier analysis of electron micrographs of the barb nanostructure revealed short-range order in the organization of fibres at the appropriate size scale needed to produce the observed colour by coherent scattering. These two-dimensional quasi-ordered penguin nanostructures are convergent with similar arrays of parallel collagen fibres in avian and mammalian skin, but constitute a novel morphology for feathers. The identification of a new class of {beta}-keratin nanostructures adds significantly to the known mechanisms of colour production in birds and suggests additional complexity in their self-assembly.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Asymptotic solutions and convergence studies of the resistive inner region equations

A recent publication presents a detailed prescription for constructing global eigenfunctions and complex growth rates for resistive instabilities in axisymmetric toroidal systems, using the method of matched asymptotic expansions. At each singular surface in the plasma, the small-x asymptotic solutions of the ideal MHD outer-region equations are matched to the large-x asymptotic solutions of the resistive inner-region equations, with x being the distance from the singular surface. The success of this method depends upon accurate asymptotic solutions of the resistive inner region equations. Extensive studies have shown that there are not only regimes of excellent agreement with straight-through methods, without asymptotic matching, but also regimes of strong disagreement. As such, in an effort to find the cause of this discrepancy, we present here a different method of constructing asymptotic solutions of the inner region equations and we study the convergence properties of the new and old solutions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

INVERSE CASCADE IN IMBALANCED ELECTRON MAGNETOHYDRODYNAMIC TURBULENCE

Electron magnetohydrodynamics (EMHD) provides a fluid-like description of small-scale magnetized plasmas. Balanced EMHD turbulence has been studied for a long time. However, driven imbalanced EMHD turbulence, in which waves moving in one direction (dominant waves) have higher amplitudes than waves moving in the other direction (sub-dominant waves), has not been well studied. In this paper, we numerically study driven three-dimensional imbalanced weak EMHD turbulence. We find the following results. First, in driven imbalanced EMHD turbulence, we clearly observe inverse cascade of magnetic helicity, as well as magnetic energy. This is because magnetic helicity is a conserved quantity and non-zero magnetic helicity is injected into the system in driven imbalanced EMHD turbulence. Second, the magnetic energy spectrum of the dominant waves on scales larger than the energy injection scale does not show a single power-law spectrum, which indicates that the inverse cascade is not a self-similar process. The peak of the spectrum roughly follows a k {sup –3/2} spectrum, which can be explained by a Kolmogorov-type argument for weak turbulence. Third, a small amount of sub-dominant waves is induced by the dominant waves on large scales and the ratio of helicity densities of the dominant and the sub-dominant waves on large scales seems to converge to a certain value.

79 ASTRONOMY AND ASTROPHYSICS↗

Numerical analysis of Weyl's method for integrating boundary layer equations

A fast method for accurate numerical integration of Blasius equation is proposed. It is based on the limit interchange in Weyl's fixed point method formulated as an iterated limit process. Each inner limit represents convergence to a discrete solution. It is shown that the error in a discrete solution admits asymptotic expansion in even powers of step size. An extrapolation process is set up to operate on a sequence of discrete solutions to reach the outer limit. Finally, this method is extended to related boundary layer equations.

Najfeld, I.↗

Improved Euler analysis of advanced turboprop propeller flows

An implicit approximate factorization scheme in conjunction with a new boundary treatment was used to compute the inviscid flow field about an advanced high-speed propeller. The method of characteristics was used to apply impermeable boundary conditions. The convergence history of numerical calculations shows substantial decrease in the residual error decay. Euler solutions were computed for SR-3 propeller geometry and the results were compared with the experimental data and previous numerical results.

Yamamoto, O.↗

On approaches to analyze the sensitivity of simulated hydrologic fluxes to model parameters in the community land model

Here, effective sensitivity analysis approaches are needed to identify important parameters or factors and their uncertainties in complex Earth system models composed of multi-phase multi-component phenomena and multiple biogeophysical-biogeochemical processes. In this study, the impacts of 10 hydrologic parameters in the Community Land Model on simulations of runoff and latent heat flux are evaluated using data from a watershed. Different metrics, including residual statistics, the Nash-Sutcliffe coefficient, and log mean square error, are used as alternative measures of the deviations between the simulated and field observed values. Four sensitivity analysis (SA) approaches, including analysis of variance based on the generalized linear model, generalized cross validation based on the multivariate adaptive regression splines model, standardized regression coefficients based on a linear regression model, and analysis of variance based on support vector machine, are investigated. Results suggest that these approaches show consistent measurement of the impacts of major hydrologic parameters on response variables, but with differences in the relative contributions, particularly for the secondary parameters. The convergence behaviors of the SA with respect to the number of sampling points are also examined with different combinations of input parameter sets and output response variables and their alternative metrics. This study helps identify the optimal SA approach, provides guidance for the calibration of the Community Land Model parameters to improve the model simulations of land surface fluxes, and approximates the magnitudes to be adjusted in the parameter values during parametric model optimization.

58 GEOSCIENCES↗

Cusps in the center of galaxies: a real conflict with observations or a numerical artefact of cosmological simulations?

Galaxy observations and N-body cosmological simulations produce conflicting dark matter halo density profiles for galaxy central regions. While simulations suggest a cuspy and universal density profile (UDP) of this region, the majority of observations favor variable profiles with a core in the center. In this paper, we investigate the convergency of standard N-body simulations, especially in the cusp region, following the approach proposed by [1]. We simulate the well known Hernquist model using the SPH code Gadget-3 and consider the full array of dynamical parameters of the particles. We find that, although the cuspy profile is stable, all integrals of motion characterizing individual particles suffer strong unphysical variations along the whole halo, revealing an effective interaction between the test bodies. This result casts doubts on the reliability of the velocity distribution function obtained in the simulations. Moreover, we find unphysical Fokker-Planck streams of particles in the cusp region. The same streams should appear in cosmological N-body simulations, being strong enough to change the shape of the cusp or even to create it. Our analysis, based on the Hernquist model and the standard SPH code, strongly suggests that the UDPs generally found by the cosmological N-body simulations may be a consequence of numerical effects. A much better understanding of the N-body simulation convergency is necessary before a 'core-cusp problem' can properly be used to question the validity of the CDM model.

79 ASTRONOMY AND ASTROPHYSICS↗

Numerical analysis of air dehumidification through electrospray-enhanced vortical flow cyclone separator

The dehumidification of air is a highly energy-intensive process, typically requiring vapor compression systems to cool the air significantly below its dew point temperature. Here, this study numerically investigates air dehumidification achieved by electrostatic and dielectrophoresis forces in a converging duct and centrifugal principles in cyclone separators. An electrospray was employed to energize water droplets, which then captured water vapor via dielectrophoresis. A novelty of the present study is the inclusion of a cyclone separator positioned downstream of the electrospray. This combination increased the dehumidification efficiency and controlled the flow of the condensate. By integrating electrospray technology and applying an electric field, this study examines the impact of particle separation efficiency, pressure drops, and cyclone performance on air dehumidification. The control parameters for air dehumidification were investigated. Three cyclones were modelled for airflow rates of 5 cubic feet per minute [cfm] (0.0024 m 3 /s), 38 cfm (0.0179 m 3 /s), and 200 cfm (0.094 m 3 /s). The results indicated air dehumidification was up to 5.7 % under environmental conditions, similar to those used to validate the CFD model. The degree of air dehumidification was primarily influenced by the electric charge acquired by the water droplets, their flight trajectory, and the vorticity intensity within the cyclone. Air dehumidification could increase by up to 13 % by increasing the number of droplets injected into the air stream, the flight time, and the electrical charge. The breakdown electric field strength threshold near the electrospray was a limiting factor for further enhancing dehumidification performance. This study examines the optimization of an electro-assisted dehumidification system that includes cyclone-based droplet separation through numerical analysis. The results provide insights into the design of energy-efficient air treatment technologies.

42 ENGINEERING↗

Efficient modeling of photonic crystals with local Hermite polynomials

Developing compact algorithms for accurate electrodynamic calculations with minimal computational cost is an active area of research given the increasing complexity in the design of electromagnetic composite structures such as photonic crystals, metamaterials, optical interconnects, and on-chip routing. We show that electric and magnetic (EM) fields can be calculated using scalar Hermite interpolation polynomials as the numerical basis functions without having to invoke edge-based vector finite elements to suppress spurious solutions or to satisfy boundary conditions. This approach offers several fundamental advantages as evidenced through band structure solutions for periodic systems and through waveguide analysis. Compared with reciprocal space (plane wave expansion) methods for periodic systems, advantages are shown in computational costs, the ability to capture spatial complexity in the dielectric distributions, the demonstration of numerical convergence with scaling, and variational eigenfunctions free of numerical artifacts that arise from mixed-order real space basis sets or the inherent aberrations from transforming reciprocal space solutions of finite expansions. The photonic band structure of a simple crystal is used as a benchmark comparison and the ability to capture the effects of spatially complex dielectric distributions is treated using a complex pattern with highly irregular features that would stress spatial transform limits. This general method is applicable to a broad class of physical systems, e.g., to semiconducting lasers which require simultaneous modeling of transitions in quantum wells or dots together with EM cavity calculations, to modeling plasmonic structures in the presence of EM field emissions, and to on-chip propagation within monolithic integrated circuits.

36 MATERIALS SCIENCE↗

Marangoni mixed convection flow with Joule heating and nonlinear radiation

Marangoni mixed convective flow of Casson fluid in a thermally stratified medium is addressed. Flow analysis has been carried out in presence of inclined magnetic field. Heat transfer analysis is discussed in the presence of viscous dissipation, Joule heating and nonlinear thermal radiation. The governing nonlinear partial differential equations are first converted into ordinary differential systems and then developed the convergent series solutions. Flow pattern with the influence of pertinent parameters namely the magnetic parameter, Casson fluid parameter, temperature ratio parameter, stratification parameter, Prandtl number, Eckert number and radiation parameter is investigated. Expression of local Nusselt number is computed and analyzed. It is found that the Nusselt number decreases by increasing magnetic parameter, temperature ratio parameter, angle of inclination and stratification parameter. Moreover the effect of buoyancy parameter on the velocity distribution is opposite in both the opposing and assisting flow phenomena. Thermal field and associated layer thickness are enhanced for larger radiation parameter.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A case study of using x-ray Thomson scattering to diagnose the in-flight plasma conditions of DT cryogenic implosions

The design of inertial confinement fusion ignition targets requires radiation-hydrodynamics simulations with accurate models of the fundamental material properties (i.e., equation of state, opacity, and conductivity). Validation of these models is required via experimentation. A feasibility study of using spatially integrated, spectrally resolved, x-ray Thomson scattering measurements to diagnose the temperature, density, and ionization of the compressed DT shell of a cryogenic DT implosion at two-thirds convergence was conducted. Synthetic scattering spectra were generated using 1D implosion simulations from the LILAC code that were post processed with the x-ray scattering model, which is incorporated within SPECT3D. Analysis of two extreme adiabat capsule conditions showed that the plasma conditions for both compressed DT shells could be resolved.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A Bayesian approach to time-domain photonic Doppler velocimetry analysis

Photonic Doppler velocimetry (PDV) is an established technique for measuring the velocities of fast-moving surfaces in high-energy-density experiments. In the standard approach to PDV analysis, the short-time Fourier transform (STFT) is used to generate a spectrogram from which the velocity history of the target is inferred. The user chooses the form, duration, and separation of the window function. Here, in this study, we present a Bayesian approach to infer the velocity directly from the PDV oscilloscope trace, without using the spectrogram for analysis. This is clearly a difficult inference problem due to the highly periodic nature of the data, but we find that with carefully chosen prior distributions for the model parameters, we can accurately recover the injected velocity from synthetic data. We validate this method using PDV data collected at the STAR two-stage light gas gun at Sandia National Laboratories, recovering shock-front velocities in quartz that are consistent with those inferred using the STFT-based approach and are interpolated across regions of low signal-to-noise data. Although this method does not rely on the same user choices as the STFT, we caution that it can be prone to misspecification if the chosen model is not sufficient to capture the velocity behavior. Analysis using posterior predictive checks can be used to establish whether a better model is required, although more complex models come with additional computational cost, often taking more than several hours to converge when sampling the Bayesian posterior. We, therefore, recommend it be viewed as a complementary method to that of the STFT-based approach.

Allison, James R. [First Light Fusion Ltd., Yarnto↗