Upscaling of the nonlinear response of composites by homogenization
Explore the source record for details and available documents.
SEARCH · Engineering Papers
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.
Explore the source record for details and available documents.
Ever since Crombie used the Doppler spectrum of RF signals scattered from the ocean surface to provide conclusive evidence for Bragg resonant scattering, scattering from the ocean surface has served as a testing ground for the identification of governing scattering mechanisms in rough surface scattering. The moving ocean surface has the advantage for rough surface scattering that, being a dispersive medium, features with different spatial scales the tendency to travel at different speeds. Thus, by looking at the temporal signature of the scattered field, one can try to isolate the scattering mechanism. This feature has been exploited by various authors. To give a recent example, we cite the work of Plant, who examined the Doppler spectrum of the return signal to infer the presence of two-scale scattering in scattering at medium incidence angles. The problem with restricting oneself to the examination of the Doppler spectrum, as has been done traditionally, is that one restricts the phenomena under investigation to temporally homogeneous phenomena: transient phenomena, which are localized in time, are not localized in frequency space. There is increasing evidence that this type of phenomena, due to specular scattering or breaking waves, may also play an important role in determining the scattering mechanism in ocean-like surfaces. To overcome this problem, we introduce the use of the Wavelet Transform to study the frequency-temporal signatures of the scattered field from moving ocean waves. We calculate this signature using various analytic scattering theories and show that the Wavelet Transform provides a useful tool for separating the different scattering mechanisms operating in scattering in ocean-like surfaces. Next we simulate realistic nonlinear moving ocean waves and calculate the temporal scattered field signature by using the method of moments and the stop-start approximation: the surface is assumed stationary during a scattering event, but moves between incident EM pulses. We compute and compare the Doppler frequency spectrum and the Wavelet Transform of the resulting signal for various types of surfaces and show how different scattering mechanisms may be isolated by this technique. Finally, we compare these numerical results against the analytical predictions in order to obtain the domain of applicability of various scattering theories for moving ocean surfaces which may include wave breaking and specular glints.
Tipping elements are components of the Earth system which may respond nonlinearly to anthropogenic climate change by transitioning toward substantially different long-term states upon passing key thresholds or “tipping points.” In some cases, such changes could produce additional greenhouse gas emissions or radiative forcing that could compound global warming. Improved understanding of tipping elements is important for predicting future climate risks and their impacts. Here we review mechanisms, predictions, impacts, and knowledge gaps associated with 10 notable Earth system components proposed to be tipping elements. We evaluate which tipping elements are approaching critical thresholds and whether shifts may manifest rapidly or over longer timescales. Some tipping elements have a higher risk of crossing tipping points under middle-of-the-road emissions pathways and will possibly affect major ecosystems, climate patterns, and/or carbon cycling within the 21st century. However, literature assessing different emissions scenarios indicates a strong potential to reduce impacts associated with many tipping elements through climate change mitigation. The studies synthesized in our review suggest most tipping elements do not possess the potential for abrupt future change within years, and some proposed tipping elements may not exhibit tipping behavior, rather responding more predictably and directly to the magnitude of forcing. Nevertheless, uncertainties remain associated with many tipping elements, highlighting an acute need for further research and modeling to better constrain risks.
A review is conducted of LF wave nonlinear properties at comets and in the earth's foreshock, engaging such compelling questions as why there are no cometary cyclotron waves, the physical mechanism responsible for 'dispersive whiskers', and the character of a general description of linear waves. Attention is given to the nonlinear properties of LF waves, whose development is illustrated by examples of waves and their features at different distances from the comet, as well as by computer simulation results. Also discussed is a curious wave mode detected from Comet Giacobini-Zinner, both at and upstream of the bow shock/wave.
Abstract Terahertz (THz) technology is critical for quantum material physics, biomedical imaging, ultrafast electronics, and next‐generation wireless communications. However, standing in the way of widespread applications is the scarcity of efficient ultrafast THz sources with on‐demand fast modulation and easy on‐chip integration capability. Here the discovery of colossal THz emission is reported from a van der Waals (vdW) ferroelectric semiconductor NbOI 2 . Using THz emission spectroscopy, a THz generation efficiency an order of magnitude higher than that of ZnTe, a standard nonlinear crystal for ultrafast THz generation is observed. The underlying generation mechanisms associated are further uncovered with its large ferroelectric polarization by studying the THz emission dependence on excitation wavelength, incident polarization, and fluence. Moreover, the ultrafast coherent amplification and annihilation of the THz emission and associated coherent phonon oscillations by employing a double‐pump scheme are demonstrated. These findings combined with first‐principles calculations, inform a new understanding of the THz light–matter interaction in emergent vdW ferroelectrics and pave the way to develop high‐performance THz devices on them for quantum materials sensing and ultrafast electronics.
We present a study of optical bi-stability in a 3.02 refractive index at 1550nm plasma enhanced chemical vapor deposition (PECVD) silicon-rich nitride (SRN) film, as it pertains to bi-stable switching, memory applications, and thermal sensing applications. In this work we utilize an SRN ring resonator device, which we first characterize at low-power and then compare thermo-optic coefficients, (2.12 ± 0.125) × 10 −4 /°C, obtained from thermal-heating induced resonance shifts to optically induced resonance shifts as well as estimated propagation loss and absorption. We then measure the time response of this nonlinearity demonstrating the relaxation time to be 18.7 us, indicating the mechanism to be thermal in nature. Finally, we demonstrate bi-stable optical switching.
A simple mechanism is proposed to explain the observed presence in the atmosphere of open or closed cellular convection. If convection is produced by cooling concentrated near the top of the cloud layer, as in radiative cooling of stratus clouds, it develops strong descending currents which are compensated by weak ascent over most of the horizontal area, and closed cells result. Conversely, heating concentrated near the bottom of a layer, as when an air mass is heated by warm water, results in strong ascending currents compensated by weak descent over most of the area, or open cells. This mechanism is similar to the one suggested by Stommel (1962) to explain the smallness of the oceans' sinking regions. The mechanism is studied numerically by means of a two-dimensional, nonlinear Boussinesq model.
This paper applies C1 cubic Hermite polynomials embedded in an orthogonal collocation scheme to the spatial discretization of the unsteady nonlinear Burgers equation as a model of the equations of fluid mechanics. The temporal discretization is carried out by means of either a noniterative finite difference or an iterative finite difference procedure. Results of this method are compared with those of a second-order finite difference scheme and a splined-cubic Taylor's series scheme. Stability limits are derived and the matrix structure of the several schemes are compared.
The development of an ultrasonic method for the nondestructive characterization of mechanical properties of engineering material is described. The method utilizes the nonlinearity parameter measurement which describes the anharmonic behavior of the solid through measurements of amplitudes of the fundamental and of the generated second harmonic ultrasonic waves. The nonlinearity parameter is also directly related to the acoustoelastic constant of the solid which can be determined by measuring the linear dependence of ultrasonic velocity on stress. A major advantage of measurements of the nonlinearity parameter over that of the acoustoelastic constant is that it may be determined without the application of stress on the material, which makes it more applicable for in-service nondestructive characterization. The relationships between the nonlinearity parameter of second-harmonic generation and the percentage of solid solution phase in engineering materials such as heat treatable aluminum alloys was established. The acoustoelastic constants are measured on these alloys for comparison and confirmation. A linear relationship between the nonlinearity parameter and the volume fraction of second phase precipitates in the alloys is indicated.
Results of an analytical and experimental study to characterize the structural response of two compression-loaded variable stiffness composite panels are presented and discussed. These variable stiffness panels are advanced composite structures, in which tows are laid down along precise curvilinear paths within each ply and the fiber orientation angle varies continuously throughout each ply. The panels are manufactured from AS4/977-3 graphite-epoxy pre-preg material using an advanced tow placement system. Both variable stiffness panels have the same layup, but one panel has overlapping tow bands and the other panel has a constant-thickness laminate. A baseline cross-ply panel is also analyzed and tested for comparative purposes. Tests performed on the variable stiffness panels show a linear prebuckling load-deflection response, followed by a nonlinear response to failure at loads between 4 and 53 percent greater than the baseline panel failure load. The structural response of the variable stiffness panels is also evaluated using finite element analyses. Nonlinear analyses of the variable stiffness panels are performed which include mechanical and thermal prestresses. Results from analyses that include thermal prestress conditions correlate well with measured variable stiffness panel results. The predicted response of the baseline panel also correlates well with measured results.
A previously developed analytical formulation has been modified in order to more accurately account for the effects of hydrostatic stresses on the nonlinear, strain rate dependent deformation of polymer matrix composites. State variable constitutive equations originally developed for metals have been modified in order to model the nonlinear, strain rate dependent deformation of polymeric materials. To account for the effects of hydrostatic stresses, which are significant in polymers, the classical J2 plasticity theory definitions of effective stress and effective inelastic strain, along with the equations used to compute the components of the inelastic strain rate tensor, are appropriately modified. To verify the revised formulation, the shear and tensile deformation of two representative polymers are computed across a wide range of strain rates. Results computed using the developed constitutive equations correlate well with experimental data. The polymer constitutive equations are implemented within a strength of materials based micromechanics method to predict the nonlinear, strain rate dependent deformation of polymer matrix composites. The composite mechanics are verified by analyzing the deformation of a representative polymer matrix composite for several fiber orientation angles across a variety of strain rates. The computed values compare well to experimentally obtained results.
We propose a model for the jetting activity that is commonly observed in the Sun's corona, especially in the open-field regions of polar coronal holes. Magnetic reconnection is the process driving the jets and a relevant magnetic configuration is the well-known null point and fan separatrix topology. The primary challenge in explaining the observations is that reconnection must occur in a short-duration energetic burst rather than quasi-continuously as is implied by the observations of long-lived structures in coronal holes, such as polar plumes, for example. The key idea underlying our model for jets is that reconnection is forbidden for an axisymmetric null-point topology. Consequently, by imposing a twisting motion that maintains the axisymmetry, magnetic stress can be built up to large levels until an ideal instability breaks the symmetry and leads to an explosive release of energy via reconnection. Using 3D MHD simulations we demonstrate that this mechanism does produce jets with high speed and mass, driven by nonlinear Alfven waves. We discuss the implications of our results for observations of the solar corona.
Incompressible magnetohydrodynamic (MHD) turbulence and magnetic dynamos, which occur in magnetofluids with large fluid and magnetic Reynolds numbers, will be discussed. When Reynolds numbers are large and energy decays slowly, the distribution of energy with respect to length scale becomes quasi-stationary and MHD turbulence can be described statistically. In the limit of infinite Reynolds numbers, viscosity and resistivity become zero and if these values are used in the MHD equations ab initio, a model system called ideal MHD turbulence results. This model system is typically confined in simple geometries with some form of homogeneous boundary conditions, allowing for velocity and magnetic field to be represented by orthogonal function expansions. One advantage to this is that the coefficients of the expansions form a set of nonlinearly interacting variables whose behavior can be described by equilibrium statistical mechanics, i.e., by a canonical ensemble theory based on the global invariants (energy, cross helicity and magnetic helicity) of ideal MHD turbulence. Another advantage is that truncated expansions provide a finite dynamical system whose time evolution can be numerically simulated to test the predictions of the associated statistical mechanics. If ensemble predictions are the same as time averages, then the system is said to be ergodic; if not, the system is nonergodic. Although it had been implicitly assumed in the early days of ideal MHD statistical theory development that these finite dynamical systems were ergodic, numerical simulations provided sufficient evidence that they were, in fact, nonergodic. Specifically, while canonical ensemble theory predicted that expansion coefficients would be (i) zero-mean random variables with (ii) energy that decreased with length scale, it was found that although (ii) was correct, (i) was not and the expected ergodicity was broken. The exact cause of this broken ergodicity was explained, after much investigation, by greatly extending the statistical theory of ideal MHD turbulence. The mathematical details of broken ergodicity, in fact, give a quantitative explanation of how coherent structure, dynamic alignment and force-free states appear in turbulent magnetofluids. The relevance of these ideal results to real MHD turbulence occurs because broken ergodicity is most manifest in the ideal case at the largest length scales and it is in these largest scales that a real magnetofluid has the least dissipation, i.e., most closely approaches the behavior of an ideal magnetofluid. Furthermore, the effects grow stronger when cross and magnetic helicities grow large with respect to energy, and this is exactly what occurs with time in a real magnetofluid, where it is called selective decay. The relevance of these results found in ideal MHD turbulence theory to the real world is that they provide at least a qualitative explanation of why confined turbulent magnetofluids, such as the liquid iron that fills the Earth's outer core, produce stationary, large-scale magnetic fields, i.e., the geomagnetic field. These results should also apply to other planets as well as to plasma confinement devices on Earth and in space, and the effects should be manifest if Reynolds numbers are high enough and there is enough time for stationarity to occur, at least approximately. In the presentation, details will be given for both theoretical and numerical results, and references will be provided.
Aviation gas-turbine combustion demands high efficiency, wide operability and minimal trace gas emissions. Performance critical design parameters include injector geometry, combustor layout, fuel-air mixing and engine cycle conditions. The present investigation explores these factors and their impact on a radially staged low-emission aviation combustor sized for a next-generation 24,000-lbf-thrust engine. By coupling multi-fidelity computational tools, a design exploration was performed using a parameterized annular combustor sector at projected 100% takeoff power conditions. Design objectives included nitrogen oxide emission indices and overall combustor pressure loss. From the design space, an optimal configuration was selected and simulated at 7.1, 30 and 85% part-power operation, corresponding to landing-takeoff cycle idle, approach and climb segments. All results were obtained by solution of the steady-state Reynolds-averaged Navier-Stokes equations. Species concentrations were solved directly using a reduced 19-step reaction mechanism for Jet-A. Turbulence closure was obtained using a nonlinear K-epsilon model. This research demonstrates revolutionary combustor design exploration enabled by multi-fidelity physics-based simulation.
Injector geometry, physical mixing, chemical processes, and engine cycle conditions together govern performance, operability and emission characteristics of aviation gas-turbine combustion systems. The present investigation explores swirl-venturi lean direct injection combustor fundamentals, characterizing the influence of key geometric injector parameters on reacting flow physics and emission production trends. In this computational study, a design space exploration was performed using a parameterized swirl-venturi lean direct injector model. From the parametric geometry, 20 three-element lean direct injection combustor sectors were produced and simulated using steady-state, Reynolds-averaged Navier-Stokes reacting computations. Species concentrations were solved directly using a reduced 18-step reaction mechanism for Jet-A. Turbulence closure was obtained using a nonlinear κ-ε model. Results demonstrate sensitivities of the geometric perturbations on axially averaged flow field responses. Output variables include axial velocity, turbulent kinetic energy, static temperature, fuel patternation and minor species mass fractions. Significant trends have been reduced to surrogate model approximations, intended to guide future injector design trade studies and advance aviation gas-turbine combustion research.
Volcanic eruptions are an important climate driver. The impact of Pinatubo-sized eruptions has been observed and is well constrained. The magnitude and duration of volcanic winter effects after super-eruptions such as Toba remain disputed due to disagreement between the strong cooling predicted by models, and much milder climate perturbations according to the paleo data. Here we present a re-evaluated climate impact of a Toba-sized super-eruption based on up-to-date GISS ModelE simulations. In this study, we account for all known primary mechanisms that govern the evolution of the volcanic plume and their nonlinear interactions. The SO2 radiative effects are evaluated for the first time in coupled climate simulations with the interactive atmospheric chemistry module. We found that SO2 effects on photochemistry, dynamics, and radiative forcing are especially prominent. Due to strong absorption in ultraviolet, SO2 feedback on photochemistry partially offsets the limiting effect associated with aerosol microphysical processes. SO2 greenhouse warming soothes the radiative cooling exerted by sulfate aerosols. SO2 absorption in the shortwave and longwave causes radiative heating and lofting of the volcanic plume, and boosts the efficiency of SO2 impact on photochemistry. Our analysis shows that SO2 lifetime and magnitude of effects scale up and increase with the amount of emitted material. For a Pinatubo-sized eruption, SO2 feedbacks on chemistry and dynamics are relevant only during the initial stage of the volcanic plume evolution, while local SO2 concentrations are high. For a Toba-sized eruption, SO2 effects are as important as sulfate aerosols, and produce a less extreme volcanic winter.
The widely used Darcy’s law specifies a linear relation between the Darcy velocity of fluid flow and the pressure gradient that drives the flow. However, studies have shown that Darcy velocity can exhibit a nonlinear dependence on the pressure gradient in low-permeability porous media such as clay and shale when the pressure gradient is adequately low. This phenomenon is referred to as low-velocity non-Darcian flow or pre-Darcy flow. This paper provides a comprehensive review of the theories, experimental data, and modeling methods for pre-Darcy flow in low-permeability porous media. The review begins by outlining the fundamental mechanisms underlying pre-Darcy flow that regulate the unique characteristics such as nonlinear dependence of the Darcy velocity on the pressure gradient and its relevance to fluid–rock interactions. The review then proceeds to present a thorough compilation of experimental investigations performed in various low-permeability geomaterials including tight sandstones, shales, and clays. Next, empirical and theoretical models and simulation methods that have been developed to fit and interpret experimental data are reviewed. Finally, the review underscores the challenges encountered in conducting and interpreting pre-Darcy flow experiments and suggests future research directions. By analyzing previous experimental investigations, this review aims to offer a valuable resource for researchers and practitioners seeking to enhance their understanding of fluid dynamics in low-permeability geomaterials. This provides insights into the application of pre-Darcy flow in numerous natural and engineered processes such as shale oil and gas recovery, contaminant transport in low-permeability aquifers, and geological disposal of nuclear waste.
Nonlinear focusing elements can enhance the stability of particle beams in high-energy colliders by means of Landau Damping, through the tune spread which is introduced. We propose an experiment at Fermilab’s Integrable Optics Test Accelerator (IOTA) to investigate the influence of nonlinear focusing elements on the transverse stability of the beam. In this experiment, we employ an anti-damper, an active transverse feedback system, as a controlled mechanism to induce coherent beam instability. By utilizing the anti-damper, we can examine the impact of the nonlinear focusing element on the beam’s transverse stability. The stability diagram, a tool used to determine the system’s stability, will be measured using a recently demonstrated method at the LHC. This measurement is carried out experimentally by selecting specific threshold gains and measuring them for a range of phases. The stability diagram is represented by gei¿ on the complex plane. The experiment at IOTA adds insight towards the stability diagram measurement method by supplying a reduced machine impedance, to investigate the impedance’s effect on the stability diagram, as well as a larger range of phase measurements.