On the Invariance of Compressible Navier–Stokes and Energy Equations Subject to Density-Weighted Filtering

The scale invariance properties of compressible Navier–Stokes and energy equations subject to density-weighted filtering are investigated. Scale or filter invariance forms of compressible moment equations require that two forms of generalized central second-order moments be defined—(1) product of two density-weighted sub-filter fluctuations and (2) product of one density-weighted and one un-weighted sub-filter fluctuation. The evolution equations for all required first and second order filtered moments are derived. These results provide the theoretical underpinning for variable-resolution calculations of reacting and compressible turbulent flows.     

The Spectrum of the Frobenius–Perron Operator and Its Discretization for Circle Diffeomorphisms

To investigate the dynamical behaviour of a discrete dynamical system given by a map f, it is nowadays a standard method to look at the discretization of the Frobenius–Perron operator _f w.r.t. a box-partition of the state space resulting in a transition matrix P _{f, N} of a finite Markov chain. We are interested in information about the dynamics of f obtained from the spectrum of P _{f, N} , especially for circle diffeomorphisms. Therefore the spectra of _f and P _{f, N} are investigated.     

Neimark–Sacker bifurcation and chaos control in discrete-time predator–prey model with parasites

In this paper, we discuss the qualitative behavior of a four-dimensional discrete-time predator–prey model with parasites. We investigate existence and uniqueness of positive steady state and find parametric conditions for local asymptotic stability of positive equilibrium point of given system. It is also proved that the system undergoes Neimark–Sacker bifurcation (NSB) at positive equilibrium point with the help of an explicit criterion for NSB. The system shows chaotic dynamics at increasing values of bifurcation parameter. Chaos control is also discussed through implementation of hybrid control strategy, which is based on feedback control methodology and parameter perturbation. Finally, numerical simulations are conducted to illustrate theoretical results.     

Hölder Continuity and Injectivity of Optimal Maps

Consider transportation of one distribution of mass onto another, chosen to optimize the total expected cost, where cost per unit mass transported from x to y is given by a smooth function c(x, y). If the source density f ⁺(x) is bounded away from zero and infinity in an open region ${U' \subset \mathbf{R}^n}$ , and the target density f ^?(y) is bounded away from zero and infinity on its support ${\overline{V} \subset \mathbf{R}^n}$ , which is strongly c-convex with respect to U′, and the transportation cost c satisfies the ${\mathbf{A3}_{\rm w}}$ condition of Trudinger and Wang (Ann Sc Norm Super Pisa Cl Sci 5, 8(1):143–174, 2009), we deduce the local Hölder continuity and injectivity of the optimal map inside U′ (so that the associated potential u belongs to ${C^{1,\alpha}_{loc}(U')}$ ). Here the exponent α > 0 depends only on the dimension and the bounds on the densities, but not on c. Our result provides a crucial step in the low/interior regularity setting: in a sequel (Figalli et al., J Eur Math Soc (JEMS), 1131–1166, 2013), we use it to establish regularity of optimal maps with respect to the Riemannian distance squared on arbitrary products of spheres. Three key tools are introduced in the present paper. Namely, we first find a transformation that under ${\mathbf{A3}_{\rm w}}$ makes c-convex functions level-set convex (as was also obtained independently from us by Liu (Calc Var Partial Diff Eq 34:435–451, 2009)). We then derive new Alexandrov type estimates for the level-set convex c-convex functions, and a topological lemma showing that optimal maps do not mix the interior with the boundary. This topological lemma, which does not require ${\mathbf{A3}_{\rm w}}$ , is needed by Figalli and Loeper (Calc Var Partial Diff Eq 35:537–550, 2009) to conclude the continuity of optimal maps in two dimensions. In higher dimensions, if the densities f ^± are Hölder continuous, our result permits continuous differentiability of the map inside U′ (in fact, ${C^{2,\alpha}_{loc}}$ regularity of the associated potential) to be deduced from the work of Liu et al. (Comm Partial Diff Eq 35(1):165–184, 2010).     

Solutions of the Dirac–Fock Equations and the Energy of the Electron-Positron Field

We consider atoms with closed shells, i.e. the electron number N is 2, 8, 10,..., and weak electron-electron interaction. Then there exists a unique solution γ of the Dirac–Fock equations with the additional property that γ is the orthogonal projector onto the first N positive eigenvalues of the Dirac–Fock operator . Moreover, γ minimizes the energy of the relativistic electron-positron field in Hartree–Fock approximation, if the splitting of into electron and positron subspace is chosen self-consistently, i.e. the projection onto the electron-subspace is given by the positive spectral projection of. For fixed electron-nucleus coupling constant g:=α Z we give quantitative estimates on the maximal value of the fine structure constant α for which the existence can be guaranteed.     

Homogenization of the nonlinear Kelvin–Voigt model of viscoelasticity and of the Prager model of plasticity

