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 469 records · Page 26

A theory of supercritical wing sections, with computer programs and examples.

Mathematical methods for the design of supercritical wings, which depend on the numerical solution of the partial differential equations of two-dimensional gas dynamics, are developed. The main contribution is a computer program for the design of shockless transonic airfoils using the hodograph transformation and analytic continuation into the complex domain. The mathematical theory is described, and a manual for users of the programs is provided. Numerical examples are given and computational results are discussed, and the computer programs themselves are listed. The analysis routine can be used to ascertain whether the profiles behave well at off-design conditions, or to smooth coordinates and obtain a desirable shape more quickly when perfectly shockless flow is not essential.

Bauer, F.↗

Numerical Simulation of Combustion and Extinction of a Solid Cylinder in Low-Speed Cross Flow

The combustion and extinction behavior of a diffusion flame around a solid fuel cylinder (PMMA) in low-speed forced flow in zero gravity was studied numerically using a quasi-steady gas phase model. This model includes two-dimensional continuity, full Navier Stokes' momentum, energy, and species equations with a one-step overall chemical reaction and second-order finite-rate Arrhenius kinetics. Surface radiation and Arrhenius pyrolysis kinetics are included on the solid fuel surface description and a parameter Phi, representing the percentage of gas-phase conductive heat flux going into the solid, is introduced into the interfacial energy balance boundary condition to complete the description for the quasi-steady gas-phase system. The model was solved numerically using a body-fitted coordinate transformation and the SIMPLE algorithm. The effects of varying freestream velocity and Phi were studied. These parameters have a significant effect on the flame structure and extinction limits. Two flame modes were identified: envelope flame and wake flame. Two kinds of flammability limits were found: quenching at low-flow speeds due to radiative loss and blow-off at high flow speeds due to insufficient gas residence time. A flammability map was constructed showing the existence of maximum Phi above which the solid is not flammable at any freestream velocity.

Tien, J. S.↗

Statistical Mechanics and Dynamics of the Outer Solar System.I. The Jupiter/Saturn Zone

We report on numerical simulations designed to understand how the solar system evolved through a winnowing of planetesimals accreeted from the early solar nebula. This sorting process is driven by the energy and angular momentum and continues to the present day. We reconsider the existence and importance of stable niches in the Jupiter/Saturn Zone using greatly improved numerical techniques based on high-order optimized multi-step integration schemes coupled to roundoff error minimizing methods.

numerical simulations planetesimals planet formati↗

Land-Atmosphere Interactions: Successes, Problems and Prospects

After two decades of active research, a much better understanding of the broader role of biospheric processes on the local climate has emerged. A surface-albedo increase, particularly in desert border regions of the subtropics (as well as the deforested tropical regions), leads to a net surface energy deficit, which in turn leads to a relative sinking and reduced rainfall. On the other hand, studies of the influence of altered ratios of evapotranspiration and sensible fluxes, in situations where the net solar income is unchanged, show that evapotranspiration is a more desirable flux for increased precipitation and vitality of the biosphere. Besides providing water vapor and convective available potential energy (CAPE) to the lower troposphere, evapotranspiration helps in building larger CAPE before "turning on" the moist-convection. Larger CAPE in the lower troposphere enables convection to reach into the deeper atmosphere thereby heating the upper troposphere; indeed, moist-convection is also accompanied by the evaporation of falling precipitation that cools and moistens the lower atmosphere. While convective, as opposed to stratiform, precipitation reduces the fractional cloud cover; it also allows more solar radiation to reach the surface thereby invigorating surface fluxes. These, together with moist convection and associated downdrafts help to maintain the characteristic upper temperature limit(s) of the moist-land as well as oceanic regions. Regardless of the above understanding, several important problems continue to hinder the accurate simulation of a realistic land atmosphere interaction in a numerical model (both GCM and/or Meso-scale models). Among the unsolved problems are parameterization of sub-grid scale land processes that include small-scale variability of soil moisture, snow-cover and snow-physics, the biodiversity of the biosphere, orography, local drainage characteristics under natural conditions, and surface flow over the natural terrain. A well-known non-linear response of surface fluxes to these variations makes the problem of parameterizing land-atmosphere interaction processes hard-to-address and simulate, particularly in a GCM. In our presentation, we will discuss how orographic, snow-cover, and water table interactions can be included into a Simple Biosphere Model such as SiB/SSiB. Figure I shows how, in the Russian region, spring snowmelt affects the soil moisture profile. Corresponding figure 2 shows how interaction with the water table decreases the natural evapotranspiration in the Sahel region simulation. While these simulations need better validation with data, the simulations reveal that surface processes are sensitive to these parameterizations. With these developments, we continue to advance our understanding of the interaction of land with the atmosphere aloft, but the intrinsic variability of the newer parameters, e. g., hydraulic properties of the soil, diminish the positive influences of these advances on the improved climate simulation with GCMs.

