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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

A Continuous-Discontinuous Galerkin Method for Electromagnetic Simulations Based on an All-Frequency Stable Formulation

In this paper, a potential-based partial-differential formulation, called the all-frequency stable formulation, is presented for the accurate and robust simulation of electromagnetic problems at all frequencies. Due to its stability from (near) dc to microwave frequencies, this formulation can be applied to simulate wide-band and multiscale problems without encountering the infamous low-frequency breakdown issue or applying basis function decompositions such as the tree-cotree splitting technique. To provide both efficient and flexible numerical solutions to the electromagnetic formulation, a mixed continuous-discontinuous Galerkin (CDG) method is proposed and implemented. In regions with homogeneous media, the continuous Galerkin method is employed to avoid the introduction of duplicated degrees of freedom (DoFs) on the elemental interfaces, while on the interfaces of two different media, the discontinuous Galerkin method is applied to permit the jump of the normal components of the electromagnetic fields. Numerical examples are provided to validate and demonstrate the proposed numerical solver for problems in a wide electromagnetic spectrum.

Yan, Su↗
NASA NTRS2004

Continuation Methods for Qualitative Analysis of Aircraft Dynamics

A class of numerical methods for constructing bifurcation curves for systems of coupled, non-linear ordinary differential equations is presented. Foundations are discussed, and several variations are outlined along with their respective capabilities. Appropriate background material from dynamical systems theory is presented.

Cummings, Peter A.↗
NASA NTRS1980

Analysis of two-dimensional incompressible flows by a subsurface panel method

A new approach to panel methods is explored for two-dimensional steady incompressible flows. The method uses linear distributions of sources and vortices on straight-line panels, but satisfies boundary conditions on the actual body surface, at nodes that are also end points of the panels. The result is continuity in body-surface velocity distribution, without recourse to numerical quadrature for the velocity influence coefficients. The method is unusually sensitive to the distribution of the nodes. For example, it almost always fails to give acceptable results when the nodes are distributed randomly. However, the continuity of the velocity distribution makes possible a unique node redistribution scheme, which may be iterated to give accurate results reliably.

Moran, J.↗
DOE OSTI2021

A solution framework for linear PDE-constrained mixed-integer problems

Abstract We present a general numerical solution method for control problems with state variables defined by a linear PDE over a finite set of binary or continuous control variables. We show empirically that a naive approach that applies a numerical discretization scheme to the PDEs to derive constraints for a mixed-integer linear program (MILP) leads to systems that are too large to be solved with state-of-the-art solvers for MILPs, especially if we desire an accurate approximation of the state variables. Our framework comprises two techniques to mitigate the rise of computation times with increasing discretization level: First, the linear system is solved for a basis of the control space in a preprocessing step. Second, certain constraints are just imposed on demand via the IBM ILOG CPLEX feature of a lazy constraint callback. These techniques are compared with an approach where the relations obtained by the discretization of the continuous constraints are directly included in the MILP. We demonstrate our approach on two examples: modeling of the spread of wildfire and the mitigation of water contamination. In both examples the computational results demonstrate that the solution time is significantly reduced by our methods. In particular, the dependence of the computation time on the size of the spatial discretization of the PDE is significantly reduced.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1996

The Prediction of Ducted Fan Engine Noise Via a Boundary Integral Equation Method

A computationally efficient Boundary Integral Equation Method (BIEM) for the prediction of ducted fan engine noise is presented. The key features of the BIEM are its versatility and the ability to compute rapidly any portion of the sound field without the need to compute the entire field. Governing equations for the BIEM are based on the assumptions that all acoustic processes are linear, generate spinning modes, and occur in a uniform flow field. An exterior boundary value problem (BVP) is defined that describes the scattering of incident sound by an engine duct with arbitrary profile. Boundary conditions on the duct walls are derived that allow for passive noise control treatment. The BVP is recast as a system of hypersingular boundary integral equations for the unknown duct surface quantities. BIEM solution methodology is demonstrated for the scattering of incident sound by a thin cylindrical duct with hard walls. Numerical studies are conducted for various engine parameters and continuous portions of the total pressure field are computed. Radiation and duct propagation results obtained are in agreement with the classical results of spinning mode theory for infinite ducts.

Boundary integral equation↗
NASA NTRS1991

Numerical simulation of the reactive flow in advanced (HSR) combustors using KIVA-2

