Human foot modeling towards footwear design

Comfort test of footwear is mainly based on subjective perception of the wearer and a large number of subjects are required to obtain a reliable result. Therefore, the subjective comfort test is expensive and time consuming. Although the foot size and shape of a subject can be obtained by using a three-dimensional (3D) foot scanner, it is still difficult to create foot motion animations of each subject suitable for computer simulation.In this paper, we propose a fast approach to model foot deformation and present its application in simulating interaction with footwear towards footwear design. The simulation determines deformation of foot and footwear models. It can also determine stress distribution in the footwear. Given an initial foot model and a captured foot motion, human foot animation is created first. Then, the footwear model is fitted to the foot to compute the deformation and stress in the footwear. In this article, the boundary element method (BEM) is adopted. We demonstrate the results by conducting simulation of a captured gait motion. Experimental results showed that the method can be used to simulate human gait motion, and can determine deformation of footwear.     

A new method for identifying and validating features from 2D sectional views

Feature identification is one of the key steps for 3D solids reconstruction from 2D vector engineering drawings using the volume-based method. In this paper, we propose a novel method to identify and validate features from sectional views. First, features are classified as explicit features (EPFs) and implicit features (IPFs), which are then identified in an order of priority using heuristic hints. We show that the problem of constructing EPFs can be formulated as a 0-1 integer linear program (ILP), and the IPFs are generated based on the understanding of semantic information of omitted projections in sectional views. Then, the Loop-Relation Graph (LRG) is introduced as a multi-connected-subgraph representation for describing the relations between loops and features. According to the LRG, a reasoning technique based on confidence is implemented to interactively validate features. This method can recover features without complete projections, and the level of understanding sectional views is improved. Full sections, partial sections, offset sections as well as revolved sections can be handled by our method. Several examples are provided to demonstrate the practicability of our approach.     

Sufficient conditions for the preservation of the boundedness in a numerical method for a physical model with transport memory and nonlinear damping

Departing from a finite-difference scheme to approximate solutions of a nonlinear, hyperbolic partial differential equation which generalizes the Burgers–Huxley equation from fluid dynamics, we investigate conditions on the model coefficients and the computational parameters under which positive and bounded initial data evolve into positive and bounded new approximations. The model under investigation includes nonlinear coefficients of damping and advection, and the reaction term extends the reaction law of the classical Fisher–Kolmogorov–Petrovsky–Piscounov equation. The method can be expressed in vector form in terms of a multiplicative matrix which, under certain parametric conditions, becomes an M-matrix. Using the fact that every M-matrix is non-singular and that the entries of its inverse are positive, real numbers, we establish sufficient conditions under which the method provides new, positive and bounded approximations from previous, positive and bounded data and boundary conditions. The numerical results confirm the fact that the conditions derived here are sufficient for the positivity and the boundedness of the approximations; moreover, computational experiments evidence the fact that the method still preserves these properties for values of the model and the numerical parameters outside of the analytic regions of positivity and boundedness. We point out that our simulations show a good agreement between the numerical approximations computed through our method and the corresponding, analytical solutions.     

NGluon: A package to calculate one-loop multi-gluon amplitudes

We present a computer library for the numerical evaluation of colour-ordered n-gluon amplitudes at one-loop order in pure Yang–Mills theory. The library uses the recently developed technique of generalised unitarity. Running in double precision the library yields reliable results for up to 14 gluons with only a small fraction of events requiring a re-evaluation using extended floating point arithmetic. We believe that the library presented here provides an important contribution to future LHC phenomenology. The program may also prove useful in cross checking results obtained by other methods. In addition, the code provides a sample implementation which may serve as a starting point for further developments.

Program summary

