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,405 records · Page 78

Pattern Search Methods for Linearly Constrained Minimization

We extend pattern search methods to linearly constrained minimization. We develop a general class of feasible point pattern search algorithms and prove global convergence to a Karush-Kuhn-Tucker point. As in the case of unconstrained minimization, pattern search methods for linearly constrained problems accomplish this without explicit recourse to the gradient or the directional derivative. Key to the analysis of the algorithms is the way in which the local search patterns conform to the geometry of the boundary of the feasible region.

Lewis, Robert Michael↗

Using TRMM Data To Understand Interannual Variations In the Tropical Water Balance

A significant element of the science rationale for TRMM centered on assembling rainfall data needed to validate climate models-- climatological estimates of precipitation, its spatial and temporal variability, and vertical modes of latent heat release. Since the launch of TRMM, a great interest in the science community has emerged for quantifying interannual variability (IAV) of precipitation and its relationship to sea-surface temperature (SST) changes. The fact that TRMM has sampled one strong warm/ cold ENSO couplet, together with the prospect for a mission lifetime approaching ten years, has bolstered this interest in these longer time scales. Variability on a regional basis as well as for the tropics as a whole is of concern. Our analysis of TRMM results so far has shown surprising lack of concordance between various algorithms in quantifying IAV of precipitation. The first objective of this talk is to quantify the sensitivity of tropical precipitation to changes in SSTs. We analyze performance of the 3A11, 3A25, and 3B31 algorithms and investigate their relationship to scattering-- based algorithms constructed from SSM/I and TRMM 85 kHz data. The physical basis for the differences (and similarities) in depicting tropical oceanic and land rainfall will be discussed. We argue that scattering-based estimates of variability constitute a useful upper bound for precipitation variations. These results lead to the second question addressed in this talk-- How do TRMM precipitation / SST sensitivities compare to estimates of oceanic evaporation and what are the implications of these uncertainties in determining interannual changes in large-scale moisture transport? We summarize results of an analysis performed using COADS data supplemented by SSM/I estimates of near-surface variables to assess evaporation sensitivity to SST. The response of near 5 W sq m/K is compared to various TRMM precipitation sensitivities. Implied moisture convergence over the tropics and its sensitivity to errors of these algorithms is discussed.

Robertson, Franklin R.↗

Real-Time Parameter Estimation in the Frequency Domain

A method for real-time estimation of parameters in a linear dynamic state space model was developed and studied. The application is aircraft dynamic model parameter estimation from measured data in flight for indirect adaptive or reconfigurable control. Equation error in the frequency domain was used with a recursive Fourier transform for the real-time data analysis. Linear and nonlinear simulation examples and flight test data from the F-18 High Alpha Research Vehicle HARV) were used to demonstrate that the technique produces accurate model parameter estimates with appropriate error bounds. Parameter estimates converged in less than 1 cycle of the dominant dynamic mode natural frequencies, using control surface inputs measured in flight during ordinary piloted maneuvers. The real-time parameter estimation method has low computational requirements, and could be implemented aboard an aircraft in real time.

Morelli, Eugene A.↗

A mixed shear flexible finite element for the analysis of laminated plates

A mixed shear flexible finite element based on the Hencky-Mindlin type shear deformation theory of laminated plates is presented and their behavior in bending is investigated. The element consists of three displacements, two rotations, and three moments as the generalized degrees of freedom per node. The numerical convergence and accuracy characteristics of the element are investigated by comparing the finite element solutions with the exact solutions. The present study shows that reduced-order integration of the stiffness coefficients due to shear is necessary to obtain accurate results for thin plates.

Putcha, N. S.↗

Adaptive remeshing method for finite-element thermal analysis

A finite-element remeshing approach that makes use of quadrilateral and triangular elements is described. The approach uses the solution on a previous mesh to create a new mesh. Meshes are completely unstructured with highly refined elements in regions of steep gradients and larger elements where gradients are smaller. Studies of convergence rates for heat conduction problems with exact solutions show that for problems with highly localized solution variations, the remeshing approach gives smaller solution errors with fewer unknowns than refinement of uniform, structured meshes.

Thornton, Earl A.↗

