Hölder stability in Type III thermoelastodynamics

This note is concerned with the linear (and linearized) Type III thermoelastodynamic theory proposed by Green and Naghdi. We here assume that the mass density is positive and the thermal conductivity tensor is positive definite. However, we do not assume the positivity of any other tensor. In this situation, we obtain Hölder continuous dependence results on the supply terms. We also sketch how to prove the continuous dependence on the initial data.     

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).     

An FMM for waveguide problems of 2-D Helmholtz’ equation and its application to eigenvalue problems

This paper discusses an FMM for solving waveguide problems and associated eigenvalue problems for Helmholtz’ equation in a two dimensional infinite strip with homogeneous Neumann boundary condition on the sides. Layer potentials with Green’s function for this problem are evaluated efficiently with the help of the method of images and FMM. We apply FMM to solve some boundary value problems in waveguides and associated resonance frequency problems using the Sakurai–Sugiura projection method after discussing the required analytic continuation of the solutions to complex frequencies. Some numerical examples show the accuracy and the efficiency of the proposed method.     

An exact solution for the three-phase thermo-electro-magneto-elastic cylinder model and its application to piezoelectric–magnetic fiber composites

A three-phase cylindrical model for analyzing fiber composite subject to in-plane mechanical load under the coupling effects of multiple physical fields (thermo, electric, magnetic and elastic) is presented. By introducing an eigenstrain corresponding to the thermo-electro-magnetic-elastic effect, the complex multi-field coupling problem can be reduced to a formal in-plane elasticity problem for which an exact closed form solution is available. The present three-phase model can be applied to fiber/interphase/matrix composites, such that a lot of interesting thermo-electro-magnetism and stress coupling phenomena induced by the interphase layer are revealed. The present model can also be applied to fiber/matrix composites, in terms of which a generalized self-consistent method (GSCM) is developed for predicting the effective properties of piezoelectric–magnetic fiber reinforced composites. The effective piezoelectric, piezomagnetic, thermoelectric and magnetoelectric moduli can be expressed in compact explicit formulae for direct references and applications. A comparison of the predictions by the GSCM with available experimental data is presented, and interesting magnification effects and peculiar product properties are discussed. As a theoretical basis for the GSCM, the equivalence of the three sets of different average field equations in predicting the effective properties are proved, and this fact provides a strong evidence of mathematical rigor and physical realism in the formulation.     

New Green’s function for stress field and a note of its application in quantum-wire structures

The elastic field caused by the lattice mismatch between the quantum wires and the host matrix can be modeled by a corresponding two-dimensional hydrostatic inclusion subjected to plane strain conditions. The stresses in such a hydrostatic inclusion can be effectively calculated by employing the Green’s functions developed by Downes and Faux, which tend to be more efficient than the conventional method based on the Green’s function for the displacement field. In this study, Downes and Faux’s paper is extended to plane inclusions subjected to arbitrarily distributed eigenstrains: an explicit Green’s function solution, which evaluates the stress field due to the excitation of a point eigenstrain source in an infinite plane directly, is obtained in a closed-form. Here it is demonstrated that both the interior and exterior stress fields to an inclusion of any shape and with arbitrarily distributed eigenstrains are represented in a unified area integral form by employing the derived Green’s functions. In the case of uniform eigenstrain, the formulae may be simplified to contour integrals by Green’s theorem. However, special care is required when Green’s theorem is applied for the interior field. The proposed Green’s function is particularly advantageous in dealing numerically or analytically with the exterior stress field and the non-uniform eigenstrain. Two examples concerning circular inclusions are investigated. A linearly distributed eigenstrain is attempted in the first example, resulting in a linear interior stress field. The second example solves a circular thermal inclusion, where both the interior and exterior stress fields are obtained simultaneously.     

A multi-scale architecture for multi-scale simulation and its application to gas–solid flows

A multi-scale hardware and software architecture implementing the EMMS (energy-minimization multi-scale) paradigm is proven to be effective in the simulation of a two-dimensional gas–solid suspension. General purpose CPUs are employed for macro-scale control and optimization, and many integrated cores (MICs) operating in multiple-instruction multiple-data mode are used for a molecular dynamics simulation of the solid particles at the meso-scale. Many cores operating in single-instruction multiple-data mode, such as general purpose graphics processing units (GPGPUs), are employed for direct numerical simulation of the fluid flow at the micro-scale using the lattice Boltzmann method. This architecture is also expected to be efficient for the multi-scale simulation of other complex systems.     

An integrated 2-D Navier–Stokes equation and its application to 3-D internal flows