Program title:NGluonCatalogue identifier: AEIZ_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEIZ_v1_0.htmlProgram obtainable from: CPC Program Library, Queen?s University, Belfast, N. IrelandLicensing provisions: GNU Public LicenseNo. of lines in distributed program, including test data, etc.: 30 677No. of bytes in distributed program, including test data, etc.: 334 896Distribution format: tar.gzProgramming language: C++Computer: Any computer platform supported by the GNU compiler suite.Operating system: No specific requirements – tested on Scientific Linux 5.2.RAM: Depending on the complexity, for realistic applications like 10 gluon production in double precision below 10 MB.Classification: 11.5External routines: QCDLoop (http://qcdloop.fnal.gov/), qd (http://crd.lbl.gov/~dhbailey/mpdist/)Nature of problem: Evaluation of next-to-leading order corrections for gluon scattering amplitudes in pure gauge theory.Solution method: Purely numerical approach based on tree amplitudes obtained via Berends–Giele recursion combined with unitarity method.Restrictions: Running in double precision the number of gluons should not exceed 14.Running time: Depending on the number of external gluons between less than a millisecond (4 gluons) up to a 1 s (14 gluons) per phase space point.     

Three dimensional discontinuous Galerkin methods for Euler equations on adaptive conforming meshes

In the numerical simulation of three dimensional fluid dynamical equations, the huge computational quantity is a main challenge. In this paper, the discontinuous Galerkin (DG) finite element method combined with the adaptive mesh refinement (AMR) is studied to solve the three dimensional Euler equations based on conforming unstructured tetrahedron meshes, that is according the equation solution variation to refine and coarsen grids so as to decrease total mesh number. The four space adaptive strategies are given and analyzed their advantages and disadvantages. The numerical examples show the validity of our methods.     

Solving differential matrix Riccati equations by a piecewise-linearized method based on diagonal Padé approximants

Differential Matrix Riccati Equations (DMREs) appear in several branches of science such as applied physics and engineering. For example, these equations play a fundamental role in control theory, optimal control, filtering and estimation, decoupling and order reduction, etc. In this paper a new method based on a theorem proved in this paper is described for solving DMREs by a piecewise-linearized approach. This method is applied for developing two block-oriented algorithms based on diagonal Padé approximants. MATLAB versions of the above algorithms are developed, comparing, under equal conditions, accuracy and computational costs with other piecewise-linearized algorithms implemented by the authors. Experimental results show the advantages of solving stiff or non-stiff DMREs by the implemented algorithms.     

A general formulation of monotonically convergent algorithms in the control of quantum dynamics beyond the linear dipole interaction

This paper presents a general method to formulate monotonically convergent algorithms for identifying optimal control fields to manipulate quantum dynamics phenomena beyond the linear dipole interaction. The method, facilitated by a field-dependent dipole moment operator, is based on an integral equation of the first kind arising from the Heisenberg equation of motion for tracking a time-dependent, dynamical invariant observable associated with a reference control field.     

An analytic algorithm for the space–time fractional advection–dispersion equation

Fractional advection–dispersion equation (FADE) is a generalization of the classical ADE in which the first order time derivative and first and second order space derivatives are replaced by Caputo derivatives of orders 0<α?1, 0<β?1 and 1<γ?2, respectively. We use Caputo definition to avoid (i) mass balance error, (ii) hyper-singular improper integral, (iii) non-zero derivative of constant, and (iv) fractional derivative involved in the initial condition which is often ill-defined. We present an analytic algorithm to solve FADE based on homotopy analysis method (HAM) which has the advantage of controlling the region and rate of convergence of the solution series via the auxiliary parameter ? over the variational iteration method (VIM) and homotopy perturbation method (HPM). We find that the proposed method converges to the numerical/exact solution of the ADE as the fractional orders α, β, γ tend to their integral values. Numerical examples are given to illustrate the proposed algorithm. Example 5 describes the intermediate process between advection and dispersion via Caputo fractional derivative.     

Accurate finite element method for atomic calculations based on density functional theory and Hartree–Fock method

An accurate finite element method is developed for atomic calculations based on density functional theory (DFT) within local density approximation (LDA) and Hartree–Fock (HF) method. The radial wave functions are expanded by cubic Hermite spline functions on a uniform mesh for , and all the associated integrals are analytically evaluated in conjunction with fitting procedures of the Hartree and the exchange–correlation potentials to the same cubic Hermite spline functions using a set of recurrence formulas. The total energy of atoms systematically converges from above, and the error algebraically decays as the mesh spacing decreases. When the mesh spacing d is taken to be , the total energy for an atom of atomic number Z can be calculated within error of ¹⁰−7 hartree for both the LDA and HF methods. The equal applicability of the method to DFT and the HF method with a similar degree of high accuracy enables the method to be a reliable platform for development of new functionals in DFT such as hybrid functionals.