Sud, Y. C.↗

Thermal Conductivity Estimation from Transient Test Data with Embedded Thermocouples using Genetic Algorithm Optimization

Inverse heat transfer methodology previously developed to estimate thermal properties of high temperature fibrous insulation from embedded thermocouples was applied to the thermoplastic polymer, polyether-ether-ketone (PEEK). A small experimental setup was utilized to test the PEEK sample between temperatures of 300 K and 525 K at atmospheric pressure. Cylindrical plugs of PEEK with three thermocouples embedded at various depths were incorporated in the test sample. The experimental PEEK thermocouple data were used as the boundary and initial conditions of a one-dimensional numerical thermal model to predict the internal temperatures of the material. The thermal conductivity of PEEK was estimated with the Continuous Genetic Algorithm optimization technique by searching for the coefficients of a functional form of thermal conductivity that minimized the difference between the experimentally measured and model predicted internal temperature values. The thermal testing, one-dimensional numerical thermal model, optimization algorithm, analysis, and results are presented.

Thermal Properties↗

Well-Balanced Second-Order Convex Limiting Technique for Solving the Serre–Green–Naghdi Equations

In this article, we introduce a numerical method for approximating the dispersive Serre–Green–Naghdi equations with topography using continuous finite elements. The method is an extension of the hyperbolic relaxation technique introduced in Guermond et al. (J Comput Phys 450:110809, 2022). It is explicit, second-order accurate in space, third-order accurate in time, and is invariant-domain preserving. It is also well balanced and parameter free. Special attention is given to the convex limiting technique when physical source terms are added in the equations. The method is verified with academic benchmarks and validated by comparison with laboratory experimental data.

97 MATHEMATICS AND COMPUTING↗

Bioinspired Diversification Approach Toward the Total Synthesis of Lycodine-Type Alkaloids

Nitrogen heterocycles (azacycles) are common structural motifs in numerous pharmaceuticals, agrochemicals, and natural products. Many powerful methods have been developed and continue to be advanced for the selective installation and modification of nitrogen heterocycles through C–H functionalization and C–C cleavage approaches, revealing new strategies for the synthesis of targets containing these structural entities. Here, we report the first total syntheses of the lycodine-type Lycopodium alkaloids casuarinine H, lycoplatyrine B, lycoplatyrine A, and lycopladine F as well as the total synthesis of 8,15-dihydrohuperzine A through bioinspired late-stage diversification of a readily accessible common precursor, N-desmethyl-β-obscurine. Key steps in the syntheses include oxidative C–C bond cleavage of a piperidine ring in the core structure of the obscurine intermediate and site-selective C–H borylation of a pyridine nucleus to enable cross-coupling reactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Stability of fractional Chern insulators with a non-Landau level continuum limit