Numerical simulation and analysis of a prefrontal squall line. I - Observations and basic simulation results. II - Propagation of the squall line as an internal gravity wave

Results of an observational and numerical study of the squall line that occurred on June 17-18, 1978, are presented. This squall line was initially triggered by the strong surface convergence along a cold front and stretched from Illinois to the Texas Panhandle. The squall line was aligned with the surface front during its initial development, but then propagated faster than the front, resulting in a separation of approximately 200 km by 0300 UTC and 300-400 km by 0600 UTC. The Colorado State University Regional Atmospheric Modeling System is used to model the squall-line development and propagation. Results are described from several experiments that tested the sensitivity to the use of the Kuo-type cumulus parameterization scheme and grid-scale microphysical processes. The movement of the squall line in the model is shown to be due to the propagation of a deep tropospheric internal gravity wave in a CISK-like process. The thermal and dynamic perturbations associated with the hypothesized wave are shown to be consistent with internal gravity wave theory.

Cram, Jennifer M.↗

Space-Split Algorithm for Sensitivity Analysis of Discrete Chaotic Systems With Multidimensional Unstable Manifolds

Accurate approximations of the change of a system's output and its statistics with respect to the input are highly desired in computational dynamics. Ruelle's linear response theory provides breakthrough mathematical machinery for computing the linear response of chaotic dynamical systems. In this paper, we propose an algorithm for sensitivity analysis of discrete chaos with an arbitrary number of positive Lyapunov exponents. We combine the concept of perturbation space-splitting, which regularizes Ruelle's original expression, together with measure-based parameterization of the expanding subspace. We use these tools to rigorously derive trajectory-following recursive relations that converge exponentially fast, and construct a memory-efficient Monte Carlo scheme for derivatives of the output statistics. Thanks to the regularization and lack of simplifying assumptions on the system's behavior, our method is immune to the common problems of other popular methods such as the exploding tangent solutions and unphysical shadowing directions. Here, we provide a ready-to-use algorithm, analyze its complexity, and demonstrate several numerical examples of sensitivity computation using physically-inspired low-dimensional systems.

97 MATHEMATICS AND COMPUTING↗

Converging shock wave focusing and interaction with a target

Converging shock waves in liquids can be used efficiently in the research of the extreme state of matter and in various applications. In this paper, the recent results related to the interaction of a shock wave with plasma preliminarily formed in the vicinity of the shock wave convergence are presented. The shock wave is produced by the underwater electrical explosion of a spherical wire array. The plasma is generated prior to the shock wave's arrival by a low-pressure gas discharge inside a quartz capillary placed at the equatorial plane of the array. Analysis of the Stark broadening of H{sub α} and H{sub β} spectral lines and line-to-continuum ratio, combined with the ratio of the relative intensities of carbon C III/C II and silicon Si III/Si II lines, were used to determine the plasma density and temperature evolution. It was found that during the first ∼200 ns with respect to the beginning of the plasma compression by the shock wave and when the spectral lines are resolved, the plasma density increases from 2 × 10{sup 17 }cm{sup −3} to 5 × 10{sup 17 }cm{sup −3}, while the temperature remains at the same value of 3–4 eV. Further, following the model of an adiabatically imploding capillary, the plasma density increases >10{sup 19 }cm{sup −3}, leading to the continuum spectra obtained experimentally, and the plasma temperature >30 eV at radii of compression of ≤20 μm. The data obtained indicate that the shock wave generated by the underwater electrical explosion of a spherical wire array retains its uniformity during the main part of its convergence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Changes in Moisture Flux over the Tibetan Plateau during 1979-2011: Insights from a High Resolution Simulation

Net precipitation (precipitation minus evapotranspiration, P-E) changes between 1979 and 2011 from a high resolution regional climate simulation and its reanalysis forcing are analyzed over the Tibet Plateau (TP) and compared to the global land data assimilation system (GLDAS) product. The high resolution simulation better resolves precipitation changes than its coarse resolution forcing, which contributes dominantly to the improved P-E change in the regional simulation compared to the global reanalysis. Hence, the former may provide better insights about the drivers of P-E changes. The mechanism behind the P-E changes is explored by decomposing the column integrated moisture flux convergence into thermodynamic, dynamic, and transient eddy components. High-resolution climate simulation improves the spatial pattern of P-E changes over the best available global reanalysis. High-resolution climate simulation also facilitates new and substantial findings regarding the role of thermodynamics and transient eddies in P-E changes reflected in observed changes in major river basins fed by runoff from the TP. The analysis revealed the contrasting convergence/divergence changes between the northwestern and southeastern TP and feedback through latent heat release as an important mechanism leading to the mean P-E changes in the TP.