An efficient two-dimensional (2-D) analytical and numerical procedure has been proposed to investigate three-dimensional (3-D) internal flows through a passage with a spatially variable depth, in which the viscous forces act significantly on both upper and lower walls. The integral 2-D version of the Navier–Stokes equation was obtained by integrating the full Navier–Stokes equation in a 3-D form over the depth of the passage. In order to examine the validity of the integrated momentum equations, fully-developed flows in straight noncircular ducts were investigated analytically prior to numerical investigations. It has been shown that the exact solutions for circular, elliptical and equilateral triangular ducts are obtainable from the integrated Navier–Stokes equation. Having confirmed its wide applicability to internal flows, numerical computations were conducted to investigate the oscillation mechanism of a fluidic oscillator. Comparison of the present prediction and experiment reveals the validity of the present treatment.     

Synthesis and surface modification of magnetic particles for application in biotechnology and biomedicine

Magnetic particles have numerous applications in biotechnology and biomedicine. In this paper we reviewed the synthesis, surface modification and some applications of magnetic particles with focus on their synthesis and surface modification. Various methods have been developed for the production of magnetic particles （magnetic nanoparticles and magnetic composite particles）. For future application magnetic particles must be modified to obtain stability and surface functional groups. Finally, the application of magnetic particles in magnetic separation, drug delivery, hyperthermia, and magnetic resonance imaging are discussed.     

The effect of a dipolar structure on the Hölder stability in Green–Naghdi thermoelasticity

This study is concerned with the mixed initial boundary value problem for a dipolar body in the context of the thermoelastic theory proposed by Green and Naghdi. For the solutions of this problem we prove a result of Hölder’s-type stability on the supply terms. We impose middle restrictions on the thermoelastic coefficients, which are common in continuum mechanics. For the same conditions we propose a continuous dependence result with regard to the initial data.     

Fhree-point explicit compact difference scheme with arbitrary order of accuracy and its application in CFD

Based on the successive iteration in the Taylor series expansion method, a three-point explicit compact difference scheme with arbitrary order of accuracy is derived in this paper. Numerical characteristics of the scheme are studied by the Fourier analysisl Unlike the conventional compact difference schemes which need to solve the equation to obtain the unknown derivatives in each node, the proposed scheme is explicit and can achieve arbitrary order of accuracy in space. Application examples for the convectiondiffusion problem with a sharp front gradient and the typical lid-driven cavity flow are given. It is found that the proposed compact scheme is not only simple to implement and economical to use, but also is effective to simulate the convection-dominated problem and obtain high-order accurate solution in coarse grid systems.     

Three-point explicit compact difference scheme with arbitrary order of accuracy and its application in CFD

Based on the successive iteration in the Taylor series expansion method,a three-point explicit compact difference scheme with arbitrary order of accuracy is derived in this paper.Numerical characteristics of the scheme are studied by the Fourier analysis. Unlike the conventional compact difference schemes which need to solve the equation to obtain the unknown derivatives in each node,the proposed scheme is explicit and can achieve arbitrary order of accuracy in space.Application examples for the convection- diffusion problem with a sharp front gradient and the typical lid-driven cavity flow are given.It is found that the proposed compact scheme is not only simple to implement and economical to use,but also is effective to simulate the convection-dominated problem and obtain high-order accurate solution in coarse grid systems.     

An inverse problem in film/substrate indentation: extracting both the Young’s modulus and thickness of films

In an indentation test, the effective Young’s modulus of a film/substrate bilayer heterostructure varies with the indentation depth, a phenomenon known as the substrate effect. In previous studies investigating this, only the Young’s modulus of the film was unknown. Once the effective Young’s modulus of a film/substrate structure is determined at a given contact depth, the Young’s modulus of the film can be uniquely determined, i.e., there is a one-to-one relation between the Young’s modulus of the film and the film/substrate effective Young’s modulus. However, at times it is extremely challenging or even impossible to measure the film thickness. Furthermore, the precise definition of the layer/film thickness for a two-dimensional material can be problematic. In the current study, therefore, the thickness of the film and its Young’s modulus are treated as two unknowns that must be determined. Unlike the case with one unknown, there are infinite combinations of film thickness and Young’s modulus which can yield the same effective Young’s modulus for the film/substrate. An inverse problem is formulated and solved to extract the Young’s modulus and thickness of the film from the indentation depth-load curve. The accuracy and robustness of the inverse problem-solving method are also demonstrated.     

A method for Lyapunov spectrum estimation using cloned dynamics and its application to the discontinuously-excited FitzHugh–Nagumo model