The stability of fractional Chern insulators is widely believed to be predicted by the resemblance of their single-particle spectra to Landau levels. Here we investigate the scope of this geometric stability hypothesis by analyzing the stability of a set of fractional Chern insulators that explicitly do not have a Landau level continuum limit. By computing the many-body spectra of Laughlin states in a generalized Hofstadter model, we analyze the relationship between single-particle metrics, such as trace inequality saturation, and many-body metrics, such as the magnitude of the many-body and entanglement gaps. We show numerically that the geometric stability hypothesis holds for Chern bands that are not continuously connected to Landau levels, as well as conventional Chern bands, albeit often requiring larger system sizes to converge for these configurations.

2-dimensional systems↗

Superdiffusion from Nonabelian Symmetries in Nearly Integrable Systems

The Heisenberg spin chain is a canonical integrable model. As such, it features stable ballistically propagating quasiparticles, but spin transport is subballistic at any nonzero temperature: An initially localized spin fluctuation spreads in time t to a width t 2/3 . This exponent as well as the functional form of the dynamical spin correlation function suggest that spin transport is in the Kardar–Parisi–Zhang (KPZ) universality class. However, the full counting statistics of magnetization is manifestly incompatible with KPZ scaling. A simple two-mode hydrodynamic description, derivable from microscopic principles, captures both the KPZ scaling of the correlation function and the coarse features of the full counting statistics, but remains to be numerically validated. These results generalize to any integrable spin chain invariant under a continuous nonabelian symmetry and are surprisingly robust against moderately strong integrability-breaking perturbations that respect the nonabelian symmetry.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Deep Learning for Ensemble Forecasting

Focal Area: (2) Predictive modeling through the use of AI techniques and AI-derived model components and the use of AI and other tools to design a prediction system comprising a hierarchy of models. Science Challenge: While both climate and weather forecast systems have continued to improve due to substantial efforts to improve computational capabilities, observations, and numerical models, the atmosphere is a chaotic system, and this puts a fundamental limit on our ability to make predictions. Forecasts made by high-resolution models initialized with only slightly different atmospheric states can quickly diverge. Quantifying uncertainty in forecasts is essential to adequately understand them and to make the best-informed policy decisions particularly when it comes to hydrology, extreme weather (including extreme precipitation events), and climate.

54 ENVIRONMENTAL SCIENCES↗

Modeling the TetraSpar Floating Offshore Wind Turbine Foundation as a Flexible Structure in OrcaFlex and OpenFAST

Floating offshore wind turbine technology has seen an increasing and continuous development in recent years. When designing the floating platforms, both experimental and numerical tools are applied, with the latter often using time-domain solvers based on hydro-load estimation from a Morison approach or a boundary element method. Commercial software packages such as OrcaFlex, or open-source software such as OpenFAST, are often used where the floater is modeled as a rigid six degree-of-freedom body with loads applied at the center of gravity. However, for final structural design, it is necessary to have information on the distribution of loads over the entire body and to know local internal loads in each component. This paper uses the TetraSpar floating offshore wind turbine design as a case study to examine new modeling approaches in OrcaFlex and OpenFAST that provide this information. The study proves the possibility of applying the approach and the extraction of internal loads, while also presenting an initial code-to-code verification between OrcaFlex and OpenFAST. As can be expected, comparing the flexible model to a rigid-body model proves how motion and loads are affected by the flexibility of the structure. OrcaFlex and OpenFAST generally agree, but there are some differences in results due to different modeling approaches. Since no experimental data are available in the study, this paper only forms a baseline for future studies but still proves and describes the possibilities of the approach and codes.

17 WIND ENERGY↗

Study of thermal effects on nickel-cadmium batteries

Isothermal continuous flow calorimeter is designed to test a nickel-cadmium battery under numerous orbital conditions. This sensitive calorimeter collects cell data such as oxygen pressure and rate of heat generation, and calculates changes in enthalpy.

Foley, R. T.↗

Effects of inclination and eccentricity on optimal trajectories between earth and Venus