Recent work has been done with the goal of establishing ultralow emission aircraft gas turbine combustors. A significant portion of the effort is the development of three dimensional computational combustor models. The KIVA-II computer code which is based on the Implicit Continuous Eulerian Difference mesh Arbitrary Lagrangian Eulerian (ICED-ALE) numerical scheme is one of the codes selected by NASA to achieve these goals. This report involves a simulation of jet injection through slanted slots within the Rich burn/Quick quench/Lean burn (RQL) baseline experimental rig. The RQL combustor distinguishes three regions of combustion. This work specifically focuses on modeling the quick quench mixer region in which secondary injection air is introduced radially through 12 equally spaced slots around the mixer circumference. Steady state solutions are achieved with modifications to the KIVA-II program. Work currently underway will evaluate thermal mixing as a function of injection air velocity and angle of inclination of the slots.

Winowich, Nicholas S.↗
DOE OSTI2020

Characterizing Variability in Lignocellulosic Biomass: A Review

Feedstock variability is a significant barrier to the scale-up and commercialization of lignocellulosic biofuel technologies. Variability in feedstock characteristics and behavior creates numerous challenges to the biorefining industry by affecting continuous operation and biofuels yields. Recently, feedstock variability is understood and explained largely on the basis of chemical composition. Physical and mechanical properties and behavior of lignocellulosic feedstock in various unit operations, studied through advanced analytical methods, can further explain variability. Such studies will enable us in developing processes and designing equipment to improve operation and conversion performance. In this perspective, we review several advanced analytical methods that measure density, moisture content, thermal properties, flowability, grindability, rheology properties, and micromorphological characteristics. We also discuss the correlations and interactions among these properties that reflect the complexity of lignocellulosic biomass as a feedstock and the associated quality metrics and logistics of supplying consistent quality feedstock to a biorefinery. We further examine methods that have not traditionally been used to characterize lignocellulosic feedstocks but have the potential to bridge the gap in our explanation of feedstock variability.

09 BIOMASS FUELS↗
DOE OSTI2022

Real-time laser ultrasonic monitoring of laser-induced thermal processes

Abstract Intra- and inter-layer integrity of components fabricated with advanced manufacturing techniques, such as laser powder bed fusion, is dependent upon rapid heating, melting, and solidification processes. There is a need for new techniques to provide in situ feedback of these processes. Here a laser-based ultrasonic technique to probe thermal effects induced by a high-power continuous wave laser in titanium samples is described. Numerical simulations were performed to show that, for a spatially uniform heating beam, laser-induced surface acoustic waves are strongly influenced by surface heating conditions, are dispersive in the case of rapid heating, and that an abrupt velocity reduction happens upon the onset of surface melting. Furthermore, laser-based ultrasound experimental results which monitor the transient change of surface wave travel time associated with high power laser surface heating are provided. A pulsed laser is used to generate high frequency surface acoustic waves that propagate through the laser-heated region and are detected using a photorefractive crystal-based interferometer. Qualitative agreement is observed between theory and experiment with both showing a rapid reduction in the surface wave velocity at the onset of illumination and further decrease in surface wave velocity associated with melting. It is demonstrated that changes in the surface wave velocity can be used to track local heating and detect the onset of surface melting in real time.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Disentangling the chemistry of Australian plant exudates from a unique historical collection

For thousands of years, the unique physicochemical properties of plant exudates have defined uses in material culture and practical applications. Native Australian plant exudates, including resins, kinos, and gums, have been used and continue to be used by Aboriginal Australians for numerous technical and cultural purposes. A historic collection of well-preserved native Australian plant exudates, assembled a century ago by plant naturalists, gives a rare window into the history and chemical composition of these materials. Here we report the full hierarchical characterization of four genera from this collection, Xanthorrhoea, Callitris, Eucalyptus, and Acacia, from the local elemental speciation, to functional groups and main molecular markers. We use high-resolution X-ray Raman spectroscopy (XRS) to achieve bulk-sensitive chemical speciation of these plant exudates, including insoluble, amorphous, and cross-linked fractions, without the limitation of invasive and/or surface specific methods. Combinatorial testing of the XRS data allows direct classification of these complex natural species as terpenoid, aromatic, phenolic, and polysaccharide materials. Differences in intragenera chemistry was evidenced by detailed interpretation of the XRS spectral features. We complement XRS with Fourier-transform infrared (FT-IR) spectroscopy, gas chromatography–mass spectrometry (GC-MS), and pyrolysis–GC-MS (Py-GC-MS). This multimodal approach provides a fundamental understanding of the chemistry of these natural materials long used by Aboriginal Australian peoples.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Rheology of concentrated fiber suspensions with a load-dependent friction coefficient