Denoting by the stress tensor, by the linearized strain tensor, by A the elasticity tensor, and assuming that is a convex potential, the inclusion accounts for nonlinear viscoelasticity, and encompasses both the linear Kelvin–Voigt model of solid-type viscoelasticity and the Prager model of rigid plasticity with linear kinematic strain-hardening. This relation is assumed to represent the constitutive behavior of a space-distributed system, and is here coupled with the dynamical equation. An initial- and boundary-value problem is formulated, and the existence and uniqueness of the solution are proved via classical techniques based on compactness and monotonicity. A composite material is then considered, in which the function and the tensor A rapidly oscillate in space. A two-scale model is derived via Nguetseng’s notion of two-scale convergence. This provides a detailed account of the mesoscopic state of the system. Any dependence on the fine-scale variable is then eliminated, and the existence of a solution of a new single-scale macroscopic model is proved. The final outcome is at variance with the nonlinear extension of the generalized Kelvin–Voigt model, which is based on an apparently unjustified mean-field-type hypothesis.     

Ill-posedness of the Hydrostatic Euler and Navier–Stokes Equations

We prove that the linearization of the hydrostatic Euler equations at certain parallel shear flows is ill-posed. The result also extends to the hydrostatic Navier–Stokes equations with a small viscosity.     

The Hydrodynamical Relevance of the Camassa–Holm and Degasperis–Procesi Equations

In recent years two nonlinear dispersive partial differential equations have attracted much attention due to their integrable structure. We prove that both equations arise in the modeling of the propagation of shallow water waves over a flat bed. The equations capture stronger nonlinear effects than the classical nonlinear dispersive Benjamin–Bona–Mahoney and Korteweg–de Vries equations. In particular, they accommodate wave breaking phenomena.     

Acoustic-Emission Monitoring of the Deformation and Fracture of Metal–Composite Pressure Vessels

This paper presents the results of experimental studies of damage accumulation in a metal–composite pressure vessel by pneumatic strength tests. The deformation and fracture of the composite structure accompanied by matrix cracking and fiber rupture are analyzed. It is shown that the cracks and fractures generate acoustic-emission signals of various types. The results of acoustic-emission monitoring were used to develop a criterion for ranking vessels according to the strength characteristics of the pressure composite shell.     

Analytical and experimental investigations of the pulsed air–water jet

This paper describes a new way of generating pulsed air–water jet by entraining and mixing air into the cavity of a pulsed water jet nozzle. Based on the theory of hydro-acoustics and fluid dynamics, a theoretical model which describes the frequency characteristic of the pulsed air–water jet is outlined aimed at gaining a better understanding of this nozzle for generating pulses. The calculated result indicates that as the air hold-up increases, the jet oscillation frequency has an abrupt decrease firstly, and then reaches a minimum gradually at α (air hold-up)=0.5, finally it gets increased slightly. Furthermore, a vibration test was conducted to validate the present theoretical result. By this way, the jet oscillation frequency can be obtained by analyzing the vibration acceleration of the equal strength beam affected by the jet impinging. Thereby, it is found that the experimental result shows similar trend with the prediction of the present model. Also, the relationship between vibration acceleration and cavity length for the pulsed water jet follows a similar tendency in accord with the pulsed air–water jet, i.e. there exists a maximum for each curve and the maximum occurs at the ratio of L/d₁ (the ratio of cavity length and upstream nozzle diameter) =2.5 and 2.2, respectively. In addition, experimental results on specimens impinged by the pulsed water jet and pulsed air–water jet show that the erosion depth increases slightly with air addition within a certain range of cavity length. Further, this behavior is very close to the vibration test results. As for erosion volume, the air entrained into the cavity significantly affects the material removal rate.     

Condensation models for the water–steam interface and the volume of fluid method

We assess the quantitative capabilities of three condensation models. These models are: (1) Numerical iteration technique; (2) heat flux balance equation; (3) phase field. The numerical iteration technique introduces a mass and energy transfer at the interface, if the temperature of the corresponding cell differs from the saturation value. The second approach solves the heat flux balance equation at the interface, hence, the resolution of the thermal boundary layer around the liquid-vapor interface is necessary to obtain an accurate value for the condensation rate. The third technique is based on a recently derived phase field model for boiling and condensation phenomena. The models were implemented in FLUENT and the interface was tracked explicitly with the volume of fluid (VOF) method. The models were tested on the LAOKOON facility, which measured direct contact condensation in a horizontal duct. The results showed that the phase field model fit best the experimental results.     

Influence of nanoprecipitates on the creep strength and ductility of a Fe–Ni–Al alloy

Creep studies of a duplex Fe–Ni–Al intermetallic alloy, in two microstructural states, have been carried out at temperatures between 725 and 800 °C (about 0.6 T_m). In the as-cast state, the alloy contains a large volume fraction of nanoprecipitates (50–100 nm) which confer a very high creep strength with a stress exponent of 3 and an activation energy of 280 kJ/mol. The different microstructure obtained in the second state of the alloy, obtained after annealing at 1000 °C for 24 h, leads to a much lower creep strength with a higher stress exponent as well as a large value of the apparent activation energy. While volume diffusion appears to control creep in the as-cast state, both thermal and athermal processes seem to contribute to the different creep rate of material in the annealed state. The latter also exhibits a much larger ductility (12%) relative to that observed in the as-cast material (3%), due to the presence of large numbers of interfaces between the two phases present where strain incompatibilities can be accommodated.     