The true optimal transfers, including the effects of the inclination and eccentricity of the planets' orbits, between earth and Venus are presented as functions of the corresponding idealized Hohmann transfers. The method of determining the optimal transfers using the calculus of variations is presented. For every possible Hohmann window, specified as a continuous function of the longitude of perihelion of the Hohmann trajectory, the corresponding numerically exact optimal two-impulse transfers are given in graphical form. The cases for which the optimal two-impulse transfer is the absolute optimal, and those for which a three-impulse transfer provides the absolute optimal transfer are indicated. This information furnishes everything necessary for quick and accurate orbit calculations for preliminary Venus mission analysis. This makes it possible to use the actual optimal transfers for advanced planning in place of the standard Hohmann transfers.

Gravier, J.-P.↗

Mars - Paleostratigraphic restoration of buried surfaces in Tharsis montes

Volcanism in the Tharsis province of Mars occurred in several different areas and was generally continuous without large time intervals between eruptive episodes. Major lava flow units are numerous and extensive, but relatively thin. In many places, impact craters on buried surfaces project above younger flows that overlie them. A new application of crater dating methods has been developed to aid in the identification of these buried surfaces and to determine their lateral extent. The technique is especially adaptable to the Tharsis region where the stratigraphic succession of major flow units has been established by detailed geologic mapping. Knowledge of the overall stratigraphy allows correlations to be made between known and unknown surfaces by comparing their crater frequencies at diameters large enough to insure their recognition on the buried unit. The method has been applied to aid in the restoration of buried rock units and to construct a series of paleostratigraphic maps showing the sequence of major eruptive events in the Tharsis region.

Scott, D. H.↗

A numerical solution for two-dimensional Fredholm integral equations of the second kind with kernels of the logarithmic potential form

Two dimensional Fredholm integral equations with logarithmic potential kernels are numerically solved. The explicit consequence of these solutions to their true solutions is demonstrated. The results are based on a previous work in which numerical solutions were obtained for Fredholm integral equations of the second kind with continuous kernels.

Gabrielsen, R. E.↗

Elasticity Solution of an Adhesively Bonded Cover Plate of Various Geometries

The plane strain of adhesively bonded structures consisting of two different isotropic adherends is considered. By expressing the x-y components of the displacements in terms of Fourier integrals and using the corresponding boundary and continuity conditions, the integral equations for the general problem are obtained and solved numerically by applying Gauss-Chebyshev integration scheme. The shear and the normal stresses in the adhesive are calculated for various geometries and material properties for a stiffened plate under uniaxial tension. Numerical results involving the stress intensity factors and the strain energy release rate are presented. The closed-form expressions for the Fredholm kernels are provided to obtain the solution for an arbitrary geometry and material properties. For the general geometry, the contribution of the normal stress is quite significant, while for symmetric geometries, the shear stress is dominant, the normal stress vanishes if the adherends are of the same material and the same thickness.

Aksel, G. N.↗

Analysis of three-dimensional solidification interface shape

The determination of the shape of the solidification interface in a cooling melt by the thermal boundary conditions is investigated analytically, extending the analysis of Siegel (1982) and Siegel and Sosoka (1982) to the three-dimensional case. The solutions are generated by continuation into the solid region from conditions at the cooled wall, and some numerical results are presented in tables and graphs. These findings and the analytical approach employed are applicable to crystallization and metal-forming processes.

Siegel, R.↗

High spatial resolution passive microwave sounding systems

Two extensive series of flights aboard the ER-2 aircraft were conducted with the MIT 118 GHz imaging spectrometer together with a 53.6 GHz nadir channel and a TV camera record of the mission. Other microwave sensors, including a 183 GHz imaging spectrometer were flown simultaneously by other research groups. Work also continued on evaluating the impact of high-resolution passive microwave soundings upon numerical weather prediction models.

Staelin, D. H.↗