Here we investigate the effects of fiber aspect ratio, roughness, flexibility, and volume fraction on the rheology of concentrated suspensions in a steady shear flow using direct numerical simulations. We model the fibers as inextensible continuous flexible slender bodies with the Euler-Bernoulli beam equation governing their dynamics suspended in an incompressible Newtonian fluid. The fiber dynamics and fluid flow coupling is achieved using the immersed boundary method. In addition, the fiber surface roughness might lead to interfiber contacts, resulting in normal and tangential forces between the fibers, which follow Coulomb’s law of friction. The surface roughness is modeled as hemispherical protrusions on the fiber surfaces. Their deformation results in a normal load-dependent friction coefficient. Our simulations accurately predict the experimentally observed shear thinning in fiber suspensions. Furthermore, we find that the suspension viscosity η increases with increasing the volume fraction, roughness, fiber rigidity, and aspect ratio. The increase in η is the macroscopic manifestation of a similar increase in the microscopic contact contribution to the total stress with these parameters. In addition, we observe positive and negative first N 1 and second N 2 normal stress differences, respectively, with |N 2 | < |N 1 |, in agreement with previous experiments. Last, we propose a modified Maron-Pierce law to quantify the reduction in the jamming volume fraction by increasing the fiber aspect ratio and roughness. Our results and analysis establish the use of fiber surface tribology to tune the suspension flow behavior.

42 ENGINEERING↗
DOE OSTI2022

Laplacian-level meta-generalized gradient approximation for solid and liquid metals

In this work, we derive and motivate a Laplacian-level, orbital-free meta-generalized-gradient approximation (LL-MGGA) for the exchange-correlation energy, targeting accurate ground-state properties of sp and sd metallic condensed matter, in which the density functional for the exchange-correlation energy is only weakly nonlocal due to perfect long-range screening. Our model for the orbital-free kinetic energy density restores the fourth-order gradient expansion for exchange to the r 2 SCAN meta-GGA [Furness et al., J. Phys. Chem. Lett. 11, 8208 (2020)], yielding a LL-MGGA we call OFR2. OFR2 matches the accuracy of SCAN for prediction of common lattice constants and improves the equilibrium properties of alkali metals, transition metals, and intermetallics that were degraded relative to the PBE GGA values by both SCAN and r 2 SCAN. We compare OFR2 to the r 2 SCAN -L LL-MGGA [Mejia-Rodriguez and Trickey, Phys. Rev. B 102, 121109(R) (2020)] and show that OFR2 tends to outperform r 2 SCAN -L for the equilibrium properties of solids, but r 2 SCAN -L much better describes the atomization energies of molecules than OFR2 does. For best accuracy in molecules and nonmetallic condensed matter, we continue to recommend SCAN and r 2 SCAN. Numerical performance is discussed in detail, and our paper provides an outlook to machine learning.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Topological Langmuir-cyclotron wave

A theory is developed to describe the topological Langmuir-cyclotron wave (TLCW), a topological excitation in magnetized plasmas recently identified by numerical simulations. As a topological wave in a continuous medium, the TLCW propagates unidirectionally without scattering in complex boundaries and can be explored as an effective mechanism to energize particles. We show that, because momentum space in continuous media is contractible in general, the topology of the wave bundles is trivial over momentum space that contains no degeneracy points. This is in stark contrast to condensed matters with periodic lattice structures that impose nontrivial topology on momentum space. In continuous media without lattice structures, nontrivial topology of the eigenmode bundles manifests over phase space, and it is the nontrivial topology over phase space that underpins the topological excitations, such as the TLCW. It is shown that the TLCW can be faithfully modeled by a generic tilted Dirac cone in phase space, whose entire spectrum, including the spectral flow, is given.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Design of quantum optical experiments with logic artificial intelligence

