Below is a list of programs I have developed and are available for download:
This program is designed for the analysis and resolution of articulated bar structures at their ends, a technique commonly used in metallic structural engineering. The program allows for solving both 2D and 3D structures, considering an approach based on first-order elasticity and using matrix methods to calculate displacements and stresses.
Articulated bar structures, also known as trusses or lattices, consist of straight elements connected at nodes that allow for the transmission of forces without generating moments. These structures are widely used in construction due to their ability to support large loads with relatively low weight.
The program is specifically designed to facilitate the analysis of this type of structure using efficient and precise routines. With it, engineers can quickly obtain internal forces and displacements of the bars, ensuring that the structure meets safety and functionality requirements.
This program is designed to solve structures formed by articulated bars at their ends, commonly used in metallic structures. It is based on first-order elasticity and matrix methods to calculate the displacements and forces in the structure’s nodes and bars. Below, we explain how to use the program and the key concepts behind its operation.
Once the above parameters are defined, the program uses the matrix method to solve the problem:
The program ensures that the displacements at the ends of the bars match the displacements at the corresponding nodes and that the equilibrium of forces is maintained at each node. Additionally, it uses first-order elasticity, meaning that the method assumes small deformations and that the structure remains within the elastic range, i.e., no plastic deformations occur.
Throughout the study of the human eye from a biomechanical perspective, a stochastic model has been developed that describes both the structure and dynamic behavior of the eyeball. This model takes into account the mechanical properties of the eye, its components, shape, and external forces acting upon it, such as the ocular muscles.
The eye is primarily modeled as an ellipsoidal system, with the cornea and the eyeball represented as prolate ellipsoids. This provides an accurate approximation of the shape of the eye. Geometric and mechanical parameters of the eye, such as corneal thickness, external radius, and asphericity, among others, are considered to create a three-dimensional model. Additionally, the normal oscillation modes of the eye, external forces, and the behavior of the eye as a rotating solid are considered, allowing for predictions of its dynamic response.
Download noteboksUse the following links to navigate through the notebooks:
Solar is an advanced program designed to calculate spectral transmittance and spectral irradiance, as well as other parameters essential for analyzing energy fluxes through the Earth's atmosphere. The program also calculates daily extraterrestrial irradiance, solar time, and various parameters critical for studying energy flows in the atmosphere.
The program takes into account complex atmospheric processes, such as Rayleigh scattering (scattering by small molecules and particles) and Mie scattering (scattering caused by larger particles like aerosols and water droplets). These phenomena are key in determining how solar radiation interacts with the Earth's atmosphere.
By allowing users to input atmospheric spectral data via detailed tables, the program can simulate a wide range of atmospheric conditions. This makes it a flexible and powerful tool for studies related to energy flux, solar energy, and climate science.
This program is designed to compute energy fluxes through the atmosphere, focusing on processes like spectral transmittance and irradiance. It incorporates both first-principles calculations and user-supplied atmospheric data to perform precise energy flow simulations. Below is a guide on how to use the program effectively:
Once the calculations are complete, the program provides output for a range of parameters, including:
The program ensures that atmospheric parameters like scattering and absorption are accurately factored into the calculation of energy fluxes. By taking into account real-time solar time and using first-order atmospheric models, it provides accurate data for solar energy applications and climate studies.
Download SolarPacumo is an advanced program designed to numerically solve the Schrödinger equation in the context of quantum mechanics, with the goal of simulating the interaction of a wave packet with a potential barrier. It uses implicit methods such as Crank-Nicholson for solving the discretized partial differential equation.
The program allows the study of quantum phenomena, such as reflection and transmission of waves, as well as free propagation in space. It is ideal for analyzing how a wave packet behaves before, during, and after interacting with the potential barrier.
Pacumo is designed to perform quantum simulations, analyzing the interaction of a wave packet with a potential barrier. Below is the usage process:
Once the simulation is complete, the program provides the following results:
Pacumo ensures that the total probability is conserved during the simulation, both for the reflected and transmitted parts. The program can handle complex numbers systems to solve differential equations accurately.
Download PacumoThe lighthouse problem is a classic exercise used to demonstrate how prior knowledge or information influences probability estimation in the context of Bayesian statistics.
Imagine a lighthouse located along the coastline at an unknown position α, along the beach, and at a fixed distance β offshore. The lighthouse emits pulses of light in random directions. A set of detectors on the coast can register the arrival of the pulses but not the exact direction from which they came. Based on the positions of the detected rays {xk} on the coast and knowing the distance β, the objective is to estimate the position α of the lighthouse.
Given that β and the positions of the detected rays {xk} are known, the mathematical relationship defining the position of a ray on the coast is given by the equation:
x_k = α + β * tan(θ_k)
Using Bayes' theorem, the posterior probability of the lighthouse's position α, given the information of the detected rays {xk}, is expressed as:
prob(α|{x_k}, β, I) ∝ prob(xk|α, β, I) * prob(α|β, I)
Where prob(α|β, I) is the prior probability of the position α, and prob(x_k|α, β, I) is the probability density function (PDF) describing the distribution of the detected rays, which follows a Cauchy distribution.
The key to the problem lies in how the probability is updated as more data from the detected rays is collected. With more information, the posterior probability becomes more precise and better approximates the true position of the lighthouse.
Bayesian statistics allow for the dynamic modeling of knowledge. The process begins with a prior probability, reflecting what is known (or unknown) about the lighthouse's position before observing the data. In the absence of information, the prior probability is usually uniform, meaning all positions along the coast are equally probable. However, if there is prior knowledge, for example, about a range where the lighthouse is likely to be, this knowledge can be incorporated into the initial estimation.
Each time a ray is detected on the coast, the prior probability is updated using the newly observed information. This updating process generates a new, more precise posterior probability. With each new observation, the Bayesian model adjusts the beliefs about the lighthouse's position, and the probability becomes increasingly concentrated around the true location.
If the prior knowledge is accurate or close to reality, the convergence process towards the correct solution is rapid. Even when the initial knowledge is incorrect or imprecise, collecting enough data allows for gradual corrections and refinement of the estimation.
Bayesian inference not only allows the incorporation of prior knowledge but also continuously improves it through the observation of data. This update cycle is iterative and dynamic: each new observation (each new detected ray) provides information that adjusts the posterior probability. As more observations are gathered, the model is refined and converges toward the true position of the lighthouse.
The use of Bayes' theorem enables continuous refinement of the lighthouse position estimate by effectively combining prior knowledge with new observations. If one starts from complete ignorance about the position of α, the collected data gradually adjust the posterior probability. If prior knowledge is reasonable, it accelerates convergence towards the correct solution.
This program is based on the doctoral thesis of María Rosario González Anera, professor at the University of Granada. Her research analyzes the Modulation Transfer Function (MTF) of the human eye, using different optical models of the eye obtained through detailed procedures.
The program allows loading several models of the human eye, compiled and provided by Pepe Fernández Dorado. These models simulate the polychromatic MTF in various configurations, and have been incorporated into the system by the author of this website. The project’s goal is to use this tool to enhance understanding of the MTF in Vision Physics practices since 2006.
This program was designed by the author of this website for students of the fifth year of the former Physics degree at the University of Granada, specifically for Vision Physics practices. With it, students can explore how human eye models influence the perceived image quality, analyzing aspects such as chromatic aberration and light dispersion.
The models implemented in this program include the classic Gullstrand model, as well as more recent models like Liou & Brennan, which take into account the refractive index variation within the lens. These simulations allow the study of image quality projected onto the retina under different optical conditions.