high resolution simulation↗

Exploiting rRNA operon copy number to investigate bacterial reproductive strategies

The potential for rapid reproduction is a hallmark of microbial life, but microbes in nature must also survive and compete when growth is constrained by resource availability. Successful reproduction requires different strategies when resources are scarce and when they are abundant but a systematic framework for predicting these reproductive strategies in bacteria has not been available. Here, we show that the number of ribosomal RNA operons (rrn) in bacterial genomes predicts two important components of reproduction—growth rate and growth efficiency—which are favoured under contrasting regimes of resource availability. We find that the maximum reproductive rate of bacteria doubles with a doubling of rrn copy number, and the efficiency of carbon use is inversely related to maximal growth rate and rrn copy number. We also identify a feasible explanation for these patterns: the rate and yield of protein synthesis mirror the overall pattern in maximum growth rate and growth efficiency. Furthermore, comparative analysis of genomes from 1,167 bacterial species reveals that rrn copy number predicts traits associated with resource availability, including chemotaxis and genome streamlining. Genome-wide patterns of orthologous gene content covary with rrn copy number, suggesting convergent evolution in response to resource availability. Our findings imply that basic cellular processes adapt in contrasting ways to long-term differences in resource availability. They also establish a basis for predicting changes in bacterial community composition in response to resource perturbations using rrn copy number measurements or inferences.

59 BASIC BIOLOGICAL SCIENCES↗

Energy Flow in Thin Shell Implosions and Explosions

Energy flow and balance in convergent systems beyond petapascal energy densities controls the fate of late-stage stars and the potential for controlling thermonuclear inertial fusion ignition. Time-resolved x-ray self-emission imaging combined with a Bayesian inference analysis is used to describe the energy flow and the potential information stored in the rebounding spherical shock at 0.22 PPa (2.2 Gbar or billions of atmospheres pressure). This analysis, together with a simple mechanical model, describes the trajectory of the shell and the time history of the pressure at the fuel-shell interface, ablation pressure, and energy partitioning including kinetic energy of the shell and internal energy of the fuel. Here, the techniques used here provide a fully self-consistent uncertainty analysis of integrated implosion data, a thermodynamic-path independent measurement of pressure in the petapascal range, and can be used to deduce the energy ow in a wide variety of implosion systems to petapascal energy densities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Subtropical Jet Regulates Arabian Winter Precipitation: A Viable Mechanism

The present study explores the mechanism governing wintertime (November–April) precipitation over the Arabian Peninsula (AP) using a 17-yr-long (2002–18) high-resolution WRF simulation. The composite analysis of strong precipitation events suggests that the equatorward extension of the upper-level jet together with the embedded upper-level trough creates a positive (cyclonic) midlevel vorticity and subsequently generates an anomalous lower-level convergence through Ekman pumping. This leads to the development of an anomalous surface low, which is further enhanced in the presence of the existing Red Sea trough over the AP. This surface low weakens the persistent anticyclone over the AP, shifting it further eastward to the Arabian Sea. The eastward shift in the lower-level anticyclone contributes to the transport of warm, moist air from the Arabian Sea and the Red Sea toward the AP. This warm, moist air converges with the cold and dry air advected by the midlatitude jet and creates a moisture convergence zone, leading to the initiation of convection. We test the proposed mechanism through numerical experiments with modified upper-level wind and demonstrate that a strong, southward intrusion of the jet can indeed lead to precipitation over the AP. The above mechanism also explains the interannual variability of precipitation over the AP. During wet years, we notice approximately 3 m s -1 stronger jet core magnitude and about a 2° equatorward shift of the jet compared to dry years. While the equatorward extension of the jet explains about 21% of the interannual variability, the jet magnitude explains around 7% of the variability during wet years.

54 ENVIRONMENTAL SCIENCES↗