This work presents a new method to calculate the Lyapunov spectrum of dynamical systems based on the time evolution of initially small disturbed copies (“clones”) of the motion equations. In this approach, it is not necessary to construct the tangent space associated with the time evolution of linearized versions of motion equations, being the Lyapunov exponents directly estimated in terms of the rate of convergence or divergence of these disturbed clones with respect to the fiducial trajectory, there being periodic correction via the Gram–Schmidt Reorthonormalization procedure. The proposed method offers the possibility of partial estimation of the Lyapunov spectrum and can also be applied to nonsmooth dynamics, since the linearization procedure is no longer required. The idea is tested for representative continuous- and discrete-time dynamical systems and validated by means of comparison with the classical method to perform this calculation. To illustrate its applicability in the nonsmooth context, the largest Lyapunov exponent of the FitzHugh–Nagumo neuronal model under discontinuous periodic excitation is calculated taking the amplitude of stimulation as control parameter. This analysis reveals some complex behaviours for this simple neuronal model, which motivates relevant discussions about the possible role of chaos in the cognitive process.     

Development of a hybrid model based on mesh and meshfree methods and its application to fluid–elastic structure interaction for free surface waves

In this paper, a hybrid scheme, Fluid–Fluid–Elastic Structure (FFES) model was developed in the time domain to address the wave breaking impact on the structure. The model is developed based on the partitioned approach with different governing equations that describe various regions of the model domain. The fluid–fluid model denotes that two different fluid models were used to describe fluid in the actual physical domain. The method is a physics-based approximation to reduce the computational time, i.e. in the far-field inviscid fluid (fully nonlinear potential flow theory model), and near to the structure, viscous fluid (Navier Stokes model) is used. The coupled model then interacts with the elastic structure (based on Euler–Bernoulli beam theory). The system of equations is strongly coupled both in space and time. The Fluid–Fluid coupling uses an implicit predictor–corrector scheme, and the fluid–structure coupling works based on an iterative scheme. This approach makes the method more robust and for future extension. Three different possibilities for introducing the coupling was identified and implemented. The model was validated against results from the analytical solution and literature. The method proposed is a reliable, robust, and efficient alternative for simulating fluid–structure interaction problems.     

A reformulated synthetic turbulence generation method for a zonal RANS–LES method and its application to zero-pressure gradient boundary layers

A synthetic turbulence generation (STG) method for subsonic and supersonic flows at low and moderate Reynolds numbers to provide inflow distributions of zonal Reynolds-averaged Navier–Stokes (RANS) – large-eddy simulation (LES) methods is presented. The STG method splits the LES inflow region into three planes where a local velocity signal is decomposed from the turbulent flow properties of the upstream RANS solution. Based on the wall-normal position and the local flow Reynolds number, specific length and velocity scales with different vorticity content are imposed at the inlet plane of the boundary layer. The quality of the STG method for incompressible and compressible zero-pressure gradient boundary layers is shown by comparing the zonal RANS–LES data with pure LES, pure RANS, and direct numerical simulation (DNS) solutions. The distributions of the time and spanwise wall-shear stress, Reynolds stress distributions, and two point correlations of the zonal RANS–LES simulations are smooth in the transition region and in good agreement with the pure LES and reference DNS findings. The STG approach reduces the RANS-to-LES transition length to less than four boundary-layer thicknesses.     

An implicit algorithm for hypoelasto-plastic and hypoelasto-viscoplastic endochronic theory in finite strain isotropic–kinematic-hardening model

This paper is concerned with objective stress update algorithm for elasto-plastic and elasto-viscoplastic endochronic theory within the framework of additive plasticity. The elastic response is stated in terms of hypoelastic model and endochronic constitutive equations are stated in unrotated frame of reference. A trivially incrementally objective integration scheme for rate constitutive equations is established. Algorithmic modulus consistent with numerical integration algorithm of constitutive equations is extracted. The implementation is validated by means of a set of simple deformation paths (simple shear, extension and rotation), two benchmark test in nonlinear mechanics (the necking of a circular bar and expansion of a thick-walled cylinder), a test which demonstrates the capabilities of the proposed model in simulation of cyclic loading and ratcheting in finite strain case (cyclically loaded notched bar) and finally, the analysis of a tensile test, which presents a shear band with a finite thickness independent of the finite element mesh using endochronic viscoplastic constitutive model.     

A micro-mechanical model of knitted fabric and its application to the analysis of buckling under tension in wale direction： buckling analysis

With the aid of the micro-mechanical model of knitted fabric proposed in Part 1 we analyze the buckling of a knitted fabric sheet when it is subjected to a tension along the wale direction. The large deformation of the fabric sheet in the critical configuration is considered and, to avoid possible deviation due to the approximation of the theory of thin plate, the three-dimensional theory of instability is used. The fabric sheet is considered as a three-dimensional body and all boundary conditions are satisfied. It is shown that the buckling of the fabric sheet is possible, two buckling modes and the corresponding buckling conditions are obtained, but only the flexural mode is physically possible as observed in experiments.The project supported by the National Natural Science Foundation of China (10272079)