Variational bound on energy dissipation in turbulent shear flow

We present numerical solutions to the extended Doering-Constantin variational principle for upper bounds on the energy dissipation rate in plane Couette flow, bridging the entire range from low to asymptotically high Reynolds numbers. Our variational bound exhibits structure, namely a pronounced minimum at intermediate Reynolds numbers, and recovers the Busse bound in the asymptotic regime. The most notable feature is a bifurcation of the minimizing wavenumbers, giving rise to simple scaling of the optimized variational parameters, and of the upper bound, with the Reynolds number.Comment: 4 pages, RevTeX, 5 postscript figures are available as one .tar.gz file from [email protected]

Variational bound on energy dissipation in plane Couette flow

We present numerical solutions to the extended Doering-Constantin variational principle for upper bounds on the energy dissipation rate in turbulent plane Couette flow. Using the compound matrix technique in order to reformulate this principle's spectral constraint, we derive a system of equations that is amenable to numerical treatment in the entire range from low to asymptotically high Reynolds numbers. Our variational bound exhibits a minimum at intermediate Reynolds numbers, and reproduces the Busse bound in the asymptotic regime. As a consequence of a bifurcation of the minimizing wavenumbers, there exist two length scales that determine the optimal upper bound: the effective width of the variational profile's boundary segments, and the extension of their flat interior part.Comment: 22 pages, RevTeX, 11 postscript figures are available as one uuencoded .tar.gz file from [email protected]

Application of internal state variable plasticity and damage models to welding

An internal state variable constitutive model coupled with a ductile void growth model is applied in two finite element simulations of welding. Shrinkage of a 304 L stainless steel pipe due to multipass gas tungsten arc welds is presented as an example of tracking distortion far from the weld. Weld solidification cracking in Al-6061 disks is presented as an application of the plasticity model coupled with the damage model

On parameters of the Levi-Civita solution

The Levi-Civita (LC) solution is matched to a cylindrical shell of an anisotropic fluid. The fluid satisfies the energy conditions when the mass parameter $\sigma$ is in the range $0 \le \sigma \le 1$. The mass per unit length of the shell is given explicitly in terms of $\sigma$, which has a finite maximum. The relevance of the results to the non-existence of horizons in the LC solution and to gauge cosmic strings is pointed out.Comment: Latex, no figure

Parameter Estimation and Quantitative Parametric Linkage Analysis with GENEHUNTER-QMOD

Objective: We present a parametric method for linkage analysis of quantitative phenotypes. The method provides a test for linkage as well as an estimate of different phenotype parameters. We have implemented our new method in the program GENEHUNTER-QMOD and evaluated its properties by performing simulations. Methods: The phenotype is modeled as a normally distributed variable, with a separate distribution for each genotype. Parameter estimates are obtained by maximizing the LOD score over the normal distribution parameters with a gradient-based optimization called PGRAD method. Results: The PGRAD method has lower power to detect linkage than the variance components analysis (VCA) in case of a normal distribution and small pedigrees. However, it outperforms the VCA and Haseman-Elston regression for extended pedigrees, nonrandomly ascertained data and non-normally distributed phenotypes. Here, the higher power even goes along with conservativeness, while the VCA has an inflated type I error. Parameter estimation tends to underestimate residual variances but performs better for expectation values of the phenotype distributions. Conclusion: With GENEHUNTER-QMOD, a powerful new tool is provided to explicitly model quantitative phenotypes in the context of linkage analysis. It is freely available at http://www.helmholtz-muenchen.de/genepi/downloads. Copyright (C) 2012 S. Karger AG, Base

Universality in fully developed turbulence

We extend the numerical simulations of She et al. [Phys.\ Rev.\ Lett.\ 70, 3251 (1993)] of highly turbulent flow with $15 \le$ Taylor-Reynolds number $Re_\lambda\le 200$ up to $Re_\lambda \approx 45000$, employing a reduced wave vector set method (introduced earlier) to approximately solve the Navier-Stokes equation. First, also for these extremely high Reynolds numbers $Re_\lambda$, the energy spectra as well as the higher moments -- when scaled by the spectral intensity at the wave number $k_p$ of peak dissipation -- can be described by {\it one universal} function of $k/k_p$ for all $Re_\lambda$. Second, the ISR scaling exponents $\zeta_m$ of this universal function are in agreement with the 1941 Kolmogorov theory (the better, the large $Re_\lambda$ is), as is the $Re_\lambda$ dependence of $k_p$. Only around $k_p$ viscous damping leads to slight energy pileup in the spectra, as in the experimental data (bottleneck phenomenon).Comment: 14 pages, Latex, 5 figures (on request), 3 tables, submitted to Phys. Rev.

Statistics of extremal intensities for Gaussian interfaces

The extremal Fourier intensities are studied for stationary Edwards-Wilkinson-type, Gaussian, interfaces with power-law dispersion. We calculate the probability distribution of the maximal intensity and find that, generically, it does not coincide with the distribution of the integrated power spectrum (i.e. roughness of the surface), nor does it obey any of the known extreme statistics limit distributions. The Fisher-Tippett-Gumbel limit distribution is, however, recovered in three cases: (i) in the non-dispersive (white noise) limit, (ii) for high dimensions, and (iii) when only short-wavelength modes are kept. In the last two cases the limit distribution emerges in novel scenarios.Comment: 15 pages, including 7 ps figure

Automatic Planning and Control of Robot Natural Motion Via Feedback

A feedback control strategy for the command of robot motion includes some limited automatic planning capabilities. These may be seen as sequential solution algorithms implemented by the robot arm interpreted as a mechanical analog computer. This perspective lends additional insight into the manner in which such control techniques may fail, and motivates a fresh look at requisite sensory capabilities. For more information: Kod*La

Fr-TM-align: a new protein structural alignment method based on fragment alignments and the TM-score

©2008 Pandit and Skolnick; licensee BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. This article is available from: http://www.biomedcentral.com/1471-2105/9/531doi:10.1186/1471-2105-9-531Background: Protein tertiary structure comparisons are employed in various fields of contemporary structural biology. Most structure comparison methods involve generation of an initial seed alignment, which is extended and/or refined to provide the best structural superposition between a pair of protein structures as assessed by a structure comparison metric. One such metric, the TM-score, was recently introduced to provide a combined structure quality measure of the coordinate root mean square deviation between a pair of structures and coverage. Using the TM-score, the TM-align structure alignment algorithm was developed that was often found to have better accuracy and coverage than the most commonly used structural alignment programs; however, there were a number of situations when this was not true. Results: To further improve structure alignment quality, the Fr-TM-align algorithm has been developed where aligned fragment pairs are used to generate the initial seed alignments that are then refined using dynamic programming to maximize the TM-score. For the assessment of the structural alignment quality from Fr-TM-align in comparison to other programs such as CE and TMalign, we examined various alignment quality assessment scores such as PSI and TM-score. The assessment showed that the structural alignment quality from Fr-TM-align is better in comparison to both CE and TM-align. On average, the structural alignments generated using Fr-TM-align have a higher TM-score (~9%) and coverage (~7%) in comparison to those generated by TM-align. Fr- TM-align uses an exhaustive procedure to generate initial seed alignments. Hence, the algorithm is computationally more expensive than TM-align. Conclusion: Fr-TM-align, a new algorithm that employs fragment alignment and assembly provides better structural alignments in comparison to TM-align. The source code and executables of Fr- TM-align are freely downloadable at: http://cssb.biology.gatech.edu/skolnick/files/FrTMalign/