Numerical Experiments and Identification of Areas for Further Consideration (B639388 Subcontract Quarter II Report)

The second quarter of the project was spent carrying out numerical experiments on a variety of test matrices, using a combination of precisions, with the intent to study the numerical properties (namely, accuracy and convergence behavior) of the Conjugate Gradient (CG) method. We summarize our findings in the remainder of the document. Other activities include attending biweekly xSDK meetings. We are also currently collaborating with Steven Thomas (NREL) on the numerical stability analysis of “low-synch” Gram-Schmidt routines, for potential application within high-performance GMRES variants. This work is in progress. The subsequent quarter will be spent delving into the numerical analysis in order to provide theoretical explanation for the behavior observed in our experiments.

97 MATHEMATICS AND COMPUTING↗

A Low Fuel Convergence Path to Inertial Confinement Fusion on the National Ignition Facility (Final Report)

After a decade of trying, the grand challenge problem of achieving hot spot ignition on NIF has not been achieved using indirect-drive of a high-convergence single shell capsule inside a hohlraum. Here, the alternate concept of obtaining volume-like ignition of liquid DT fuel via multi-shell implosion using laser polar direct drive (PDD) on NIF was assessed through theoretical analysis, computer simulation, fabrication development and implosion experiments on both the Laboratory for Laser Energetics’ (LLE) Omega laser and the National Ignition Facility (NIF). This research has shown that multi-shell ignition has merit and should be further advanced by the National ICF Program to determine its ultimate limitations for achieving ignition on NIF.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A semi-implicit dynamic relaxation algorithm for static nonlinear structural analysis

A semiimplicit dynamic relaxation technique for solution of the nonlinear structural equlibrium equation is presented. A previously presented basic transient response analysis algorithm is employed, permitting use of one solution method and one software module for both static and dynamic analyses. A theoretical comparison of the method with explicit dynamic relaxation techniques shows that it offers a substantially improved convergence property without additional computational overhead.

Park, K. C.↗

Investigation of mesoscale pre-convective forcing mechanisms

The installation and testing of a numerical mesoscale boundary layer model is discussed. A data analysis is first made using the rawinsonde network which shows a strong baroclinic zone developing in the region where convection was initiated. The observed wind field shows an easterly component developing in opposition to the synoptic flow producing convergence. The numerical boundary layer model is used to diagnose the physical development of the thermal and flow fields. Two physical mechanisms are examined; topography and cloud shading. The results shows that while heated sloping topography alone does create a baroclinic zone and convergence values, cloud shading is a stronger mechanism for producing a thermal field representative of observed conditions

Mcnider, R. T.↗

Analysis models for the estimation of oceanic fields

A general model for statistically optimal estimates is presented for dealing with scalar, vector and multivariate datasets. The method deals with anisotropic fields and treats space and time dependence equivalently. Problems addressed include the analysis, or the production of synoptic time series of regularly gridded fields from irregular and gappy datasets, and the estimate of fields by compositing observations from several different instruments and sampling schemes. Technical issues are discussed, including the convergence of statistical estimates, the choice of representation of the correlations, the influential domain of an observation, and the efficiency of numerical computations.

Carter, E. F.↗

Numerical study of shock-wave/boundary layer interactions in premixed hydrogen-air hypersonic flows

A computational study of shock wave/boundary layer interactions involving premixed combustible gases, and the resulting combustion processes is presented. The analysis is carried out using a new fully implicit, total variation diminishing (TVD) code developed for solving the fully coupled Reynolds-averaged Navier-Stokes equations and species continuity equations in an efficient manner. To accelerate the convergence of the basic iterative procedure, this code is combined with vector extrapolation methods. The chemical nonequilibrium processes are simulated by means of a finite-rate chemistry model for hydrogen-air combustion. Several validation test cases are presented and the results compared with experimental data or with other computational results. The code is then applied to study shock wave/boundary layer interactions in a ram accelerator configuration. Results indicate a new combustion mechanism in which a shock wave induces combustion in the boundary layer, which then propagates outwards and downstream. At higher Mach numbers, spontaneous ignition in part of the boundary layer is observed, which eventually extends along the entire boundary layer at still higher values of the Mach number.

Yungster, Shaye↗