A computational study of the wing–wing and wing–body interactions of a model insect

The aerodynamic interaction between the contralateral wings and between the body and wings of a model insect are studied, by using the method of numerically solving the Navier-Stokes equations over moving overset grids, under typical hovering and forward flight conditions. Both the interaction between the contralateral wings and the interaction between the body and wings are very weak, e.g. at hovering, changes in aerodynamic forces of a wing due to the present of the other wing are less than 3% and changes in aerodynamic forces of the wings due to presence of the body are less than 2%. The reason for this is as following. During each down- or up-stroke, a wing produces a vortex ring, which induces a relatively large jet-like flow inside the ring but very small flow outside the ring. The vortex rings of the left and right wings are on the two sides of the body. Thus one wing is outside vortex ring of the other wing and the body is outside the vortex rings of the left and right wings, resulting in the weak interactions.     

Setting and Analysis of the Multi-configuration Time-dependent Hartree–Fock Equations

In this paper, we formulate and analyze the multi-configuration time-dependent Hartree–Fock (MCTDHF) equations for molecular systems with pairwise interaction. This set of coupled nonlinear PDEs and ODEs is an approximation of the N-particle time-dependent Schrödinger equation based on (time-dependent) linear combinations of (time-dependent) Slater determinants. The “one-electron” wave-functions satisfy nonlinear Schrödinger-type equations coupled to a linear system of ordinary differential equations for the expansion coefficients. The invertibility of the one-body density matrix (full-rank hypothesis) plays a crucial rôle in the analysis. Under the full-rank assumption a fiber bundle structure emerges and produces unitary equivalence between different useful representations of the MCTDHF approximation. For a large class of interactions (including Coulomb potential), we establish existence and uniqueness of maximal solutions to the Cauchy problem in the energy space as long as the density matrix is not singular. A sufficient condition in terms of the energy of the initial data ensuring the global-in-time invertibility is provided (first result in this direction). Regularizing the density matrix violates energy conservation. However, global well-posedness for this system in L ² is obtained with Strichartz estimates. Eventually, solutions to this regularized system are shown to converge to the original one on the time interval when the density matrix is invertible.     

Existence and Regularity of Very Weak Solutions of the Stationary Navier–Stokes Equations

We consider the stationary Navier–Stokes equations in a bounded domain Ω in R ⁿ with smooth connected boundary, where n = 2, 3 or 4. In case that n = 3 or 4, existence of very weak solutions in L ⁿ (Ω) is proved for the data belonging to some Sobolev spaces of negative order. Moreover we obtain complete L ^q -regularity results on very weak solutions in L ⁿ (Ω). If n = 2, then similar results are also proved for very weak solutions in with any q ₀ > 2. We impose neither smallness conditions on the external force nor boundary data for our existence and regularity results.     

Estimate of the Latent Free Energy and Damage at the Tip of an Opening–Mode Crack

A method of estimating the latent elastic energy associated with the microinhomogeneity of the stress and plastic–strain fields inside the plastic zone localized near the tip of an opening–mode crack (Dugdale zone) under conditions of plane stresses is proposed. The microinhomogeneity of plastic flow upon small strain hardening is taken into account only in the form of considerable distortion of the geometry of the free surfaces of the plastic zone. The damage that developes because of release of the latent free energy is estimated depending on the magnitude of the crack opening.     

Simulation of the Emission Spectrum of SiH (A2Δ → X2Π) and Measurement of the Rotational Temperature of the A2Δ State in an Electron-Beam Plasma

The 0–0 band emission spectrum of the A² X² transition of the SiH molecule was modeled numerically. The results obtained agree well with known calculated and experimental data. The rotational temperature of the A² state of SiH in a free stream of pure monosilane (SiH₄) and in a mixture with helium (He+SiH₄) activated by an electron beam is determined by comparing calculated and experimental spectra. The assumption that the emission of SiH results from dissociative excitation of SiH₄ by electron impact is confirmed. Rotational temperatures for various monosilane concentrations and distances from the nozzle are given. The spectra obtained exhibit the emission of silicon ions at wavelengths of 412.807 and 413.089 nm.     

Spectral Analysis of the Neumann–Poincaré Operator and Characterization of the Stress Concentration in Anti-Plane Elasticity

When holes or hard elastic inclusions are closely located, stress which is the gradient of the solution to the anti-plane elasticity equation can be arbitrarily large as the distance between two inclusions tends to zero. It is important to precisely characterize the blow-up of the gradient of such an equation. In this paper we show that the blow-up of the gradient can be characterized by a singular function defined by the single layer potential of an eigenfunction corresponding to the eigenvalue 1/2 of a Neumann–Poincaré type operator defined on the boundaries of the inclusions. By comparing the singular function with the one corresponding to two disks osculating to the inclusions, we quantitatively characterize the blow-up of the gradient in terms of explicit functions. In electrostatics, our results apply to the electric field, which is the gradient of the solution to the conductivity equation, in the case where perfectly conducting or insulating inclusions are closely located.