Logic Artificial Intelligence (AI) is a subfield of AI where variables can take two defined arguments, True or False, and are arranged in clauses that follow the rules of formal logic. Several problems that span from physical systems to mathematical conjectures can be encoded into these clauses and solved by checking their satisfiability (SAT). In contrast to machine learning approaches where the results can be approximations or local minima, Logic AI delivers formal and mathematically exact solutions to those problems. In this work, we propose the use of logic AI for the design of optical quantum experiments. We show how to map into a SAT problem the experimental preparation of an arbitrary quantum state and propose a logic-based algorithm, called Klaus, to find an interpretable representation of the photonic setup that generates it. We compare the performance of Klaus with the state-of-the-art algorithm for this purpose based on continuous optimization. We also combine both logic and numeric strategies to find that the use of logic AI significantly improves the resolution of this problem, paving the path to developing more formal-based approaches in the context of quantum physics experiments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Experimental Investigation of Low-Frequency Distributed Acoustic Sensor Responses to Two Parallel Propagating Fractures

Low-frequency distributed acoustic sensing (LF-DAS) is a diagnostic tool for hydraulic fracture propagation with far-field monitoring using fiber optic sensors. LF-DAS senses strain rate variation caused by stress field change due to fracture propagation. Fiber optic sensors are installed in the monitoring wells in the vicinity of a fractured well. From the strain responses, fracture propagation can be evaluated. To understand subsurface conditions with multiple propagating fractures, a laboratory-scale hydraulic fracture experiment was performed simulating the LF-DAS response to fracture propagation with embedded distributed optical fiber strain sensors under these conditions. The experiment was performed using a transparent cube of epoxy with two parallel radial initial flaws centered in the cube. Fluid was injected into the sample to generate fractures along the initial flaws. The experiment used distributed high-definition fiber optic strain sensors with tight spatial resolutions. The sensors were embedded at two different locations on opposite sides of the initial flaws, serving as observation/monitoring locations. We also employed finite element modeling to numerically solve the linear elastic equations of equilibrium continuity and stress–strain relationships. The measured strains from the experiment were compared to simulation results from the finite element model. The experimentally derived strain and strain-rate waterfall plots from this study show the responses to both fractures propagating, while the fracture at the lower position took most of the fluid during the experiment. Interestingly, a fracture first began propagating from the upper flaw of the two flaws, but once the lower fracture was initiated, it grew much faster than the upper fracture. Both fibers were intercepted by the lower fracture, further verifying the strain signature as a fracture is approaching and intersecting an offset fiber.

Chemistry↗
DOE OSTI2020

Microwave Radiometer: Meteorological Temperature Profiler (raw)

The microwave radiometer profiler (MWRP) provides vertical profiles of temperature, humidity, and cloud liquid water content as a function of height or pressure at approximately 5-minute intervals for nearly all weather conditions. The profiles are derived from measurements of absolute microwave radiances (expressed as "brightness temperatures") obtained at 12 frequencies in the range of 22-30 GHz and 51-59 GHz. The data are useful for input to numerical weather forecast models and others that require continuous, high-temporal-resolution profiles. The microwave radiance measurements are useful for testing models of microwave radiation transfer used to derive the profiles.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1975

Perturbation theory for charged-particle transport in one dimension

Perturbation theory, when applied to charged-particle transport, generates a series solution that requires a double quadrature per term. The continuity of higher-order terms leads to numerical evaluation of the series. The high rate of convergence of the series makes the method a practical tool for charged-particle transport problems. The coupling of the neutron component in the case of proton transport in tissue does not greatly alter the rate of convergence. The method holds promise for a practical high-energy proton transport theory.

Wilson, J. W.↗
NASA NTRS1980

Composition of the nighttime ionospheric F 1 region near the magnetic equator

The effects of vertical E x B transport on NO(+), O2(+) and O(+) densities in the nighttime equatorial ionospheric F 1 region are investigated. Ion densities are calculated as functions of altitude, latitude and local time by the numerical solution of coupled, time-dependent ion continuity equations, taking into account production, loss by charge exchange and dissociative recombination and transport by diffusion and E x B drift. The results of the calculations are compared with measurements of NO(+), O2(+) and O(+) ion densities obtained at low altitudes by a mass spectrometer on board the Atmospheric Explorer C satellite, and are found to be consistent with the observations, suggesting that in the equatorial region, vertical transport by E x B drift is primarily responsible for producing the observed NO(+), O2(+) and O(+) density profiles. In addition, the reaction of O2(+) with N(4S) is found to be an important sink for O2(+) and a source of NO(+) ions. Implications of the observed and calculated near constancy of electron and ion densities with altitude when NO(+) is the dominant ion on the growth of large-scale irregularities are also considered

Anderson, D. N.↗